Catastrophe of Coronal Magnetic Flux Ropes Caused by Photospheric Motions

Using a 2.5-D, time-dependent ideal MHD model in Cartesian coordinates, we carried out numerical simulations to investigate the equilibrium and evolution properties of a magnetic configuration that consists of a coronal magnetic flux rope and a partly open photospheric background field, which is equivalent to that produced by a two-patch magnetic source on the photospheric surface. The axial and annular magnetic fluxes of the flux rope are given and fixed. The global magnetic configuration evolves in response to three types of changes of the background field: decreasing of the distance between the two sources, shrinking of the size of each source, and increasing of the shear in the closed component of the background field. As a result, the geometrical parameters of the flux rope, i.e. the height of the rope axis, the half-width of the rope and the length of the vertical current sheet below the rope, change due to the variation of the background field. It is shown that for a given coronal magnetic flux rope in a partly open background field, the variation of the geometrical parameters of the flux rope displays a catastrophic behavior, namely, there exists a critical point for each case, at which an infinitesimal change of the background field leads to a loss of equilibrium, and thus a jump of the flux rope. The implication of such a catastrophe in solar active phenomena is briefly discussed.     

Initiation of cmes by magnetic flux emergence

The initiation of solar Coronal Mass Ejections (CMEs) is studied in the framework of numerical magnetohydrodynamics (MHD). The initial CME model includes a magnetic flux rope in spherical, axi-symmetric geometry. The initial configuration consists of a magnetic flux rope embedded in a gravitationally stratified solar atmosphere with a background dipole magnetic field. The flux rope is in equilibrium due to an image current below the photosphere. An emerging flux triggering mechanism is used to make this equilibrium system unstable. When the magnetic flux emerges within the filament below the flux rope, this results in a catastrophic behavior similar to previous models. As a result, the flux rope rises and a current sheet forms below it. It is shown that the magnetic reconnection in the current sheet below the flux rope in combination with the outward curvature forces results in a fast ejection of the flux rope as observed for solar CMEs. We have done a parametric study of the emerging flux rate.     

Predictions of Energy and Helicity in Four Major Eruptive Solar Flares

In order to better understand the solar genesis of interplanetary magnetic clouds (MCs), we model the magnetic and topological properties of four large eruptive solar flares and relate them to observations. We use the three-dimensional Minimum Current Corona model (Longcope, 1996, Solar Phys. 169, 91) and observations of pre-flare photospheric magnetic field and flare ribbons to derive values of reconnected magnetic flux, flare energy, flux rope helicity, and orientation of the flux-rope poloidal field. We compare model predictions of those quantities to flare and MC observations, and within the estimated uncertainties of the methods used find the following: The predicted model reconnection fluxes are equal to or lower than the reconnection fluxes inferred from the observed ribbon motions. Both observed and model reconnection fluxes match the MC poloidal fluxes. The predicted flux-rope helicities match the MC helicities. The predicted free energies lie between the observed energies and the estimated total flare luminosities. The direction of the leading edge of the MC’s poloidal field is aligned with the poloidal field of the flux rope in the AR rather than the global dipole field. These findings compel us to believe that magnetic clouds associated with these four solar flares are formed by low-corona magnetic reconnection during the eruption, rather than eruption of pre-existing structures in the corona or formation in the upper corona with participation of the global magnetic field. We also note that since all four flares occurred in active regions without significant pre-flare flux emergence and cancelation, the energy and helicity that we find are stored by shearing and rotating motions, which are sufficient to account for the observed radiative flare energy and MC helicity.     

On current sheet approximations in models of eruptive flares

We consider an approximation sometimes used for current sheets in flux-rope models of eruptive flares. This approximation is based on a linear expansion of the background field in the vicinity of the current sheet, and it is valid when the length of the current sheet is small compared to the scale length of the coronal magnetic field. However, we find that flux-rope models which use this approximation predict the occurrence of an eruption due to a loss of ideal-MHD equilibrium even when the corresponding exact solution shows that no such eruption occurs. Determination of whether a loss of equilibrium exists can only be obtained by including higher order terms in the expansion of the field or by using the exact solution.     

Coronal mass ejections: Some possibilities for prediction

The onset stage of coronal mass ejections (CMEs) is difficult to observe and is poorly studied. In spite of their practical importance, methods for CME predictions with sufficient lead times are only in the nascent stages of development. The most probable CME mechanism is a catastrophic loss of equilibrium of a large-scale current system in the corona (a flux rope). A twisted magnetic rope is maintained by the tension of field lines of photospheric sources until parameters of the system reach critical values and the equilibrium is lost. Unfortunately, there is low-density plasma (coronal cavity) in most of the rope volume; thus, it is difficult to observe a rope. However, the lower parts of the helical field lines of a rope are fine traps for the dense cold plasma of prominences. Thus, prominences are the best tracers of flux ropes in the corona. The maximal height up to which the rope is in stable equilibrium can be found by analyzing the distribution of the magnetic field generated by photospheric sources in the corona. Comparing this critical height with the actually observed prominence height, one can estimate the probability of the loss of equilibrium by a magnetic rope with a following eruption of prominences and coronal mass ejections.     

Preflare current sheets in the solar atmosphere

Neutral current sheets are expected to form in the solar atmosphere when photospheric motions or the emergence of new magnetic flux causes oppositely directed magnetic fields to be pressed together. Magnetic energy may thus be stored slowly in excess of the minimum energy associated with a purely potential field and released suddenly during a solar flare. For simplicity, we investigate the neutral sheet which forms between two parallel line dipoles when either the distance between them decreases or their dipole moments increase. It is found that, when the dipoles have approached by an amount equal to a tenth of their original separation distance, the stored energy is comparable with that released in a major flare. In addition, a similarity solution for one-dimensional magnetohydro-dynamic flow within such a neutral sheet is presented; it demonstrates that rapid conversion of magnetic energy into heat is possible provided conditions at the edge of the neutral sheet are changing sufficiently quickly.     

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.     

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.     

Does the region of flare-energy release work as a vacuum-cleaner?

We aim to explore the unusual flare event which took place in the solar atmosphere on September 22, 2011 and propose its theoretical interpretation. We analyze the process of energy release in the twisted magnetic flux-rope associated with the event, assuming the excitation of anomalous resistivity of turbulent plasma in the rope, and solve numerically nonlinear two-dimensional (2D) magnetohydrodynamic (MHD) equations. The analytical approach to the problem of flare-energy release show that the conditions of excitation of anomalous resistivity can be satisfied in the twisted magnetic flux-rope whose parameters fits well the SDO observational findings. One of the most remarkable properties of the flare phenomenon under the present consideration was the permanent sucking of the coronal/chromospheric gas from the very remote points to the flare filament, i.e. into the low-lying hot region of the flare energy release. This unusual phenomenon has been simulated by numerical methods in terms of ideal MHD. The numerical results reveal that siphon back-flow exhibits characteristic spatial signatures which mimic the observational findings. The flare-energy release region, as a part of strongly twisted magnetic flux-rope, is able to work as a vacuum-cleaner.     

Numerical Experiments Base on the Catastrophe Model of Solar Eruptionstwo

On the basis of the catastrophe model developed by Isenberg et al., we have used the NIRVANA code to perform the magnetohydrodynamics (MHD) numerical experiments to look into the various behaviors of the coronal magnetic configuration that includes a current-carrying flux rope for modelling the prominence levitation in the corona. These behaviors include the evolution of the equilibrium height of magnetic flux rope with the background magnetic field, the corresponding interior equilibrium of magnetic flux rope, the dynamic properties of magnetic flux rope after the system loses equilibrium, as well as the impact of the reference radius on the equilibrium height of magnetic flux rope. In our calculations, an empirical model of the coronal density distribution given by Sittler & Guhathakurta is used, and the physical dissipation is included. Our experiments show that a deviation between the simulated equilibrium height of magnetic flux rope and the theoretical result of Isenberg et al. exists, but it is not apparent, and the evolutionary features of the two results are similar. If the magnetic flux rope is initially located at the stable branch of the theoretical equilibrium curve, the magnetic flux rope will quickly reach the equilibrium position after several rounds of oscillations as a result of the self-adjustment of the system; when the system is located at the critical point it will quickly lose equilibrium and evolve to the eruptive state; the impact of the variation of reference radius on the equilibrium height of magnetic flux rope is consistent with the prediction of the theory; in the eruptive state, the kinetic properties of magnetic flux rope are consistent with the results given by the Lin-Forbes model and observation, and the fast-mode shock in front of the magnetic flux rope is observed in our experiments; furthermore, because that the dissipation is included in our numerical experiments, the energy conversion from the magnetic energy to other forms of energy is very apparent in the eruptive process.     

Three-Dimensional Evolution of Flux-Rope CMEs and Its Relation to the Local Orientation of the Heliospheric Current Sheet

Flux ropes ejected from the Sun may change their geometrical orientation during their evolution, which directly affects their geoeffectiveness. Therefore, it is crucial to understand how solar flux ropes evolve in the heliosphere to improve our space-weather forecasting tools. We present a follow-up study of the concepts described by Isavnin, Vourlidas, and Kilpua (Solar Phys. 284, 203, 2013). We analyze 14 coronal mass ejections (CMEs), with clear flux-rope signatures, observed during the decay of Solar Cycle 23 and rise of Solar Cycle 24. First, we estimate initial orientations of the flux ropes at the origin using extreme-ultraviolet observations of post-eruption arcades and/or eruptive prominences. Then we reconstruct multi-viewpoint coronagraph observations of the CMEs from ≈?2 to 30 R_⊙ with a three-dimensional geometric representation of a flux rope to determine their geometrical parameters. Finally, we propagate the flux ropes from ≈?30 R_⊙ to 1 AU through MHD-simulated background solar wind while using in-situ measurements at 1 AU of the associated magnetic cloud as a constraint for the propagation technique. This methodology allows us to estimate the flux-rope orientation all the way from the Sun to 1 AU. We find that while the flux-ropes’ deflection occurs predominantly below 30 R_⊙, a significant amount of deflection and rotation happens between 30 R_⊙ and 1 AU. We compare the flux-rope orientation to the local orientation of the heliospheric current sheet (HCS). We find that slow flux ropes tend to align with the streams of slow solar wind in the inner heliosphere. During the solar-cycle minimum the slow solar-wind channel as well as the HCS usually occupy the area in the vicinity of the solar equatorial plane, which in the past led researchers to the hypothesis that flux ropes align with the HCS. Our results show that exceptions from this rule are explained by interaction with the Parker-spiraled background magnetic field, which dominates over the magnetic interaction with the HCS in the inner heliosphere at least during solar-minimum conditions.     

Catastrophe of coronal magnetic rope in partly open multipolar magnetic field

Catastrophe of coronal magnetic rope embedded in a partly open multipolar background magnetic field is studied by using a 2-dimensional, 3-component ideal MHD model in spherical coordinates. The background field is composed of three closed bipolar fields of a coronal streamer and an open field with an equatorial current sheet. The magnetic rope lies below the central bipolar field, and it is characterized by its annular and axial magnetic fluxes. For a given annual flux, there is a critical value of the axial flux, and for a given axial flux, there is a critical value of annual flux such that, below the critical value, the magnetic rope is attached to the solar surface and the system stays in equilibrium, but when the critical value is exceeded, the magnetic rope breaks free and erupts upward. This implies that catastrophe can occur in a coronal magnetic rope embedded in a partly open multipolar background magnetic field. Our computation gives a threshold value of magnetic energy that is about 15% greater than the energy of the partly open magnetic field (the central bipolar field open and the fields on either side closed). The excess energy may serve as source for solar explosions such as coronal mass ejections.     

Properties of Magnetic Tongues over a Solar Cycle

The photospheric spatial distribution of the main magnetic polarities of bipolar active regions (ARs) present during their emergence deformations are known as magnetic tongues. They are attributed to the presence of twist in the toroidal magnetic-flux tubes that form the ARs. The aim of this article is to study the twist of newly emerged ARs from the evolution of magnetic tongues observed in photospheric line-of-sight magnetograms. We apply the procedure described by Poisson et al. (Solar Phys. 290, 727, 2015a) to ARs observed over the full Solar Cycle 23 and the beginning of Cycle 24. Our results show that the hemispherical rule obtained using the tongues as a proxy of the twist has a weak sign dominance (53 % in the southern hemisphere and 58 % in the northern hemisphere). By defining the variation of the tongue angle, we characterize the strength of the magnetic tongues during different phases of the AR emergence. We find that there is a tendency of the tongues to be stronger during the beginning of the emergence and to become weaker as the AR reaches its maximum magnetic flux. We compare this evolution with the emergence of a toroidal flux-rope model with non-uniform twist. The variety of evolution of the tongues in the analyzed ARs can only be reproduced when using a broad range of twist profiles, in particular having a large variety of twist gradients in the direction vertical to the photosphere. Although the analytical model used is a special case, selected to minimize the complexity of the problem, the results obtained set new observational constraints to theoretical models of flux-rope emergence that form bipolar ARs.     

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

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

Helical structures around quiescent solar prominences computed from observable magnetic fields

We analyse the magnetic support of solar prominences in two-dimensional linear force-free fields. A line current is added to model a helical configuration, well suited to trap dense plasma in its bottom part. The prominence is modeled as a vertical mass-loaded current sheet in equilibrium between gravity and magnetic forces.We use a finite difference numerical technique which incorporates both vertical photospheric and horizontal prominence magnetic field measurements. The solution of this mixed boundary problem generally presents singularities at both the bottom and top of the model prominence. The removal of the singularities is achieved by superposition of solutions. Together with the line current equilibrium, these three conditions determine the amplitude of the magnetic field in the prominence, the flux below the prominence and the current intensity, for a given height of the line current. A numerical check of accuracy in the removal of singularities, is done by using known analytical solutions in the potential limit.We have investigated both bipolar and quadrupolar photospheric regions. In this mixed boundary problem the polarity of the field component orthogonal to the prominence is mainly fixed by the imposed height of the line current. For bipolar regions above (respectively below) a critical height the configuration is inverse (respectively normal). For quadrupolar regions the polarity is reversed if we refer the prominence polarity to the closest photospheric polarities. We introduce the polarity of the component parallel to the prominence axis with reference to a sheared arcade. Increasing the shear with fixed boundary conditions can increase or decrease the mass supported depending on the configuration.     

The Importance of Photospheric Intense Flux Tubes for Coronal Heating

Dissipation of magnetic energy in the corona requires the creation of very fine scale-lengths because of the high magnetic Reynolds number of the plasma. The formation of current sheets is a natural possible solution to this problem and it is now known that a magnetic field that is stressed by continous photospheric motions through a series of equilibria can easily form such sheets. Furthermore, in a large class of 3D magnetic fields without null points there are locations, called quasi-separatrix layers (QSLs), where the field-line linkage changes drastically. They are the relevant generalisation of normal separatrices to configurations without nulls: along them concentrated electric currents are formed by smooth boundary motions and 3D magnetic reconnection takes place when the layers are thin enough. With a homogenous normal magnetic field component at the boundaries, the existence of thin enough QSL to dissipate magnetic energy rapidly requires that the field is formed by flux tubes that are twisted by a few turns. However, the photospheric field is not homogeneous but is fragmented into a large number of thin flux tubes. We show that such thin tubes imply the presence of a large number of very thin QSLs in the corona. The main parameter on which their presence depends is the ratio between the magnetic flux located outside the flux tubes to the flux inside. The thickness of the QSLs is approximately given by the distance between neighbouring flux tubes multiplied by the ratio of fluxes to a power between two and three (depending on the density of flux tubes). Because most of the photospheric magnetic flux is confined in thin flux tubes, very thin QSLs are present in the corona with a thickness much smaller than the flux tube size. We suggest that a turbulent resistivity is triggered in a QSL, which then rapidly evolves into a dynamic current sheet that releases energy by fast reconnection at a rate that we estimate to be sufficient to heat the corona. We conclude that the fragmentation of the photospheric magnetic field stimulates the dissipation of magnetic energy in the corona.     

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.     

A Twisted Flux Rope Model for Coronal Mass Ejections and Two-Ribbon Flares

We present a new approach to the theory of large-scale solar eruptive phenomena such as coronal mass ejections and two-ribbon flares, in which twisted flux tubes play a crucial role. We show that it is possible to create a highly nonlinear three-dimensional force-free configuration consisting of a twisted magnetic flux rope representing the magnetic structure of a prominence (surrounded by an overlaying, almost potential, arcade) and exhibiting an S-shaped structure, as observed in soft X-ray sigmoid structures. We also show that this magnetic configuration cannot stay in equilibrium and that a considerable amount of magnetic energy is released during its disruption. Unlike most previous models, the amount of magnetic energy stored in the configuration prior to its disruption is so large that it may become comparable to the energy of the open field.     

Catastrophe of coronal magnetic flux ropes in partially open magnetic fields

Using a 2.5-D, time-dependent ideal MHD model in Cartesian coordinates, a numerical study is carried out to find equilibrium solutions associated with a magnetic flux rope in the corona. The ambient magnetic field is partially open, consisting of a closed arcade in the center and an open field at the flank. The coronal magnetic flux rope is characterized by its magnetic properties, including the axial and annular magnetic fluxes and the magnetic helicity, and its geometrical features, including the height of the rope axis, the halfwidth of the rope and the length of the vertical current sheet below the rope. It is shown that for a given partially open ambient magnetic field, the dependence of the geometrical features on the magnetic properties displays a catastrophic behavior, namely, there exists a certain critical point, across which an infinitesimal enhancement of the magnetic parameters causes a finite jump of the geometrical parameters for the rope. The amplitude of the jump depends on the extent to which the ambient magnetic field in open, and approaches to zero when the ambient magnetic field becomes completely closed. The implication of such a catastrophe in solar active phenomena is briefly discussed.     

Coronal Magnetic Flux Rope Equilibria and Magnetic Helicity

1 INTRODUCTIONObservations show that the magnetic helicity of solar magnetic structures has a predominantsign in each hemisphere of the Sun, positive in the southern hemisphere and negative in thenorthern, regardless of the solar cycle (Rust, 1994). The magnetic helicity is strictly conservedin the frame of ideal MHD (WOltjer, 1958), and approximately conserved in the presence ofresistive dissipation and magnetic reconnection in a highly conductive plajsma (Taylor, 1974;Berger, 1984; H…