The magnetic non-equilibrium of buoyant flux tubes in the solar corona

The equilibrium shape of a slender flux tube in the stratified solar atmosphere is studied. The path is determined by a balance between the downwards magnetic tension, which depends on the curvature of the loop, and the upwards buoyancy force. Previous results for untwisted slender tubes are extended to include twisted tubes embedded in an external magnetic field.The path of an untwisted tube in an atmosphere with an ambient magnetic field is calculated. For a given footpoint separation, the height of the tube is lowered by increasing the strength of the external magnetic field. If the footpoints are slowly moved apart, the tube rises, until a threshold separation is reached beyond which there is no possible equilibrium height. This threshold width does not depend on the strength of the external field.The effects of twisting up a curved loop are studied, using an extension of results obtained for slender curved tubes with a straight axis. It is shown that for a twisted tube of given width, there can be two possible values of the equilibrium height. If, however, the tube is twisted more than a certain amount or if the footpoints are too widely separated there is no equilibrium. The critical footpoint separation for non-equilibrium is smaller for a twisted tube that an untwisted one.Twisting a tube or moving its feet apart is thus likely to result in non-equilibrium, causing the tube to rise indefinitely under the influence of the unbalanced buoyant force. It is suggested that this phenomenon could be important in the preflare stage of a large two-ribbon solar flare, by causing the initial slow rise of an active region filament. As well as being involved in the onset of an erupting prominence, this non-equilibrium may also be relevant to the formation of coronal loop transients.     

Numerical simulation of reconnection in an emerging magnetic flux region

The resistive MHD equations are numerically solved in two dimensions for an initial-boundary-value problem which simulates reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. The emerging region is modelled by a cylindrical flux tube with a poloidal magnetic field lying in the same plane as the external, coronal field. The plasma betas of the emerging and coronal regions are 1.0 and 0.1, respectively, and the magnetic Reynolds number for the system is 2 × 10³. At the beginning of the simulation the tube starts to emerge through the base of the rectangular computational domain, and, when the tube is halfway into the computational domain, its position is held fixed so that no more flux of plasma enters through the base. Because the time-scale of the emergence is slower than the Alfvén time-scale, but faster than the reconnection time-scale, a region of closed loops forms at the base. These loops are gradually opened and reconnected with the overlying, external magnetic field as time proceeds.The evolution of the plasma can be divided into four phases as follows: First, an initial, quasi-steady phase during which most of the emergence is completed. During this phase, reconnection initially occurs at the slow rate predicted by the Sweet model of diffusive reconnection, but increases steadily until the fast rate predicted by the Petschek model of slow-shock reconnection is approached. Second, an impulsive phase with large-scale, super-magnetosonic flows. This phase appears to be triggered when the internal mechanical equilibrium inside the emerging flux tube is upset by reconnection acting on the outer layers of the flux tube. During the impulsive phase most of the flux tube pinches off from the base to form a cylindrical magnetic island, and temporarily the reconnection rate exceeds the steady-state Petschek rate. (At the time of the peak reconnection rate, the diffusion region at the X-line is not fully resolved, and so this may be a numerical artifact.) Third, a second quasi-steady phase during which the magnetic island created in the impulsive phase is slowly dissipated by continuing, but low-level, reconnection. And fourth, a static, non-evolving phase containing a potential, current-free field and virtually no flow.During the short time in the impulsive phase when the reconnection rate exceeds the steady-state Petschek rate, a pile-up of magnetic flux at the neutral line occurs. At the same time the existing Petschek-slow-mode shocks are shed and replaced by new ones; and, for a while, both new and old sets of slow shocks coexist.     

Evolution of emerging force-free flux tube to be a flare loop

Evolution of a filamentary magnetic flux tube emerging from the photosphere is investigated in the assumption that the magnetic field is force-free and unchanged during the evolution.If a characteristic radius of the flux tube is 3 km or less setting the field to 1000G, the temperature increases at first due to Joule heating up to about one million degree keeping the plasma density almost constant, and then the density decreases down to a critical value at which a current instability may occur. Thus, a.strong field-aligned electric field of 200 million volts or more is expected to be produced during the following anomalous Ohmic decay of the magnetic field as already shown by a numerical simulation.     

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.     

The eruption of a prominence and coronal mass ejection which drive reconnection

Two possible limiting scenarios are proposed for the production of a coronal mass ejection. In the first the magnetic field around a prominence evolves until it loses equilibrium and erupts, which drives reconnection below the prominence and an eruption of the overlying magnetic arcade. In the second a large-scale magnetic arcade evolves until it loses equilibrium and erupts, thereby causing a prominence to erupt. In general it is likely to be the non-equilibrium of the coupled system which creates the eruption. Furthermore, large quiescent prominences are expected to be centred within the magnetic bubble of a coronal mass ejection whereas when active-region prominences erupt they are likely to be located initially to one side of the bubble.A model is set up for the eruption of a magnetically coupled prominence and coronal mass ejection. This represents a development of the Anzer and Pneuman (1982) model by overcoming two limitations of it, namely that: it is not globally stable initially and so one wonders how it can be set up in a stable way before the eruption; it has reconnection driving the CME whereas recent observations suggest that the reverse may be happening. In our model we assume that magnetic reconnection below the prominence is driven by the eruption and the driver is magnetic non-equilibrium in the coupled prominence-mass ejection system. The prominence is modelled as a twisted flux tube and the mass ejection as an overlying void and magnetic bubble. Two different models of the prominence are considered. In one a globally stable equilibrium becomes unstable when a threshold magnetic flux below the prominence is exceeded and, in the other, equilibrium ceases to exist. In both cases, the prominence and mass-ejection accelerate upwards before reaching constant velocities in a manner that is consistent with observations. It is found that the greater the reconnection that is driven by the eruption, the higher is the final speed.     

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.     

Magnetic Energy of Force-Free Fields with Detached Field Lines

Using an axisymmetrical ideal MHD model in spherical coordinates, we present a numerical study of magnetic configurations characterized by a levitating flux rope embedded in a bipolar background field whose normal field at the solar surface is the same or very close to that of a central dipole. The characteristic plasma β (the ratio between gas pressure and magnetic pressure) is taken to be sosmall (β= 10^-4) that the magnetic field is close to being force-free. The system as a whole is then let evolve quasi-statically with a slow increase of either the annular magnetic flux or the axial magnetic flux of the rope, and the total magneticenergy of the system grows accordingly. It is found that there exists an energy threshold: the flux rope sticks to the solar surface in equilibrium if the magneticenergy of the system is below the threshold, whereas it loses equilibrium if the threshold is exceeded. The energy threshold is found to be larger than that of thecorresponding fully-open magnetic field by a factor of nearly 1.08 irrespective as towhether the background field is completely closed or partly open, or whether the magnetic energy is enhanced by an increase of annular or axial flux of the rope.This gives an example showing that a force-free magnetic field may have an energy larger than the corresponding open field energy if part of the field lines is allowed tobe detached from the solar surface. The implication of such a conclusion in coronal mass ejections is briefly discussed and some comments are made on the maximum energy of force-free magnetic fields.     

Non-ideal MHD Properties of Magnetic Flux Tubes in the Solar Photosphere

Magnetic flux tubes reaching from the solar convectivezone into the chromosphere have to pass through the relatively cool, and therefore non-ideal (i.e. resistive) photospheric region enclosed between the highly ideal sub-photospheric and chromospheric plasma. It is shown that stationary MHD equilibria of magnetic flux tubes which pass through this region require an inflow of photospheric material into the flux tube and a deviation from iso-rotation along the tube axis. This means that there is a difference in angular velocity of the plasma flow inside the tube below and above the non-ideal region. Both effects increase with decreasing cross section of the tube. Although for characteristic parameters of thick flux tubes the effect is negligible, a scaling law indicates its importance for small-scale structures. The relevance of this `inflow effect' for the expansion of flux tubes above the photosphere is discussed.     

Twist transport in strongly torsioned astrophysical flux tubes

Riemannian geometrical effects on the expansion of the electron magnetohydrodynamical (EMH) superconductivity modeled twisted nonplanar thin magnetic flux tubes are considered. A solution is found which represents almost incompressible plasma flows, where the twist of flux tube is computed in terms of the continuous variation of its cross-section. It is shown that the twist increases in regions where twisted flux tube expands as in Parker’s conjecture. From computation of compression along the tube we show that when the torsion is weak a centrifugal or vorticity effect on the longitudinal direction of the tube enhances the screening effect on the “superconductor”. Throughout the paper we consider helical flux tubes where torsion and curvature of the tube are constants. Thus we show that the Parker’s conjecture is valid in a continuos manner for these type II superconducting twisted flux tubes. Throughout the paper we adopt the approximation that the radial component of the magnetic field varies so slowly along the tube axis that it can be approximated to zero along the tube. It is suggested that the models discussed here may also be applied to DNA and nanotubes.     

Model development of supersonic trough wind with shocks

The time dependent one dimensional hydrodynamic equations describe the evolution of the thermal plasma flow along closed magnetic field lines outside of the plasmasphere. The convection of the supersonic polar wind onto a closed field line results in the assumed formation of collisionless plasma shocks. These shocks move earthward as the field line with its ‘frozen-in’ plasma remains fixed or contracts with time to smaller L coordinates. The high equatorial plasma temperature (of the order of electron volts) produced by the shock process decreases with time if the flow is isothermal but it will increase if the contraction is under adiabatic conditions. Assuming adiabaticity a peak in the temperature forms at the equator in conjunction with a depression in the ion density. After an initial contraction, if the flux tube drifts to higher L coordinates the direction of the shock motion can be reversed so that the supersonic region will expand along the field line towards the state characterizing the supersonic polar wind. A rapid expansion will lower the equatorial density while the temperature decreases with time under adiabatic but not isothermal conditions.     

Thermal evolution of current sheets and flash phase of solar flares

The physical conditions in a stationary flow of the Petchek type, allowing reconnection between flux emerging from below the solar photosphere and a preexisting magnetic field, are discussed. It is shown that, when rising in the solar atmosphere, the reconnection region has at first a rather low temperature as compared with its environment. Above a certain critical height, however, this low temperature thermal equilibrium often ceases to be possible, and the sheet rapidly heats, seeking a new thermal equilibrium. During this dynamical process, current-driven microinstabilities may be triggered in the current sheet, giving rise to an enhanced resistivity. High energy particles might be produced by the induced electric field developed during the rapid readjustment of MHD flows that results from this change in the transport properties of the plasma.     

Sausage MHD Waves in Incompressible Flux Tubes with Twisted Magnetic Fields

Twisted magnetic flux tubes are of considerable interest because of their natural occurrence from the Sun’s interior, throughout the solar atmosphere and interplanetary space up to a wide range of applicabilities to astrophysical plasmas. The aim of the present work is to obtain analytically a dispersion equation of linear wave propagation in twisted incompressible cylindrical magnetic waveguides and find appropriate solutions for surface, body and pseudobody sausage modes (i.e. m = 0) of a twisted magnetic flux tube embedded in an incompressible but also magnetically twisted plasma. Asymptotic solutions are derived in long- and short-wavelength approximations. General solutions of the dispersion equation for intermediate wavelengths are obtained numerically. We found, that in case of a constant, but non-zero azimuthal component of the equilibrium magnetic field outside the flux tube the index ν of Bessel functions in the dispersion relation is not integer any more in general. This gives rise to a rich mode-structure of degenerated magneto-acoustic waves in solar flux tubes. In a particular case of a uniform magnetic twist the total pressure is found to be constant across the boundary of the flux tube. Finally, the effect of magnetic twist on oscillation periods is estimated under solar atmospheric conditions. It was found that a magnetic twist will increase, in general, the periods of waves approximately by a few percent when compared to their untwisted counterparts.     

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.     

A New Model for Propagating Parts of EIT Waves: A Current Shell in a CME

EIT waves are observed in EUV as bright fronts. Some of these bright fronts propagate across the solar disk. EIT waves are all associated with a flare and a CME and are commonly interpreted as fast-mode magnetosonic waves. Propagating EIT waves could also be the direct signature of the gradual opening of magnetic field lines during a CME. We quantitatively addressed this alternative interpretation. Using two independent 3D MHD codes, we performed nondimensional numerical simulations of a slowly rotating magnetic bipole, which progressively result in the formation of a twisted magnetic flux tube and its fast expansion, as during a CME. We analyse the origins, the development, and the observability in EUV of the narrow electric currents sheets that appear in the simulations. Both codes give similar results, which we confront with two well-known SOHO/EIT observations of propagating EIT waves (7 April and 12 May 1997), by scaling the vertical magnetic field components of the simulated bipole to the line of sight magnetic field observed by SOHO/MDI and the sign of helicity to the orientation of the soft X-ray sigmoids observed by Yohkoh/SXT. A large-scale and narrow current shell appears around the twisted flux tube in the dynamic phase of its expansion. This current shell is formed by the return currents of the system, which separate the twisted flux tube from the surrounding fields. It intensifies as the flux tube accelerates and it is co-spatial with weak plasma compression. The current density integrated over the altitude has the shape of an ellipse, which expands and rotates when viewed from above, reproducing the generic properties of propagating EIT waves. The timing, orientation, and location of bright and faint patches observed in the two EIT waves are remarkably well reproduced. We conjecture that propagating EIT waves are the observational signature of Joule heating in electric current shells, which separate expanding flux tubes from their surrounding fields during CMEs or plasma compression inside this current shell. We also conjecture that the bright edges of halo CMEs show the plasma compression in these current shells.     

Evidence of Magnetic Helicity in Emerging Flux and Associated Flare

The aim of this paper is to look at the magnetic helicity structure of an emerging active region and show that both emergence and flaring signatures are consistent with a same sign for magnetic helicity. We present a multiwavelength analysis of an M1.6 flare occurring in the NOAA active region 10365 on 27 May 2003, in which a large new bipole emerges in a decaying active region. The diverging flow pattern and the “tongue” shape of the magnetic field in the photosphere with elongated polarities are highly suggestive of the emergence of a twisted flux tube. The orientation of these tongues indicates the emergence of a flux tube with a right-hand twist (i.e., positive magnetic helicity). The flare signatures in the chromosphere are ribbons observed in Hα by the MSDP spectrograph in the Meudon solar tower and in 1600 Å by TRACE. These ribbons have a J shape and are shifted along the inversion line. The pattern of these ribbons suggests that the flare was triggered by magnetic reconnection at coronal heights below a twisted flux tube of positive helicity, corresponding to that of the observed emergence. It is the first time that such a consistency between the signatures of the emerging flux through the photosphere and flare ribbons has been clearly identified in observations. Another type of ribbons observed during the flare at the periphery of the active region by the MSDP and SOHO/EIT is related to the existence of a null point, which is found high in the corona in a potential field extrapolation. We discuss the interpretation of these secondary brightenings in terms of the “breakout” model and in terms of plasma compression/heating within large-scale separatrices.     

Topological Departures from Translational Invariance along a Filament Observed by THEMIS

We study the topology of the 3D magnetic field in a filament channel to address the following questions: Is a filament always formed in a single flux tube? How does the photospheric magnetic field lead to filament interruptions and to feet formation? What is the relation between feet-related field lines and the parasitic polarities? What can topological analyses teach us about EUV filament channels? To do so, we consider a filament observed on 6 October 2004 with THEMIS/MTR, in Hα with the full line profile simultaneously and cospatially with its photospheric vector magnetic field. The coronal magnetic field was calculated from a “linear magnetohydrostatic” extrapolation of a composite THEMIS-MDI magnetogram. Its free parameters were adjusted to get the best match possible between the distribution of modeled plasma-supporting dips and the Hα filament morphology. The model results in moderate plasma β≤1 at low altitudes in the filament, in conjunction with non-negligible departures from force-freeness measured by various metrics. The filament here is formed by a split flux tube. One part of the flux tube is rooted in the photosphere aside an observed interruption in the filament. This splitted topology is due to strong network polarities on the edge of the filament channel, not to flux concentrations closer to the filament. We focus our study to the northwest portion of the filament. The related flux tube is highly fragmented at low altitudes. This fragmentation is due to small flux concentrations of two types. First, some locally distort the tube, leading to noticeable thickness variations along the filament body. Second, parasitic polarities, associated with filament feet, result in secondary dips above the related local inversion line. These dips belong to long field lines that pass below the flux tube. Many of these field lines are not rooted near the related foot. Finally, the present model shows that the coronal void interpretation cannot be ruled out to interpret the wideness of EUV filament channels.     

Dynamical behavior of thermal protons in the mid-latitude ionosphere and magnetosphere

Satellite and other observations have shown that H⁺ densities in the mid-latitude topside ionosphere are greatly reduced during magnetic storms when the plasmapause and magnetic field convection move to relatively low L-values. In the recovery phase of the magnetic storm the convection region moves to higher L-values and replenishment of H⁺ in the empty magnetospheric field tubes begins. The upwards flow of H⁺, which arises from O⁺—H charge exchange, is initially supersonic. However, as the field tubes fill with plasma, a shock front moves downwards towards the ionosphere, eventually converting the upwards flow to subsonic speeds. The duration of this supersonic recovery depends strongly on the volume of the field tube; for example calculations indicate that for L = 5 the time is approximately 22 hours. The subsonic flow continues until diffusive equilibrium is reached or a new magnetic storm begins. Calculations of the density and flux profiles expected during the subsonic phase of the recovery show that diffusive equilibrium is still not reached after an elapsed time of 10 days and correspondingly there is still a net loss of plasma from the ionosphere to the magnetosphere at that time. This slow recovery of the H⁺ density and flux patterns, following magnetic storms, indicates that the mid-latitude topside ionosphere may be in a continual dynamic state if the storms occur sufficiently often.     

Basic properties of magnetic flux tubes and restrictions on theories of solar activity

It is shown that the mean longitudinal field in a magnetic flux tube is reduced, rather than enhanced, by twisting the tube to form a rope. It is shown that there is no magnetohydrostatic equilibrium when one twisted rope is wound around another. Instead there is rapid line cutting (neutral point annihilation). It is shown that the twisting increases, and the field strength decreases, along a flux tube extending upward through a stratified atmosphere.These facts are at variance with Piddington's recent suggestion that solar activity is to be understood as the result of flux tubes which are enormously concentrated by twisting, which consist of several twisted ropes wound around each other, and which came untwisted where they emerge through the photosphere.This work was supported in part by the National Aeronautics and Space Administration under Grant NGL 14-001-001.     

The equilibrium of coronal flux tubes under toroidal forces

The properties of slender isolated flux tubes, taking into account curvature effects, were investigated by Parker (1975, 1979) and Spruit (1981), and many studies have been made concerning the equilibrium of slender flux tubes in the solar corona. In this paper we use a different approach considering the coronal loop as a part of a circular torus and studying the position of its top when the loop is in equilibrium under toroidal forces. Toroidal forces were considered by Shafranov (1966) for toroidal pinches and the equilibrium can be studied for different values of the toroidal current intensity and external magnetic field. The results show that it is possible to have a coronal flux tube in equilibrium without considering gravity and external magnetic field. Furthermore, the total twist of the flux tube and its variation with the toroidal intensity has been studied.     

Evolution of emerging force-free flux tube to be a flare loop

Solar Physics - Evolution of a filamentary magnetic flux tube emerging from the photosphere is investigated in the assumption that the magnetic field is force-free and unchanged during the...