Current Fragmentation and Particle Acceleration in Solar Flares

Particle acceleration in solar flares remains an outstanding problem in plasma physics and space science. While the observed particle energies and timescales can perhaps be understood in terms of acceleration at a simple current sheet or turbulence site, the vast number of accelerated particles, and the fraction of flare energy in them, defies any simple explanation. The nature of energy storage and dissipation in the global coronal magnetic field is essential for understanding flare acceleration. Scenarios where the coronal field is stressed by complex photospheric motions lead to the formation of multiple current sheets, rather than the single monolithic current sheet proposed by some. The currents sheets in turn can fragment into multiple, smaller dissipation sites. MHD, kinetic and cellular automata models are used to demonstrate this feature. Particle acceleration in this environment thus involves interaction with many distributed accelerators. A series of examples demonstrate how acceleration works in such an environment. As required, acceleration is fast, and relativistic energies are readily attained. It is also shown that accelerated particles do indeed interact with multiple acceleration sites. Test particle models also demonstrate that a large number of particles can be accelerated, with a significant fraction of the flare energy associated with them. However, in the absence of feedback, and with limited numerical resolution, these results need to be viewed with caution. Particle in cell models can incorporate feedback and in one scenario suggest that acceleration can be limited by the energetic particles reaching the condition for firehose marginal stability. Contemporary issues such as footpoint particle acceleration are also discussed. It is also noted that the idea of a “standard flare model” is ill-conceived when the entire distribution of flare energies is considered.     

Towards a new paradigm: from a quasi-steady description to a dynamical description of the magnetosphere

A great deal of the research done on the dynamical process of the solar wind- magnetosphere interaction is based on large-scale, quasi-steady theoretical models, such as the classical reconnection model. However, it can be argued that the theoretical and observational foundations of these commonly believed paradigms are not always strong, and support for these models is sometimes weak, controversial or inconsistent. This paper discusses the need for a transition from an oversimplified quasi-steady paradigm towards a more realistic one including the dynamics of MHD waves and wave packets. The effects of localized wave packets may be most important in active plasma regions, where ideal MHD breaks down and localized, time-dependent processes become dominant. New insights into the theories of field-aligned current generation, auroral particle acceleration and the concept of reconnection may be found by including MHD wave propagation and wave packet dynamics.     

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.     

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.     

Kinetic Characterization of Plasma Sheet Dynamics

The Wind spacecraft made 26 perigee passes through the near-earth plasma sheet region during 1994 to 1997. Nearly all of these passes obtained plasma data from substorms and bursty bulk flow (BBF) events. New features of ion distributions have been observed in both the plasma sheet boundary layer (PSBL) and the central plasma sheet (CPS) in the vicinity of the current sheet that are relevant for understanding the structure of the PSBL and the mechanisms of particle acceleration to MeV energies associated with the BBF events. Kinetic processes are key to understanding these new observations that are not adequately explained by existing magnetohydrodynamics (MHD) models and theories. This article will feature the phase space distribution functions as the primary data product. The main purpose of this article is to establish an observational framework for new improved models and theories. The new observations should challenge modelers and theorists.     

Coronal streamers' theories

Some theoretical aspects of solar coronal streamers are discussed with emphasis on the current sheet and reconnection processes going on along the axis of the streamer. The dynamics of the streamer is a combination of MHD and transport, with acceleration of particles due to reconnection and leakage of plasma outwards as a slow solar wind as the observable results. The presence of the almost-closed magnetic bottles of streamers that can store high-energy particles for significant times provides the birdcage for solar cosmic rays, the reconnection in the sheet feeds medium-energy protons into the corona for the large-scale storage needed for certain flare models, and the build-up of excess density sets the stage for coronal mass ejections.     

Computer modeling of test particle acceleration at oblique shocks

We review the basic techniques and results of numerical codes used to model the acceleration of charged particles at oblique, fast-mode, collisionless shocks. The emphasis is upon models in which accelerated particles (ions) are treated as test particles, and particle dynamics is calculated by numerically integrating along exact phase-space orbits. We first review the case where ions are sufficiently energetic so that the shock can be approximated by a planar discontinuity, and where the electromagnetic fields on both sides of the shock are defined at the outset of each computer run. When the fields are uniform and static, particles are accelerated by the scatter-free drift acceleration process at a single shock encounter. We review the characteristics of scatter-free drift acceleration by considering how an incident particle distribution is modified by interacting with a shock. Next we discuss drift acceleration when magnetic fluctuations are introduced on both sides of the shock, and compare these results with those obtained under scatter-free conditions. We describe the modeling of multiple shock encounters, discuss specific applications, and compare the model predictions with theory. Finally, we review some recent numerical simulations that illustrate the importance of shock structure to both the ion injection process and to the acceleration of ions to high energies at quasi-perpendicular shocks.     

What auroral electron and ion beams tell us about magnetosphere-ionosphere coupling

The electron and ion beams which have been detected on many rockets and satellites are of particular interest because beam particles carry information about both the ionosphere and the magnetosphere out to the distant tail. Stability analyses have shown that even the most dramatic beams have evolved until the particle distribution functions are only weakly unstable. The shortest plasma wave growth lengths in the auroral region are usually comparable to the size of an arc. The resulting clearest electron beams generally are relatively minor features of distribution functions which are dominated by plateaus, loss cones, broad or stretched out field aligned features, and hot or cold isotropic components. The true electron beams therefore represent a small fraction of the total electron number density. Ion beams carry a much larger fraction of all ions, but also are only weakly unstable. The electron beams seen at low altitudes can drive whistlers (both electromagnetic and electrostatic, including lower hybrid waves) and upper hybrid waves, which may be particularly intense near electron gyroharmonics. Ion beams can drive low frequency electromagnetic waves that are related to gyrofrequencies of several ion species as well as ion acoustic and electrostatic ion cyclotron waves. These latter waves can be driven both by the drift of ion beams relative to cold stationary ions and by the drift of electrons relative to either stationary or drifting ions. Abrupt changes or boundaries in the electron and ion velocity space distribution functions (e.g. beams and loss cones) have been analyzed to provide information about the plasma source, acceleration process, and regions of strong wave-particle interactions. Fluid analyses have shown that upgoing ion beams carry a great deal of momentum flux from the ionosphere. This aspect of ion beams is analyzed by treating the entire acceleration region as a black box, and determining the forces that must be applied to support the upgoing beams. This force could be provided by moderate energy (10's of eV) electrons which are heated near the lower border of the acceleration region. It is difficult to use standard particle detectors to measure the particles which carry electric current in much of the magnetosphere. Such measurements may be relatively easy within upgoing ion beams because there is some evidence that few of the hard-to-measure cold plasma particles are present. Therefore, ion beam regions may be good places to study fluid or MHD properties of magnetospheric plasmas, including the identification of current carriers, a study of current continuity, and some aspects of the substorm and particle energization processes. Finally, some of the experimental results which would be helpful in an analysis of several magnetospheric problems are summarized.     

MHD Waves and Heating in Coronal Holes

Coronal holes have been identified as source regions of the fast solar wind, and MHD wave activity has been detected in coronal holes by remote sensing, and in situ in fast solar wind streams. I review some of the most suggestive wave observations, and discuss the theoretical aspects of MHD wave heating and solar wind acceleration in coronal holes. I review the results of single fluid 2.5D MHD, as well as multi-fluid 2.5D MHD models of waves in coronal holes, the heating, and the acceleration of the solar wind be these waves.     

Working Group 1 Report: Solar Wind Models from the Sun to 1 AU: Constraints by "in situ" and Remote Sensing Measurements

The goal of Working Group 1 was to discuss constraints on solar wind models. The topics for discussion, outlined by Eckart Marsch in his introduction, were: (1) what heats the corona, (2) what is the role of waves, (3) what determines the solar wind mass flux, (4) can stationary, multi-fluid models describe the fast and slow solar wind, or (5) do we need time dependent fluid models, kinetic models, and/or MHD models to describe solar wind acceleration. The discussion in the working group focused on observations of "temperatures" in the corona, mainly in coronal holes, and whether the observations of line broadening should be interpreted as thermal broadening or wave broadening. Observations of the coronal electron density and the flow speed in coronal holes were also discussed. There was only one contribution on observations of the distant solar wind, but we can place firm constraints on the solar wind particle fluxes and asymptotic flow speeds from observations with Ulysses and other spacecraft. Theoretical work on multi-fluid models, higher-order moment fluid models, and MHD models of the solar wind were also presented. This revised version was published online in June 2006 with corrections to the Cover Date.     

Particle Acceleration in Impulsive Solar Flares

We present the major observationally-derived requirements for a solar flare particle acceleration mechanism, briefly discuss some general electrodynamic constraints that also need to be considered, and suggest a unified electron and ion acceleration theory. This theory consists of two elements: cascading MHD turbulence generated at large scales during the primary flare energy release, which is responsible for the energization of electrons and all ions except 3He, and an electron beam, which excites the waves necessary for 3He acceleration. An issue of special importance for understanding ion acceleration is the convincing measurement of the charge state of Fe, which can be accomplished by the Advanced Composition Explorer in the upcoming solar maximum. This revised version was published online in June 2006 with corrections to the Cover Date.     

Particle Acceleration in the Heliosphere: Implications for Astrophysics

There has been a remarkable discovery concerning particles that are accelerated in the solar wind. At low energies, in the region where the particles are being accelerated, the spectrum of the accelerated particles is always the same: when expressed as a distribution function, the spectrum is a power law in particle speed with a spectral index of ?5, and a rollover at higher particle speeds that can often be described as exponential. This common spectral shape cannot be accounted for by any conventional acceleration mechanism, such as diffusive shock acceleration or traditional stochastic acceleration. It has thus been necessary to invent a new acceleration mechanism to account for these observations, a pump mechanism in which particles are pumped up in energy through a series of adiabatic compressions and expansions. The conditions under which the pump acceleration is the dominant acceleration mechanism are quite general and are likely to occur in other astrophysical plasmas. In this paper, the most compelling observations of the ?5 spectra are reviewed; the governing equation of the pump acceleration mechanism is derived in detail; the pump acceleration mechanism is applied to acceleration at shocks; and, as an illustration of the potential applicability of the pump acceleration mechanism to other astrophysical plasmas, the pump mechanism is applied to the acceleration of galactic cosmic rays in the interstellar medium.     

CME Theory and Models

This chapter provides an overview of current efforts in the theory and modeling of CMEs. Five key areas are discussed: (1) CME initiation; (2) CME evolution and propagation; (3) the structure of interplanetary CMEs derived from flux rope modeling; (4) CME shock formation in the inner corona; and (5) particle acceleration and transport at CME driven shocks. In the section on CME initiation three contemporary models are highlighted. Two of these focus on how energy stored in the coronal magnetic field can be released violently to drive CMEs. The third model assumes that CMEs can be directly driven by currents from below the photosphere. CMEs evolve considerably as they expand from the magnetically dominated lower corona into the advectively dominated solar wind. The section on evolution and propagation presents two approaches to the problem. One is primarily analytical and focuses on the key physical processes involved. The other is primarily numerical and illustrates the complexity of possible interactions between the CME and the ambient medium. The section on flux rope fitting reviews the accuracy and reliability of various methods. The section on shock formation considers the effect of the rapid decrease in the magnetic field and plasma density with height. Finally, in the section on particle acceleration and transport, some recent developments in the theory of diffusive particle acceleration at CME shocks are discussed. These include efforts to combine self-consistently the process of particle acceleration in the vicinity of the shock with the subsequent escape and transport of particles to distant regions.     

The Cusp as a Source of Magnetospheric Energetic Particles, Currents, and Electric Fields: A New Paradigm

Energetic particle instrumentation on the Polar satellite has discovered that significant fluxes of energetic particles are continuously present in the region of the dayside magnetosphere where they cannot be stably trapped. This region is associated with either open magnetic field lines or a magnetic topology associated with pseudo-trapping. Two distinct features [Time-Energy Dispersion (TED) signatures and Cusp Energetic Particle (CEP) events] are observed in these energetic particle fluxes that strongly suggest a local acceleration of mostly shocked solar wind particles. As the solar wind particles ram themselves into the cusp geometry, they form diamagnetic cavities with strong turbulence that are capable of accelerating particles to energies of 100s and 1000s of kiloelectronvolts. This process forms a layer of energetic particles on the magnetopause as well as permits such particles to enter via drift the equatorial nightside magnetosphere to distances as close as six Earth radii under the influence of gradient and curvature effects in the local magnetic field. The fluxes of these particles have all of the properties associated with the ring current and can supply the magnitude of the cross tail current required. ISEE-1 energetic particle data and their pitch angle distributions [PAD] are examined at the magnetic equatorial plane on the night side to investigate and possibly validate the insights gains from the Polar data and energetic particle trajectory tracing in a realistic magnetic field. The existence and properties of butterfly-type PADs strongly supports the concept of a dayside high latitude source of energetic particle fluxes. Because the CEP process is impulsive and time variable the charge separation produced by the drifting electrons (eastward) and ions (westward) on the magnetospheric nightside may be responsible for the cross tail electric field that has been ascribed to the reconnection/convection process.     

Modelling mesoscale processes in the global geospace system

Intermediate or mesoscale processes mediate the transfer of mass, momentum, and energy across the dynamic solar wind-magnetosphere interface, and the propagation of this input through the system to the ionosphere and atmosphere. The Dartmouth-Berkeley-Minnesota theory team has identified a number of mesoscale phenomena to be investigated as part of the GGS program, including: (1) effects of upstream density fluctuations on magnetopause dynamics, (2) three-dimensional reconnection, (3) magnetopause depletion layer studies, (4) ring current interaction with Pc 1 and Pc 5 waves, (5) generation of ion Larmor-scale current layers in the near Earth plasmasheet, (6) test particle studies in the magnetotail, (7) simulation of magnetosphere- ionosphere coupling including effects of kinetic Alfvén waves and (8) auroral acceleration region studies of the effects of kinetic Alfvén waves on particle distribution functions. A broad range of techniques will be implemented including ideal and reduced MHD, two fluid, hybrid, particle-in-cell and test particle simulations. Detailed comparison of simulation results with GGS satellite and ground based data will be undertaken.     

Particle Acceleration at Interplanetary Shocks

This paper briefly reviews proton acceleration at interplanetary shocks. This is key to describing the acceleration of heavy ions at interplanetary shocks because wave excitation—and hence particle scattering—at oblique shocks is controlled by the protons and not the heavy ions. Heavy ions behave as test particles, and their acceleration characteristics are controlled by the properties of proton-excited turbulence. As a result, the resonance condition for heavy ions introduces distinctly different signatures in abundance, spectra, and intensity profiles, depending on ion mass and charge. Self-consistent models of heavy-ion acceleration and the resulting fractionation are discussed. This includes discussion of the injection problem and the acceleration characteristics of quasi-parallel and quasi-perpendicular shocks.     

Electron Acceleration in the Heliosphere

We review the evidence for electron acceleration in the heliosphere putting emphasis on the acceleration processes. There are essentially four classes of such processes: shock acceleration, reconnection, wave particle interaction, and direct acceleration by electric fields. We believe that only shock and electric field acceleration can in principle accelerate electrons to very high energies. The shocks known in the heliosphere are coronal shocks, traveling interplanetary shocks, CME shocks related to solar type II radio bursts, planetary bow shocks, and the termination shock of the heliosphere. Even in shocks the acceleration of electrons requires the action of wave particle resonances of which beam driven whistlers are the most probable. Other mechanisms of acceleration make use of current driven instabilities which lead to electron and ion hole formation. In reconnection acceleration is in the current sheet itself where the particles perform Speiser orbits. Otherwise, acceleration takes place in the slow shocks which are generated in the reconnection process and emanate from the diffusion region in the Petschek reconnection model and its variants. Electric field acceleration is found in the auroral zones of the planetary magnetospheres and may also exist on the sun and other stars including neutron stars. The electric potentials are caused by field aligned currents and are concentrated in narrow double layers which physically are phase space holes in the ion and electron distributions. Many of them add up to a large scale electric field in which the electrons may be impulsively accelerated to high energies and heated to large temperatures.     

Acceleration of Particles to High Energies in Earth’s Radiation Belts

Discovered in 1958, Earth’s radiation belts persist in being mysterious and unpredictable. This highly dynamic region of near-Earth space provides an important natural laboratory for studying the physics of particle acceleration. Despite the proximity of the radiation belts to Earth, many questions remain about the mechanisms responsible for rapidly energizing particles to relativistic energies there. The importance of understanding the radiation belts continues to grow as society becomes increasingly dependent on spacecraft for navigation, weather forecasting, and more. We review the historical underpinning and observational basis for our current understanding of particle acceleration in the radiation belts.     

Stochastic Acceleration of Energetic Particles in the Heliosphere

In this paper we assess possible roles of stochastic acceleration by random electric field and plasma motion in the production and transport of energetic particles in the heliosphere. Stochastic acceleration can occur in the presence of multiple small-scale magnetohydrodynamic waves propagating in different directions. Usually, this type of stochastic acceleration is closely related to particle pitch angle scattering or parallel diffusion. Given the values of the parallel diffusion coefficient inferred from the observations of cosmic ray modulation or other energetic particle phenomena in the heliosphere, stochastic acceleration by small-scale waves is much slower than acceleration by shock waves and it is also much slower than adiabatic cooling by the expansion of the solar wind; thus it is considered as inefficient for producing heliospheric energetic particles or for the modulation of cosmic rays. Another type of stochastic acceleration occurs when particles go through random compressions or expansions due to large-scale plasma motion. This acceleration mechanism could be very fast when the correlation time of the fluctuations in plasma compression is short compared to the diffusion time. Particle acceleration by an ensemble of small shock waves or intermittent long wavelength compressible turbulence belongs to this category. It tends to establish an asymptotic p ^?3 universal distribution function quickly if there is no or little large-scale adiabatic cooling. Such a particle distribution will contain an infinite amount of pressure. Back reaction from the pressure is expected to modify the amplitude of plasma waves to an equilibrium state. At that point, the pressure of accelerated particles must remain finite and the accelerated particles could approach a p ^?5 distribution function.     

The plasma physics of shock acceleration

The notion that plasma shocks in astrophysical settings can and do accelerate charged particles to high energies is not a new one. However, in recent years considerable progress has been achieved in understanding the role particle acceleration plays both in astrophysics and in the shock process itself. In this paper we briefly review the history and theory of shock acceleration, paying particular attention to theories of parallel shocks which include the backreaction of accelerated particles on the shock structure. We discuss in detail the work that computer simulations, both plasma and Monte Carlo, are playing in revealing how thermal ions interact with shocks and how particle acceleration appears to be an inevitable and necessary part of the basic plasma physics that governs collisionless shocks. We briefly describe some of the outstanding problems that still confront theorists and observers in this field.