Nonlinear Force-Free Modeling of Coronal Magnetic Fields Part I: A Quantitative Comparison of Methods

We compare six algorithms for the computation of nonlinear force-free (NLFF) magnetic fields (including optimization, magnetofrictional, Grad–Rubin based, and Green's function-based methods) by evaluating their performance in blind tests on analytical force-free-field models for which boundary conditions are specified either for the entire surface area of a cubic volume or for an extended lower boundary only. Figures of merit are used to compare the input vector field to the resulting model fields. Based on these merit functions, we argue that all algorithms yield NLFF fields that agree best with the input field in the lower central region of the volume, where the field and electrical currents are strongest and the effects of boundary conditions weakest. The NLFF vector fields in the outer domains of the volume depend sensitively on the details of the specified boundary conditions; best agreement is found if the field outside of the model volume is incorporated as part of the model boundary, either as potential field boundaries on the side and top surfaces, or as a potential field in a skirt around the main volume of interest. For input field (B) and modeled field (b), the best method included in our study yields an average relative vector error E_n = 〈 |B−b|〉/〈 |B|〉 of only 0.02 when all sides are specified and 0.14 for the case where only the lower boundary is specified, while the total energy in the magnetic field is approximated to within 2%. The models converge towards the central, strong input field at speeds that differ by a factor of one million per iteration step. The fastest-converging, best-performing model for these analytical test cases is the Wheatland, Sturrock, and Roumeliotis (2000) optimization algorithm as implemented by Wiegelmann (2004).     

Nonlinear Force-Free Modeling of the Corona in Spherical Coordinates

We present a code for solving the nonlinear force-free equations in spherical polar geometry, with the motivation of modeling the magnetic field in the corona. The code is an implementation of the Grad–Rubin method. Our method is applicable to a spherical domain of arbitrary angular size. The implementation is based on a global spectral representation for the magnetic field that makes no explicit assumptions about the form of the magnetic field at the transverse boundaries of the domain. We apply the code to a bipolar test case with analytic boundary conditions, and demonstrate the convergence of the Grad–Rubin method and the self-consistency of the resulting numerical solution.     

Calculating and Testing Nonlinear Force-Free Fields

Improvements to an existing method for calculating nonlinear force-free magnetic fields (Wheatland, Solar Phys. 238, 29, 2006) are described. In particular a solution of the 3-D Poisson equation using 2-D Fourier transforms is presented. The improved nonlinear force-free method is demonstrated in application to linear force-free test cases with localized nonzero values of the normal component of the field in the boundary. These fields provide suitable test cases for nonlinear force-free calculations because the boundary conditions involve localized nonzero values of the normal components of the field and of the current density, and because (being linear force-free fields) they have more direct numerical solutions. Despite their simplicity, fields of this kind have not been recognized as test cases for nonlinear methods before. The examples illustrate the treatment of the boundary conditions on current in the nonlinear force-free method, and in particular the limitations imposed by field lines that connect outside of the boundary region.     

Helical structures around quiescent solar prominences computed from observable magnetic fields

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

Topological differences and similarities between force-free and potential models of coronal magnetic fields

Potential fields and linear force-free fields are often used as models for the magnetic field of the Sun's corona. They can be written as analytical expressions in terms of boundary values at the photosphere. Because of their relative simplicity compared with nonlinear force-free fields, these two models are of particular importance in topological analysis of solar phenomena. However, it has been suggested by Hudson and Wheatland (1999) that the topologies of potential and force-free models are in general not even qualitatively equivalent. In this paper, their example is re-examined and it is found that the opposite conclusions hold. In general, potential and force-free fields are topologically similar sufficiently close to localized sources. The exception to this are structurally unstable states, such as bifurcation states, where a small change of current can produce a significant change of topology.     

Nonlinear Force-Free Reconstruction of the Global Solar Magnetic Field: Methodology

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday’s equation, give rise to a respective normal-field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modeled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal-field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition – the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time that a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but it depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.     

A Fast Current-Field Iteration Method for Calculating Nonlinear Force-Free Fields

Existing methods for calculating nonlinear force-free magnetic fields are slow, and are likely to be inadequate for reconstructing coronal magnetic fields based on high-resolution vector magnetic field data from a new generation of spectro-polarimetric instruments. In this paper a new implementation of the current-field iteration method is presented, which is simple, fast, and accurate. The time taken by the method scales as N ⁴, for a three-dimensional grid with N ³ points. The method solves the field-updating part of the iteration by exploiting a three-dimensional Fast Fourier Transform solution of Ampere’s law with a current density field constructed to satisfy the required boundary conditions, and uses field line tracing to solve the current-updating part of the iteration. The method is demonstrated in application to a known nonlinear force-free field and to a bipolar test case.     

Preprocessing of Vector Magnetograph Data for a Nonlinear Force-Free Magnetic Field Reconstruction

Knowledge regarding the coronal magnetic field is important for the understanding of many phenomena, like flares and coronal mass ejections. Because of the low plasma beta in the solar corona, the coronal magnetic field is often assumed to be force-free and we use photospheric vector magnetograph data to extrapolate the magnetic field into the corona with the help of a nonlinear force-free optimization code. Unfortunately, the measurements of the photospheric magnetic field contain inconsistencies and noise. In particular, the transversal components (say B _x and B _y) of current vector magnetographs have their uncertainties. Furthermore, the magnetic field in the photosphere is not necessarily force free and often not consistent with the assumption of a force-free field above the magnetogram. We develop a preprocessing procedure to drive the observed non–force-free data towards suitable boundary conditions for a force-free extrapolation. As a result, we get a data set which is as close as possible to the measured data and consistent with the force-free assumption.     

A Magnetostatic Grad�CRubin Code for Coronal Magnetic Field Extrapolations

The coronal magnetic field cannot be directly observed, but, in principle, it can be reconstructed from the comparatively well observed photospheric magnetic field. A?popular approach uses a nonlinear force-free model. Non-magnetic forces at the photosphere are significant, meaning the photospheric data are inconsistent with the force-free model, and this causes problems with the modeling (De Rosa et?al., Astrophys. J. 696, 1780, 2009). In this paper we present a numerical implementation of the Grad?CRubin method for reconstructing the coronal magnetic field using a magnetostatic model. This model includes a pressure force and a non-zero magnetic Lorentz force. We demonstrate our implementation on a simple analytic test case and obtain the speed and numerical error scaling as a function of the grid size.     

Force-free magnetic fields in the solar atmosphere

Reliable measurements of the solar magnetic field are restricted to the level of the photosphere. For about half a century attempts have been made to calculate the field in the layers above the photosphere, i.e. in the chromosphere and in the corona, from the measured photospheric field. The procedure is known as magnetic field extrapolation. In the superphotospheric parts of active regions the magnetic field is approximately force-free, i.e. electric currents are aligned with the magnetic field. The practical application to solar active regions has been largely confined to constant-α or linear force-free fields, with a spatially constant ratio, α, between the electric current and the magnetic field. We review results obtained from extrapolations with constant-α force-free fields, in particular on magnetic topologies favourable for flares and on magnetic and current helicities. Presently, different methods are being developed to calculate non-constant-α or nonlinear force-free fields from photospheric vector magnetograms. We also briefly discuss these methods and present a comparison of a linear and a nonlinear force-free magnetic field extrapolation applied to the same photospheric boundary data. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Topology of Magnetic Field and Coronal Heating in Solar Active Regions

Force-free magnetic fields can be computed by making use of a new numerical technique, in which the fields are represented by a boundary integral equation based on a specific Green's function. Vector magnetic fields observed on the photospheric surface can be taken as the boundary conditions of this equation. In this numerical computation, the following two points are emphasized: (1) A new method for data reduction is proposed, for removing uncertainties in boundary data and determining the parameter in this Green's function, which is important for solving the boundary integral equation. In this method, the transverse components of the observed boundary field are calibrated with a linear force-free field model without changing their azimuth. (2) The computed 3-D fields satisfy the divergence-free and force-free conditions with high precision. The alignment of these field lines is mostly in agreement with structures in Hα and Yohkoh soft X-ray images. Since the boundary data are calibrated with a linear force-free field model, the computed 3-D magnetic field can be regarded as a quasi-linear force-free field approximation. The reconstruction of 3-D magnetic field in active region NOAA 7321 was taken as an example to quantitatively exhibit the capability of our new numerical technique.     

Topological Differences Between Force-Free Field Models

The potential and linear force-free field models for the magnetic field in the solar corona are often used in the analysis of flares. The field is calculated using boundary values measured in the low solar atmosphere. The topology of the field calculated using these models is then compared to the position of flare emissions. We demonstrate that the topology of the field according to each of these models, with the same boundary conditions in place, is not in general even qualitatively equivalent. An argument is given for a similar discrepancy between a linear force-free field solution and a nonlinear force-free field solution.     

How Should One Optimize Nonlinear Force-Free Coronal Magnetic Field Extrapolations from SDO/HMI Vector Magnetograms?

The Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) provides photospheric vector magnetograms with a high spatial and temporal resolution. Our intention is to model the coronal magnetic field above active regions with the help of a nonlinear force-free extrapolation code. Our code is based on an optimization principle and has been tested extensively with semianalytic and numeric equilibria and applied to vector magnetograms from Hinode and ground-based observations. Recently we implemented a new version which takes into account measurement errors in photospheric vector magnetograms. Photospheric field measurements are often affected by measurement errors and finite nonmagnetic forces inconsistent for use as a boundary for a force-free field in the corona. To deal with these uncertainties, we developed two improvements: i) preprocessing of the surface measurements to make them compatible with a force-free field, and ii) new code which keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which must be optimized for use with data from SDO/HMI. In this work we describe the corresponding analysis method and evaluate the force-free equilibria by how well force-freeness and solenoidal conditions are fulfilled, by the angle between magnetic field and electric current, and by comparing projections of magnetic field lines with coronal images from the Atmospheric Imaging Assembly (SDO/AIA). We also compute the available free magnetic energy and discuss the potential influence of control parameters.     

A Fluid Dynamics Approach for the Computation of Non—linear Force—Free Magnetic Field

Inspired by the analogy between the magnetic field and velocity field of incompressible fluid flow, we propose a fluid dynamics approach for computing nonlinear force-free magnetic fields. This method has the advantage that the divergence-free condition is automatically satisfied, which is a sticky issue for many other algorithms, and we can take advantage of modern high resolution algorithms to process the force-free magnetic field. Several tests have been made based on the well-known analytic solution proposed by Low & Lou. The numerical results are in satisfactory agreement with the analytic ones. It is suggested that the newly proposed method is promising in extrapolating the active region or the whole sun magnetic fields in the solar atmosphere based on the observed vector magnetic field on the photosphere.     

Optimization code with weighting function for the reconstruction of coronal magnetic fields

We developed a code for the reconstruction of nonlinear force-free and non-force-free coronal magnetic fields. The 3D magnetic field is computed numerically with the help of an optimization principle. The force-free and non-force-free codes are compiled in one program. The force-free approach needs photospheric vector magnetograms as input. The non-force-free code additionally requires the line-of-sight integrated coronal density distribution in combination with a tomographic inversion code. Previously the optimization approach has been used to compute magnetic fields using all six boundaries of a computational box. Here we extend this method and show how the coronal magnetic field can be reconstructed only from the bottom boundary, where the boundary conditions are measured with vector magnetographs. The program is planed for use within the Stereo mission.     

Development of force-free current in astrophysical circumstances

The physical conditions needed for the development of field-aligned force-free current in astrophysical circumstances are considered. It is shown that a large-scale differential motion of magnetic regions can lead to the development of magnetic field with the preferential enhancement of force-free current. Other physical consequences of force-free current in evolving magnetic field are also discussed.     

RECONSTRUCTION OF THE THREE-DIMENSIONAL CORONAL MAGNETIC FIELD

Studies of solar flares indicate that the mechanism of flares is magnetic in character and that the coronal magnetic field is a key to understanding solar high-energy phenomena. In our ongoing research we are conducting a systematic study of a large database of observations which includes both coronal structure (from the Soft X-ray Telescope on the Yohkoh spacecraft) and photospheric vector magnetic fields (from the Haleakala Stokes Polarimeter at Mees Solar Observatory). We compare the three-dimensional nonlinear force-free coronal magnetic field, computed from photospheric boundary data, to images of coronal structure. In this paper we outline our techniques and present results for active region AR 7220/7222. We show that the computed force-free coronal magnetic field agrees well with Yohkoh X-ray coronal loops, and we discuss the properties of the coronal magnetic field and the soft X-ray loops.     

Bifurcation of force-free solar magnetic fields: A numerical approach

Numerical calculations of two-dimensional force-free fields as models of solar active regions are presented. For a given toroidal component of the photospheric magnetic field two branches of solutions are numerically obtained which merge at the critical point of maximum allowed toroidal magnetic field. Depending on boundary conditions magnetic islands may or may not form. The results are discussed with respect to their relevance to the flare process.     

Singular lines of one-dimensional force-free magnetic field

The behavior of adiabatically slow deformations of the force-free field is investigated. Using the linear approximation it is shown that for a rather wide class of boundary perturbations of one-dimensional force-free field there appear singular magnetic force lines or surfaces. Hence the problem of quasi-steady deformation of frozen-in magnetic field has no solution. Relating to the problem of magnetic field in the solar corona it means that there will appear discontinuities (current sheets), when the magnetic field is deformed, for example, due to photospheric motion.     

The Free Energy of NOAA Solar Active Region AR 11029

The NOAA active region (AR) 11029 was a small but highly active sunspot region which produced 73 GOES soft X-ray flares during its transit of the disk in late October 2009. The flares appear to show a departure from the well-known power law frequency-size distribution. Specifically, too few GOES C-class and no M-class flares were observed by comparison with a power law distribution (Wheatland, Astrophys. J. 710, 1324, 2010). This was conjectured to be due to the region having insufficient magnetic energy to power the missing large events. We construct nonlinear force-free extrapolations of the coronal magnetic field of AR 11029 using data taken on 24 October by the SOLIS Vector SpectroMagnetograph (SOLIS/VSM) and data taken on 27 October by the Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). Force-free modeling with photospheric magnetogram data encounters problems, because the magnetogram data are inconsistent with a force-free model. We employ a recently developed “self-consistency” procedure which addresses this problem and accommodates uncertainties in the boundary data (Wheatland and Régnier, Astrophys. J. 700, L88, 2009). We calculate the total energy and free energy of the self-consistent solution, which provides a model for the coronal magnetic field of the active region. The free energy of the region was found to be ≈?4×10²⁹?erg on 24 October and ≈?7×10³¹?erg on 27 October. An order of magnitude scaling between RHESSI non-thermal energy and GOES peak X-ray flux is established from a sample of flares from the literature and is used to estimate flare energies from the observed GOES peak X-ray flux. Based on the scaling, we conclude that the estimated free energy of AR 11029 on 27 October when the flaring rate peaked was sufficient to power M-class or X-class flares; hence, the modeling does not appear to support the hypothesis that the absence of large flares is due to the region having limited energy.