Analysis of Reconstruction Methods for Nonlinear Force-Free Fields

We focus on quantitative evaluation of several methods for the reconstruction of force-free magnetic fields in the solar corona. We have studied two topics. The first is to carry out a comparison test of two approaches to implement the optimization method: (1) using a boundary layer with a weighting function and (2) involving an optimization algorithm for unknown boundary field vector values. The second approach is shown to provide a better approximation to a desired true solution that is finite in an unlimited semispace. The approximation obtained by the second approach is reasonably close to the true solution up to the boundaries of the reconstruction region. Among the applications of the optimization method to real data, we demonstrate its ability to obtain accurate enough energy estimates and find that the pre- to post-flare energy difference is sufficient in powering the flare.     

Magnetic Connectivity Between Active Regions 10987, 10988, and 10989 by Means of Nonlinear Force-Free Field Extrapolation

Extrapolation codes for modelling the magnetic field in the corona in Cartesian geometry do not take the curvature of the Sun’s surface into account and can only be applied to relatively small areas, e.g., a single active region. We apply a method for nonlinear force-free coronal magnetic field modelling of photospheric vector magnetograms in spherical geometry which allows us to study the connectivity between multi-active regions. We use Vector Spectromagnetograph (VSM) data from the Synoptic Optical Long-term Investigations of the Sun (SOLIS) survey to model the coronal magnetic field, where we study three neighbouring magnetically connected active regions (ARs 10987, 10988, 10989) observed on 28, 29, and 30 March 2008, respectively. We compare the magnetic field topologies and the magnetic energy densities and study the connectivities between the active regions. We have studied the time evolution of the magnetic field over the period of three days and found no major changes in topologies, as there was no major eruption event. From this study we have concluded that active regions are much more connected magnetically than the electric current.     

Nonlinear Force-Free Extrapolation of Emerging Flux with a Global Twist and Serpentine Fine Structures

We study the flux emergence process in NOAA active region 11024, between 29 June and 7 July 2009, by means of multi-wavelength observations and nonlinear force-free extrapolation. The main aim is to extend previous investigations by combining, as much as possible, high spatial resolution observations to test our present understanding of small-scale (undulatory) flux emergence, whilst putting these small-scale events in the context of the global evolution of the active region. The combination of these techniques allows us to follow the whole process, from the first appearance of the bipolar axial field on the east limb, until the buoyancy instability could set in and raise the main body of the twisted flux tube through the photosphere, forming magnetic tongues and signatures of serpentine field, until the simplification of the magnetic structure into a main bipole by the time the active region reaches the west limb. At the crucial time of the main emergence phase high spatial resolution spectropolarimetric measurements of the photospheric field are employed to reconstruct the three-dimensional structure of the nonlinear force-free coronal field, which is then used to test the current understanding of flux emergence processes. In particular, knowledge of the coronal connectivity confirms the identity of the magnetic tongues as seen in their photospheric signatures, and it exemplifies how the twisted flux, which is emerging on small scales in the form of a sea-serpent, is subsequently rearranged by reconnection into the large-scale field of the active region. In this way, the multi-wavelength observations combined with a nonlinear force-free extrapolation provide a coherent picture of the emergence process of small-scale magnetic bipoles, which subsequently reconnect to form a large-scale structure in the corona.     

Comparison of Nonlinear Force-Free Field and Potential Field in the Quiet Sun

In this paper, a potential field extrapolation and three nonlinear force-free (NLFF) field extrapolations (optimization, direct boundary integral (DBIE), and approximate vertical integration (AVI) methods) are used to study the spatial configuration of magnetic field in the quiet Sun. It is found that differences in the computed field strengths among the three NLFF and potential fields exist in the low layers. However, they tend to disappear as the height increases, and the differences are of the order of 0.1 gauss when the height exceeds ≈ 2000 km above the photosphere. The difference in azimuth angles between each NLFF field model and the potential field is as follows: for the optimization field, it decreases evidently as the height increases; for the DBIE field, it almost stays constant and shows no significant change as the height increases; for the AVI field, it increases slowly as the height increases. Our analysis shows that the reconstructed NLFF fields deviate significantly from the potential field in the quiet Sun.     

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

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

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

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

A Comparison Between Nonlinear Force-Free Field and Potential Field Models Using Full-Disk SDO/HMI Magnetogram

Measurements of magnetic fields and electric currents in the pre-eruptive corona are crucial to the study of solar eruptive phenomena, like flares and coronal mass ejections (CMEs). However, spectro-polarimetric measurements of certain photospheric lines permit a determination of the vector magnetic field only at the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field above multiple active regions with the help of a potential field and a nonlinear force-free field (NLFFF) extrapolation code over the full solar disk using Helioseismic and Magnetic Imager (SDO/HMI) data as boundary conditions. We compare projections of the resulting magnetic field lines with full-disk coronal images from the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study has found that the NLFFF model reconstructs the magnetic configuration closer to observation than the potential field model for full-disk magnetic field extrapolation. We conclude that many of the trans-equatorial loops connecting the two solar hemispheres are current-free.     

The Force-Free Electrodynamics Method for the Extrapolation of Coronal Magnetic Fields from Vector Magnetograms

We present a new improved version of our force-free electrodynamics (FFE) numerical code in spherical coordinates that extrapolates the magnetic field in the inner solar corona from a photospheric vector magnetogram. The code satisfies the photospheric boundary condition and the condition ??B=0 to machine accuracy. The performance of our method is evaluated with standard convergence parameters, and is found to be comparable to that of other nonlinear force-free extrapolations.     

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. I. Theory

We derive an analytical approximation of nonlinear force-free magnetic field solutions (NLFFF) that can efficiently be used for fast forward-fitting to solar magnetic data, constrained either by observed line-of-sight magnetograms and stereoscopically triangulated coronal loops, or by 3D vector-magnetograph data. The derived NLFFF solutions provide the magnetic field components B _x (x), B _y (x), B _z (x), the force-free parameter α(x), the electric current density j(x), and are accurate to second-order (of the nonlinear force-free α-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.     

无力场模型黑子上空的纵向磁场梯度

本文按常α无力场模型计算了1980年10月23日Boulder 2744活动区前导黑子的纵向磁场随高度的变化,并与用CIV 1548谱线观测得到的色球一日冕过渡区的磁场资料相结合,求得CIV 1548发射区的有效高度。这些结果与文献[4]中对同一黑子用势场模型推求的结果有很大差别。从而表明,势场和无力场在某些方面导致的结果是极不相同的。鉴于观测已表明活动区上空存在电流的事实,在活动区磁场的模拟中,特别是在强扭曲活动区磁场的计算中,应当避免采用势场,而尽可能采用无力场模型。     

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

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

Three-Dimensional Nonlinear Force-Free Field Reconstruction of Solar Active Region 11158 by Direct Boundary Integral Equation

A three-dimensional coronal magnetic field is reconstructed for the NOAA active region 11158 on 14 February 2011. A GPU-accelerated direct boundary integral equation (DBIE) method is implemented which is approximately 1000 times faster than the original DBIE used on solar non-linear force-free field modeling. Using the SDO/HMI vector magnetogram as the bottom boundary condition, the reconstructed magnetic field lines are compared with the projected EUV loop structures as observed in the front-view (SDO/AIA) and the side-view (STEREO-A/B) images for the first time; they show very good agreement three-dimensionally. A quantitative comparison with some stereoscopically reconstructed coronal loops shows that the average misalignment angles in our model are at the same order as the state-of-the-art results obtained from reconstructed coronal loops. It is found that the observed coronal loop structures can be grouped into a number of closed and open field structures with some central bright coronal loop features around the polarity inversion line. The reconstructed highly sheared magnetic field lines agree very well with the low-lying sigmoidal filament along the polarity inversion line. This central low-lying magnetic field loop system must have played a key role in powering the flare. It should be noted that while a strand-like coronal feature along the polarity inversion line may be related to the filament, one cannot simply interpret all the coronal bright features along the polarity inversion line as manifestation of the filament without any stereoscopic information.     

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

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

Nonlinear Force-Free Magnetic Field Extrapolations: Comparison of the Grad Rubin and Wheatland Sturrock Roumeliotis Algorithm

We compare the performance of two alternative algorithms which aim to construct a force-free magnetic field given suitable boundary conditions. For this comparison, we have implemented both algorithms on the same finite element grid which uses Whitney forms to describe the fields within the grid cells. The additional use of conjugate gradient and multigrid iterations result in quite effective codes. The Grad Rubin and Wheatland Sturrock Roumeliotis algorithms both perform well for the reconstruction of a known analytic force-free field. For more arbitrary boundary conditions the Wheatland Sturrock Roumeliotis approach has some difficulties because it requires overdetermined boundary information which may include inconsistencies. The Grad Rubin code on the other hand loses convergence for strong current densities. For the example we have investigated, however, the maximum possible current density seems to be not far from the limit beyond which a force-free field cannot exist anymore for a given normal magnetic field intensity on the boundary.     

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. III. The Free Energy

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time–elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east–west asymmetry in its brightness and degree of polarisation. Therefore, measurement of the interplanetary polarised brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.     

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests

Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-fitting to six potential field and six nonpotential field cases simulated with our analytical model, as well as by forward-fitting to an exactly force-free solution of the Low and Lou (Astrophys. J. 352, 343, 1990) model. The forward-fitting tests demonstrate: i) a satisfactory convergence behavior (with typical misalignment angles of μ≈1^°?–?10^°), ii) relatively fast computation times (from seconds to a few minutes), and iii) the high fidelity of retrieved force-free α-parameters (α _fit/α _model≈0.9?–?1.0 for simulations and α _fit/α _model≈0.7±0.3 for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-force-free magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing nonlinear force-free field (NLFFF) codes based on photospheric magnetograms without coronal constraints.     

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

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

太阳线性无力场评注和快速傅氏分析法的应用

首先对30年来用无力场计算太阳活动区磁场方法作了一些评注，特别研究了线性无力场解中的几个误区，提出了准线性无力场方法并探讨了快速傅氏变换在线性无力场计算中的应用，以2000年9月15日NOAA 9165活动区磁场为例，应用3种不同的方法计算出磁力线走向，然后同紫金山天文台赣榆观测站观测的精细Hα图形对比。     

A Non-Linear Force-Free Field Model for the Evolving Magnetic Structure of Solar Filaments

In this paper the effect of a small magnetic element approaching the main body of a solar filament is considered through non-linear force-free field modeling. The filament is represented by a series of magnetic dips. Once the dips are calculated, a simple hydrostatic atmosphere model is applied to determine which structures have sufficient column mass depth to be visible in Hα. Two orientations of the bipole are considered, either parallel or anti-parallel to the overlying arcade. The magnetic polarity that lies closest to the filament is then advected towards the filament. Initially for both the dominant and minority polarity advected elements, right/left bearing barbs are produced for dextral/sinsitral filaments. The production of barbs due to dominant polarity elements is a new feature. In later stages the filament breaks into two dipped sections and takes a highly irregular, non-symmetrical form with multiple pillars. The two sections are connected by field lines with double dips even though the twist of the field is less than one turn. Reconnection is not found to play a key role in the break up of the filament. The non-linear force-free fields produce very different results to extrapolated linear-force free fields. For the cases considered here the linear force-free field does not produce the break up of the filament nor the production of barbs as a result of dominant polarity elements.     

Comments on “Comments on Solar Linear Force-Free Field and Application of FFT Analysis”

The paper “Comments on Solar Linear Force-free Field and Application of FFT Analysis”^[19] made some comments on certain solar linear force-free field works and accordingly proposed an improved method of fast Fourier transform (FFT). Some viewpoints of this paper are worthy of reconsideration and it is necessary to clarify certain questions existing in it. We would also like to make some comments on its proposed method from the theoretical point of view.