Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, forcefree vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.

Solar flares and coronal mass ejections are associated with rapid changes in field connectivity and are powered by the partial dissipation of electrical currents in the solar atmosphere. A critical unanswered question is whether the currents involved are induced by the motion of preexisting atmospheric magnetic flux subject to surface plasma flows or whether these currents are associated with the emergence of flux from within the solar convective zone. We address this problem by applying state-of-the-art nonlinear force-free field (NLFFF) modeling to the highest resolution and quality vector-magnetographic data observed by the recently launched Hinode satellite on NOAA AR 10930 around the time of a powerful X3.4 flare. We compute 14 NLFFF models with four different codes and a variety of boundary conditions. We find that the model fields differ markedly in geometry, energy content, and force-freeness. We discuss the relative merits of these models in a general critique of present abilities to model the coronal magnetic field based on surface vector field measurements. For our application in particular, we find a fair agreement of the best-fit model field with the observed coronal configuration, and argue (1) that strong electrical currents emerge together with magnetic flux preceding the flare, (2) that these currents are carried in an ensemble of thin strands, (3) that the global pattern of these currents and of field lines are compatible with a large-scale twisted flux rope topology, and (4) that the $10 32 erg change in energy associated with the coronal electrical currents suffices to power the flare and its associated coronal mass ejection.

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

No abstract

No abstract

No abstract

This Letter is devoted to the still open problem of the evolution of a three-dimensional coronal flux tube embedded in a low-beta ideal plasma and having its footpoints twisted by slow photospheric motions. Such a process has been simulated with a recently developed magnetohydrodynamic code. In the particular calculation reported here, the system occupies a large cubic box. The field is initially potential, being generated by an underlying horizontal dipole, and it is twisted by two vortices located on the lower face { z ϭ 0} of the box, on both sides of the neutral line. In a first phase, the field roughly evolves quasi-statically through a sequence of force-free configurations. Thus, it enters a dynamical phase during which it suffers a very fast expansion, closely approaching after some finite time a semiopen configuration. The energy increases monotonically during all the evolution, and it tends to a limit, which is equal to about 80% of the energy of the totally open field associated with B z .

We present and compare two methods for the reconstruction of the solar coronal magnetic field, assumed to be force-free, from photospheric boundary data. Both methods rely on a well posed mathematical boundary value problem and are of the Grad-Rubin type, i.e., the couple (B, α) is computed iteratively. They do differ from each other on the one hand by the way they address the zero-divergence of B issue, and on the other hand by the scheme they use for computing α at each iteration. The comparison of the two methods is done by numerically computing two examples of nonlinear force-free fields associated to large scale strong electric current distributions, whose exact forms can be otherwise determined semi-analytically. In particular, the second solution has a large nonlinearity even in the weak field region -a feature which is not present in the actual magnetograms, but is interesting to consider as it does allow to push the methods to the limits of their range of validity. The best results obtained with those methods give a relative vector error smaller than 0.01. For the latter extreme case, our results show that higher resolution reconstructions with bounded convergence improve the approximated solution, which may be of some interest for the treatment of the data of future magnetographs.