Nonlinear forcefree modeling of the solar coronal magnetic field
Abstract
The coronal magnetic field is an important quantity because the magnetic field dominates the structure of the solar corona. Unfortunately direct measurements of coronal magnetic fields are usually not available. The photospheric magnetic field is measured routinely with vector magnetographs. These photospheric measurements are extrapolated into the solar corona. The extrapolated coronal magnetic field depends on assumptions regarding the coronal plasma, e.g. forcefreeness. Forcefree means that all nonmagnetic forces like pressure gradients and gravity are neglected. This approach is well justified in the solar corona due to the low plasma beta. One has to take care, however, about ambiguities, noise and nonmagnetic forces in the photosphere, where the magnetic field vector is measured. Here we review different numerical methods for a nonlinear forcefree coronal magnetic field extrapolation: GradRubin codes, upward integration method, MHDrelaxation, optimization and the boundary element approach. We briefly discuss the main features of the different methods and concentrate mainly on recently developed new codes.
WIEGELMANN \titlerunningheadCoronal magnetic fields. \authoraddrMaxPlanckInstitut für Sonnensystemforschung, MaxPlanckStrasse 2, 37191 KatlenburgLindau, Germany, (wiegelmann@mps.mpg.de)
Contents

Introduction

Why do we need nonlinear forcefree fields?


Nonlinear forcefree codes

GradRubin method

Upward integration method

MHD relaxation

Optimization approach

Boundary element or Greens function like method


How to deal with nonforcefree boundaries and noise ?

Consistency check of vector magnetograms

Preprocessing


Code testing and code comparisons

The NLFFFconsortium


Conclusions and Outlook
1 Introduction
Information regarding the coronal magnetic field are important for space weather
application like the onset of flares and coronal mass ejections (CMEs).
Unfortunately we usually cannot measure the coronal magnetic field directly,
although recently some progress has been made
(see e.g., Judge, 1998; Solanki et al., 2003; Lin et al., 2004). Due to the optically thin coronal
plasma direct measurements of the coronal magnetic field have a lineofsight
integrated character and to derive the accurate 3D structure of the coronal magnetic
field a vector tomographic inversion is required. Corresponding feasibility studies
based on coronal Zeeman and Hanle effect measurements have recently been done by
Kramar et al. (2006) and Kramar and Inhester (2006). These direct measurements are only
available for a few individual cases and usually one has to extrapolate the coronal
magnetic field from photospheric magnetic measurements. To do so, one has to make
assumptions regarding the coronal plasma. It is helpful that the low solar corona is
strongly dominated by the coronal magnetic field and the magnetic pressure is orders
of magnitudes higher than the plasma pressure. The quotient of plasma pressure
and magnetic pressure, is small compared to unity
. In lowest order nonmagnetic forces like pressure
gradient and gravity can
be neglected which leads to the forcefree assumption. Force free fields are
characterized by the equations:
(1)  
(2)  
(3) 
where is the magnetic field, the electric current density and the permeability of vacuum. Equation (1) implies that for forcefree fields the current density and the magnetic field are parallel, i.e.
(4) 
or by replacing with Eq. (2)
(5) 
where is called the forcefree function. To get some insights in the structure of the space dependent function , we take the divergence of Eq. (4) and make use of Eqs. (2) and (3):
(6) 
which tells us that the forcefree function is constant on every field line, but will usually change from one field line to another. This generic case is called nonlinear forcefree approach.
Popular simplifications are (current free potential fields, see e.g., Schmidt (1964); Semel (1967); Schatten et al. (1969); Sakurai (1982)) and (linear forcefree approach, see e.g., Nakagawa and Raadu (1972); Chiu and Hilton (1977); Seehafer (1978); Alissandrakis (1981); Seehafer (1982); Semel (1988)). These simplified models have been in particular popular due to their relative mathematical simplicity and because only lineofsight photospheric magnetic field measurements are required. Linear forcefree fields still contain one free global parameter , which can be derived by comparing coronal images with projections of magnetic field lines (e.g., Carcedo et al. (2003)). It is also possible to derive an averaged value of from transverse photospheric magnetic field measurements (e.g. Pevtsov et al. (1994); Wheatland (1999); Leka and Skumanich (1999)). Despite the popularity and frequent use of these simplified models in the past, there are several limitations in these models (see below) which ask for considering the more sophisticated nonlinear forcefree approach.
Our aim is to review recent developments of the extrapolation of nonlinear forcefree fields (NLFFF). For earlier reviews on forcefree fields we refer to (Sakurai, 1989; Aly, 1989; Amari et al., 1997; McClymont et al., 1997) and chapter 5 of Aschwanden (2005). Here we will concentrate mainly on new developments which took place after these earlier reviews. Our main emphasis is to study methods which extrapolate the coronal magnetic field from photospheric vector magnetograms. Several vector magnetographs are currently operating or planed for the nearest future, e.g., ground based: the solar flare telescope/NAOJ (Sakurai et al., 1995), the imaging vector magnetograph/MEES Observatory (Mickey et al., 1996), Big Bear Solar Observatory, Infrared Polarimeter VTT, SOLIS/NSO (Henney et al., 2006) and space born: Hinode/SOT (Shimizu, 2004), SDO/HMI (Borrero et al., 2006). Measurements from these vector magnetograms will provide us eventually with the magnetic field vector on the photosphere, say for the normal and and for the transverse field. Deriving these quantities from the measurements is an involved physical process, which includes measurements based on the Zeeman and Hanle effect, the inversion of Stokes profiles (e.g., LaBonte et al. (1999)) and removing the ambiguity (e.g., Metcalf (1994); Metcalf et al. (2006)) of the horizontal magnetic field component. Special care has to be taken for vector magnetograph measurements which are not close to the solar disk, when the lineofsight and normal magnetic field component are far apart (e.g. Gary and Hagyard (1990)). For the purpose of this paper we do not address the observational methods and recent developments and problems related to deriving the photospheric magnetic field vector. We rather will concentrate on how to use the photospheric and to derive the coronal magnetic field.
The transverse photospheric magnetic field can be used to approximate the normal electric current distribution by
(7) 
and from this one gets the distribution of on the photosphere by
(8) 
By using Eq. (8) one has to keep in mind that rather large uncertainties in the transverse field component and numerical derivations used in (7) can cumulate in significant errors for the current density. The problem becomes even more severe by using (8) to compute in regions with a low normal magnetic field strength . Special care has to be taken at photospheric polarity inversion lines, i.e. lines along which (see e.g., Cuperman et al., 1991). The nonlinear forcefree coronal magnetic field extrapolation is a boundary value problem. As we will see later some of the NLFFFcodes make use of (8) to specify the boundary conditions while other methods use the photospheric magnetic field vector more directly to extrapolate the field into the corona.
Pure mathematical investigations of the nonlinear forcefree equations (see e.g. Aly, 1984; Boulmezaoud and Amari, 2000; Aly, 2005) and modelling approaches not based on vector magnetograms are important and occasionally mentioned in this paper. A detailed review of these topics is well outside the scope of this paper, however. Some of the modelapproaches not based on vector magnetograms are occasionally used to test the nonlinear forcefree extrapolation codes described here.
1.1 Why do we need nonlinear forcefree fields?

A comparison of global potential magnetic field models with TRACEimages by Schrijver et al. (2005) revealed that significant nonpotentially occurs regulary in active regions, in particular when new flux has emerged in or close to the regions.

Usually changes in space, even inside one active region. This can be seen, if we try to fit for the optimal linear forcefree parameter by comparing field lines with coronal plasma structures. An example can be seen in Wiegelmann and Neukirch (2002) where stereoscopic reconstructed loops by Aschwanden et al. (1999) have been compared with a linear forcefree field model. The optimal value of changes even sign within the investigated active regions, which is a contradiction to the linear forcefree approach (see Fig. 1).

Potential and linear forcefree fields are too simple to estimate the free magnetic energy and magnetic topology accurately. The magnetic energy of linear forcefree fields is unbounded in a halfspace (Seehafer, 1978) which makes this approach unsuitable for energy approximations of the coronal magnetic field. Potential fields have a minimum energy for an observed lineofsight photospheric magnetic field. An estimate of the excess of energy a configuration has above that of a potential field is an important quantity which might help to understand the onset of flares and coronal mass ejections.

A direct comparison of measured fields in a newly developed active region by Solanki et al. (2003) with extrapolations from the photosphere with a potential, linear and nonlinear forcefree model by Wiegelmann et al. (2005b) showed that nonlinear fields are more accurate than simpler models. Fig. 2 shows some selected magnetic field lines for the original measured field and extrapolations from the photosphere with the help of a potential, linear and nonlinear forcefree model.
These points tell us that nonlinear forcefree modelling is required for an accurate reconstruction of the coronal magnetic field. Simpler models have been used frequently in the past. Global potential fields provide some information of the coronal magnetic field structure already, e.g. the location of coronal holes. The generic case of forcefree coronal magnetic field models are nonlinear forcefree fields, however. Under generic we understand that can (and usually will) change in space, but this approach also includes the special cases and . Some active regions just happen to be more potential (or linear forcefree) and if this is the case they can be described with simpler models. Linear forcefree models might provide a rough estimate of the true 3D magnetic field structure if the nonlinearity is weak. The use of simpler models was often justified due to limited observational data, in particular if only the lineofsight photospheric magnetic field has been measured.
While the assumption of nonlinear forcefree fields is well accepted for the coronal magnetic fields in active regions, this is not true for the photosphere. The photospheric plasma is a finite plasma and nonmagnetic forces like pressure gradient and gravity cannot be neglected here. As a result electric currents have a component perpendicular to the magnetic field, which contradicts the forcefree assumption. We will discuss later, how these difficulties can be overcome.
2 Nonlinear forcefree codes
Different methods have been proposed to extrapolate nonlinear forcefree fields from photospheric vector magnetic field measurements.
2.1 GradRubin method
The GradRubin method reformulates the nonlinear forcefree equations in such a way, that one has to solve a well posed boundary value problem. This makes this approach also interesting for a mathematical investigation of the structure of the nonlinear forcefree equations. Bineau (1972) demonstrated that the used boundary conditions (vertical magnetic field on the photosphere and distribution at one polarity) ensure, at least for small values of and weak nonlinearities the existence of a unique nonlinear forcefree solution. A detailed analysis of the mathematical problem of existence and uniqueness of nonlinear forcefree fields is outside the scope of this review and can be found e.g. in Amari et al. (1997, 2006).
The method first computes a potential field, which can be obtained from the observed lineofsight photospheric magnetic field (say in Cartesian geometry) by different methods, e.g., a Greens function method as described in Aly (1989). It is also popular to use linear forcefree solvers, e.g. as implemented by Seehafer (1978); Alissandrakis (1981) with the linear forcefree parameter to compute the initial potential field. The transverse component of the measured magnetic field is then used to compute the distribution of on the photosphere by Eq. (8). While is described this way on the entire photosphere, for both polarities, a well posed boundary value problem requires that the distribution becomes only described for one polarity. The basic idea is to iteratively calculate for a given field from (6), then calculate the current via (4) and finally update from the BiotSavart problem (5). These processes are repeated until the full current as prescribed by the distribution has been injected into the magnetic field and the updated magnetic field configuration becomes stationary in the sense that eventually the recalculation of the magnetic field with Amperes law does not change the configuration anymore. To our knowledge the GradRubin approach has been first implemented by Sakurai (1981). has been prescribed on several nodal points along a number of magnetic field lines of the initial potential field. The method used a finiteelementlike discretization of current tubes associated with magnetic field lines. Each current tube was divided into elementary current tubes of cylindrical shape. The magnetic field is updated with Ampere’s law using a superposition of the elementary current tubes. The method was limited by the number of currentcarrying field lines, nodalpoints and the corresponding number of nonlinear equations to solve with the available computer resources more than a quarter century ago.
Computer resources have increased rapidly since the first NLFFFimplementation by Sakurai (1981) and about a decade ago the GradRubin method has been implemented on a finite difference grid by Amari et al. (1997, 1999). This approach decomposes the equations (13) into a hyperbolic part for evolving along the magnetic field lines and an elliptic one to iterate the updated magnetic field from Amperes law. For every iteration step one has to solve iteratively for:
(9)  
(10) 
which evolves in the volume and
(12)  
(13)  
(14)  
(15) 
where corresponds to the photospheric distribution of for either on the positive or the negative polarity. The GradRubin method as described in Amari et al. (1997, 1999) has been applied to investigate particular active regions in Bleybel et al. (2002) and a comparison of the extrapolated field with 2D projections of plasma structures as seen in H, EUV and Xray has been done in Régnier et al. (2002); Régnier and Amari (2004). The code has also been used to investigate mutual and self helicity in active regions by Régnier et al. (2005) and to flaring active regions by Régnier and Canfield (2006).
A similar approach as done by Sakurai (1981) has been implemented by Wheatland (2004). The implemented method computes the magnetic field directly on the numerical grid from Ampere’s law. This is somewhat simpler and faster as Sakurai’s approach which required solving a large system of nonlinear equations for this aim. The implementation by Wheatland (2004) has, in particular, been developed with the aim of parallelization. The parallelization approach seems to be effective due to a limited number of interprocess communications. This is possible because as the result of the linearity of Ampere’s law the contributions of the different current carrying field lines are basically independent from each other. In the original paper Wheatland reported problems for large currents on the field lines. These problems have been related to an error in current representation of the code and the corrected code worked significantly better, see also Schrijver et al. (2006). The method has been further developed in Wheatland (2006). This newest Wheatlandimplementation scales with the number of grid points for a volume, rather than for the earlier Wheatland (2004) implementation. The main new development is a faster implementation of the currentfield iteration. To do so the magnetic field has been separated into a currentfree and a current carrying part at each iteration step. Both parts are solved using a discrete Fast Fourier Transformation, which imposes the required boundary conditions implicitly. The code has been parallelized on shared memory distributions with OpenMP.
Amari et al. (2006) developed two new versions of their GradRubin code. The first version is a finite difference method and the code was called ’XTRAPOL’. This code prescribe the coronal magnetic field with the help of a vector potential . The code has obviously it’s heritage from the earlier implementation of Amari et al. (1999), but with several remarkable differences:

The code includes a divergence cleaning routine, which takes care about . The condition is fulfilled with high accuracy in the new code compared to in the earlier implementation.

The lateral and top boundaries are more flexible compared to the earlier implementation and allow a finite and non zero values for one polarity on all boundaries. This treats the whole boundary (all six faces) as a whole.

The slow current input as reported for the earlier implementation, which lead to a two level iteration, has been replaced. Now the whole current is injected at once and only the inner iteration loop of the earlier code remained in the new version.

The computation of the characteristics has been improved with an adaptive AdamsBashforth integration scheme (see Press, 2002, for details).

The fixed number of iteration loops have been replaced by a quantitative convergence criterium.
In the same paper Amari et al. (2006) introduced another GradRubin approach based on finite elements, which they called ’FEMQ’. Different from alternative implementation this code does not use a vector potential but iterates the coupled divergence and curl system, which is solved with the help of a finite element discretization. The method transforms the nonlinear forcefree equations into a global linear algebraic system.
Inhester and Wiegelmann (2006) implemented a GradRubin code on a finite element grid with staggered field components (see Yee (1966)) which uses discrete Whitney forms (Bossavit (1988)). Whitney forms allow to transform standard vector analysis (as the differential operators gradient, curl and divergence) consistently into the discrete space used for numerical computations. Whitney forms contain four types of finite elements (form 03). They can be considered as a discrete approximation of differential forms. The finite element base may consist of polynominals of any order. In its simplest form, the 0forms have as parameters the function values at the vertices of the cells and are linearly interpolated within each cell. 1forms are a discrete representation of a vector field defined on the cell edges. 2forms are defined as the field component normal to the surfaces of the cells. The 3forms are finite volume elements for a scalar function approximation, which represents the average of a scalar over the entire cell. The four forms are related to each other by GRAD (0 to 1 form), CURL (1 to 2 form) and DIV (2 to 3 form). As for continous differential forms, double differentiation (CURL of GRAD, DIV of CURL) give exactly zero, independent of the numerical precision. A dual grid, shifted by half a grid size in each axis, was introduced in order to allow for Laplacians. Whitney forms on the dual grid are related to forms on the primary grid in a consistent way.
The GradRubin implementation uses a vector potential representation of the magnetic field, where the vector potential is updated with a Poisson equation in each iteration step. The Poisson equation is effectively solved with the help of a multigrid solver. The main computing time is spend to distribute along the field lines with (6). This seems to be a general property of GradRubin implementations. One can estimate the scaling of (6) by , where the number of field lines to compute is and the length of a field line . The BiotSavart step (5) solved with FFT or multigrid methods scales only with . Empirical tests show that the number of iteration steps until a stationary state is reached does not depend on the number of grid points for GradRubin solvers. We have explained before, that the GradRubin implementation requires the prescription of only for one polarity to have a well posed mathematical problem. The Inhester and Wiegelmann (2006) implementation allows these choice of boundary conditions as a special case. In general one does not need to make the distinction between and in the new implementation. A well posed mathematical problem is still ensured, however, in the following way. Each boundary value of is attached with a weight. The final version of on each field line is then determined by a weighted average of the values on both endpoints of a field lines. By this way the influence of uncertain boundary values, e.g. on the side walls and imprecise photospheric measurements can be suppressed.
2.2 Upward integration method
The basic equations for the upward integration method (or progressive extension method) have been published already by Nakagawa (1974) and a corresponding code has been developed by Wu et al. (1985, 1990a). The upward integration method is a straight forward approach to use the nonlinear forcefree equations directly to extrapolate the photospheric magnetic field into the corona. To do so one reformulates the forcefree equations (13) in order to extrapolate the measured photospheric magnetic field vector into the solar corona.
As a first step the magnetic field vector on the lower boundary is used to compute the component of the electric current with Eq. (7) and the photospheric distribution (say ) by Eq. (8). With the help of Eq. (4) we calculate the and component of the current density
(16)  
(17) 
We now use Eq. (3) and the and component of Eq. (2) to obtain expressions for the derivatives of all three magnetic field components in the form
(18)  
(19)  
(20) 
The idea is to integrate this set of equations numerically upwards in by repeating the previous steps at each height. As a result we get in principle the 3D magnetic field vector in the corona. While this approach is straight forward, easy to implement and computational fast (no iteration is required), a serious drawback is that it is unstable. Several authors (e.g., Cuperman et al., 1990; Amari et al., 1997) pointed out that the formulation of the forcefree equations in this way is unstable because it is based on an illposed mathematical problem. In particular one finds that exponential growth of the magnetic field with increasing height is a typical behaviour. What makes this boundary value problem illposed is that the solution does not depend continuously on the boundary data. Small changes or inaccuracies in the measured boundary data lead to a divergent extrapolated field (see Low and Lou, 1990, for a more detailed discussion). As pointed out by Low and Lou meaningful boundary conditions are required also on the outer boundaries of the computational domain. It is also possible to prescribe open boundaries in the sense that the magnetic field vanishes at infinity. This causes an additional problem for the upward integration method, because the method transports information only from the photosphere upwards and does not incorporate boundary information on other boundaries or at infinity. Attempts have been made to regularize the method (e.g., Cuperman et al., 1991; Demoulin and Priest, 1992), but cannot be considered as fully successful.
Wu et al. (1990b) compared the GradRubin method in the implementation of Sakurai (1981) with the upward integration method in the implementation of Wu et al. (1990a) ^{1}^{1}1The authors used a somewhat different nomenclature: The upward integration method was called ’progressive extension method’ and the GradRubin method ’iterative method’. That time the term ’iterative method’ was reasonable because GradRubin was the only iterative approach available, but now, 17 years later, several other iterative methods are available to compute nonlinear forcefree fields.. The comparison showed qualitatively similar results for extrapolations from an observed magnetogram, but quantitatively differences. The NLFFFcomputations have been very similar to potential field extrapolations, however, too. One reason for this behaviour was, that the method of Sakurai (1981) is limited to small values of and an ’by eye’ comparison shows that the corresponding NLFFF field is very close to a potential field configuration. The field computed with the upward integration method deteriorated if the height of the extrapolation exceeded a typical horizontal scale length.
The upward integration method has been recently reexamined by Song et al. (2006) who developed a new formulation of this approach. The new implementation uses smooth continuous functions and the equations are solved in asymptotic manner iteratively. The original upward integration equations are reformulated into a set of ordinary differential equations and uniqueness of the solution seems to be guarantied at least locally. While Demoulin and Priest (1992) stated that ’no further improvement has been obtained with other types of smoothing functions’ the authors of Song et al. (2006) point out that the transformation of the original partial differential equations into ordinary ones eliminates the growing modes in the upward integration method, which have been reported before in Wu et al. (1990a) and subsequent papers. The problem that all three components of the photospheric magnetic field and the photospheric distribution has to be prescribed in a consistent way remains in principle, but some compatibility conditions to compute a slowly varying have been provided by Song et al. (2006). These compatibility conditions are slightly different for real photospheric observations and tests with smooth boundaries extracted from semianalytic equilibria. For the latter kind of problems the new formulation provided reasonable results with the standard test equilibrium found by Low and Lou (1990). The method seems to be also reasonable fast. Of course further tests with more sophisticated equilibria and real data are necessary to evaluate this approach in more detail.
2.3 MHD relaxation
MHD relaxation codes (e.g., Chodura and Schlueter, 1981) can be applied to solve nonlinear forcefree fields as well. The idea is to start with a suitable magnetic field which is not in equilibrium and to relax it into a forcefree state. This is done by using the MHD equations in the following form:
(21)  
(22)  
(23)  
(24) 
where is a viscosity and the electric field. As the MHDrelaxation aims for a quasi physical temporal evolution of the magnetic field from a nonequilibrium towards a (nonlinear forcefree) equilibrium this method is also called ’evolutionary method’ or ’magnetofrictional method’. The basic idea is that the velocity field in the equation of motion (21) is reduced during the relaxation process. Ideal Ohm’s law (22) ensures that the magnetic connectivity remains unchanged during the relaxation. The artificial viscosity plays the role of a relaxation coefficient which can be chosen in such way that it accelerates the approach to the equilibrium state. A typical choice is
(25) 
with constant. Combining Eqs. (21), (22), (23) and (25) we get an equation for the evolution of the magnetic field during the relaxation process:
(26) 
with
(27) 
This equation is then solved numerically starting with a given initial condition for , usually a potential field. Equation (26) ensures that Eq. (24) is satisfied during the relaxation if the initial magnetic field satisfies it. ^{2}^{2}2As we will see below the ’optimization’ approach leads to a similar iteration equations for the magnetic field, but a different artificial driving force . The difficulty with this method is that it cannot be guaranteed that for given boundary conditions and initial magnetic field (i.e. given connectivity), a smooth forcefree equilibrium exists to which the system can relax. If such a smooth equilibrium does not exist the formation of current sheets is to be expected which will lead to numerical difficulties. Therefore, care has to be taken when choosing an initial magnetic field.
Yang et al. (1986) developed a magneto frictional method which represent the magnetic field with the help of Euler (or Clebsch) potentials.
(28) 
where the potentials and are scalar functions. The general method has been developed for three dimensional fields and iterative equations for and have been derived. The Clebsch representation automatically ensures . The method has been explicitly tested in the paper by Yang et al. (1986) with the help of an equilibrium with one invariant coordinate. In principle it should be possible to use this representation for the extrapolation of nonlinear forcefree fields, but we are not aware of a corresponding implementation. Due to the discussion in Yang et al. (1986) a difficulty seems to be that one needs to specify boundary conditions for the potentials, rather than for the magnetic fields. It seems in particular to be difficult to find boundaries conditions for potentials which correspond to the transverse component of the photospheric magnetic field vector. One problem is that boundary conditions for and prescribe the connectivity. Every field line can be labelled by its values. Hence boundary values for and establish foot point relations although the field is not known yet.
The MHDrelaxation (or evolutionary) method has been implemented by Mikic and McClymont (1994); McClymont et al. (1997) based on the time dependent MHDcode by Mikic et al. (1988). The code uses a nonuniform mesh and the region of interested is embedded in a large computational domain to reduce the influence of the lateral boundaries. The method has been applied to extrapolate the magnetic field above an active region by Jiao et al. (1997). The computations have been carried out with a resolution of the order of points. A supercomputer was required for these computations that time (10 years ago), but due to the rapid increase of computer speed and memory within the last decade this restriction is very probably not valid anymore.
Roumeliotis (1996) developed the so called stress and relax method. In this approach the initial potential field becomes disturbed by the observed transverse field component on the photosphere. The boundary conditions are replaced in subsequently in several small steps and always relaxed with a similar MHDrelaxation scheme as described above towards a forcefree equilibrium. The code by Roumeliotis (1996) has implemented a function which allows to give a lower weight to regions where the transverse photospheric field has been measured with lower accuracy. Additional to the iterative equations as discussed above, the method includes a resistivity (or diffusivity) by adding a term on the right hand site of Ohms law (22). This relaxes somewhat the topological constrains of ideal MHD relaxation, because a finite resistivity allows a kind of artificial reconnection and corresponding changes of the initial potential field topology. The method has been tested with a forcefree equilibrium found by Klimchuk and Sturrock (1992) and applied to an active region measured with the MSFC vectormagnetograph.
The stress and relax method has been revisited by Valori et al. (2005). Different from the earlier implementation by Roumeliotis (1996) the new implementation uses directly the magnetic field, rather than the vector potential in order to keep errors from taking numerical deviations from noisy magnetograms minimal. The solenoidal condition is controlled by a diffusive approach by Dedner et al. (2002) which removes effectively a numerically created finite divergence of the relaxed magnetic field. The new implementation uses a single stress step, rather than the multiple small stress used by Roumeliotis (1996) to speed up the computation.
The single step stress and relax method is connected with a suitable control of artificial plasma flows by the Courant criterium. The authors reported that a multistep and singlestep implementation do not reveal significant differences. The numerical implementation is based on the timedepended full MHDcode ’AMRVAC’ by Keppens et al. (2003). Valori et al. (2005) tested their nonlinear forcefree implementation with a numerically constructed nonlinear forcefree twisted loop computed by Török and Kliem (2003).
2.4 Optimization approach
The optimization approach has been developed in Wheatland et al. (2000). The solution is found by minimizing the functional
(29) 
Obviously, is bound from below by . This bound is attained if the magnetic field satisfies the forcefree equations (1)(3).
By taking the functional derivatives with respect to some iteration parameter we get:
(30) 
with
(31)  
(32) 
The surface term vanishes if the magnetic field vector is kept constant on the surface, e.g. prescribed from photospheric measurements. In this case decreases monotonically if the magnetic field is iterated by
(33) 
Let us remark that as defined in Equation (27) and used for MHDrelaxation is identical with the first term on the righthandside of Eq. (31), but Eq. (31) contains additional terms.
For this method the vector field is not necessarily solenoidal during the computation, but will be divergencefree if the optimal state with is reached. A disadvantage of the method is that it cannot be guaranteed that this optimal state is indeed reached for a given initial field and boundary conditions. If this is not the case then the resulting will either be not forcefree or not solenoidal or both.
McTiernan has implemented the optimization approach basically as described in Wheatland et al. (2000) in IDL (see Schrijver et al. (2006) for a brief description of the McTiernan implementation.) This code allows the use of a nonuniform computational grid. In a code intercomparison by Schrijver et al. (2006) the IDL optimization code by McTiernan was about a factor of slower compared to an implementation in parallelized C by Wiegelmann (2004). To our knowledge McTiernan has translated his IDLcode into FORTRAN in the meantime for faster computation (personal communication on the NLFFFworkshop PaloAlto, june 2006. See also Metcalf et al. (2007).)
Several tests have been performed with the optimization approach in Wiegelmann and Neukirch (2003). It has been investigated how the unknown lateral and top boundary influence the solution. The original optimization approach by Wheatland et al. (2000) has been extended towards more flexible boundaryconditions, which allow on the lateral and top boundaries. This has been made with the help of the surface integral term in (30) and led to an additional term on the boundaries. This approached improved the performance of the code for cases, where only the bottom boundary was prescribed. No improvement was found for a slow multistep replacement of the boundary and this possibility has been abandoned in favour of a single step method.
It has been also investigated how noise influences the optimization code and this study revealed that noise in the vector magnetograms leads to less accurate nonlinear forcefree fields.
Wiegelmann (2004) has reformulated the optimization principle by introducing weighting functions One defines the functional
(34) 
where is a weighting function. It is obvious that (for ) the forcefree equations (13) are fulfilled when L is equal zero. Minimization of the functional (34) lead to:
(35) 
(36)  
(37)  
(38) 
with as defined in (31). With this approach reduces to the Wheatland et al. (2000) method as described above. The weighting function is useful if only the bottom boundary data are known. In this case we a buffer boundary of several grid points towards the lateral and top boundary of the computational box is introduced. The weighting function is chosen constant in the inner, physical domain and drop to with a cosine profile in the buffer boundary towards the lateral and top boundary of the computational box. In Schrijver et al. (2006) some tests have been made with different weighting functions for the forcefree and solenoidal part of the functional , but the best results have been obtained if both terms got the same weight. The computational implementation involves the following steps.

Compute start equilibrium (e.g. a potential field) in the computational box.

Replace the bottom boundary with the vector magnetogramm.

Minimize the functional (34) with the help of Eq. (35). The continuous form of (35) guaranties a monotonically decreasing . This is as well ensured in the discretized form if the iteration step is sufficiently small. The code checks if after each time step. If the condition is not fulfilled, the iteration step is repeated with reduced by a factor of 2. After each successful iteration step we increase slowly by a factor of to allow the time step to become as large as possible with respect to the stability condition.

The iteration stops if becomes stationary. Stationarity is assumed if for consecutive iteration steps.
The program has been tested with the semianalytic nonlinear forcefree configuration by Low and Lou (1990) and Titov and Démoulin (1999) in Wiegelmann et al. (2006a). The code has been applied to extrapolate the coronal magnetic field in active regions in Wiegelmann et al. (2005b, a).
A finite element optimization approach has been implemented by Inhester and Wiegelmann (2006) using the Whitney elements as for the GradRubin code (which has been described above). The optimization method uses exactly the same staggered finite element grid as described above, which is different from the finite difference grids used in the earlier implementations by Wheatland et al. (2000); Wiegelmann and Neukirch (2003); Wiegelmann (2004). Another difference is that earlier implementations discretized the analytical derivative of the functional (31), while the new code takes the numerical more consistent derivative of the discretized function . All other implementations used a simple Landweber scheme for updating the magnetic field, which is replaced here by an unpreconditioned conjugate gradient iteration, which at every time step performs an exact line search to the minimum of in the current search direction and additional selects an improve search direction instead of the gradient of the functional . To do so the Hessian matrix of the functional is computed during every iteration step. An effective computation of the Hessian matrix is possible, because the reformulated function is a fourth order polynomial in and all five polynomial coefficients can be computed in one go. The code has been tested with Low and Lou (1990) and the result of twisted loop computations of the GradRubin implementation on the same grid.
The optimization code in the implementation of Wiegelmann (2004) has recently be extended towards using a multiscale implementation. The main difference from the original code are (see also Metcalf et al. (2007)):

The method is not full multigrid, but computes the solution on different grids only once, e.g., something like .

The main idea is to get a better (than potential field) start equilibrium on the full resolution box.

Solution of smaller grids are interpolated onto larger grids as initial state for the magnetic field in the computational domain of the next larger box.
The multiscale implementation has been tested as part of a codeintercomparison test in Metcalf et al. (2007) with the help of solarlike reference model computed by van Ballegooijen (2004); van Ballegooijen et al. (2007).
The optimization approach has recently be implemented in spherical geometry by Wiegelmann (2007) and tested with Low and Lou (1990). The original longitudinal symmetric Low and Lou solution has been shifted by of a solar radius to test the code without any symmetry with respect to the Suns surface. The numerical implementation is very similar as the Cartesian implementation described in Wiegelmann (2004). The spherical implementation converged fast for low latitude regions, but the computing time increased significantly if polar regions have been included. It has been suggested to implement the code on a so called ’Yin and Yang’ grid as developed by Kageyama and Sato (2004) to reduce the computing time. The ’Yin and Yang’ grid is suitable for massive parallelization, which is necessary for full sphere high resolution NLFFFcomputations.
2.5 Boundary element or Greens function like method
The boundary integral method has been developed by Yan and Sakurai (2000). The method relates the measured boundary values with the nonlinear forcefree field in the entire volume by:
(39) 
where for points in the volume and for boundary points and is the measured vector magnetic field on the photosphere. The auxiliary vector function is defined as
(40) 
and the are computed in the original approach by Yan and Sakurai (2000) with integrals over the whole volume, which define the implicitly:
(41) 
This volume integration, which has to be carried out for every point in the volume is certainly very time consuming (a sixth order process). The have the same dimension as the magnetic field. The existence of the has been confirmed for the semianalytic field of Low and Lou (1990) by Li et al. (2004). While the work of Li et al. (2004) showed that one can find the auxiliary function for a given forcefree field in 3D, the difficulty is that is apriori unknown if only the photospheric magnetic field vector is given. Yan and Sakurai (2000) proposed an iterative scheme to compute the auxiliary functions and the nonlinear forcefree magnetic field selfconsistently. They use the approximate solution on the right hand side of Eq. (39) to compute a better solution by:
(42) 
where the initial guess for the magnetic field in the volume is and also the initial . In principle it would be also possible to compute a potential field first and derive the auxiliary functions for this field as done in Li et al. (2004) and iterate subsequently for the nonlinear forcefree fields and the associated auxiliary functions with Eq. (42). This possibility has not been tried out to our knowledge until now, however. The method iterates the magnetic field until and converge. In an inter code comparison by Schrijver et al. (2006) one iteration step of (42) took about for this method and only this one step was carried out without further iteration. This seems, however, not to be sufficient to derive an accurate nonlinear forcefree solution. The method has been applied for the comparison with soft Xray loops observed with YOHKOH by Wang et al. (2000); Liu et al. (2002) and to model a magnetic flux robe by Yan et al. (2001a, b).
In a new implementation of the boundary element method by Yan and Li (2006) the auxiliary functions are computed iteratively with the help of a simplex method. This avoids the numerical expensive computation of the volume integral (41). The boundary element method is still rather slow if a magnetic field has to be computed in an entire 3D domain. Different from other method, it allows, however, to evaluate the NLFFFfield at every arbitrary point within the domain from the boundary data, without the requirement to compute the field in an entire domain. This is in particular useful if one is interested to compute the NLFFFfield only along a given loop.
He and Wang (2006) investigated the validity of the boundary integral representation for a spherical implementation. The method has been tested with the longitudinal invariant Low and Lou (1990) solution. The spherical implementation method of this method revealed reasonable results for smooth modestly nonlinear fields, but a poor convergence for complex magnetic field structures and large values of .
3 How to deal with nonforcefree boundaries and noise ?
Given arbitrary boundary conditions of the magnetic field vector on the photosphere, the solution to the forcefree equations in 3D may not exist. Nonlinear forcefree coronal magnetic field models assume, however, that the solution exists. It is certainly possible and necessary to check after or during the computation if a solution has been found. In the following we will discuss what we can do if the measured photospheric data are incompatible with the assumption of a forcefree coronal magnetic field.
3.1 Consistency check of vector magnetograms
We reexamine some necessary conditions with the photospheric field (or bottom boundary of a computational box). These conditions have to be fulfilled in order to be suitable boundary conditions for a nonlinear forcefree coronal magnetic field extrapolation. An apriori assumption about the photospheric data is that the magnetic flux from the photosphere is sufficiently distant from the lateral boundaries of the observational domain and the net flux is in balance, i.e.,
(43) 
Molodensky (1969, 1974); Aly (1989); Sakurai (1989) used the virial theorem to define which conditions a vector magnetogram has to fulfill to be consistent with the assumption of a forcefree field in the corona above the boundary. These conditions are:

The total force on the boundary vanishes
(44) (45) 
The total torque on the boundary vanishes
(46) (47) (48)
In an earlier review Aly (1989) has mentioned already that the magnetic field is probably not forcefree in the photosphere, where is measured because the plasma in the photosphere is of the order of one and pressure and gravity forces are not negligible. The integral relations (44)(48) are not satisfied in this case in the photosphere and the measured photospheric field is not a suitable boundary condition for a forcefree extrapolation. Metcalf et al. (1995) concluded that the solar magnetic field is not forcefree in the photosphere, but becomes forcefree only at about above the photosphere. Gary (2001) pointed out that care has to be taken when extrapolating the coronal magnetic field as a forcefree field from photospheric measurements, because the forcefree low corona is sandwiched between two regions (photosphere and higher corona) with a plasma , where the forcefree assumption might break down. An additional problem is that measurements of the photospheric magnetic vector field contain inconsistencies and noise. In particular the transverse components (say and ) of current vector magnetographs include uncertainties.
The forcefree field in a domain requires the Maxwell stress (44)(48) to sum to zero over the boundary. If these conditions are not fulfilled a forcefree field cannot be found in the volume. A faithful algorithm should therefore have the capability of rejecting a prescription of the vector field at the boundary that fails to produce zero net Maxwell stress. A simple way to incorporate these conditions would be to evaluate the integrals (44)(48) within or prior to the NLFFFcomputation and to refuse the vector field if the conditions are not fulfilled with sufficient accuracy. Current codes do run, however, although if feeded with inconsistent boundary data, but they certainly cannot find a forcefree solution in this case (because it does not exist). This property of current codes does, however, not challenge the trustworthiness of the algorithms, because the forcefree and solenoidal conditions are checked in 3D, e.g., with the help of the functional as defined in (29). A non zero value of (within numerical accuracy) tells the user that a forcefree state has not been reached. In principle it would be possible that the codes do refuse to output the magnetic field in this case. For current codes this is not automatically controlled but responsibility of the user.
Unfortunately current measurements of the magnetic field vector are only available routinely in the photosphere, where we have a finite plasma and nonmagnetic forces might become important. The forcefree compatibility conditions (44)(48) are not fulfilled in the photosphere, but they should be fulfilled in the low chromospheric and coronal plasma above. The question is if we still can use the photospheric measurements to find suitable consistent boundary conditions for a nonlinear forcefree modelling. Such an approach has been called preprocessing of vector magnetograms.
3.2 Preprocessing
The preprocessing routine has been developed by Wiegelmann et al. (2006b). The integral relations (44)(48) have been used to define a 2D functional of quadratic forms:
(49) 
where
(50)  
(51)  
(52)  
(53) 
The surface integrals are here replaced by a summation over all grid nodes of the bottom surface grid and the differentiation in the smoothing term is achieved by the usual 5point stencil for the 2DLaplace operator. Each constraint is weighted by a yet undetermined factor . The first term (=1) corresponds to the forcebalance conditions (44)(45), the next (=2) to the torquefree condition (46)(48). The following term (=3) ensures that the optimized boundary condition agrees with the measured photospheric data and the last terms (=4) controls the smoothing. The 2DLaplace operator is designated by . The aim of the preprocessing procedure is to minimize so that all terms if possible are made small simultaneously. A strategy on how to find the optimal yet undefined parameters is described in Wiegelmann et al. (2006b). As result of the preprocessing we get a data set which is consistent with the assumption of a forcefree magnetic field in the corona but also as close as possible to the measured data within the noise level.
4 Code testing and code comparisons
Newly developed codes for the extrapolation of nonlinear forcefree fields from boundary data have to be tested before they are applied to measurements. In principle any analytical or numerically created solution of the forcefree equations (1)(3) can be used as a reference case. One cuts a plane (artificial photosphere, bottom boundary) out of the 3D reference solution ^{3}^{3}3For the pure task of code testing it is also acceptable to use all six boundaries of the reference solution. These kind of data are not available for real solar cases of course. and uses the above described extrapolation codes to reconstruct the magnetic field. The result of this extrapolation is then compared with the reference to rate the quality of the reconstruction. Unfortunately, it is very hard to find a truly nonlinear 3D solution of (1)(3) analytically and very few solutions are known. Low and Lou (1990) (LL) found a class of solutions which have become a standard reference for testing NLFFFextrapolation codes. LL found axisymmetric equilibria which are separable in spherical coordinates. They are selfsimilar in the radial coordinate, and the polar angle dependence is determined from a nonlinear eigenvalue equation. The symmetry is broken by cutting out a rectangular chunk of the solution by using a Cartesian coordinate system which is shifted and rotated with respect to the original coordinate system in which the LL equilibria are calculated. The parameters of the LL solutions and the parameters of the new Cartesian coordinate system allow for a large number of different situations which can be used for tests. The original axisymmetric spherical LL solution has also been used (with and without symmetry breaking by shifting the origin of the coordinate system) to test spherical NLFFF programs. To our knowledge all recent implementations of the described NLFFF approaches have been tested with LL, either immediately in the original codedescribing papers or in subsequent works, e.g., in a blindalgorithm test within the NLFFFconsortium, as described below.
The MHDrelaxation method and the optimization approach have been compared by Wiegelmann and Neukirch (2003). Both methods have been applied to the Low and Lou (1990) equilibrium with exactly the same finite difference grid. The iterative equations for MHDrelaxation and optimization have both the form but the structure of is more complicated for optimization than for MHDrelaxation. The MHDrelaxation term is indeed identical with the first term of the optimization approach. While MHDrelaxation minimizes only the Lorentzforce, the optimization does additional minimize , while a decreasing magnetic field divergence during MHDrelaxation (as seen in Wiegelmann and Neukirch (2003)) is the result of numerical diffusion. Despite the numerical overhead in computing for the optimization code, optimization provided more accurate results and faster convergence.
A practical advantage of the MHDapproach is that several timedependent MHDcodes are well known and established and can be used for the forcefree relaxation discussed here. The inclusion of nonmagnetic forces like pressure gradients and gravity looks straight forward for the MHD approach. Other methods are usually developed with the only task of computing nonlinear forcefree coronal magnetic fields, also a generalization towards magnetohydrostatic and stationary MHDequilibria is possible and has been done for the optimization approach (see Wiegelmann and Inhester, 2003; Wiegelmann and Neukirch, 2006). Another advantage of using timedependent MHDcodes for relaxation is that the computed forcefree equilibrium can be used on the same grid and with the same code as initial state for timedependent MHDsimulations. One can, in principle, use the forcefree equilibria computed with any of the described method as initial state for timedependent MHDsimulation, but having the initial equilibrium state already directly on the MHDgrid might be very handy, because no further adjustments are needed.
4.1 The NLFFFconsortium
Since the year 2004 activities are ongoing to bring NLFFFmodelers together and to compare the different existing codes. A workshop series has been organized for this aim by Karel Schrijver and three workshops took place so far from 20042006. The next workshop is planed for june 2007. As we have been asked to summarize the workshop results on the CSWMmeeting, we give also a very brief overview in the corresponding special issuepaper here. The main results of the first two workshops have been published in Schrijver et al. (2006). In this paper six different NLFFFimplementations (GradRubin codes of Amari et al. (1999) and Wheatland (2004), MHDrelaxation code of Valori et al. (2005), optimization codes by McTiernan and Wiegelmann (2004), boundary element method by Yan and Sakurai (2000)) have been compared. The codes have been tested in a blind algorithm test with the help of the semianalytic equilibrium by Low and Lou (1990) in two cases. In case I all six boundaries of a computational box have been described and in case II only the bottom boundary. The comparison of the extrapolation results with the reference solution has been done qualitatively by magnetic field line plots (Shown here in Fig. 3 for the central region of case II) and quantitatively by a number of sophisticated comparison metrices. All NLFFFfields agreed best with the reference field for the low lying central magnetic field region, where the magnetic field and electric currents are strongest and the influence of the boundaries lowest. The code converged with speeds that differed by a factor of one million per iteration steps ^{4}^{4}4The codes run on different machines, have been written in different programming languages and used different compilers. A real test of the exact computing time would comprise a proper operation count, e.g., the number of fixed point additions and multiplications per iteration step.. The fastestconverging and bestperforming code was the Wheatland et al. (2000) optimization code as implemented by Wiegelmann (2004). Recent implementations of the GradRubin code by Amari et al. (2006); Inhester and Wiegelmann (2006) and a new implementation of the upward integration method by Song et al. (2006) did not participate in the blindalgorithm intercomparison by Schrijver et al. (2006), but these three new codes have been tested by the authors with similar measures and revealed similar accuracy as the best performing codes in the blind algorithm test. It seems that the somewhat more flexible boundary conditions used in the GradRubin approaches of Amari et al. (2006) and Inhester and Wiegelmann (2006) are responsible for the better performance compared to the earlier implementation by Amari et al. (1999), which has been used in the blind algorithm test.
The widely used LLequilibrium contains a very smooth photospheric magnetic field and an extended current distribution. It is therefore also desirable to test NLFFFcodes also with other, more challenging boundary fields, which are less smooth, have localized current distribution and to investigate also the effects of noise and effects from non forcefree boundaries. A somewhat more challenging reference case is the equilibrium found by Titov and Démoulin (1999) (TD). Similar as LL, the TD equilibrium is an axisymmetric equilibrium. The TDmodel contains a potential field which is disturbed by a toroidal nonlinear forcefree current. This equilibrium has been used for testing the MHDrelaxation code (Valori and Kliem, private communication) and the optimization code in Wiegelmann et al. (2006a).
Any numerically created NLFFFmodel might be suitable for code testing, too. It is in particular interesting to use models, which are partly related on observational data. Very recently van Ballegooijen et al. (2007) used lineofsight photospheric measurements from SOHO/MDI to compute a potential field, which was then disturbed by inserting a twisted flux robe and relaxed towards a nonlinear forcefree state with a magnetofrictional method as described in van Ballegooijen (2004). The van Ballegooijen et al. (2007) model is not forcefree in the entire computational domain, but only above a certain height above the bottom boundary (artificial chromosphere). On the lowest boundary (photosphere) the model contains significant nonmagnetic forces. Both the chromospheric as well as the photospheric magnetic field vector from the van Ballegooijen et al. (2007)model have been used to test four of the recently developed extrapolation codes (One GradRubin method, one MHDrelaxation code and two optimization approaches) in a second blind algorithm test by Metcalf et al. (2007). While the I. NLFFFconsortium paper (Schrijver et al. (2006)) used a domain of just pixel, the II. paper used a computational domain of pixel and modern NLFFFcodes where able to compute the nonlinear forcefree field in such relatively large boxes within a few hours for a moderate parallelization on only 14 processors and a memory requirement of GB of Ram. This very recent codecomparison shows a major improvement regarding computing time and suitable grid sizes within less than three years. On the first NLFFFconsortium meeting in 2004, box sizes of some have been a kind of standard or computing times of some two weeks have been reported for boxes. We briefly summarize the results of Metcalf et al. (2007) as:

NLFFFextrapolations from chromospheric data recover the original reference field with high accuracy.

When the extrapolations are applied to the photospheric data, the reference field is not well recovered.

Preprocessing of the photospheric data improve the result, but the accuracy is still lower as for extrapolations from the chromosphere.
5 Conclusions and Outlook
Within the last few years the scientific community showed a growing interest into coronal magnetic fields. ^{5}^{5}5Publications containing the phrase ’coronal magnetic fields’ in title or abstract have been cited less than about times per year until the early 1990th and this number increased to about citations per year in 2004. A peak year was 2006 (last year) with more than citations. (Source: ISI Web of Knowledge, march 2007) The development of new ground based and space born vector magnetographs provide us measurements of the magnetic field vector on the suns photosphere. Accompanied from these hardware development, software has been developed to extrapolate the photospheric measurements into the corona. Special attention has recently been given to nonlinear forcefree codes. Five different numerical approaches (GradRubin, upward integration, MHDrelaxation, optimization, boundary elements) have been developed for this aim. It is remarkable that new codes or major updates of existing codes have been published for all five methods within the last two years, mainly in the last year (2006). A workshop series (NLFFFconsortium) since 2004 on nonlinear forcefree fields has recently released synergy effects, by bringing modelers of the different numerical implementations together to compare, evaluate and improve the programs. Several of the most recent new codes and utility programs (e.g. preprocessing) have at least been partly inspired by these workshops. The new implementations have been tested with the smooth semianalytic LowLouequilibrium and showed reasonable agreement with this reference field. While all methods aim for a reconstruction of the coronal magnetic field from the photospheric magnetic field vector, the way how these measurements are used to prescribe the boundaries of the codes is different.

MHDrelaxation and optimization use on the bottom boundary. This overdetermines the boundary value problem. Both methods are closely related and compute the magnetic field in a computational box with
(54) where the structure of F is somewhat different (the optimization approach has more terms) for both methods. Usually a potential field is used as initial state for both approaches, also the use of a linear forcefree initial state is possible. Recently a multiscale version of optimization has been installed, which uses a low resolution NLFFFfield as input for higher resolution computations. Specifying the entire magnetic field vector on the bottom boundary is an overimposed problem and a unique NLFFFfield (or a solution at all) requires that the boundary data fulfill certain consistency criteria. A recently developed preprocessingroutine helps to find suitable consistent boundary data from inconsistent photospheric measurements. Earlier and current comparisons showed a somewhat higher accuracy for the optimization approach. A practical advantage of the MHDapproach is that in principle any available timedependent MHDcode can be adjusted to compute the NLFFFfield.

The GradRubin approach uses and the distribution of computed with Eq. (8) for one polarity, which corresponds to well posed mathematical problem. A practical problem is that the computation of requires numerical differences of the noisy and forced transverse photospheric field with (7) leading to inaccuracies in the normal electric current distribution and in . For smooth semianalytic test cases this is certainly not a problem, but real data require special attention (smoothing, preprocessing, limiting to regions where is above a certain limit) to derive a meaningful distribution of . While the method requires only for one polarity, the computation from photospheric data provide for both polarities. We are not aware of any tests on how well NLFFFsolutions computed from prescribed on the positive and negative polarity coincide. It is also unclear how well the computed transverse field components on the bottom boundary agree with the measured values ^{6}^{6}6In principle may have an additional field without making a difference for and hence for the GradRubin result. of . More tests on this topics are necessary, including the recently installed possibility to prescribe for both polarities and adjust the boundary by a weighed average of on both polarities to fulfill Eq. (6). As initial state the GradRubin method uses a potential field, which is also true for MHDrelaxation and optimization.

The upward integration and the boundary element method prescribe both all components of the bottom boundary magnetic field vector and the distribution computed with Eq. (8). This approach over imposes the boundary and and have to be consistent which each other and the forcefree assumption. This is certainly not a problem at all for smooth semianalytic test equilibria and strategies to derive consistent boundary data from measured data have been developed recently. Different from the three approaches discussed above, upward integration and boundary element methods do not require to compute first an initial potential field in the computational domain. It is well known that the upward integration method is based on an illposed problem and the method has not been considered for several years, but a recent implementation with smooth analytic functions might help to regularize this method. First tests showed a reasonable results for computations with the smooth semianalytic LowLou solution.
The boundary element method has the problem to be very slow and an earlier implementation of this method could not reach a converged state for a boxed used in the I. NLFFFconsortium paper due to this problem. A new ’direct boundary method’ has been developed, which seems to be faster than the original ’boundary element method’, but still slower compared with the four other NLFFFapproaches if the task is to compute a 3D magnetic field in an entire 3Ddomain. Different from all other described methods the boundary element approach allows to compute the nonlinear forcefree field vector at any arbitrary point above the boundary and it is not necessary to compute the entire 3Dfield above the photosphere. This might be a very useful feature if one is interested in computing the magnetic field only along a single loop and not interested in an entire active region.
The new implementations of upward integration and boundary element method show both reasonable results for first tests with the smooth semianalytic Low and Lou equilibrium. Further tests with more sophisticated equilibria, e.g. a solarlike test case as used in the II. NLFFFconsortium paper would be useful to come to more sound conclusions regarding the feasibility of these methods.
Most of the efforts done in nonlinear forcefree modelling until now concentrated mainly on developing these models and testing their accuracy and speed with the help of well known test configuration. Not too many applications of nonlinear forcefree models to real data are currently available, from which we learned new physics. One reason was the insufficient access to high accuracy photospheric vector magnetograms and a second one were limitations of the models. Forcefree field extrapolation is a mere tool, if properly employed on vector magnetograms, it can help to understand physical, magnetic field dominated processes in the corona. Both the computational methods as well as the accuracy of required measurements (e.g. with Hinode, SDO) are rapidly improving. Within the NLFFFconsortium we just started (since april 2007) to apply the different codes to compute nonlinear forcefree coronal magnetic fields from Hinode vectormagnetograms. This project might provide us already some new insights about coronal physics.
To conclude, we can say that the capability of Cartesian nonlinear forcefree extrapolation codes has rapidly increased in recent years. Only three years ago most codes run usually on grids of about pixel. Recently developed or updated codes (GradRubin by Wheatland, MHDrelaxation by Valori, optimization by Wiegelmann, optimization by McTiernan) have been applied to grids of about pixel. Although this increase of traceable grid sizes is certainly encouraging, the resolution of current and near future vector magnetographs (which of course measure only data in 2D!) is significantly higher. We should keep in mind, however, that the currently implemented NLFFFcodes have been only moderately parallelized using only a few processors. The CSWMconference, where this paper has been presented, took place at the ’Earth simulator’ in Yokohama, which contains several thousands of processors used for Earthscience computer simulations. An installation of NLFFFcodes on such massive parallel computers (which has been briefly addressed on NLFFFconsortium meetings) combined with adaptive mesh refinements might enable drastically improved grid sizes. One should not underestimate the time and effort necessary to program and install such massive parallelized versions of existing codes. As full disk vectormagnetograms will become available soon (SOLIS, SDO/HMI) it is also an important task to take a spherical geometry into account. First steps in this direction have been carried out with the optimization and boundary element methods. Spherical NLFFFgeometries are currently still in it’s infancy and have been tested until now only with smooth semianalytic Low and Lou equilibria and require further developments.
Attention has also recently been drawn to the problem that the coronal magnetic field is forcefree, but the photospheric one is not. Tests with extrapolations from solarlike artificial photospheric and chromospheric measurements within the II. NLFFFconsortium paper revealed that extrapolations from the (forcefree) chromospheric field provide significantly better results as extrapolations using directly the (forced) photospheric field. Applying a preprocessing program on the photospheric data, which effectively removes the nonmagnetic forces, leads to significantly better results, but they are not as good as by using the chromospheric magnetic field vector as boundary condition. An area of current research is the possibility to use chromospheric images to improve the preprocessing of photospheric magnetic field measurements. Improvements in measuring the chromospheric magnetic field directly (e.g. Lagg et al., 2004) might further improve to find suitable boundary conditions for NLFFFextrapolations. Forcefree extrapolations are not suitable, however, to understand the details of physical processes on how the magnetic field evolves from the forced photosphere into the chromosphere, because nonmagnetic forces are important in the photosphere. For a better understanding of these phenomena more sophisticated models which take pressure gradients and gravity (and maybe also plasma flow) into account are required. Some first steps have been done with a generalization of the optimization method by Wiegelmann and Neukirch (2006), but such approaches are still in their infancy and have been tested so far only with smooth MHDequilibria. It is also not entirely clear how well necessary information regarding the plasma (density, pressure, temperature, flow) can be derived from measurements. Nonmagnetic forces become important also in quiet sun regions (Schrijver and van Ballegooijen (2005)) and in the higher layers of the corona, where the plasma is of the order of unity. Coronagraph measurements, preferably from two viewpoints as provided by the STEREOmission, combined with a tomographic inversion might help here to get insights in the required 3D structure of the plasma density. One should also pay attention to the combination of extrapolation methods, as described here, with measurements of the Hanle and Zeeman effects in coronal lines which allows the reconstruction of the coronal magnetic field as proposed in feasibility studies of vector tomography by Kramar et al. (2006); Kramar and Inhester (2006). Other measurements of coronal features, e.g., coronal plasma images from two STEREOviewpoints, can be used for observational tests of coronal magnetic field models. Using two viewpoints provide a much more restrictive test of models as images from only one view direction. While a nonlinear forcefree coronal magnetic field model helps us to derive the topology, magnetic field and electric current strength in coronal loops, they do not provide plasma parameters. One way to get insights regarding the coronal plasma is the use of scaling laws to model the plasma along the reconstructed 3D field lines and compare correspondent artificial plasma images with real coronal images. Schrijver et al. (2004) applied such an approach to global potential coronal magnetic fields and compared simulated and real coronal images from one viewpoint. A generalization of such methods towards the use of more sophisticated magnetic field models and coronal images from two STEREOviewpoints will probably provide many insights regarding the structure and physics of the coronal plasma. An important challenge is for example the coronal heating problem. The dominating coronal magnetic field is assumed to play an important role here, because magnetic field configuration containing free energy can under certain circumstances reconnect (Priest (1996, 1999)) and supply energy for coronal heating. Priest et al. (2005) pointed out that magnetic reconnection at separators and separatrices plays an important role for coronal heating. Nonlinear forcefree models can help here to identify the magnetic field topology, magnetic null points, separatrices and localized strong current concentration. While magnetic reconnection (see e.g. Priest and Schrijver, 1999) is a dynamical phenomenon, the static magnetic field models discussed here can help to identify the locations favourable for reconnection. Time sequences of nonlinear forcefree models computed from corresponding vector magnetograms will also tell wether the topology of the coronal magnetic field has changed due to reconnection, even if the physics of reconnection is not described by forcefree models. Sophisticated 3D coronal magnetic field models and plasma images from two viewpoints might help to constrain the coronal heating function further, which has been done so far with plasma images from one viewpoint (Aschwanden, 2001a, b, by using data from Yokoh, Soho and Trace).
Acknowledgements.
I would like to thank the organizers of CSWM for inviting me to give a review talk on nonlinear forcefree magnetic field modelling. The author thanks Bernd Inhester for useful comments and acknowledge inspiriting discussions with many other NLFFFmodelers on three workshops organized by Karel Schrijver (NLFFFconsortium) in Palo Alto between 2004 and 2006 and at various other occasions. This work was supported by DLRgrant 50 OC 0501.References
 Alissandrakis (1981) Alissandrakis, C. E. (1981), On the computation of constant alpha forcefree magnetic field, A&A, 100, 197–200.
 Aly (1984) Aly, J. J. (1984), On some properties of forcefree magnetic fields in infinite regions of space, apj, 283, 349–362, doi:10.1086/162313.
 Aly (1989) Aly, J. J. (1989), On the reconstruction of the nonlinear forcefree coronal magnetic field from boundary data, Sol. Phys., 120, 19–48.
 Aly (2005) Aly, J. J. (2005), A uniqueness result for a simple forcefree magnetic field submitted to a topological constraint, aap, 429, 15–18, doi:10.1051/00046361:20041547.
 Amari et al. (1997) Amari, T., J. J. Aly, J. F. Luciani, T. Z. Boulmezaoud, and Z. Mikic (1997), Reconstructing the Solar Coronal Magnetic Field as a ForceFree Magnetic Field, Sol. Phys., 174, 129–149.
 Amari et al. (1999) Amari, T., T. Z. Boulmezaoud, and Z. Mikic (1999), An iterative method for the reconstructionbreak of the solar coronal magnetic field. I. Method for regular solutions, A&A, 350, 1051–1059.
 Amari et al. (2006) Amari, T., T. Z. Boulmezaoud, and J. J. Aly (2006), Well posed reconstruction of the solar coronal magnetic field, A&A, 446, 691–705, doi:10.1051/00046361:20054076.
 Aschwanden (2001a) Aschwanden, M. J. (2001a), Revisiting the Determination of the Coronal Heating Function from Yohkoh Data, apjl, 559, L171–L174, doi:10.1086/323788.
 Aschwanden (2001b) Aschwanden, M. J. (2001b), An Evaluation of Coronal Heating Models for Active Regions Based on Yohkoh, SOHO, and TRACE Observations, apj, 560, 1035–1044, doi:10.1086/323064.
 Aschwanden (2005) Aschwanden, M. J. (2005), Physics of the Solar Corona. An Introduction with Problems and Solutions (2nd edition), Pour la Science.
 Aschwanden et al. (1999) Aschwanden, M. J., J. S. Newmark, J.P. Delaboudinière, W. M. Neupert, J. A. Klimchuk, G. A. Gary, F. PortierFozzani, and A. Zucker (1999), Threedimensional Stereoscopic Analysis of Solar Active Region Loops. I. SOHO/EIT Observations at Temperatures of (1.01.5) X 10^6 K, ApJ, 515, 842–867, doi:10.1086/307036.
 Bineau (1972) Bineau, M. (1972), Existence of forcefree magneticfields, Comm. Pure Appl. Math., 25, 77–.
 Bleybel et al. (2002) Bleybel, A., T. Amari, L. van DrielGesztelyi, and K. D. Leka (2002), Global budget for an eruptive active region . I. Equilibrium reconstruction approach, A&A, 395, 685–695, doi:10.1051/00046361:20021332.
 Borrero et al. (2006) Borrero, J. M., S. Tomczyk, A. Norton, T. Darnell, J. Schou, P. Scherrer, R. Bush, and Y. Lui (2006), Magnetic field vector retrieval with HMI, ArXiv Astrophysics eprints.
 Bossavit (1988) Bossavit, A. (1988), Mixed finite elements and the complex of Whitney forms, in The mathemetics of finite elements and applications VI, edited by J. Whiteman, pp. 137–144, Academic Press, London.
 Boulmezaoud and Amari (2000) Boulmezaoud, T. Z., and T. Amari (2000), On the existence of nonlinear forcefree fields in threedimensional domains, Zeitschrift Angewandte Mathematik und Physik, 51, 942–967.
 Carcedo et al. (2003) Carcedo, L., D. S. Brown, A. W. Hood, T. Neukirch, and T. Wiegelmann (2003), A Quantitative Method to Optimise Magnetic Field Line Fitting of Observed Coronal Loops, Sol. Phys., 218, 29–40.
 Chiu and Hilton (1977) Chiu, Y. T., and H. H. Hilton (1977), Exact Green’s function method of solar forcefree magneticfield computations with constant alpha. I  Theory and basic test cases, ApJ, 212, 873–885.
 Chodura and Schlueter (1981) Chodura, R., and A. Schlueter (1981), A 3D code for MHD equilibrium and stability, Journal of Computational Physics, 41, 68–88.
 Cuperman et al. (1990) Cuperman, S., L. Ofman, and M. Semel (1990), Extrapolation of photospheric potential magnetic fields using oblique boundary values  A simplified approach, A&A, 227, 583–590.
 Cuperman et al. (1991) Cuperman, S., P. Demoulin, and M. Semel (1991), Removal of singularities in the Cauchy problem for the extrapolation of solar forcefree magnetic fields, A&A, 245, 285–288.
 Dedner et al. (2002) Dedner, A., F. Kemm, D. Kröner, C.D. Munz, T. Schnitzer, and M. Wesenberg (2002), Hyperbolic Divergence Cleaning for the MHD Equations, Journal of Computational Physics, 175, 645–673.
 Demoulin and Priest (1992) Demoulin, P., and E. R. Priest (1992), The properties of sources and sinks of a linear forcefree field, A&A, 258, 535–541.
 Gary (2001) Gary, G. A. (2001), Plasma Beta above a Solar Active Region: Rethinking the Paradigm, Sol. Phys., 203, 71–86.
 Gary and Hagyard (1990) Gary, G. A., and M. J. Hagyard (1990), Transformation of vector magnetograms and the problems associated with the effects of perspective and the azimuthal ambiguity, solphys, 126, 21–36.
 Grad and Rubin (1958) Grad, H., and H. Rubin (1958), , in Proc. 2nd Intern. Conf. on Peaceful Uses of Atomic Energy, Vol 31, pp. 190–+.
 He and Wang (2006) He, H., and H. Wang (2006), The validity of the boundary integral equation for magnetic field extrapolation in open space above a spherical surface, mnras, 369, 207–215, doi:10.1111/j.13652966.2006.10288.x.
 Henney et al. (2006) Henney, C. J., C. U. Keller, and J. W. Harvey (2006), SOLISVSM Solar Vector Magnetograms, ArXiv Astrophysics eprints.
 Inhester and Wiegelmann (2006) Inhester, B., and T. Wiegelmann (2006), Nonlinear ForceFree Magnetic Field Extrapolations: Comparison of the Grad Rubin and Wheatland Sturrock Roumeliotis Algorithm, Sol. Phys., 235, 201–221, doi:10.1007/s112070060065x.
 Jiao et al. (1997) Jiao, L., A. N. McClymont, and Z. Mikic (1997), Reconstruction of the ThreeDimensional Coronal Magnetic Field, solphys, 174, 311–327.
 Judge (1998) Judge, P. G. (1998), Spectral Lines for Polarization Measurements of the Coronal Magnetic Field. I. Theoretical Intensities, ApJ, 500, 1009–+.
 Kageyama and Sato (2004) Kageyama, A., and T. Sato (2004), “YinYang grid”: An overset grid in spherical geometry, Geochemistry, Geophysics, Geosystems, 5, 9005–+, doi:10.1029/2004GC000734.
 Keppens et al. (2003) Keppens, R., M. Nool, G. Tóth, and J. P. Goedbloed (2003), Adaptive Mesh Refinement for conservative systems: multidimensional efficiency evaluation, Computer Physics Communications, 153, 317–339.
 Klimchuk and Sturrock (1992) Klimchuk, J. A., and P. A. Sturrock (1992), Threedimensional forcefree magnetic fields and flare energy buildup, apj, 385, 344–353, doi:10.1086/170943.
 Kramar and Inhester (2006) Kramar, M., and B. Inhester (2006), Inversion of coronal Zeeman and Hanle Observations to reconstruct the coronal magnetic field, ArXiv Astrophysics eprints.
 Kramar et al. (2006) Kramar, M., B. Inhester, and S. K. Solanki (2006), Vector tomography for the coronal magnetic field. I. Longitudinal Zeeman effect measurements, A&A, 456, 665–673, doi:10.1051/00046361:20064865.
 LaBonte et al. (1999) LaBonte, B. J., D. L. Mickey, and K. D. Leka (1999), The Imaging Vector Magnetograph at Haleakalā  II. Reconstruction of Stokes Spectra, solphys, 189, 1–24.
 Lagg et al. (2004) Lagg, A., J. Woch, N. Krupp, and S. K. Solanki (2004), Retrieval of the full magnetic vector with the He I multiplet at 1083 nm. Maps of an emerging flux region, A&A, 414, 1109–1120.
 Leka and Skumanich (1999) Leka, K. D., and A. Skumanich (1999), On the value of ‘AR’ from vector magnetograph data  I. Methods and Caveats, Sol. Phys., 188, 3–19.
 Li et al. (2004) Li, Z., Y. Yan, and G. Song (2004), Properties of the boundary integral equation for solar nonconstant forcefree magnetic fields, MNRAS, 347, 1255–1265, doi:10.1111/j.13652966.2004.07309.x.
 Lin et al. (2004) Lin, H., J. R. Kuhn, and R. Coulter (2004), Coronal Magnetic Field Measurements, ApJL, 613, L177–L180, doi:10.1086/425217.
 Liu et al. (2002) Liu, Y., X. P. Zhao, J. T. Hoeksema, P. H. Scherrer, J. Wang, and Y. Yan (2002), On Formation of the Sigmoidal Structure in Solar Active Region NOAA 8100, Sol. Phys., 206, 333–346.
 Low and Lou (1990) Low, B. C., and Y. Q. Lou (1990), Modeling solar forcefree magnetic fields, ApJ, 352, 343–352.
 McClymont et al. (1997) McClymont, A. N., L. Jiao, and Z. Mikic (1997), Problems and Progress in Computing ThreeDimensional Coronal Active Region Magnetic Fields from Boundary Data, Sol. Phys., 174, 191–218.
 Metcalf (1994) Metcalf, T. R. (1994), Resolving the 180degree ambiguity in vector magnetic field measurements: The ’minimum’ energy solution, Sol. Phys., 155, 235–242.
 Metcalf et al. (1995) Metcalf, T. R., L. Jiao, A. N. McClymont, R. C. Canfield, and H. Uitenbroek (1995), Is the solar chromospheric magnetic field forcefree?, apj, 439, 474–481, doi:10.1086/175188.
 Metcalf et al. (2007) Metcalf, T. R., M. L. DeRosa, C. J. Schrijver, G. Barnes, A. VanBallegooijen, T. Wiegelmann, M. S. Wheatland, G. Valori, and J. M. McTiernan (2007), Nonlinear forcefree modeling of coronal magnetic fields. II. Modeling a filament arcade from simulated chromospheric and photospheric vector fields., solphys, submitted.
 Metcalf et al. (2006) Metcalf, T. R., et al. (2006), An Overview of Existing Algorithms for Resolving the 180degree Ambiguity in Vector Magnetic Fields: Quantitative Tests with Synthetic Data, solphys, 237, 267–296, doi:10.1007/s112070060170x.
 Mickey et al. (1996) Mickey, D. L., R. C. Canfield, B. J. Labonte, K. D. Leka, M. F. Waterson, and H. M. Weber (1996), The Imaging Vector Magnetograph at Haleakala, solphys, 168, 229–250.
 Mikic and McClymont (1994) Mikic, Z., and A. N. McClymont (1994), Deducing Coronal Magnetic Fields from Vector Magnetograms, in ASP Conf. Ser. 68: Solar Active Region Evolution: Comparing Models with Observations, edited by K. S. Balasubramaniam and G. W. Simon, pp. 225–232.
 Mikic et al. (1988) Mikic, Z., D. C. Barnes, and D. D. Schnack (1988), Dynamical evolution of a solar coronal magnetic field arcade, apj, 328, 830–847, doi:10.1086/166341.
 Molodensky (1969) Molodensky, M. M. (1969), Integral properties of forcefree fields, Soviet Astron.AJ, 12, 585–588.
 Molodensky (1974) Molodensky, M. M. (1974), Equilibrium and stability of forcefree magnetic field, Sol. Phys., 39, 393–404.
 Nakagawa (1974) Nakagawa, Y. (1974), Dynamics of the Solar Magnetic Field. I. Method of Examination of ForceFree Magnetic Fields, ApJ, 190, 437–440.
 Nakagawa and Raadu (1972) Nakagawa, Y., and M. A. Raadu (1972), On Practical Representation of Magnetic Field, solphys, 25, 127–+.
 Pevtsov et al. (1994) Pevtsov, A. A., R. C. Canfield, and T. R. Metcalf (1994), Patterns of helicity in solar active regions, ApJL, 425, L117–L119, doi:10.1086/187324.
 Press (2002) Press, W. H. (2002), Numerical recipes in C++ : the art of scientific computing, Numerical recipes in C++ : the art of scientific computing by William H. Press. xxviii, 1,002 p. : ill. ; 26 cm. Includes bibliographical references and index. ISBN : 0521750334.
 Priest (1996) Priest, E. (1996), Coronal Heating by Magnetic Reconnection, apss, 237, 49–73.
 Priest (1999) Priest, E. R. (1999), Heating the Solar Corona by Magnetic Reconnection, apss, 264, 77–100.
 Priest and Schrijver (1999) Priest, E. R., and C. J. Schrijver (1999), Aspects of ThreeDimensional Magnetic Reconnection (Invited Review), solphys, 190, 1–24, doi:10.1023/A:1005248007615.
 Priest et al. (2005) Priest, E. R., D. W. Longcope, and J. Heyvaerts (2005), Coronal Heating at Separators and Separatrices, apj, 624, 1057–1071, doi:10.1086/429312.
 Régnier and Amari (2004) Régnier, S., and T. Amari (2004), 3D magnetic configuration of the H filament and Xray sigmoid in NOAA AR 8151, A&A, 425, 345–352.
 Régnier et al. (2002) Régnier, S., T. Amari, and E. Kersalé (2002), 3D Coronal magnetic field from vector magnetograms: nonconstantalpha forcefree configuration of the active region NOAA 8151, A&A, 392, 1119–1127.
 Régnier and Canfield (2006) Régnier, S., and R. C. Canfield (2006), Evolution of magnetic fields and energetics of flares in active region 8210, \aap, 451, 319–330, doi:10.1051/00046361:20054171.
 Régnier et al. (2005) Régnier, S., T. Amari, and R. C. Canfield (2005), Self and mutual magnetic helicities in coronal magnetic configurations, \aap, 442, 345–349, doi:10.1051/00046361:20053509.
 Roumeliotis (1996) Roumeliotis, G. (1996), The “StressandRelax” Method for Reconstructing the Coronal Magnetic Field from Vector Magnetograph Data, ApJ, 473, 1095–+.
 Sakurai (1981) Sakurai, T. (1981), Calculation of ForceFree Magnetic Field with Non Constant Alpha, Sol. Phys., 69, 343–+.
 Sakurai (1982) Sakurai, T. (1982), Green’s Function Methods for Potential Magnetic Fields, Sol. Phys., 76, 301–+.
 Sakurai (1989) Sakurai, T. (1989), Computational modeling of magnetic fields in solar active regions, Space Science Reviews, 51, 11–48.
 Sakurai et al. (1995) Sakurai, T., et al. (1995), Solar flare telescope at Mitaka, Publ. Astron. Soc. Japan, 47, 81–92.
 Schatten et al. (1969) Schatten, K. H., J. M. Wilcox, and N. F. Ness (1969), A model of interplanetary and coronal magnetic fields, Sol. Phys., 6, 442–455.
 Schmidt (1964) Schmidt, H. U. (1964), On the Observable Effects of Magnetic Energy Storage and Release Connected With Solar Flares, in The Physics of Solar Flares, pp. 107–+.
 Schrijver and van Ballegooijen (2005) Schrijver, C. J., and A. A. van Ballegooijen (2005), Is the QuietSun Corona a Quasisteady, Forcefree Environment?, ApJ, 630, 552–560, doi:10.1086/431754.
 Schrijver et al. (2004) Schrijver, C. J., A. W. Sandman, M. J. Aschwanden, and M. L. DeRosa (2004), The Coronal Heating Mechanism as Identified by FullSun Visualizations, ApJ, 615, 512–525, doi:10.1086/424028.
 Schrijver et al. (2005) Schrijver, C. J., M. L. DeRosa, A. M. Title, and T. R. Metcalf (2005), The Nonpotentiality of ActiveRegion Coronae and the Dynamics of the Photospheric Magnetic Field, ApJ, 628, 501–513, doi:10.1086/430733.
 Schrijver et al. (2006) Schrijver, C. J., M. L. Derosa, T. R. Metcalf, Y. Liu, J. McTiernan, S. Régnier, G. Valori, M. S. Wheatland, and T. Wiegelmann (2006), Nonlinear ForceFree Modeling of Coronal Magnetic Fields Part I: A Quantitative Comparison of Methods, Sol. Phys., 235, 161–190, doi:10.1007/s1120700600687.
 Seehafer (1978) Seehafer, N. (1978), Determination of constant alpha forcefree solar magnetic fields from magnetograph data, Sol. Phys., 58, 215–223.
 Seehafer (1982) Seehafer, N. (1982), A comparison of different solar magnetic field extrapolation procedures, Sol. Phys., 81, 69–80.
 Semel (1967) Semel, M. (1967), Contribution à létude des champs magnétiques dans les régions actives solaires, Annales d’Astrophysique, 30, 513–513.
 Semel (1988) Semel, M. (1988), Extrapolation functions for constantalpha forcefree fields  Green’s method for the oblique boundary value, A&A, 198, 293–299.
 Shimizu (2004) Shimizu, T. (2004), SolarB Solar Optical Telescope (SOT), in ASP Conf. Ser. 325: The SolarB Mission and the Forefront of Solar Physics, edited by T. Sakurai and T. Sekii, pp. 3–+.
 Solanki et al. (2003) Solanki, S. K., A. Lagg, J. Woch, N. Krupp, and M. Collados (2003), Threedimensional magnetic field topology in a region of solar coronal heating, Nature, 425, 692–695.
 Song et al. (2006) Song, M. T., C. Fang, Y. H. Tang, S. T. Wu, and Y. A. Zhang (2006), A New and Fast Way to Reconstruct a Nonlinear Forcefree Field in the Solar Corona, \apj, 649, 1084–1092, doi:10.1086/506249.
 Titov and Démoulin (1999) Titov, V. S., and P. Démoulin (1999), Basic topology of twisted magnetic configurations in solar flares, A&A, 351, 707–720.
 Török and Kliem (2003) Török, T., and B. Kliem (2003), The evolution of twisting coronal magnetic flux tubes, aap, 406, 1043–1059, doi:10.1051/00046361:20030692.
 Valori et al. (2005) Valori, G., B. Kliem, and R. Keppens (2005), Extrapolation of a nonlinear forcefree field containing a highly twisted magnetic loop, A&A, 433, 335–347.
 van Ballegooijen (2004) van Ballegooijen, A. A. (2004), Observations and Modeling of a Filament on the Sun, apj, 612, 519–529, doi:10.1086/422512.
 van Ballegooijen et al. (2007) van Ballegooijen, A. A., E. E. Deluca, K. Squires, and D. H. Mackay (2007), Modeling magnetic flux ropes in the solar atmosphere, Journal of Atmospheric and Terrestrial Physics, 69, 24–31, doi:10.1016/j.jastp.2006.06.007.
 Wang et al. (2000) Wang, H., Y. Yan, T. Sakurai, and M. Zhang (2000), Topology of Magnetic Field and Coronal Heating in Solar Active Regions  II. The Role of QuasiSeparatrix Layers, Sol. Phys., 197, 263–273.
 Wheatland (1999) Wheatland, M. S. (1999), A Better Linear Forcefree Field, apj, 518, 948–953, doi:10.1086/307301.
 Wheatland (2004) Wheatland, M. S. (2004), Parallel Construction of Nonlinear ForceFree Fields, Sol. Phys., 222, 247–264.
 Wheatland (2006) Wheatland, M. S. (2006), A Fast CurrentField Iteration Method for Calculating Nonlinear ForceFree Fields, Sol. Phys., 238, 29–39, doi:10.1007/s1120700602320.
 Wheatland et al. (2000) Wheatland, M. S., P. A. Sturrock, and G. Roumeliotis (2000), An Optimization Approach to Reconstructing Forcefree Fields, ApJ, 540, 1150–1155.
 Wiegelmann (2004) Wiegelmann, T. (2004), Optimization code with weighting function for the reconstruction of coronal magnetic fields, Sol. Phys., 219, 87–108.
 Wiegelmann (2007) Wiegelmann, T. (2007), Computing Nonlinear ForceFree Coronal Magnetic Fields in Spherical Geometry, Sol. Phys., 240, 227–239, doi:10.1007/s1120700602663.
 Wiegelmann and Inhester (2003) Wiegelmann, T., and B. Inhester (2003), Magnetic modeling and tomography: First steps towards a consistent reconstruction of the solar corona, Sol. Phys., 214, 287–312.
 Wiegelmann and Neukirch (2002) Wiegelmann, T., and T. Neukirch (2002), Including stereoscopic information in the reconstruction of coronal magnetic fields, Sol. Phys., 208, 233–251.
 Wiegelmann and Neukirch (2003) Wiegelmann, T., and T. Neukirch (2003), Computing nonlinear force free coronal magnetic fields, Nonlinear Processes in Geophysics, 10, 313–322.
 Wiegelmann and Neukirch (2006) Wiegelmann, T., and T. Neukirch (2006), An optimization principle for the computation of MHD equilibria in the solar corona, A&A, 457, 1053–1058, doi:10.1051/00046361:20065281.
 Wiegelmann et al. (2005a) Wiegelmann, T., B. Inhester, A. Lagg, and S. K. Solanki (2005a), How To Use Magnetic Field Information For Coronal Loop Identification, Sol. Phys., 228, 67–78, doi:10.1007/s1120700525116.
 Wiegelmann et al. (2005b) Wiegelmann, T., A. Lagg, S. K. Solanki, B. Inhester, and J. Woch (2005b), Comparing magnetic field extrapolations with measurements of magnetic loops, A&A, 433, 701–705.
 Wiegelmann et al. (2006a) Wiegelmann, T., B. Inhester, B. Kliem, G. Valori, and T. Neukirch (2006a), Testing nonlinear forcefree coronal magnetic field extrapolations with the TitovDémoulin equilibrium, A&A, 453, 737–741, doi:10.1051/00046361:20054751.
 Wiegelmann et al. (2006b) Wiegelmann, T., B. Inhester, and T. Sakurai (2006b), Preprocessing of vector magnetograph data for a nonlinear forcefree magnetic field reconstruction., Sol. Phys., 233, 215–232.
 Wu et al. (1985) Wu, S. T., H. M. Chang, and M. J. Hagyard (1985), On the numerical computation of nonlinear forcefree magnetic fields, in Measurements of Solar Vector Magnetic Fields, NASA Conf. Pub.,2374, pp. 17–40.
 Wu et al. (1990a) Wu, S. T., M. T. Sun, H. M. Chang, M. J. Hagyard, and G. A. Gary (1990a), On the numerical computation of nonlinear forcefree magnetic fields, ApJ, 362, 698–708.
 Wu et al. (1990b) Wu, S. T., M. T. Sun, and T. Sakurai (1990b), A comparison between progressive extension method (PEM) and iterative method (IM) for magnetic field extrapolations in the solar atmosphere, Memorie della Societa Astronomica Italiana, 61, 477–484.
 Yan and Li (2006) Yan, Y., and Z. Li (2006), Direct Boundary Integral Formulation for Solar Nonconstant Forcefree Magnetic Fields, ApJ, 638, 1162–1168, doi:10.1086/499064.
 Yan and Sakurai (2000) Yan, Y., and T. Sakurai (2000), New Boundary Integral Equation Representation for Finite Energy ForceFree Magnetic Fields in Open Space above the Sun, Sol. Phys., 195, 89–109.
 Yan et al. (2001a) Yan, Y., M. J. Aschwanden, S. Wang, and Y. Deng (2001a), Evolution of Magnetic Flux Rope in the Active Region NOAA 9077 on 14 July 2000, Sol. Phys., 204, 27–40, doi:10.1023/A:1014265123601.
 Yan et al. (2001b) Yan, Y., Y. Deng, M. Karlický, Q. Fu, S. Wang, and Y. Liu (2001b), The Magnetic Rope Structure and Associated Energetic Processes in the 2000 July 14 Solar Flare, ApJL, 551, L115–L119, doi:10.1086/319829.
 Yang et al. (1986) Yang, W. H., P. A. Sturrock, and S. K. Antiochos (1986), Forcefree magnetic fields  The magnetofrictional method, ApJ, 309, 383–391, doi:10.1086/164610.
 Yee (1966) Yee, K. S. (1966), Numerical solition of initial boundary value problems involving Maxwell’s equations in isotropic media, IEEE Transactions on Antennas and Propagation, AP14, No.4, 302–307.