HYDRODYNAMIC AND MHD EQUATIONS ACROSS THE BOW SHOCK AND ALONG THE SURFACES OF PLANETARY OBSTACLES

Examinations of the magnetohydrodynamic (MHD) equations across a bow shock are presented. These equations are written in the familiar Rankine–Hugoniot set, and an exact solution to this set is given which involves the upstream magnetosonic Mach number, plasma , polytropic index, and _B-v , as a function of position along the shock surface. The asymptotic Mach cone angle of the shock surface is also given as a function of the upstream parameters, as a set of transcendental equations. The standoff position of a detached bow shock from an obstacle is also reviewed. In addition, a detailed examination of the hydrodynamic equations along the boundary of the obstacle is performed. Lastly, the MHD relations along the obstacle surface are examined, for specific orientations of the upstream interplanetary magnetic field (IMF) in relation to the upstream flow velocity vector.     

Finite Larmor radius effects in the magnetosphere

This article reviews theories and observations related to effects produced by finite (and large) Larmor radii of charged particles in the magnetosphere. The FLR effects depend on =r _H /L, wherer _H is the Larmor radius andL is the spatial scale for field/plasma inhomogeneity. The parameter is a basic expansion parameter for most equations describing plasma dynamics in the magnetosphere. The FLR effects enter naturally the drift approximation for particle motion and represent also non-ideal MHD terms in the fluid formalism. The linear and higher order terms in lead to charge separation, energization of particles, and produce viscosity without collisions. The FLR effects introduce also important corrections to the dispersion relations for MHD waves and drift instabilities. Expansion of plasma into magnetic field leads to filamentation of the plasma boundary and to creation of structures with thickness less than an ion gyroradius. Large Larmor radius effects (1) in curved magnetic field geometry lead to stochastic behaviour of particle trajectories and to deterministic chaos. The tiny scale of the electron and ion gyroradii does not necessarily mean that FLR/LLR phenomena have negligible effect on the macroscopic dynamics and energetics of the whole magnetosphere. On the contrary, the small scale gyro-effects may provide the physical mechanism for gyroviscous coupling between the solar wind and the magnetosphere, the mechanism for triggering disruption of the magnetotail current layer, and the mechanism for parallel electric field that accelerate auroral particles.     

Possible shock wave in the local interstellar plasma, very close to the heliosphere

We show, using the HST — GHRS data on velocity and temperature in the nearby interstellar medium, that the observed 3 – 4 km s^–1 relative velocity between the Local Interstellar Cloud (LIC) and the so-called G-cloud located in the Galactic Center hemisphere can be quite naturally explained assuming that the two clouds do interact with each other. In the proposed interpretation the two media are separated by a (quasiperpendicular) MHD shock front propagating from the LIC into the G-cloud. The LIC plasma is then nothing else but the shocked (compression 1.3 – 1.4) gas of the G-cloud. A 1-D single-fluid solution of the Rankine — Hugoniot equations can fit the most probable observed values of the relative velocity (3.75 km/s), LIC (6700 K) and G-cloud (5400 K) kinetic temperatures, if the plasma-beta of the LIC plasma is in the range 1.3 – 1.5 (Table 1). This corresponds to a super — fast magnetosonic motion of the heliosphere through the LIC, independently of LIC density. The LIC magnetic field strength is 1.9 (3.1) G for the LIC electron density n_e = 0.04 (0.10) cm^–3. In this case the shock is less than 30 000 AU away and moves at about 10 km s^–1 relative to the LIC plasma. The Sun is chasing the shock and should catch up with it in about 10⁴ years. If the heliospheric VLP emissions cutoff at 1.8 kHz is indicative of n_e (LIC) = 0.04 cm^–3 (Gurnett et al., 1993), the (pure plasma) bowshock ahead of the heliopause could be the source of quasi-continuous heliospheric 2-kHz emission band. We believe that with the expected increase in the performance of modern spectroscopic instrumentation the proposed method of magnetic field evaluation may in the future find wider application in the studies of the interstellar medium.     

Solar Wind Turbulence and the Role of Ion Instabilities

Solar wind is probably the best laboratory to study turbulence in astrophysical plasmas. In addition to the presence of magnetic field, the differences with neutral fluid isotropic turbulence are: (i) weakness of collisional dissipation and (ii) presence of several characteristic space and time scales. In this paper we discuss observational properties of solar wind turbulence in a large range from the MHD to the electron scales. At MHD scales, within the inertial range, turbulence cascade of magnetic fluctuations develops mostly in the plane perpendicular to the mean field, with the Kolmogorov scaling $k_{\perp}^{-5/3}$ for the perpendicular cascade and $k_{\|}^{-2}$ for the parallel one. Solar wind turbulence is compressible in nature: density fluctuations at MHD scales have the Kolmogorov spectrum. Velocity fluctuations do not follow magnetic field ones: their spectrum is a power-law with a ?3/2 spectral index. Probability distribution functions of different plasma parameters are not Gaussian, indicating presence of intermittency. At the moment there is no global model taking into account all these observed properties of the inertial range. At ion scales, turbulent spectra have a break, compressibility increases and the density fluctuation spectrum has a local flattening. Around ion scales, magnetic spectra are variable and ion instabilities occur as a function of the local plasma parameters. Between ion and electron scales, a small scale turbulent cascade seems to be established. It is characterized by a well defined power-law spectrum in magnetic and density fluctuations with a spectral index close to ?2.8. Approaching electron scales, the fluctuations are no more self-similar: an exponential cut-off is usually observed (for time intervals without quasi-parallel whistlers) indicating an onset of dissipation. The small scale inertial range between ion and electron scales and the electron dissipation range can be together described by $\sim k_{\perp}^{-\alpha}\exp(-k_{\perp}\ell_{d})$ , with α?8/3 and the dissipation scale ? _d close to the electron Larmor radius ? _d ?ρ _e . The nature of this small scale cascade and a possible dissipation mechanism are still under debate.     

Quasi-Periodic Pulsations in Solar Flares

Quasi-periodic pulsations (QPP) are a common feature of flaring energy releases in the solar atmosphere, observed in all bands, from radio to hard X-ray. In this review we concentrate on QPP with the periods longer than one second. Physical mechanisms responsible for the generation of long QPP split into two groups: “load/unload” mechanisms and MHD oscillations. Load/unload mechanisms are repetitive regimes of flaring energy releases by magnetic reconnection or by other means. MHD oscillations can affect all elements of the flaring emission generation: triggering of reconnection and modulation of its rate, acceleration and dynamics of non-thermal electrons, and physical conditions in the emitting plasmas. In the case of MHD oscillations, the periodicity of QPP is determined either by the presence of some resonances, e.g. standing modes of plasma structures, or by wave dispersion. Periods and other parameters of QPP are linked with properties of flaring plasmas and their morphology. Observational investigation of the QPP generation mechanisms based upon the use of spatial information, broadband spectral coverage and multi-periodicity is discussed.     

Magnetic field reconnection theory and the solar wind — Magnetosphere interaction: A review

This review considers the theory of the magnetic field line reconnection and its application to the problem of the interaction between the solar wind and the Earth's magnetosphere. In particular, we discuss the reconnection models by Sonnerup and by Petschek (for both incompressible and compressible plasmas, for the asymmetric and nonsteady-state cases), the magnetic field annihilation model by Parker; Syrovatsky's model of the current sheet; and Birn's and Schindler's solution for the plasma sheet structure. A review of laboratory and numerical modelling experiments is given.Results concerning the field line reconnection, combined with the peculiarities of the MHD flow, were used in investigating the solar wind flow around the magnetosphere. We found that in the presence of a frozen-in magnetic field, the flow differs significantly from that in a pure gas dynamic case; in particular, at the subsolar. part of the magnetopause a stagnation line appears (i.e., a line along which the stream lines are branching) instead of a stagnation point. The length and location of the stagnation line determine the character of the interaction of the solar wind with the Earth's magnetosphere. We have developed the theory of that interaction for a steady-state case, and compare the results of the calculations with the experimental data.In the last section of the review, we propose a qualitative model of the solar wind — the Earth's magnetosphere interaction in the nonsteady-state case on the basis of the solution of the problem of the spontaneous magnetic field line reconnection.     

Why space physics needs to go beyond the MHD box

Magnetohydrodynamic (MHD) theory has been used in space physics for more than forty years, yet many important questions about space plasmas remain unanswered. We still do not understand how the solar wind is accelerated, how mass, momentum and energy are transported into the magnetosphere and what mechanisms initiate substorms. Questions have been raised from the beginning of the space era whether MHD theory can describe correctly space plasmas that are collisionless and rarely in thermal equilibrium. Ideal MHD fluids do not induce electromotive force, hence they lose the capability to interact electromagnetically. No currents and magnetic fields are generated, rendering ideal MHD theory not very useful for space plasmas. Observations from the plasma sheet are used as examples to show how collisionless plasmas behave. Interpreting these observations using MHD and ideal MHD concepts can lead to misleading conclusions. Notably, the bursty bulk flows (BBF) with large mean velocities left(〈 v 〉≥400 km s right) that have been interpreted previously as E×B flows are shown to involve much more complicated physics. The sources of these nonvanishing 〈 v 〉 events, while still not known, are intimately related to mechanisms that create large phase space gradients that include beams and acceleration of ions to MeV energies. The distributions of these nonvanishing 〈 v 〉 events are associated with large amplitude variations of the magnetic field at frequencies up to and exceeding the local Larmor frequency where MHD theory is not valid. Understanding collisionless plasma dynamics such as substorms in the plasma sheet requires the self-consistency that only kinetic theory can provide. Kinetic modeling is still undergoing continual development with many studies limited to one and two dimensions, but there is urgent need to improve these models as more and more data show kinetic physics is fundamentally important. Only then will we be able to make progress and obtain a correct picture of how collisionless plasmas work in space.     

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.     

Self-consistent Simulations of Alfv��n Wave Driven Winds from the Sun and Stars

We review our recent results of Alfvén wave-driven winds. First, we present the result of self-consistent 1D MHD simulations for solar winds from the photosphere to interplanetary region. Here, we emphasize the importance of the reflection of Alfvén waves in the density stratified corona and solar winds. We also introduce the recent Hinode observation that might detect the reflection signature of transverse (Alfvénic) waves by Fujimura and Tsuneta (Astrophys. J. 702:1443, 2009). Then, we show the results of Alfvén wave-driven winds from red giant stars. As a star evolves to the red giant branch, the properties of stellar winds drastically change from steady coronal winds to intermittent chromospheric winds. We also discuss how the stellar evolution affects the wave reflection in the stellar atmosphere and similarities and differences of accretion disk winds by MHD turbulence.     

Cosmic-ray acceleration and transport,and diffuse galactic γ-ray emission

Cosmic-ray acceleration and transport is considered from the point of view of application to diffuse galactic -ray sources. As an introduction we review several source models, in particular supernovae exploding inside or near large interstellar clouds. The complex problem of cosmic ray transport in random electromagnetic fields is reduced to three cases which should be sufficient for practical purposes. As far as diffusive acceleration is concerned, apart from reviewing the basic physical principles, we point out the relation between shock acceleration and 2nd order Fermi acceleration, and the relative importance of the two processes around interstellar shock waves. For -ray source models the interaction of cosmic rays with dense clouds assumes great importance. Past discussions had been confined to static interactions of clouds with the ambient medium in the sense that no large scale mass motions in the ambient interstellar medium were considered. The well-known result then is that down to some tens of MeV or less, cosmic-ray nucleons should freely penetrate molecular clouds of typical masses and sizes. The self-exclusion of very low energy nucleons however may affect electron transport with consequences for the Bremsstrahlung -luminosity of such clouds.In this paper we consider also the dynamical interaction of dense clouds with a surrounding hot interstellar medium. Through cloud evaporation and accretion there exist mass flows in the cloud surroundings. We argue that in the case of (small) cloud evaporation the galactic cosmic rays will be essentially excluded from the clouds. The dynamic effects of cosmic rays on the flow should be minor in this case. For the opposite case of gas accretion onto (large) clouds, cosmic-ray effects on the flow will in general be large, limiting the cosmic-ray compression inside the cloud to dynamic pressure equilibrium. This should have a number of interesting and new consequences for -ray astronomy. A first, qualitative discussion is given in the last section.Proceedings of the XVIII General Assembly of the IAU: Galactic Astrophysics and Gamma-Ray Astronomy, held at Patras, Greece, 19 August 1982.     

Turbulent transition in solar surges

It has been suggested that a surge can be modelled as a jet travelling in a sheared magnetic field, and that the transition to turbulence of this MHD tearing jet can explain several key observed features. In this paper we present our preliminary results of the transition to turbulencevia secondary instabilities of the MHD tearing jet. Our results confirm that turbulent transition can decelerate the surge, with decay times which compare well with surge data. Furthermore, we find that the turbulent MHD tearing jet forms magnetic field-aligned velocity filaments similar to those often observed in the surge flow field.     

Cosmological Distance Indicators

We review three distance measurement techniques beyond the local universe: (1) gravitational lens time delays, (2) baryon acoustic oscillation (BAO), and (3) HI intensity mapping. We describe the principles and theory behind each method, the ingredients needed for measuring such distances, the current observational results, and future prospects. Time-delays from strongly lensed quasars currently provide constraints on \(H_{0}\) with \(<4\%\) uncertainty, and with \(1\%\) within reach from ongoing surveys and efforts. Recent exciting discoveries of strongly lensed supernovae hold great promise for time-delay cosmography. BAO features have been detected in redshift surveys up to \(z\lesssim0.8\) with galaxies and \(z\sim2\) with Ly-\(\alpha\) forest, providing precise distance measurements and \(H_{0}\) with \(<2\%\) uncertainty in flat \(\Lambda\)CDM. Future BAO surveys will probe the distance scale with percent-level precision. HI intensity mapping has great potential to map BAO distances at \(z\sim0.8\) and beyond with precisions of a few percent. The next years ahead will be exciting as various cosmological probes reach \(1\%\) uncertainty in determining \(H_{0}\), to assess the current tension in \(H_{0}\) measurements that could indicate new physics.     

Theory and Modeling of Planetary Dynamos

Numerical dynamo models are increasingly successful in modeling many features of the geomagnetic field. Moreover, they have proven to be a useful tool for understanding how the observations connect to the dynamo mechanism. More recently, dynamo simulations have also ventured to explain the surprising diversity of planetary fields found in our solar system. Here, we describe the underlying model equations, concentrating on the Boussinesq approximations, briefly discuss the numerical methods, and give an overview of existing model variations. We explain how the solutions depend on the model parameters and introduce the primary dynamo regimes. Of particular interest is the dependence on the Ekman number which is many orders of magnitude too large in the models for numerical reasons. We show that a minor change in the solution seems to happen at $\mbox {E}=3\mbox {$\times 10^{-6}$}$ whose significance, however, needs to be explored in the future. We also review three topics that have been a focus of recent research: field reversal mechanisms, torsional oscillations, and the influence of Earth’s thermal mantle structure on the dynamo. Finally we discuss the possibility of tidally or precession driven planetary dynamos.     

Kinetic-Gasdynamic Modeling of the Heliospheric Interface

Heliospheric energetic neutral atoms (ENAs) that will be measured by the Interstellar Boundary Explorer (IBEX) originate from the heliosheath. The heliosheath is formed as a result of the interaction of the solar wind (SW) with the circum-heliospheric interstellar medium (CHISM). The expected fluxes of ENAs are strongly dependent on the nature of this interaction. In turn, the interaction of the solar wind with the local interstellar cloud has a complex and multi-component nature. Detailed theoretical modeling of the interaction between the SW and the local interstellar medium is required to understand the physics of the heliosheath and to predict and explain the heliospheric ENAs. This paper summarizes current state-of-art kinetic-gasdynamic models of the SW/CHISM interaction. We shall restrict our discussion to the kinetic-gasdynamic and kinetic-magnetohydrodynamic (MHD) models developed by the Moscow group. This paper summarizes briefly the main results of the first self-consistent, two-component, kinetic-gasdynamic model by Baranov and Malama (J. Geophys. Res. 98:15157–15163, 1993), presents new results obtained from the 3D kinetic-MHD model by Izmodenov et al. (Astron. Astrophys. 437:L35–L38, 2005a), describes the basic formulation and results of the model by Malama et al. (Astron. Astrophys. 445:693–701, 2006) as well as reports current developments in the model. This self-consistent model considers pickup protons as a separate non-equilibrium component. Then we discuss a stochastic acceleration model for pickup protons in the supersonic solar wind and in the heliosheath. We also present the expected heliospheric ENA fluxes obtained in the framework of the models.     

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.     

Order out of Randomness: Self-Organization Processes in Astrophysics

Self-organization is a property of dissipative nonlinear processes that are governed by a global driving force and a local positive feedback mechanism, which creates regular geometric and/or temporal patterns, and decreases the entropy locally, in contrast to random processes. Here we investigate for the first time a comprehensive number of (17) self-organization processes that operate in planetary physics, solar physics, stellar physics, galactic physics, and cosmology. Self-organizing systems create spontaneous “order out of randomness”, during the evolution from an initially disordered system to an ordered quasi-stationary system, mostly by quasi-periodic limit-cycle dynamics, but also by harmonic (mechanical or gyromagnetic) resonances. The global driving force can be due to gravity, electromagnetic forces, mechanical forces (e.g., rotation or differential rotation), thermal pressure, or acceleration of nonthermal particles, while the positive feedback mechanism is often an instability, such as the magneto-rotational (Balbus-Hawley) instability, the convective (Rayleigh-Bénard) instability, turbulence, vortex attraction, magnetic reconnection, plasma condensation, or a loss-cone instability. Physical models of astrophysical self-organization processes require hydrodynamic, magneto-hydrodynamic (MHD), plasma, or N-body simulations. Analytical formulations of self-organizing systems generally involve coupled differential equations with limit-cycle solutions of the Lotka-Volterra or Hopf-bifurcation type.     

The Effects of Transverse Magnetic Field and Density Variations on the Particle Energy Spectra in a Reconnecting 3D Current Sheet

Electron and proton acceleration by a super-Dreicer electric field is further investigated in a non-neutral reconnecting current sheet (RCS) with a variable plasma density. The tangential B _z and transverse magnetic field components B _x are assumed to vary with the distances x and z from the X nullpoint linearly and exponentially, respectively; the longitudinal component (a ‘guiding field’) is accepted constant. Particles are found to gain a bulk of their energy in a thin region close to the X nullpoint where the RCS density increases with z exponentially with the index λ and the tangential magnetic field B _x also increases with z exponentially with the index α. For the RCS with a constant density (λ = 0), the variations of the tangential magnetic field lead to particle power-law energy spectra with the spectral indices γ₁ being dependent on the exponent α as: for protons and for electrons in a strong guiding field (β > 10⁻²) and for electrons in a moderate or weak guiding field (β > 10⁻⁴). For the RCS with an exponential density increase in the vicinity of the X nullpoint (λ≥ 0) there is a further increase of the resulting spectral indices γ that depends on the density exponent index λ as for protons and for electrons in weaker guiding fields and as for electrons in stronger guiding fields. These dependencies can explain a wide variety (1.5–10) of particle spectral indices observed in solar flares by the variations of a magnetic field topology and physical conditions in a reconnecting region. This can be used as a diagnostic tool for the investigation of the RCS dynamics from the accelerated particle spectra found from hard X-ray and microwave emission.     

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 by magnetic reconnection and shocks during current loop coalescence in solar flares

This article reviews recent development of the theory of current loop coalescence and shock waves, giving particular attention to particle acceleration caused by these processes. First, explosive reconnection driven by the current loop coalescence and associated particle acceleration are studied by theoretical and magnetohydrodynamic simulation methods and the results are compared with observations of solar flares; this model gives a good explanation for the quasi-periodic structure of some solar flare bursts. Next follows a discussion of particle acceleration in association with fast magnetosonic shock waves. It is shown theoretically and by relativistic particle simulation that a quasi-perpendicular shock wave can accelerate trapped ions in the direction perpendicular to the ambient magnetic field up to speeds much greater than the Alfvén speed, . When the ambient magnetic field is rather strong ( _{ce pe} ), both ions and electrons can be accelerated to relativistic energies. For both the nonrelativistic and relativistic cases, the time needed for the acceleration is very short; it is for the ions. These results are compared with the rapid and simultaneous acceleration of ions and electrons in the impulsive phase of solar flares.