The topology of the coronal magnetic field produces a strong impact on the properties of the solar corona and presumably on the origin of the slow solar wind. To advance our understanding of this impact, we revisit the concept of so-called slip-back mapping (Titov et al. 2009) and adapt it for determining open, closed, and disconnected flux systems that are formed in the solar corona by magnetic reconnection during a given time interval. The developed method allows us, in particular, to describe the magnetic flux transfer between open and closed flux regions via so-called interchange reconnection with an unprecedented level of details. We illustrate the application of this method to the analysis of the global MHD evolution of the solar corona that is driven by an idealized differential rotation of the photospheric plasma.

GPU accelerators have had a notable impact on high-performance computing across many disciplines. They provide high performance with low cost/power, and therefore have become a primary compute resource on many of the largest supercomputers. Here, we implement multi-GPU acceleration into our Solar MHD code (MAS) using OpenACC in a fully portable, single-source manner. Our preliminary implementation is focused on MAS running in a reduced physics "zero-beta" mode. While valuable on its own, our main goal is to pave the way for a full physics, thermodynamic MHD implementation. We describe the OpenACC implementation methodology and challenges. "Time-to-solution" performance results of a production-level flux rope eruption simulation on multi-CPU and multi-GPU systems are shown. We find that the GPU-accelerated MAS code has the ability to run "zero-beta" simulations on a single multi-GPU server at speeds previously requiring multiple CPU server-nodes of a supercomputer.

Solar eruptions are the main driver of space-weather disturbances at the Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 14 July 2000 Bastille Day eruption, which produced a very strong geomagnetic storm. After constructing a thermodynamic MHD model of the corona and solar wind, we insert a magnetically stable flux rope along the polarity inversion line of the eruption's source region and initiate the eruption by boundary flows. More than 10^33 ergs of magnetic energy are released in the eruption within a few minutes, driving a flare, an EUV wave, and a coronal mass ejection (CME) that travels in the outer corona at about 1500 km/s, close to the observed speed. We then propagate the CME to Earth, using a heliospheric MHD code. Our simulation thus provides the opportunity to test how well in situ observations of extreme events are matched if the eruption is initiated from a stable magnetic-equilibrium state. We find that the flux-rope center is very similar in character to the observed magnetic cloud, but arrives about 8.5 hours later and about 15 degrees too far to the North, with field strengths that are too weak by a factor of about 1.6. The front of the flux rope is highly distorted, exhibiting localized magnetic-field concentrations as it passes 1 AU. We discuss these properties with regard to the development of space-weather predictions based on MHD simulations of solar eruptions.

Many existing models assume that magnetic flux ropes play a key role in solar flares and coronal mass ejections (CMEs). It is therefore important to develop efficient methods for constructing flux-rope configurations constrained by observed magnetic data and the morphology of the pre-eruptive source region. For this purpose, we have derived and implemented a compact analytical form that represents the magnetic field of a thin flux rope with an axis of arbitrary shape and circular cross-sections. This form implies that the flux rope carries axial current $I$ and axial flux $F$, so that the respective magnetic field is the curl of the sum of axial and azimuthal vector potentials proportional to $I$ and $F$, respectively. We expressed the vector potentials in terms of modified Biot-Savart laws whose kernels are regularized at the axis in such a way that, when the axis is straight, these laws define a cylindrical force-free flux rope with a parabolic profile for the axial current density. For the cases we have studied so far, we determined the shape of the rope axis by following the polarity inversion line of the eruptions' source region, using observed magnetograms. The height variation along the axis and other flux-rope parameters are estimated by means of potential field extrapolations. Using this heuristic approach, we were able to construct pre-eruption configurations for the 2009 February 13 and 2011 October 1 CME events. These applications demonstrate the flexibility and efficiency of our new method for energizing pre-eruptive configurations in simulations of CMEs.

The physical conditions that determine whether or not solar active regions (ARs) produce strong flares and coronal mass ejections (CMEs) are not yet well understood. Here we investigate the association between electric-current neutralization, magnetic shear along polarity inversion lines (PILs), and eruptive activity in four ARs; two emerging and two well-developed ones. We find that the CME-producing ARs are characterized by a strongly non-neutralized total current, while the total current in the ARs that did not produce CMEs is almost perfectly neutralized. The difference in the PIL-shear between these two groups is much less pronounced, which suggests that the degree of current-neutralization may serve as a better proxy for assessing the ability of ARs to produce CMEs.

Using a potential field source surface (PFSS) model, we recently analyzed the global topology of the background coronal magnetic field for a sequence of coronal mass ejections (CMEs) that occurred on 2010 August 1-2. Here we repeat this analysis for the background field reproduced by a magnetohydrodynamic (MHD) model that incorporates plasma thermodynamics. As for the PFSS model, we find that all three CME source regions contain a coronal hole that is separated from neighboring coronal holes by topologically very similar pseudo-streamer structures. However, the two models yield very different results for the size, shape, and flux of the coronal holes. We find that the helmet-streamer cusp line, which corresponds to a source-surface null line in the PFSS model, is structurally unstable and does not form in the MHD model. Our analysis indicates that generally, in MHD configurations, this line rather consists of a multiple-null separator passing along the edge of disconnected flux regions. Some of these regions are transient and may be the origin of so-called streamer blobs. We show that the core topological structure of such blobs is a three-dimensional "plasmoid", consisting of two conjoined flux ropes of opposite handedness, which connect at a spiral null point of the magnetic field. Our analysis reveals that such plasmoids appear also in pseudo-streamers on much smaller scales. These new insights into the coronal magnetic topology provide some intriguing implications for solar energetic particle events and for the properties of the slow solar wind.

Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal jets using a three-dimensional (3D) magnetohydrodynamic (MHD) code with thermodynamics in a large spherical domain that includes the solar wind. Our model is coupled to 3D MHD flux-emergence simulations, i.e, we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. The model includes parametric coronal heating, radiative losses, and thermal conduction, which enables us to simulate the dynamics and plasma properties of coronal jets in a more realistic manner than done so far. Here we employ these simulations to calculate the amount of mass and energy transported by coronal jets into the outer corona and inner heliosphere. Based on observed jet-occurrence rates, we then estimate the total contribution of coronal jets to the mass and energy content of the solar wind to (0.4-3.0) % and (0.3-1.0) %, respectively. Our results are largely consistent with the few previous rough estimates obtained from observations, supporting the conjecture that coronal jets provide only a small amount of mass and energy to the solar wind. We emphasize, however, that more advanced observations and simulations are needed to substantiate this conjecture.

We investigate the evolution of NOAA Active Region 11817 during 2013 August 10--12, when it developed a complex field configuration and produced four confined, followed by two eruptive, flares. These C-and-above flares are all associated with a magnetic flux rope (MFR) located along the major polarity inversion line, where shearing and converging photospheric flows are present. Aided by the nonlinear force-free field modeling, we identify the MFR through mapping magnetic connectivities and computing the twist number $\mathcal{T}_w$ for each individual field line. The MFR is moderately twisted ($|\mathcal{T}_w| < 2$) and has a well-defined boundary of high squashing factor $Q$. We found that the field line with the extremum $|\mathcal{T}_w|$ is a reliable proxy of the rope axis, and that the MFR's peak $|\mathcal{T}_w|$ temporarily increases within half an hour before each flare while it decreases after the flare peak for both confined and eruptive flares. This pre-flare increase in $|\mathcal{T}_w|$ has little effect on the active region's free magnetic energy or any other parameters derived for the whole region, due to its moderate amount and the MFR's relatively small volume, while its decrease after flares is clearly associated with the stepwise decrease in free magnetic energy due to the flare. We suggest that $\mathcal{T}_w$ may serve as a useful parameter in forewarning the onset of eruption, and therefore, the consequent space weather effects. The helical kink instability is identified as the prime candidate onset mechanism for the considered flares.

Coronal jets are transient, collimated eruptions that occur in regions of predominantly open magnetic field in the solar corona. Our understanding of these events has greatly evolved in recent years but several open questions, such as the contribution of coronal jets to the solar wind, remain. Here we present an overview of the observations and numerical modeling of coronal jets, followed by a brief description of "next-generation" simulations that include an advanced description of the energy transfer in the corona ("thermodynamic MHD"), large spherical computational domains, and the solar wind. These new models will allow us to address some of the open questions.

Force-free equilibria containing two vertically arranged magnetic flux ropes of like chirality and current direction are considered as a model for split filaments/prominences and filament-sigmoid systems. Such equilibria are constructed analytically through an extension of the methods developed in Titov & D\'emoulin (1999) and numerically through an evolutionary sequence including shear flows, flux emergence, and flux cancellation in the photospheric boundary. It is demonstrated that the analytical equilibria are stable if an external toroidal (shear) field component exceeding a threshold value is included. If this component decreases sufficiently, then both flux ropes turn unstable for conditions typical of solar active regions, with the lower rope typically being unstable first. Either both flux ropes erupt upward, or only the upper rope erupts while the lower rope reconnects with the ambient flux low in the corona and is destroyed. However, for shear field strengths staying somewhat above the threshold value, the configuration also admits evolutions which lead to partial eruptions with only the upper flux rope becoming unstable and the lower one remaining in place. This can be triggered by a transfer of flux and current from the lower to the upper rope, as suggested by the observations of a split filament in Paper~I (Liu et al. 2012). It can also result from tether-cutting reconnection with the ambient flux at the X-type structure between the flux ropes, which similarly influences their stability properties in opposite ways. This is demonstrated for the numerically constructed equilibrium.

We report the observation of an X-class long-duration flare which is clearly confined. It appears as a compact-loop flare in the traditional EUV passbands (171 and 195 {\AA}), but in the passbands sensitive to flare plasmas (94 and 131 {\AA}), it exhibits a cusp-shaped structure above an arcade of loops like other long-duration events. Inspecting images in a running difference approach, we find that the seemingly diffuse, quasi-static cusp-shaped structure consists of multiple nested loops that repeatedly rise upward and disappear approaching the cusp edge. Over the gradual phase, we detect numerous episodes of loop rising, each lasting minutes. A differential emission measure analysis reveals that the temperature is highest at the top of the arcade and becomes cooler at higher altitudes within the cusp-shaped structure, contrary to typical long-duration flares. With a nonlinear force-free model, our analysis shows that the event mainly involves two adjacent sheared arcades separated by a T-type hyperbolic flux tube (HFT). One of the arcades harbors a magnetic flux rope, which is identified with a filament that survives the flare owing to the strong confining field. We conclude that a new emergence of magnetic flux in the other arcade triggers the flare, while the preexisting HFT and flux rope dictate the structure and dynamics of the flare loops and ribbons during the long-lasting decay phase, and that a quasi-separatrix layer high above the HFT could account for the cusp-shaped structure.

There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this question, and numerical simulations have, surprisingly, barely been used to address it. Here we quantify the evolution of electric currents during the formation of a bipolar active region by considering a three-dimensional magnetohydrodynamic simulation of the emergence of a sub-photospheric, current-neutralized magnetic flux rope into the solar atmosphere. We find that a strong deviation from current neutralization develops simultaneously with the onset of significant flux emergence into the corona, accompanied by the development of substantial magnetic shear along the active region's polarity inversion line. After the region has formed and flux emergence has ceased, the strong magnetic fields in the region's center are connected solely by direct currents, and the total direct current is several times larger than the total return current. These results suggest that active regions, the main sources of coronal mass ejections and flares, are born with substantial net currents, in agreement with recent observations. Furthermore, they support eruption models that employ pre-eruption magnetic fields containing such currents.

Using solar vector magnetograms of the highest available spatial resolution and signal-to-noise ratio we perform a detailed study of electric current patterns in two solar active regions: a flaring/eruptive, and a flare-quiet one. We aim to determine whether active regions inject non-neutralized (net) electric currents in the solar atmosphere, responding to a debate initiated nearly two decades ago that remains inconclusive. We find that well-formed, intense magnetic polarity inversion lines (PILs) within active regions are the only photospheric magnetic structures that support significant net current. More intense PILs seem to imply stronger non-neutralized current patterns per polarity. This finding revises previous works that claim frequent injections of intense non-neutralized currents by most active regions appearing in the solar disk but also works that altogether rule out injection of non-neutralized currents. In agreement with previous studies, we also find that magnetically isolated active regions remain globally current-balanced. In addition, we confirm and quantify the preference of a given magnetic polarity to follow a given sense of electric currents, indicating a dominant sense of twist in active regions. This coherence effect is more pronounced in more compact active regions with stronger PILs and must be of sub-photospheric origin. Our results yield a natural explanation of the Lorentz force, invariably generating velocity and magnetic shear along strong PILs, thus setting a physical context for the observed pre-eruption evolution in solar active regions.

A sequence of apparently coupled eruptions was observed on 2010 August 1-2 by SDO and STEREO. The eruptions were closely synchronized with one another, even though some of them occurred at widely separated locations. In an attempt to identify a plausible reason for such synchronization, we study the large-scale structure of the background magnetic configuration. The coronal field was computed from the photospheric magnetic field observed at the appropriate time period by using the potential field source-surface model. We investigate the resulting field structure by analyzing the so-called squashing factor calculated at the photospheric and source-surface boundaries, as well as at different coronal cross-sections. Using this information as a guide, we determine the underlying structural skeleton of the configuration, including separatrix and quasi-separatrix surfaces. Our analysis reveals, in particular, several pseudo-streamers in the regions where the eruptions occurred. Of special interest to us are the magnetic null points and separators associated with the pseudo-streamers. We propose that magnetic reconnection triggered along these separators by the first eruption likely played a key role in establishing the assumed link between the sequential eruptions. The present work substantiates our recent simplified magnetohydrodynamic model of sympathetic eruptions and provides a guide for further deeper study of these phenomena. Several important implications of our results for the S-web model of the slow solar wind are also addressed.

We study an active-region dextral filament which was composed of two branches separated in height by about 13 Mm. This "double-decker" configuration sustained for days before the upper branch erupted with a GOES-class M1.0 flare on 2010 August 7. Analyzing this evolution, we obtain the following main results. 1) During hours before the eruption, filament threads within the lower branch were observed to intermittently brighten up, lift upward, and then merge with the upper branch. The merging process contributed magnetic flux and current to the upper branch, resulting in its quasi-static ascent. 2) This transfer might serve as the key mechanism for the upper branch to lose equilibrium by reaching the limiting flux that can be stably held down by the overlying field or by reaching the threshold of the torus instability. 3) The erupting branch first straightened from a reverse S shape that followed the polarity inversion line and then writhed into a forward S shape. This shows a transfer of left-handed helicity in a sequence of writhe-twist-writhe. The fact that the initial writhe is converted into the twist of the flux rope excludes the helical kink instability as the trigger process of the eruption, but supports the occurrence of the instability in the main phase, which is indeed indicated by the very strong writhing motion. 4) A hard X-ray sigmoid, likely of coronal origin, formed in the gap between the two original filament branches in the impulsive phase of the associated flare. This supports a model of transient sigmoids forming in the vertical flare current sheet. 5) Left-handed magnetic helicity is inferred for both branches of the dextral filament. 6) Two types of force-free magnetic configurations are compatible with the data, a double flux rope equilibrium and a single flux rope situated above a loop arcade.

Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one eruption can trigger another one remain poorly understood. We present a 3D MHD simulation that suggests two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that contains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global sympathetic event on 2010 August 1 and provides a scenario for so-called twin filament eruptions. The suggested mechanisms are applicable also for sympathetic eruptions occurring in other magnetic configurations.

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind also has large angular width, up to ~ 60{\circ}, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far from the heliospheric current sheet. We then use an MHD code and MDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

In recent work, Antiochos and coworkers argued that the boundary between the open and closed field regions on the Sun can be extremely complex with narrow corridors of open flux connecting seemingly disconnected coronal holes from the main polar holes, and that these corridors may be the sources of the slow solar wind. We examine, in detail, the topology of such magnetic configurations using an analytical source surface model that allows for analysis of the field with arbitrary resolution. Our analysis reveals three important new results: First, a coronal hole boundary can join stably to the separatrix boundary of a parasitic polarity region. Second, a single parasitic polarity region can produce multiple null points in the corona and, more important, separator lines connecting these points. It is known that such topologies are extremely favorable for magnetic reconnection, because they allow this process to occur over the entire length of the separators rather than being confined to a small region around the nulls. Finally, the coronal holes are not connected by an open-field corridor of finite width, but instead are linked by a singular line that coincides with the separatrix footprint of the parasitic polarity. We investigate how the topological features described above evolve in response to the motion of the parasitic polarity region. The implications of our results for the sources of the slow solar wind and for coronal and heliospheric observations are discussed.

We examine the nonlinear magnetofrictional extrapolation scheme using the solar active region model by Titov and D\'emoulin as test field. This model consists of an arched, line-tied current channel held in force-free equilibrium by the potential field of a bipolar flux distribution in the bottom boundary. A modified version, having a parabolic current density profile, is employed here. We find that the equilibrium is reconstructed with very high accuracy in a representative range of parameter space, using only the vector field in the bottom boundary as input. Structural features formed in the interface between the flux rope and the surrounding arcade-"hyperbolic flux tube" and "bald patch separatrix surface"-are reliably reproduced, as are the flux rope twist and the energy and helicity of the configuration. This demonstrates that force-free fields containing these basic structural elements of solar active regions can be obtained by extrapolation. The influence of the chosen initial condition on the accuracy of reconstruction is also addressed, confirming that the initial field that best matches the external potential field of the model quite naturally leads to the best reconstruction. Extrapolating the magnetogram of a Titov-D\'emoulin equilibrium in the unstable range of parameter space yields a sequence of two opposing evolutionary phases which clearly indicate the unstable nature of the configuration: a partial buildup of the flux rope with rising free energy is followed by destruction of the rope, losing most of the free energy.

A general method for describing magnetic reconnection in arbitrary three-dimensional magnetic configurations is proposed. The method is based on the field-line mapping technique previously used only for the analysis of magnetic structure at a given time. This technique is extended here so as to analyze the evolution of magnetic structure. Such a generalization is made with the help of new dimensionless quantities called "slip-squashing factors". Their large values define the surfaces that border the reconnected or to-be-reconnected magnetic flux tubes for a given period of time during the magnetic evolution. The proposed method is universal, since it assumes only that the time sequence of evolving magnetic field and the tangential boundary flows are known. The application of the method is illustrated for simple examples, one of which was considered previously by Hesse and coworkers in the framework of the general magnetic reconnection theory. The examples help us to compare these two approaches; they reveal also that, just as for magnetic null points, hyperbolic and cusp minimum points of a magnetic field may serve as favorable sites for magnetic reconnection. The new method admits a straightforward numerical implementation and provides a powerful tool for the diagnostics of magnetic reconnection in numerical models of solar-flare-like phenomena in space and laboratory plasmas.

A simple model of the coronal magnetic field prior to the CME eruption on May 12 1997 is developed. First, the magnetic field is constructed by superimposing a large-scale background field and a localized bipolar field to model the active region (AR) in the current-free approximation. Second, this potential configuration is quasi-statically sheared by photospheric vortex motions applied to two flux concentrations of the AR. Third, the resulting force-free field is then evolved by canceling the photospheric magnetic flux with the help of an appropriate tangential electric field applied to the central part of the AR. To understand the structure of the modeled configuration, we use the field line mapping technique by generalizing it to spherical geometry. It is demonstrated that the initial potential configuration contains a hyperbolic flux tube (HFT) which is a union of two intersecting quasi-separatrix layers. This HFT provides a partition of the closed magnetic flux between the AR and the global solar magnetic field. The vortex motions applied to the AR interlock the field lines in the coronal volume to form additionally two new HFTs pinched into thin current layers. Reconnection in these current layers helps to redistribute the magnetic flux and current within the AR in the flux-cancellation phase. In this phase, a magnetic flux rope is formed together with a bald patch separatrix surface wrapping around the rope. Other important implications of the identified structural features of the modeled configuration are also discussed.

The study of magnetic connectivity in the solar corona reveals a need to generalize the field line mapping technique to arbitrary geometry of the boundaries and systems of coordinates. Indeed, the global description of the connectivity in the corona requires the use of the photospheric and solar wind boundaries. Both are closed surfaces and therefore do not admit a global regular system of coordinates. At least two overlapping regular systems of coordinates for each of the boundaries are necessary in this case to avoid spherical-pole-like singularities in the coordinates of the footpoints. This implies that the basic characteristic of magnetic connectivity - the squashing degree or factor $Q$ of elemental flux tubes (Titov et al., 2002) - must be rewritten in covariant form. Such a covariant expression of $Q$ is derived in this work. The derived expression is very flexible and highly efficient for describing the global magnetic connectivity in the solar corona. In addition, a general expression for a new characteristic $Q_\perp$ which defines a squashing of the flux tubes in the directions perpendicular to the field lines is determined. This new quantity makes it possible to filter out the quasi-separatrix layers whose large values of $Q$ are caused by a projection effect at the field lines nearly touching the photosphere. Thus, the value $Q_\perp$ provides a much more precise description of the volumetric properties of the magnetic field structure. The difference between $Q$ and $Q_\perp$ is illustrated by comparing their distributions for two configurations, one of which is the Titov-Demoulin (1999) model of a twisted magnetic field.

We study dynamical consequences of the kink instability of a twisted coronal flux rope, using the force-free coronal loop model by Titov & D\'emoulin (1999) as the initial condition in ideal-MHD simulations. When a critical value of the twist is exceeded, the long-wavelength ($m=1$) kink mode develops. Analogous to the well-known cylindrical approximation, a helical current sheet is then formed at the interface with the surrounding medium. In contrast to the cylindrical case, upward-kinking loops form a second, vertical current sheet below the loop apex at the position of the hyperbolic flux tube (generalized X line) in the model. The current density is steepened in both sheets and eventually exceeds the current density in the loop (although the kink perturbation starts to saturate in our simulations without leading to a global eruption). The projection of the field lines that pass through the vertical current sheet shows an S shape whose sense agrees with the typical sense of transient sigmoidal (forward or reverse S-shaped) structures that brighten in soft X rays prior to coronal eruptions. The upward-kinked loop has the opposite S shape, leading to the conclusion that such sigmoids do not generally show the erupting loops themselves but indicate the formation of the vertical current sheet below them that is the central element of the standard flare model.

The force-free coronal loop model by Titov & D\'emoulin (1999} is found to be unstable with respect to the ideal kink mode, which suggests this instability as a mechanism for the initiation of flares. The long-wavelength ($m=1$) mode grows for average twists $\Phi\ga3.5\pi$ (at a loop aspect ratio of $\approx$ 5). The threshold of instability increases with increasing major loop radius, primarily because the aspect ratio then also increases. Numerically obtained equilibria at subcritical twist are very close to the approximate analytical equilibrium; they do not show indications of sigmoidal shape. The growth of kink perturbations is eventually slowed down by the surrounding potential field, which varies only slowly with radius in the model. With this field a global eruption is not obtained in the ideal MHD limit. Kink perturbations with a rising loop apex lead to the formation of a vertical current sheet below the apex, which does not occur in the cylindrical approximation.

The concept of hyperbolic flux tubes (HFTs) is a generalization of the concept of separator field lines for coronal magnetic fields with a trivial magnetic topology. An effective mechanism of a current layer formation in HFTs is proposed. This mechanism is called magnetic pinching and it is caused by large-scale shearing motions applied to the photospheric feet of HFTs in a way as if trying to twist the HFT. It is shown that in the middle of an HFT such motions produce a hyperbolic flow that causes an exponentially fast growth of the current density in a thin force-free current layer. The magnetic energy associated with the current layer that is built up over a few hours is sufficient for a large flare. Other implications of HFT pinching for solar flares are discussed as well.

Although the analysis of observational data indicates that quasi-separatrix layers (QSLs) have to play an important role in magnetic configurations of solar flares, the corresponding theory is only at an initial stage so far. In particular, there is still a need in a proper definition of QSL. This problem is analyzed here on the basis of geometrical properties of the mapping produced by the field lines which connect photospheric areas of positive and negative magnetic polarities of active regions. In general, one can find on the photosphere a unique pair of locally perpendicular line elements at one footpoint of a given field line, so that this pair is mapped into a similar pair of elements at the other footpoint. Along the directions of these elements the field line mapping only stretches and compresses the corresponding displacements of the neighboring field-line footpoints. This fact enables us to get a remarkably clear and concise definition of QSL as a volume filled by coronal field lines for which the ratio of the corresponding stretching and compressing coefficients is anomalously large. The new definition is also compared with the ones previously introduced by other authors. The theory is applied to flare events of the so-called sigmoid-type, using an analytical model of a twisted force-free configuration. It is shown that such a configuration may contain a QSL even if genuine separatrix surfaces are absent. Identifying QSLs is an essential prerequisite for understanding the mechanism of magnetic energy release in this type of flares. It is also demonstrated that the magnetic field under study has a Hamiltonian structure, which makes it possible to reveal a geometrical reason for the appearance of QSL in this case.

We have collected the largest existing set of radio source lists in machine-readable form: 320 tables with 1.75 million records. Only a minor fraction of these is accessible via public databases. We describe our plans to make this huge amount of heterogeneous data accessible in a homogeneous way via the World Wide Web, with reliable cross-identifications, and searchable by various observables.