In the low corona of a star like the Sun, the magnetic pressure dominates over both the gas pressure and the kinetic energy density of plasma flows. Thus to first approximation an equilibrium magnetic field must have a vanishing Lorentz force,force-free. It follows that
Unfortunately, for non-zero , the extrapolation of coronal force-free magnetic fields is not straightforward. If is constant everywhere, then taking the curl of (15) leads to the vector Helmholtz equationlinear force-free fields (Nakagawa, 1973; Levine and Altschuler, 1974), their use is limited for two main reasons. Firstly, the constant in such a solution scales as , where is the horizontal size of the area under consideration. So for global solutions, only rather small values of , close to potential, can be used. Secondly, observations of twist in magnetic structures indicate that both the sign and magnitude of ought to vary significantly between different regions on the Sun. In addition, a mathematical problem arises if one attempts to apply the same boundary conditions as in the potential field model (namely a given distribution of on , and on ). A strictly force-free field satisfying and on an outer boundary must have on that boundary (Aly and Seehafer, 1993), so that on all open magnetic field lines. Unless , this is incompatible with the linear force-free field.5 From this, it is clear that the use of linear force-free fields is restrictive and will not be discussed further. Instead, we will consider nonlinear force-free fields where is a function of position.
The problem of extrapolating nonlinear force-free fields from given photospheric data is mathematically challenging, with open questions about the existence and uniqueness of solutions. Nevertheless, several numerical techniques have been developed in recent years, though they have largely been applied to single active regions in Cartesian geometry (Schrijver et al., 2006; DeRosa et al., 2009). Applications to global solutions in spherical geometry are in their infancy: we describe here the three main approaches tried.
Wiegelmann (2007) has developed a method for numerically computing nonlinear force-free fields in spherical geometry, based on the optimisation procedure of Wheatland et al. (2000) and using the vector magnetic field in the photosphere as input. The idea is to minimise the functionalet al. (2000) show that
Wiegelmann (2007) demonstrates that this method recovers a known analytical force-free field solution (Figure 8), although the most accurate solution is obtained only if the boundary conditions at both the photosphere and upper source surface are matched to the analytical field. The method is yet to be applied to observational data on a global scale, largely due to the limitation that vector magnetogram data are required for the lower boundary input. Such data are not yet reliably measured on a global scale, though SDO or SOLIS should be a key step forward in obtaining this data. Another consequence is that even where they are measured there is the complication that the photospheric magnetic field is not force-free. Pre-processing techniques to mitigate these problems are under development (Tadesse et al., 2011).
Contopoulos et al. (2011) have recently proposed an alternative method to compute global nonlinear force-free field extrapolations, adapted from a technique applied to pulsar magnetospheres. It has the advantage that only the radial component is required as observational input at the photosphere. The computation is initialised with an arbitrary 3D magnetic field such as or that of a dipole, and also with zero electric field . The fields and are then evolved through the equations of force-free electrodynamics,
The authors have tested their method with an observed synoptic magnetogram: an example force-free field produced is shown in Figure 9. They find that active region magnetic fields are reproduced quite rapidly, but convergence to the weaker fields in polar regions is much slower. Slow convergence is undesirable because numerical diffusion was found to erode the coronal currents before the photosphere reached the target configuration. However, the convergence rate could be improved by choosing the initial condition to be a dipole, approximating the average polar field in the magnetogram. This highlights an important feature of the method: non-uniqueness. The force-free field produced is not defined solely by the photospheric boundary condition, but depends both on (i) the choice of initialisation and (ii) the path followed to reach the final state. Contopoulos et al. (2011) suggest that one way to choose between possible solutions would be to incorporate measurements from vector magnetograms. The model described in the next section takes a different approach in treating the construction of the coronal magnetic field as an explicitly time-dependent problem.
Recently, van Ballegooijen et al. (2000) and Mackay and van Ballegooijen (2006) have developed a new technique to study the long-term evolution of coronal magnetic fields, which has now been applied to model the global solar corona (Yeates et al., 2007, 2008b). The technique follows the build-up of free magnetic energy and electric currents in the corona by coupling together two distinct models. The first is a data driven surface flux transport model (Yeates et al., 2007). This uses observations of newly emerging magnetic bipoles to produce a continuous evolution of the observed photospheric magnetic flux over long periods of time. Coupled to this is a quasi-static coronal evolution model (Mackay and van Ballegooijen, 2006; Yeates et al., 2008b) which evolves the coronal magnetic field through a sequence of nonlinear force-free fields in response to the observed photospheric evolution and flux emergence. The model follows the long-term continuous build-up of free magnetic energy and electric currents in the corona. It differs significantly from the extrapolation approaches which retain no memory of magnetic flux or connectivity from one extrapolation to the next.
The photospheric component of the model evolves the radial magnetic field on with a standard flux transport model (Section 2.2.1), except that newly emerging bipolar active regions are inserted not just in the photosphere but also in the 3D corona. These regions take an analytical form (Yeates et al., 2007), with parameters (location, size, flux, tilt angle) chosen to match observed active regions. An additional twist parameter allows the emerging 3D regions to be given a non-zero helicity: in principle this could be determined from vector magnetogram observations, but these are not yet routinely available.
The coronal part of the model evolves the large-scale mean field (van Ballegooijen et al., 2000) according to the induction equationet al. (2008b) assume the form
Rather than solving the full MHD equations for the velocity , which is not computationally feasible over long timescales, the quasi-static evolution of the coronal magnetic field is approximated using magneto-frictional relaxation (Yang et al., 1986), setting
Figure 10b shows an example nonlinear force-free field using this model, seen after 100 days of evolution. Figure 10a shows a PFSS extrapolation from the same photospheric distribution. The coronal field in the non-potential model is significantly different, comprising highly twisted flux ropes, slightly sheared coronal arcades, and near potential open field lines. Since this model can be run for extended periods without resetting the coronal field, it can be used to study long-term helicity transport across the solar surface from low to high latitudes (Yeates et al., 2008a).
There are several parameters in the model. As in the PFSS model the location of the upper boundary is arbitrary, although it has less influence since the non-potential magnetic field strength falls off more slowly with radius than that of the potential field. In the magneto-frictional evolution, the turbulent diffusivity and friction coefficient are arbitrary, and must be calibrated by comparison with observed structures or timescales for flux ropes to form or lose equilibrium. Finally a 3D model is required for newly-emerging active regions: the simple analytical bipoles of the existing simulations could in fact be replaced with more detailed extrapolations from observed photospheric fields in active regions.
Living Rev. Solar Phys. 9, (2012), 6
This work is licensed under a Creative Commons License.