On the maximum energy release in flux-rope models of Eruptive Flares

We determine the photospheric boundary conditions which maximize the magnetic energy released by a loss of ideal-MHD equilibrium in two-dimensional flux-rope models. In these models a loss of equilibrium causes a transition of the flux rope to a lower magnetic energy state at a higher altitude. During the transition a vertical current sheet forms below the flux rope, and reconnection in this current sheet releases additional energy. Here we compute how much energy is released by the loss of equilibrium relative to the total energy release. When the flux-rope radius is small compared to its height, it is possible to obtain general solutions of the Grad-Shafranov equation for a wide range of boundary conditions. Variational principles can then be used to find the particular boundary condition which maximizes the magnetic energy released for a given class of conditions. We apply this procedure to a class of models known as cusp-type catastrophes, and we find that the maximum energy released by the loss of equilibrium is 20.8% of the total energy release for any model in this class. If the additional restriction is imposed that the photospheric magnetic field forms a simple arcade in the absence of coronal currents, then the maximum energy release reduces to 8.6%.     

Magnetic field evolution during prominence eruptions and two-ribbon flares

Simple models for the MHD eruption of a solar prominence are presented, in which the prominence is treated as a twisted magnetic flux tube that is being repelled from the solar surface by magnetic pressure forces. The effects of different physical assumptions to deal with this magneto-hydrodynamically complex phenomenon are evaluated, such as holding constant the prominence current, radius, flux or twist or modelling the prominence as a current sheet. Including a background magnetic field allows the prominence to be in equilibrium initially with an Inverse Polarity and then to erupt due to magnetic non-equilibrium when the background magnetic field is too small or the prominence twist is too great. The electric field at the neutral point below the prominence rapidly increases to a maximum value and then declines. Including the effect of gravity also allows an equilibrium with Normal Polarity to exist. Finally, an ideal MHD solution is found which incorporates self-consistently a current sheet below the prominence and which implies that a prominence will still erupt and form a current sheet even if no reconnection occurs. When reconnection is allowed it is, therefore, driven by the eruption.     

Coronal Mass Ejection of 15 May 2001: II. Coupling of the Cme Acceleration and the Flare Energy Release

We analyze the relationship between the dynamics of the coronal mass ejection (CME) of 15 May 2001 and the energy release in the associated flare. The flare took place behind the east limb and was disclosed by a growing system of hot soft X-ray (SXR) loops that appeared from behind the limb around the onset of the rapid acceleration of the CME. The highly correlated behavior of the SXR light-curve derivative and the time profile of the CME acceleration reveals an intrinsic relationship between the CME dynamics and the flare energy release. Furthermore, we found that the CME acceleration peak occurs simultaneously with the fastest growth (100 km s^-1) of X-ray loops, indicating that the reconnection plays an essential role in the eruption. Inspecting the CME/flare morphology we recognized in the Yohkoh-SXT images an oval feature that formed within the rising structure at the onset of the rapid acceleration phase, simultaneously with the appearance of the X-ray loops. The eruptive prominence was imbedded within the lower half of the oval, suggestive of a flux-rope/prominence magnetic configuration. We interpret the observed morphological evolution in terms of a reconnection process in the current sheet that presumably formed below the erupting flux-rope at the onset of the CME acceleration. Measurements of the tip-height of the cusped X-ray loop system and the height of the lower edge of the oval, enable us to trace the stretching of the current sheet. The initial distance between the oval and the loops amounted to 35 – 40 Mm. In about 1 h the inferred length of the current sheet increased to 150 – 200 Mm, which corresponds to a mean elongation speed of 35 – 45 km s^-1. The results are discussed in the framework of CME models that include the magnetic reconnection below the erupting flux-rope.     

Initiation of Coronal Mass Ejections by Sunspot Rotation

We study a filament eruption, two-ribbon flare, and coronal mass ejection (CME) that occurred in NOAA Active Region 10898 on 6 July 2006. The filament was located South of a strong sunspot that dominated the region. In the evolution leading up to the eruption, and for some time after it, a counter-clockwise rotation of the sunspot of about 30 degrees was observed. We suggest that the rotation triggered the eruption by progressively expanding the magnetic field above the filament. To test this scenario, we study the effect of twisting the initially potential field overlying a pre-existing flux-rope, using three-dimensional zero-β MHD simulations. We first consider a relatively simple and symmetric system, and then study a more complex and asymmetric magnetic configuration, whose photospheric-flux distribution and coronal structure are guided by the observations and a potential field extrapolation. In both cases, we find that the twisting leads to the expansion of the overlying field. As a consequence of the progressively reduced magnetic tension, the flux-rope quasi-statically adapts to the changed environmental field, rising slowly. Once the tension is sufficiently reduced, a distinct second phase of evolution occurs where the flux-rope enters an unstable regime characterised by a strong acceleration. Our simulations thus suggest a new mechanism for the triggering of eruptions in the vicinity of rotating sunspots.     

Generalized analytical models of Syrovatskii’s current sheet

Two-dimensional stationary magnetic reconnection models that include a thin Syrovatskii-type current sheet and four discontinuous magnetohydrodynamic flows of finite length attached to its endpoints are considered. The flow pattern is not specified but is determined from a self-consistent solution of the problem in the approximation of a strong magnetic field. Generalized analytical solutions that take into account the possibility of a current sheet discontinuity in the region of anomalous plasma resistivity have been found. The global structure of the magnetic field in the reconnection region and its local properties near the current sheet and attached discontinuities are studied. In the reconnection regime in which reverse currents are present in the current sheet, the attached discontinuities are trans-Alfvénic shock waves near the current sheet endpoints. Two types of transitions from nonevolutionary shocks to evolutionary ones along discontinuous flows are shown to be possible, depending on the geometrical model parameters. The relationship between the results obtained and numerical magnetic reconnection experiments is discussed.     

Understanding the Role of Mass-Unloading in a Filament Eruption

We describe a partial filament eruption on 11 December 2011 that demonstrates that the inclusion of mass is an important next step for understanding solar eruptions. Observations from the Solar Terrestrial Relations Observatory-Behind (STEREO-B) and the Solar Dynamics Observatory (SDO) spacecraft were used to remove line-of-sight projection effects in filament motion and correlate the effect of plasma dynamics with the evolution of the filament height. Flux cancellation and nearby flux emergence are shown to have played a role in increasing the height of the filament prior to eruption. The two viewpoints allow the quantitative estimation of a large mass-unloading, the subsequent radial expansion, and the eruption of the filament to be investigated. A 1.8 to 4.1 lower-limit ratio between gravitational and magnetic-tension forces was found. We therefore conclude that following the loss-of-equilibrium of the flux-rope, the radial expansion of the flux-rope was restrained by the filamentary material until 70% of the mass had evacuated the structure through mass-unloading.     

The formation of solar prominences by thermal instability in a current sheet

The energy balance equation for the upper chromosphere or lower corona contains a radiative loss term which is destabilizing, because a slight decrease in temperature from the equilibrium value causes more radiation and hence a cooling of the plasma; also a slight increase in temperature has the effect of heating the plasma. In spite of this tendency towards thermal instability, most of the solar atmosphere is remarkably stable, since thermal conduction is very efficient at equalizing any temperature irregularity which may arise. However, the effectiveness of thermal conduction in transporting heat is decreased considerably in a current sheet or a magnetic flux tube, since heat can be conducted quickly only along the magnetic field lines. This paper presents a simple model for the thermal equilibrium and stability of a current sheet. It is found that, when its length exceeds a certain maximum value, no equilibrium is possible and the plasma in the sheet cools. The results may be relevant for the formation of a quiescent prominence.     

灾变理论在太阳物理和天体物理其他研究领域中的应用

叙述和介绍了太阳爆发的磁通量绳灾变理论和模型的发展过程,强调了建立这样的模型所需要的观测基础。讨论了由模型所预言的爆发磁结构的几个重要特征以及观测结果对这种预言的证实。在此模型的基础上,讨论了一个典型的爆发过程中所出现的不同现象及它们之间的相互关系。最后,介绍了作者的一项最新尝试:将太阳爆发的灾变理论和模型应用到对黑洞吸积盘间歇性喷流的理论研究当中,以及研究所取得的初步结果。     

Catastrophic Behavior of Multiple Coronal Flux Rope System

A major solar active event called Bastille Day Event occurred in AR 9077 on July 14, 2000. Simultaneous occurrence of a filament eruption, a flare and a coronal mass ejection was observed in this event. Previous analyses of this event show that before the event, there existed an activation and eruption of a huge trans-equatorial filament, which might play a crucial role in triggering the Bastille Day event. This implies that independent flux systems are closely related to and affect each other, which has encouraged us to investigate the catastrophic behavior of a multiple coronal flux rope system with the use of a 2.5-D time-dependent MHD model. A force-free field that contains three separate coronal flux ropes is taken to be the initial state. Starting from this state, we increase either the annular or the axial flux of a certain flux rope to examine the catastrophic behavior of the system in two regimes, the ideal MHD regime and the resistive MHD regime. It is found that a catastrophe occurs if the flux exceeds a certain critical value, or the magnetic energy of the system exceeds a certain threshold: the rope of interest breaks away from the base and escapes to infinity, leaving a current sheet below. Moreover, the destiny of the remainder flux ropes relies on whether reconnection takes place across the current sheet. In the ideal MHD regime, i.e., in the absence of reconnection, these ropes remain to be attached to the base in equilibrium, whereas in the resistive MHD regime they abruptly erupt upward during reconnection and escape to infinity. Reconnection causes the field lines to close back to the base and thus changes the background field outside the attached flux ropes in such a way that the constraint on these ropes is substantially relaxed and the corresponding catastrophic energy threshold is reduced accordingly, leading to a catastrophic eruption of these ropes. Since magnetic reconnection is generally inevitable when a current sheet forms and develops through an eruption of one flux rope, the eruption of this flux rope must lead to an eruption of the others. This provides an example to demonstrate the interaction between several independent magnetic flux systems in different regions, as implied by the Bastille Day event, and may serve as a possible mechanism for sympathetic events occurring on the Sun.     

Evidence of Magnetic Flux Ropes in the Solar Wind From Sigmoidal and non-Sigmoidal Active Regions

We have examined WIND magnetic field and plasma data during the first half of 1998 in order to find encounters of this spacecraft with magnetic clouds. From the events obtained through this search, we have selected four of them taking into account their solar origin. The four magnetic clouds are related to halo or partial halo CMEs, but the morphology of the active region before the eruption is sigmoidal for three of them and non-sigmoidal for the other one. We have analyzed these events in the solar wind by fitting the experimental data to a non-force-free flux-rope model. We conclude that both kinds of active regions develop in the solar wind an ejection with a flux-rope topology.     

Eruptive instability of cylindrical prominences

The stability of prominences and the dynamics of an eruption are studied. The prominence is represented by an uniformly twisted, curved, magnetic tube, anchored at both ends in the photosphere. Several stages of the eruption are analyzed, from the pre-eruptive phase and the onset of the instability, up to the late phases of the process. Before the eruption, the prominence evolves through a series of equilibrium states, slowly ascending either due to an increase of the electric current or to mass loss. The eruption starts when the ratio of the current to the total mass attains a critical value after which no neighbouring equilibrium exists. The linearized equation of motion was used to obtain the instability threshold, which is presented in a form enabling comparison with the observations. The height at which the prominence erupts depends on the twist, and is typically comparable with the footpoint half-separation. Low-lying prominences are stable even for large twists. The importance of the external field reconnection below the filament, and the mass loss through the legs in the early phases of the eruption is stressed. The oscillations of stable prominences with periods on the Alfvén time-scale are discussed. The results are compared with the observations.     

Non-linear tearing structures in equilibrium current sheet

A family of exact analytic solutions of the time-independent Vlasov-Maxwell equations is presented. The solutions describe two-dimensional equilibrium current sheet with magnetic field structures resembling that produced by the tearing instability. In particular, the solutions presented here do not restrict the field in the magnetic island to small magnitude. It is shown that as the scale length of the magnetic island increases, the thickness of the current sheet increases while the average current and the average magnetic energy decrease. The tearing structures described by the solutions may exist in the magnetotail current sheet, the magnetopause current layer and the field-aligned auroral sheet current.     

Electric Current Equilibrium in the Corona

A hyperbolic flux-tube configuration containing a null point below the flux rope is considered as a pre-eruptive state of coronal mass ejections that start simultaneously with flares. We demonstrate that this configuration is unstable and cannot exist for a long time in the solar corona. The inference follows from general equilibrium conditions and from analyzing simple models of the flux-rope equilibrium. A direct consequence of the stable flux-rope equilibrium in the corona are separatrices in the horizontal-field distribution in the chromosphere. They can be recognized as specific “herring-bone structures” in a chromospheric fibril pattern.     

Analytical description of charged particle motion in a reconnecting current sheet

We have obtained an analytical solution to the equation of motion in the guiding center approximation for nonrelativistic charged particles in a reconnecting current sheet with a three-component magnetic field. Given the electric field attributable to magnetic reconnection, the solution describes stable and unstable three-dimensional particle orbits. We have found the domain of input parameters at which the motion is stable. A physical interpretation of the processes affecting the stability of the motion is given. Charge separation is shown to take place in the sheet during the motion: oppositely charged particles are localized mostly in different regions of the current sheet. A formula is derived for the particle energy in stable and unstable orbits. The results obtained by numerical and analytical methods are compared.     

A 2-D model for the support of a polar-crown solar prominence

We present a 2-D potential-field model for the magnetic structure in the environment of a typical quiescent polar-crown prominence. The field is computed using the general method of Titov (1992) in which a curved current sheet, representing the prominence, is supported in equilibrium by upwardly directed Lorentz forces to balance the prominence weight. The mass density of the prominence sheet is computed in this solution using a simple force balance and observed values of the photospheric and prominence magnetic field. This calculation gives a mass density of the correct order of magnitude. The prominence sheet is surrounded by an inverse-polarity field configuration adjacent to a region of vertical, open polar field in agreement with observations.A perturbation analysis provides a method for studying the evolution of the current sheet as the parameters of the system are varied together with an examination of the splitting of an X-type neutral point into a current sheet.Program Systems Institute of the Russian Academy of Sciences, Pereslavl-Zalessky 152140, Russia.     

On the Flux-Rope Topology of Ejecta Observed in the Period 1997 – 2006

In the following study our aim is to analyse the magnetic flux-rope topology of some events observed in the interplanetary medium related to ejecta. The magnetic field structures associated with interplanetary coronal mass ejections are globally classified in magnetic clouds and ejecta. One of the main questions regarding these phenomena concerns their flux-rope or non-flux-rope magnetic field line configuration. From the experimental measurements the only way to elucidate such a question is analysing the corresponding data by means of a flux-rope physical model. After selecting the ejecta events observed during the period 1997?–?2006, we have analysed them in light of an analytical model with that topology for the magnetic field components, initially developed for magnetic clouds, and with a non-force-free character; then, incorporating the expansion of the magnetic structure during their evolution in the interplanetary medium. Different parameters obtained from the fitting of the model are related to the orientation of the axis of the magnetic flux-rope structure and, additionally, the closest distance approach of the spacecraft to its axis. One of the main conclusions achieved concerns the fact that the axes of most of those structures are close to the Sun–Earth line, which implies that the passage of the spacecraft through the corresponding ejecta event is by its flank. In general, we show a rough procedure for the analysis and classification of ejecta in terms of their magnetic field topology.     

On the stationary configuration of the heliospheric sheet

The expansion of solar coronal plasma is considered for the model described in Koutchmy et al. (1999). In addition to a spherical solar surface, the initial configuration represents a heliospheric sheet of dense plasma in the dipole equatorial plane. The heliospheric-sheet current decreases with distance as 1/r ², with its sign being opposite to the sign of the initial-dipole current. The latter follows from the fact that the plasma sheet is denser than the surrounding corona and that the equilibrium condition for the sheet in the gravitational and magnetic fields is satisfied. The field lines of this configuration are nearly straight. We have obtained a general solution of the steady-state MHD equations, which depends not only on distance r but also on latitude θ. Applicability of the solution to interpreting observational data, in particular, those obtained from the Ulysses spacecraft, is discussed.     

A circuit model for filament eruptions and two-ribbon flares

We derive a circuit model for solar filament eruptions and two-ribbon flares which reproduces the slow energy build up and eruption of the filament, and the energy dissipation in a current sheet at the top of post-flare loops during the two-ribbon flare. In our model the free magnetic energy is concentrated in a current through the filament, another current through an underlying current sheet, and surface return currents. The magnetic field configuration, generated by these currents and a general photospheric background field, has a topology similar to the field topology derived from observations.We consider two circuits, that of the filament and its return current, and that of the current sheet and its return current. These circuits are inductively coupled and free energy stored in the filament in the pre-flare phase is found to be transferred to the sheet during the impulsive phase, and rapidly dissipated there. A comparable amount of magnetic energy is converted into kinetic energy of the ejected filament. The basic equations of the model are the momentum equations for the filament and the current sheet, and the induction equations for the filament and sheet circuits. The derivation of the equations is an extension of previous models by Kuperus and Raadu, Van Tend and Kuperus, Syrovatskii, and Kaastra. The set of equations is closed in the sense that only the initial conditions and a number of parameters, all related to pre-flare observables, are needed to calculate the evolution of the system. The pre-flare observations we need to determine these parameters, are: (1) a magnetogram, (2) an picture, (3) a measurement of the coronal density in the region, and (4) estimates of the photospheric velocity fields in the region.In the solutions for the evolution of the filament current sheet system we distinghuish 4 phases: (1) a slow energy build up, lasting for almost two days, during which the filament evolves quasi-statically, (2) a metastable state, lasting for about three hours, during which the filament is susceptible to flare triggers, and during which a current sheet emerges, (3) the eruptive phase, with strong acceleration of the filament, during which a large current is induced and dissipated in the current sheet, and energy is injected in the post-flare loops, and finally (4) a post-flare phase, in which the filament acceleration declines and the current sheet vanishes.From further numerical work we derive the following conclusions: (1) The magnetic flux input into the filament circuit has to surpass a certain threshold for an eruption to occur. Below that threshold we find solutions representing quiescent filaments. (2)Flare triggers are neither necessary nor sufficient for an eruption, but may set off the eruption during the metastable state. (3) The model reproduces the increase in shear in the filament prior to the eruption, through adecline of the filament current, in contrast to most models for filament eruptions. (4) The ratio of energy lost as kinetic energy of ejecta to the energy radiated away in the post-flare loops is sensitively dependent on the resistance of the current sheet. (5) Flare prediction is possible with this model, but the potential for triggering during the metastable state complicates the prediction of the exact moment of eruption.Former NAS/NRC Resident Research Associate.ST Systems Corporation.     

Breakout and Tether-Cutting Eruption Models Are Both Catastrophic (Sometimes)

We present a simplified analytic model of a quadrupolar magnetic field and flux rope to model coronal mass ejections. The model magnetic field is two-dimensional, force-free and has current only on the axis of the flux rope and within two current sheets. It is a generalization of previous models containing a single current sheet anchored to a bipolar flux distribution. Our new model can undergo quasi-static evolution either due to changes at the boundary or due to magnetic reconnection at either current sheet. We find that all three kinds of evolution can lead to a catastrophe, known as loss of equilibrium. Some equilibria can be driven to catastrophic instability either through reconnection at the lower current sheet, known as tether cutting, or through reconnection at the upper current sheet, known as breakout. Other equilibria can be destabilized through only one and not the other. Still others undergo no instability, but they evolve increasingly rapidly in response to slow steady driving (ideal or reconnective). One key feature of every case is a response to reconnection different from that found in simpler systems. In our two-current-sheet model a reconnection electric field in one current sheet causes the current in that sheet to increase rather than decrease. This suggests the possibility for the microscopic reconnection mechanism to run away.     

Dispersion relations in the electrostatic approximation for waves in a magnetic neutral sheet

Dispersion equations in the electrostatic approximation are derived for waves propagating near the centre of a magnetic neutral sheet system. The unperturbed equilibrium is based on the Alfvén-Cowley neutral sheet model, in which the sheet current is carried by accelerated non-adiabatic electrons oscillating about the field reversal, and moving through a cold neutralizing ion background. Detailed account is taken of the non-adiabatic nature of the electron motion. It is also recognized that the zeroth order electron distribution may differ significantly from a convecting isotropic Maxwellian. A companion paper presents a detailed numerical study of the dispersion relations derived here.