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.     

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…     

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.     

Coronal Flux Rope Equilibria in Closed Magnetic Fields

Using a 2.5-dimensional ideal MHD model in Cartesian coordinates,we investigate the equilibrium properties of coronal magnetic flux ropes in background magnetic fields that are completely closed.The background fields are produced by a dipole,a quadrupole,and an octapole,respectively,located below the photosphere at the same depth.A magnetic flux rope is then launched from below the photo-sphere,and its magnetic properties,i.e,the annular magnetic fluxφp and the axial magnetic fluxφz,are controlled by a single emergence parameter.The whole sys-tem eventually evolves into equilibrium,and the resultant flux rope is characterized by three geometrical parameters:the height of the rope axis,the half-width of the rope,and the length of the vertical current sheet below the rope.It is found that the geometrical parameters increase monotonically and continuously with increasing φp and φz:no catastrophe occurs.Moreover,there exists a steep segment in the profiles of the geometrical parameters versus either φp or φz,and the faster the background field decays with height,the larger both the gradient and the growth amplitude within the steep segment will be.     

磁云传播的动力学演化过程研究进展

磁云因其独特的磁场结构经常是重大灾害性空间天气的驱动源. 近来从磁云的边界层结构、环向通量、大尺度结构等方面关于磁云传播的动力学演化过程的研究取得了一些进展. 在磁云边界存在一个由于磁场重联而形成的边界层结构. 在磁云传播过程中, 这种发生在边界处的磁场重联可能会把磁云的磁场剥蚀掉, 进而引起其磁通量绳结构环向通量的减少以及不对称. 在磁云内部, 经常会观测到多个子通量绳结构. 这些特性各异的子通量绳可以通过磁场重联而合并, 进而引起磁云磁结构的改变. 关于磁云大尺度磁场拓扑位形的演化机制, 除了较早提出的交换重联外, 目前的研究表明在行星际空间中, 磁云边界处的重联过程也可以将磁云闭合或半开放的磁场线打开或断开. 尽管在相关研究中已经取得了较大进展, 但关于磁云传播的动力学演化过程还有许多问题尚不清楚. 在行星际小尺度磁通量绳边界也发现了边界层结构, 那么磁云是否会因剥蚀而成为小尺度通量绳? 磁云内子通量绳结构在相互作用中会不会引起某些不稳定性而导致整个通量绳系统的崩溃? 这些问题的解决还有待于进一步的理论、观测和数值模拟研究.     

A magnetohydrodynamical model for the formation of episodic jets

Episodic ejection of plasma blobs has been observed in many black hole systems. While steady, continuous jets are believed to be associated with large-scale open magnetic fields, what causes the episodic ejection of blobs remains unclear. Here by analogy with the coronal mass ejection on the Sun, we propose a magnetohydrodynamical model for episodic ejections from black holes associated with the closed magnetic fields in an accretion flow. Shear and turbulence of the accretion flow deform the field and result in the formation of a flux rope in the disc corona. Energy and helicity are accumulated and stored until a threshold is reached. The system then loses its equilibrium and the flux rope is thrust outward by the magnetic compression force in a catastrophic way. Our calculations show that for parameters appropriate for the black hole in our Galactic centre, the plasmoid can attain relativistic speeds in about 35 min.     

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.     

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.     

Coronal Flux Rope Catastrophe in Octapole Magnetic Fields

As demonstrated by many previous studies, a system consisting of an isolated coronal flux rope and a surrounding background magnetic field exhibits a catastrophic behavior. In particular, if the magnetic field of the system is force-free and axisymmetric in spherical geometry, the magnetic energy at the catastrophic point, referred to as the catastrophic energy threshold, is found to be larger than the corresponding partly or fully open field energy. This paper takes an octapole field as the background and introduces a flux rope within the central arcade of the octapole field. A relaxation method based on time-dependent ideal magnetohydrodynamic (MHD) simulations is used to find axisymmetric force-free field solutions in spherical geometry associated with the flux rope system. With respect to an increase of either the annular flux Φ_p or the axial flux Φ_ϕ of the rope, the system exhibits a catastrophic behavior as expected, and the catastrophic energy threshold is larger than that of the corresponding partly open field, in which the central arcade is opened up, but the remainder remains closed. For a given octapole field, the energy threshold depends on either Φ_p or Φ_ϕ at the catastrophic point, and it increases with increasing Φ_p or decreasing Φ_ϕ. On the other hand, the extent to which the central bipolar component of the octapole field is open also affects the energy threshold. These results differ from those for the bipolar background field case, in which the catastrophic energy threshold is almost independent of the magnetic properties of the flux rope at the catastrophic points and the extent to which the background field is open. The reason for such a difference is briefly discussed.     

Numerical experiments on the evolution in coronal magnetic configurations including a filament in response to the change in the photosphere

We investigate equilibrium height of a flux rope, and its internal equilibrium in a realistic plasma environment by carrying out numerical simulations of the evolution of systems including a current-carrying flux rope. We find that the equilibrium height of a flux rope is approximately described by a power-law function of the relative strength of the background field. Our simulations indicate that the flux rope can escape more easily from a weaker background field. This further confirms that a catastrophe in the magnetic configuration of interest can be triggered by a decrease in strength of the background field. Our results show that it takes some time to reach internal equilibrium depending on the initial state of the flux rope. The plasma flow inside the flux rope due to the adjustment for the internal equilibrium of the flux rope remains small and does not last very long when the initial state of the flux rope commences from the stable branch of the theoretical equilibrium curve. This work also confirms the influence of the initial radius of the flux rope in its evolution; the results indicate that a flux rope with a larger initial radius erupts more easily. In addition, by using a realistic plasma environment and a much higher resolution in our simulations,we notice some different characteristics compared to previous studies in Forbes.     

Basic Properties of Mutual Magnetic Helicity

We derive the magnetic helicity for configurations formed by flux tubes contained fully or only partially in the spatial domain considered (called closed and open configurations, respectively). In both cases, magnetic helicity is computed as the sum of mutual helicity over all possible pairs of magnetic flux tubes weighted by their magnetic fluxes. We emphasize that these mutual helicities have properties which are not those of mutual inductances in classical circuit theory. For closed configurations, the mutual helicity of two closed flux tubes is their relative winding around each other (known as the Gauss linkage number). For open configurations, the magnetic helicity is derived directly from the geometry of the interlaced flux tubes so it can be computed without reference to a ground state (such as a potential field). We derive the explicit expression in the case of a planar and spherical boundary. The magnetic helicity has two parts. The first one is given only by the relative positions of the flux tubes on the boundary. It is the only part if all flux tubes are arch-shaped. The second part counts the integer number of turns each pair of flux tubes wind about each other. This provides a general method to compute the magnetic helicity with discrete or continuous distributions of magnetic field. The method sets closed and open configurations on an equal level within the same theoretical framework.     

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.     

Reconnection of a Kinking Flux Rope Triggering the Ejection of a Microwave and Hard X-Ray Source II. Numerical Modeling

Numerical simulations of the helical (m=1) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range Φ≳6π. Such a reconnection is complex, involving also the ambient field. In addition to breaking up the original rope, it can form a new, low-lying, less twisted flux rope. The new flux rope is pushed downward by the reconnection outflow, which typically forces it to break as well by reconnecting with the ambient field. The top part of the original rope, largely rooted in the sources of the ambient flux after the break-up, can fully erupt or be halted at low heights, producing a “failed eruption.” The helical current sheet associated with the instability is squeezed between the approaching legs, temporarily forming a double current sheet. The leg – leg reconnection proceeds at a high rate, producing sufficiently strong electric fields that it would be able to accelerate particles. It may also form plasmoids, or plasmoid-like structures, which trap energetic particles and propagate out of the reconnection region up to the top of the erupting flux rope along the helical current sheet. The kinking of a highly twisted flux rope involving leg – leg reconnection can explain key features of an eruptive but partially occulted solar flare on 18 April 2001, which ejected a relatively compact hard X-ray and microwave source and was associated with a fast coronal mass ejection.     

Multiple magnetic clouds in interplanetary space

An interplanetary magnetic cloud (MC) is usually considered the byproduct of a coronal mass ejection (CME). Due to the frequent occurrence of CMEs, multiple magnetic clouds (multi-MCs), in which one MC catches up with another, should be a relatively common phenomenon. A simple flux rope model is used to get the primary magnetic field features of multi-MCs. Results indicate that the magnetic field configuration of multi-MCs mainly depends on the magnetic field characteristics of each member of multi-MCs. It may be entirely different in another situation. Moreover, we fit the data from the Wind spacecraft by using this model. Comparing the model with the observations, we verify the existence of multi-MCs, and propose some suggestions for further work.     

Magnetic braking in magnetic binary stars

The role of an external magnetic field in the magnetic braking of a star with a dipolar field is investigated. In a magnetic cataclysmic variable system (i.e. the primary compact star has a strong magnetic field), the field external to the braking star (a late-type main-sequence star with a dynamo-generated field) originates from the compact star. A closed field region — the system dead zone — is formed within the binary system, and it does not take part in magnetic braking. The overall braking rate depends on the extent of this region and of the open flux, and is dependent on centrifugal effects. In the case of two interacting dipoles, the dipole orientations relative to the spin axes and to each other are found to be important, leading to different amounts of open flux and therefore of magnetic braking, owing to different centrifugal effects on closed field regions. However, in circumstances consistent with observations and dynamo theory, the white dwarf's field reduces the magnetic braking of the secondary significantly, a finding qualitatively similar to the results previously obtained for two anti-aligned dipoles perpendicular to the orbital plane. In the cases where the two dipole axes are not perpendicular to the orbital plane, but are inclined in the plane that links them, the 'cut-off' in magnetic braking is less abrupt, and this effect is more obvious as the inclinations increase. Only in the extreme cases when the two dipole axes are aligned in the orbital plane does the braking increase with white dwarf field strength. We conclude that detailed evolutionary modelling of AM Herculis systems needs to take account of the inclination effect.     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

Convective intensification of magnetic fields in the quiet Sun

Kilogauss-strength magnetic fields are often observed in intergranular lanes at the photosphere in the quiet Sun. Such fields are stronger than the equipartition field B _e, corresponding to a magnetic energy density that matches the kinetic energy density of photospheric convection, and comparable with the field B _p that exerts a magnetic pressure equal to the ambient gas pressure. We present an idealized numerical model of three-dimensional compressible magnetoconvection at the photosphere, for a range of values of the magnetic Reynolds number. In the absence of a magnetic field, the convection is highly supercritical and characterized by a pattern of vigorous, time-dependent, 'granular' motions. When a weak magnetic field is imposed upon the convection, magnetic flux is swept into the convective downflows where it forms localized concentrations. Unless this process is significantly inhibited by magnetic diffusion, the resulting fields are often much greater than B _e and the high magnetic pressure in these flux elements leads to their being partially evacuated. Some of these flux elements contains ultraintense magnetic fields that are significantly greater than B _p. Such fields are contained by a combination of the thermal pressure of the gas and the dynamic pressure of the convective motion, and they are constantly evolving. These ultraintense fields develop owing to non-linear interactions between magnetic fields and convection; they cannot be explained in terms of 'convective collapse' within a thin flux tube that remains in overall pressure equilibrium with its surroundings.     

A Numerical Study of the Response of the Coronal Magnetic Field to Flux Emergence

Large-scale solar eruptions, known as coronal mass ejections (CMEs), are regarded as the main drivers of space weather. The exact trigger mechanism of these violent events is still not completely clear; however, the solar magnetic field indisputably plays a crucial role in the onset of CMEs. The strength and morphology of the solar magnetic field are expected to have a decisive effect on CME properties, such as size and speed. This study aims to investigate the evolution of a magnetic configuration when driven by the emergence of new magnetic flux in order to get a better insight into the onset of CMEs and their magnetic structure. The three-dimensional, time-dependent equations for ideal magnetohydrodynamics are numerically solved on a spherical mesh. New flux emergence in a bipolar active region causes destabilisation of the initial stationary structure, finally resulting in an eruption. The initial magnetic topology is suitable for the ??breakout?? CME scenario to work. Although no magnetic flux rope structure is present in the initial condition, highly twisted magnetic field lines are formed during the evolution of the system as a result of internal reconnection due to the interaction of the active region magnetic field with the ambient field. The magnetic energy built up in the system and the final speed of the CME depend on the strength of the overlying magnetic field, the flux emergence rate, and the total amount of emerged flux. The interaction with the global coronal field makes the eruption a large-scale event, involving distant parts of the solar surface.     

On-Disc Observations of Flux Rope Formation Prior to Its Eruption

Coronal mass ejections (CMEs) are one of the primary manifestations of solar activity and can drive severe space weather effects. Therefore, it is vital to work towards being able to predict their occurrence. However, many aspects of CME formation and eruption remain unclear, including whether magnetic flux ropes are present before the onset of eruption and the key mechanisms that cause CMEs to occur. In this work, the pre-eruptive coronal configuration of an active region that produced an interplanetary CME with a clear magnetic flux rope structure at 1 AU is studied. A forward-S sigmoid appears in extreme-ultraviolet (EUV) data two hours before the onset of the eruption (SOL2012-06-14), which is interpreted as a signature of a right-handed flux rope that formed prior to the eruption. Flare ribbons and EUV dimmings are used to infer the locations of the flux rope footpoints. These locations, together with observations of the global magnetic flux distribution, indicate that an interaction between newly emerged magnetic flux and pre-existing sunspot field in the days prior to the eruption may have enabled the coronal flux rope to form via tether-cutting-like reconnection. Composition analysis suggests that the flux rope had a coronal plasma composition, supporting our interpretation that the flux rope formed via magnetic reconnection in the corona. Once formed, the flux rope remained stable for two hours before erupting as a CME.     

A Comparative Study of Coronal Mass Ejections with and Without Magnetic Cloud Structure near the Earth: Are All Interplanetary CMEs Flux Ropes?

An outstanding question concerning interplanetary coronal mass ejections (ICMEs) is whether all ICMEs have a magnetic flux rope structure. We test this question by studying two different ICMEs, one having a magnetic cloud (MC) showing smooth rotation of magnetic field lines and the other not. The two ICMEs are chosen in such a way that their progenitor CMEs are very similar in remote sensing observations. Both CMEs originated from close to the central meridian directly facing the Earth. Both CMEs were associated with a long-lasting post-eruption loop arcade and appeared as an elliptical halo in coronagraph images, indicating a flux rope origin. We conclude that the difference in the in-situ observation is caused by the geometric selection effect, contributed by the deflection of flux ropes in the inner corona and interplanetary space. The first event had its nose pass through the observing spacecraft; thus, the intrinsic flux rope structure of the CME appeared as a magnetic cloud. On the other hand, the second event had the flank of the flux rope intercept the spacecraft, and it thus did not appear as a magnetic cloud. We further argue that a conspicuous long period of weak magnetic field, low plasma temperature, and density in the second event should correspond to the extended leg portion of the embedded magnetic flux rope, thus validating the scenario of the flank-passing. These observations support the idea that all CMEs arriving at the Earth include flux rope drivers.