Towards the circuit theory of solar flares

It has been shown that the main problems of the circuit theory of solar flares - unlikely huge current growth time and the origin of the current interruption - have been resolved considering the case of magnetic loop emergence and the correct application of Ohm's law. The generalized Ohm's law for solar flares is obtained. The conditions for flare energy release are as follows: large current value, > 10¹¹ A, nonsteady-state character of the process, and the existence of a neutral component in a flare plasma. As an example, the coalescence of a flare loop and a filament is considered. It has been shown that the current dissipation has increased drastically as compared with that in a completely ionized plasma. The current dissipation provides effective Joule heating of the plasma and particle acceleration in a solar flare. The ion-atom collisions play the decisive role in the energy release process. As a result the flare loop resistance can grow by 8–10 orders of magnitude. For this we do not need the anomalous resistivity driven by small-scale plasma turbulence. The energy release emerging from the upper part of a flare loop stimulates powerful energy release from the chromospheric level.     

Small-scale Langmuir wave instability in preflare chromosphere of solar active region

Necessary conditions have been investigated for the appearance of instability of high-frequency electron Langmuir waves in plasma of solar chromosphere near the foot-point of loop structure. We have considered the earliest stage of a flare process in solar active region. At the chromospheric part of current circuit of a flare loop such instability can appear and develop as the result of combined action of large-scale electric field, Landau damping and collisional processes in preflare plasma. We have investigated the process of instability development for two possible scenarios: (a) when preflare loop plasma has a classical Coulomb conductivity and (b) when anomalous resistance appears due to saturation of Bernstein turbulence. The growth rates of instability have been obtained and analyzed in detail. It has been assumed in the process of calculation that preflare plasma can be described by the FAL model of the solar atmosphere, which takes into account the process of helium diffusion. It has been shown that Langmuir wave instability can appear in its marginal form in the area under investigation either in the presence of Coulomb conductivity or in the presence of saturated Bernstein turbulence. Existence of instability with the growth rate, which changes its sign, proves the principal possibility of generation of nondamping Langmuir waves with small amplitudes.     

Plasma waves caused by transient heat conduction in a coronal loop and anomalous resistivity

Numerical simulation is made of the transient heat conduction during local heating in a model coronal magnetic loop with an axial electric current. It is assumed that a segment near the top of the normal coronal loop is heated to above 10⁷ K by a sufficiently small heat input as compared with the total flare energy. A hump appears in the velocity distribution of electrons moving down the temperature gradient with speeds slightly below the thermal one. Consequently, electron plasma waves are excited. The high intensity of the waves persists in the upper region of the loop for more than a second until the termination of the simulation. The energy density of the plasma waves normalized with respect to thermal density is 10^–3.5 at maximum. A theoretical estimate gives an anomalous resistivity 5 orders of magnitude greater than an initial value. Based on the above result, we propose a model for impulsive loop flares.     

太阳爆发过程中的大尺度磁重联电流片:理论和观测

从理论和观测两个方面来介绍和讨论出现在太阳爆发过程中的磁重联电流片及其物理本质和动力学特征。首先介绍在理论研究和理论模型中,磁重联电流片是如何在爆发磁结构当中形成并发展的,对观测研究有什么指导意义。然后介绍观测工作是从哪几个方面对理论模型预测的电流片进行证认和研究的。第三,将介绍观测研究给出了哪些过去所没有能够预期的结果,这些结果对深入研究耀斑一CME电流片以及其中的磁重联过程的理论工作有什么重要的、挑战性的意义。第四,讨论最新的与此有关的理论研究和数值实验。最后,对未来的研究方向和重要课题进行综述和展望。     

Current sheet models of solar flares

Current sheets have been suggested as the site for flare energy release because they can convert magnetic energy very rapidly into both heat and directed plasma energy. Also they contain electric fields with the potential of accelerating particles to high energies.The basic properties of current sheets are first reviewed. For instance, magnetic flux may be carried into a current sheet and annihilated. An exact solution for such a process in an infinitely long sheet has been found; it describes the annihilation of fields which are inclined at any angle, not just 180°. Moreover, field lines which are expelled from the ends of a current sheet can be described as having been reconnected. The only workable model for fast reconnection in the solar atmosphere, namely Petschek's mechanism, has recently been put on a firm foundation; it gives a reconnection rate which depends on the electrical conductivity but is typically a tenth or a hundredth of the Alfvén speed. A current sheet may be formed when the sources of an initially potential field start to move; a simple analytic technique for finding the position and shape of such a sheet in two dimensions now exists. Finally, a sheet with no transverse magnetic field component is subject to the tearing-mode instability, which rapidly produces a series of loops in the field.The main ways in which current sheets have been used for solar flare models is described. Syrovatskii's mechanism relies on the increase of the electric current density during the formation of a sheet, to a value in excess of the critical value j ^* for the onset of microinstabilities. But Anzer has recently demonstrated that the critical value is most unlikely to be reached during the initial formation process. Sturrock, on the other hand, has advocated the occurrence of the tearing-mode instability in an open streamer-like configuration (which may result from the eruption of a force-free field). But recent observations do not point to that as the relevant configuration. Rather, they suggest that flares are triggered by the emergence of new magnetic flux from below the solar photosphere. This has led Heyvaerts, Priest, and Rust (1976) to propose a new emerging flux model, according to which, as more and more flux emerges, so reconnection occurs, producing some preflare heating. When the current sheet reaches such a height (around the transition region) that its current density exceeds j ^*, then the impulsive phase of the flare is triggered. The main phase is caused by an enhanced level of magnetic energy conversion in a turbulent current sheet. The type of flare depends on the magnetic environment in which the emerging flux finds itself. A surge flare results if the flux appears near a strong unipolar region such as a simple sunspot, whereas a two ribbon flare may be produced by flux emergence near an active region filament, in which case the main phase energy is released from the field that surrounds the filament.     

Dynamics of cosmic current layers with localized lower-hybrid-drift turbulence

A two-dimensional magnetohydrodynamic model of the dynamics of tail-like current layers caused by anomalous electrical resistivity in a plasma with lower-hybrid-drift (LHD) turbulence is considered. Additionally to the LHD-resistivity, a resistivity pulse in the magnetic neutral sheet is given initiating a magnetic reconnection process. Then the temporal and spatial evolution of the magnetic and electric fields, the plasma convection and the anomalous resistivity are obtained numerically. Taking into account more exact expressions for the LHD-resistivity in the current layer as done in former works, the LHD-turbulence is found to be excited farther from the neutral sheet, and thus, with the time, secondary current sheets are obtained in the plasma-magnetic field system. It is shown that the inductive electric field moving from the magnetic neutral sheet to the current layer periphery during the reconnection process may be considered as indicator of the plasma disturbances.     

Particle acceleration during the explosive phase of a solar flare

Starting with the quasi-linear equation of the distribution function of particles in a regular electric field, a combined diffusion coefficient in the momentum space conbining the effects of the regular field and a turbulent field is obtained and a combined mechanism of acceleration by the regular and turbulent fields in the neutral sheet of solar proton flares is proposed. It is shown by calculation that conditions in solar proton flares are such that the charged particles can be effectively accelerated to tens of MeV, even ~1 GeV. It is shown that the combined acceleration by a regular electric field and ion-acoustic turbulence pumps the protons and other heavy ions into ranges of energy where they can be accelerated by Langmuir turbulence. By considering the combined acceleration by Langmuir turbulence and the regular electric field, the observed spectrum of energetic protons and the power-law spectrum of energetic electrons can be reproduced.     

Laboratory simulation of energy release in solar flares

A laboratory simulation method is proposed for energy release processes occurring in a fragment of the flare current sheet on the Sun. The method relies on the assumption that the spatial scale of such processes is represented by the current sheet's thickness whose values can be close for both the solar and laboratory conditions. It is shown that in an extended current sheet, current dissipation on anomalous resistivity that ensures the specific power of energy release close to that observed in a flare, is the main energy release mechanism. A rapid compression of the sheet by external magnetic fields can provide the condition for switching on a powerful energy release. The tearing instability developing in a homogeneous neutral sheet, leads to the formation of magnetic islands in which the energy release is localized.     

Reconnection driven by an erupting filament in the May 14, 1981 flare

We present observations of the flare of May 14, 1981, which can be classified as a three-ribbon flare. After a detailed analysis in metric, decimetric, microwave, optical, and X-ray ranges we propose that the event was caused by a reconnection process driven by erupting filament. The energy was liberated in the current sheet above the filament in the region between the erupting flux and the overlying field. It is shown that plasma microinstabilities develop as the plasma enters the current sheet. The observations indicate that during the precursor phase a certain low-frequency turbulence, such as ion-accoustic turbulence had to be present.The reconnection rate was growing due to the increasing tension of the stretched overlying field. It is shown that the reconnection proceeded in the Sonnerup-Petschek regime during the precursor, and changed to the pile-up regime in the fast reconnection phase, when the maximal lateral expansion (50 km s^–1) of the H ribbons was observed. The proposed process of reconnection driven by an erupting filament can be applied to three- and four-ribbon flares.     

Energies of Electrons Accelerated in Turbulent Reconnecting Current Sheets in Solar Flares

Yohkoh observations strongly suggest that electron acceleration in solar flares occurs in magnetic reconnection regions in the corona above the soft X-ray flare loops. Unfortunately, models for particle acceleration in reconnecting current sheets predict electron energy gains in terms of the reconnection electric field and the thickness of the sheet, both of which are extremely difficult to measure. It can be shown, however, that application of Ohm's law in a turbulent current sheet, combined with energy and Maxwell's equations, leads to a formula for the electron energy gain in terms of the flare power output, the magnetic field strength, the plasma density and temperature in the sheet, and its area. Typical flare parameters correspond to electron energies between a few tens of keV and a few MeV. The calculation supports the viewpoint that electrons that generate the continuum gamma-ray and hard X-ray emissions in impulsive solar flares are accelerated in a large-scale turbulent current sheet above the soft X-ray flare loops.     

First phase acceleration mechanisms and implications for hard X-ray burst models in solar flares

Requirements for the number of nonthermal electrons which must be accelerated in the impulsive phase of a flare are reviewed. These are uncertain by two orders of magnitude depending on whether hard X-rays above 25 keV are produced primarily by hot thermal electrons which contain a small fraction of the flare energy or by nonthermal streaming electrons which contain > 50% of the flare energy. Possible acceleration mechanisms are considered to see to what extent either X-ray production scenario can be considered viable. Direct electric field acceleration is shown to involve significant heating. In addition, candidate primary energy release mechanisms to convert stored magnetic energy into flare energy, steady reconnection and the tearing mode instability, transfer at least half of the stored energy into heat and most of the remaining energy to ions. Acceleration by electron plasma waves requires that the waves be driven to large amplitude by electrons with large streaming velocities or by anisotropic ion-acoustic waves which also require streaming electrons for their production. These in turn can only come from direct electric field acceleration since it is shown that ion-acoustic waves excited by the primary current cannot amplify electron plasma waves. Thus, wave acceleration is subject to the same limitations as direct electric field acceleration. It is concluded that at most 0.1% of the flare energy can be deposited into nonthermal streaming electrons with the energy conversion mechanisms as they have been proposed and known acceleration mechanisms. Thus, hard X-ray production above 10 keV primarily by hot thermal electrons is the only choice compatible with models for the primary energy release as they presently exist.     

Turbulence in fast-mode shocks as a triggering mechanism in a solar flare

If a solar flare originates from the dissipation of magnetic energy, available in abundance in a larger region, this dissipation must take place very rapidly. A local topological change in the magnetic field structure may be sufficient to start the dissipation process. Such a change in topology might be obtained by fast reconnection in a smaller region, such as e.g. in the Sweet-Parker model, as a result of current-driven microinstabilities.Among the candidates satisfying the requirements to obtain large enough currents, such as magnetically neutral or current sheets and MHD shocks, the latter are shown to be most probable. In a fast MHD shock the (thermal) results of turbulence do in fact destroy the conditions for turbulence. However, in this work we show numerically that the nonlinear steepening mechanism of such a shock is able to restore the driving current for a large range of parameters and over a long time. This is still true if the most difficult threshold for turbulence, being that for Langmuir turbulence, is to be achieved. The critical parameter, not only for the occurrence of turbulence but also for the restoration of the driving current, is the shock thickness.     

2 1/2-Dimensional Reconnection Model and Energy Release in Solar Flares

We employ a 2 1/2-dimensional reconnection model to analyse different aspects of the energy release in two-ribbon flares. In particular, we investigate in which way the systematic change of inflow region variables, associated with the vertical elongation of current sheet, affects the flare evolution. It is assumed that as the transversal magnetic field decreases, the ambient plasma-to-magnetic pressure ratio increases, and the reconnection rate diminishes. As the transversal field decreases due to the arcade stretching, the energy release enhances and the temperature rises. Furthermore, the magnetosonic Mach number of the reconnection outflow increases, providing the formation of fast mode standing shocks above the flare loops and below the erupting flux rope. Eventually, in the limit of a very small transversal field the reconnection becomes turbulent due to a highly non-linear response of the system to small fluctuations of the transversal field. The turbulence results in the energy release fragmentation which increases the release efficiency, and is likely to be responsible for the impulsive phase of the flare. On the other hand, as the current sheet stretches to larger heights, the ambient plasma-to-magnetic pressure ratio increases which causes a gradual decrease of the reconnection rate, energy release rate, and temperature in the late phase of flare. The described magnetohydrodynamical changes affect also the electron distribution function in space and time. At large reconnection rates (impulsive phase of the flare) the ratio of the inflow-to-outflow magnetic field strength is much smaller than at lower reconnection rates (late phase of the flare), i.e., the corresponding loss-cone angle becomes narrower. Consequently, in the impulsive phase a larger fraction of energized electrons can escape from the current sheet downwards to the chromosphere and upwards into the corona – the dominant flare features are the foot-point hard X-ray sources and type III radio bursts. On the other hand, at low reconnection rates, more particles stay trapped in the outflow region, and the thermal conduction flux becomes strongly reduced. As a result, a superhot loop-top, and above-the-loop plasma appears, as sometimes observed, to be a dominant feature of the gradual phase.     

Possible spectral diagnostics for turbulent electric fields in solar flares

A high level of electrostatic turbulence can be generated by plasma instabilities in solar flares. The turbulence, in general, has both quasistatic (ionic) and dynamic (electronic) components. Hydrogenic line profiles develop distinct features under the simultaneous effect of the quasistatic and dynamic turbulent electric fields: these features are discussed and the possible inferences that can be drawn as to the nature of the turbulent electric fields through their observation are pointed out. Solar flare spectral diagnostics along these lines can provide an observational test of the existence of turbulent electric fields in the flare region. When a detailed determination of the characteristics of these fields is feasible, it would, in turn, help in identifying the underlying plasma instability mechanisms that give rise to these fields during the flare build-up and associated processes.     

The structure of high-temperature solar flare plasma in non-thermal flare models

We solve the energy equation for the high-temperature (coronal) component of flare plasma for two models of energy input: (i) direct collisional heating by a beam of suprathermal electrons, and (ii) ohmic heating by the beam-neutralizing reverse current. We discuss the regimes where each case is applicable, and solve for the differential emission measure distribution of the coronal plasma in each case. Scaling laws between loop temperatures and injected electron fluxes are derived for both models; these are testable observationally through coordinated soft X-ray and hard X-ray observations, thus providing a method of discriminating between the two cases. We also readdress the question of the energetic importance of a return current which is below the instability threshold for generation of ion-acoustic plasma turbulence. We find that unless the ambient coronal density is very low ( 10⁹ cm ^–3), collisional heating will always dominate there, in agreement with the findings of previous authors. However, in the chromosphere/corona transition region, the relatively low temperature and correspondingly high plasma resistivity imply that reverse current ohmic heating can predominate the flare energetics, by up to an order of magnitude.Presidential Young Investigator.     

Robust Scalings in Compressible Flux Pile-Up Reconnection

Flux pile-up magnetic reconnection is traditionally considered only for incompressible plasmas. The question addressed in this paper is whether the pile-up scalings with resistivity are robust when plasma compressibility is taken into account. A simple analytical argument makes it possible to understand why the transition from a highly compressible limit to the incompressible one is difficult to discern in typical simulations spanning a few decades in resistivity. From a practical standpoint, however, flux pile-up reconnection in a compressible plasma can lead to anomalous electric resistivity in the current sheet and flare-like energy release of magnetic energy in the solar corona.     

On the possibility of the development of longitudinal wave instabilities on the background of the small-scale Bernstein turbulence in preflare chromosphere of a solar active region

The process of origination and development of instabilities of the longitudinal waves of two types, namely, low-frequency ion-acoustic and high-frequency (“electronic”) Langmuir waves, in the preflare atmosphere of an active solar region are studied. The area under study is located at the chromospheric part of the flare loop near its footpoint. A weak large-scale electric field of flaring loop is the main source of these instabilities. The velocity of an electronic flow in the preflare plasma is supposed to be much lower than thermal electron velocity. Instability development is considered against the background of small-scale Bernstein wave turbulence, which exists in the preflare plasma and has an extremely low threshold of excitation. The necessary conditions for the instability origination and development, as well as the boundary values of the main plasma and wave perturbation parameters, are calculated.     

Magnetic field rearrangement by accelerated particles near a shock front

We solve the self-consistent problem of the generation of a static magnetic field by the electric current of accelerated particles near a strong plane MHD shock front. We take into account the back reaction of the field on the particle diffusion tensors and the background plasma parameters near the front. Various states that differ significantly in static magnetic-field strength are shown to be possible near a strong front. If the initial field has a component normal to the front, then its components parallel to the front are suppressed by accelerated particles by several orders of magnitude. Only the component perpendicular to the front remains. This field configuration for uniform particle injection at the front does not lead to the generation of an additional field, and, in this sense, it is stable. If the initial field is parallel to the front, then either its significant enhancement by two or three orders of magnitude or its suppression by several orders of magnitude is possible. The phenomenon under consideration is an example of the self-organization of plasma with a magnetic field in a strongly nonequilibrium system. It can significantly affect the efficiency of particle acceleration by the shock front and the magnetobremsstrahlung of the accelerated particles.     

Spontaneous magnetic reconnection

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

A model of impulsive loop flares caused by anomalous heat conduction

Numerical simulation is made of the impulsive loop flare caused by transient heat conduction along the loop with an applied axial electric current.It is assumed that a segment near the top of the coronal loop is heated to above 10⁷ K by a heat input that is small compared with the total flare energy, which is given by the magnetic energy of the initial current. Due to the heat conduction, a hump appears in the velocity distribution of electrons, which may excite electron plasma waves with a sufficiently high intensity to cause an anomalous resistivity, as shown theoretically in a previous paper. In that paper, an effect of the plasma waves on the dynamics of electrons was taken into account consistently, but an anomalous heating due to an ohmic dissipation of the initial current under the anomalous resistivity was not taken into account.The aim of the present study is to study the subsequent dynamics of the heated gas caused by the anomalous heating, but in order to avoid an unpractically long computation time, the energy density of the plasma waves is estimated by the energy density of electrons in the velocity hump, without taking into account the effect of the plasma waves consistently in the dynamics of the electrons.The initial current starts to decay gradually by an ohmic dissipation under the anomalous resistivity occurring near the top of the loop to heat this region more. The enhanced heat conduction causes the velocity humps in a wider location. Consequently, the anomalous heating continues and spreads in a self-generating way even after the end of the initial minor heating. Thus the temperature near the loop top becomes above 10⁸ K and the high-temperature region spreads in both directions along the loop with such a high speed as (2–3) × 10⁴ km s^–1, which is nearly equal to the speed of flux-limited heat conduction. On the other hand, induced electric field estimated from the anomalous resistivity is 3.3 × 10⁷ V at the termination of the present simulation, under the modest initial current of 1.5 A m^–2.X-ray emissions expected from the present model loop, show three sources, two footpoints with unequal brightness and a coronal source expanding along the loop in both directions.