Force-free magnetic fields in the solar atmosphere

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

RECONSTRUCTION OF THE THREE-DIMENSIONAL CORONAL MAGNETIC FIELD

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

Computation of inner coronal magnetic fields from longitudinal field components on a spherical photosphere

A method is developed which for a certain day permits the approximate calculation of closed small and large scale magnetic field lines. From the photospheric longitudinal components of the magnetic field measured at this day normal components are derived taking into account the curvature of the solar surface. The magnetic fields are assumed to be potential or force-free fields.The method is applied to observations of September 5 and September 7, 1973. The projected magnetic field lines are compared with the loop structures which are visible in XUV pictures taken on these days. In the cases where no good agreement could be obtained for potential fields, force-free fields are calculated and fitted to the observed structures.     

Force-free magnetic arcades relevant to two-ribbon solar flares

Simple analytic models for the passive evolution of arcade-like magnetic fields through a series of force-free equilibria are presented. At the photospheric boundary, the normal magnetic field component is prescribed together with either the longitudinal field component or the photospheric shear. Analytic progress is made by considering either cylindrically symmetric solutions or using the separation of variables technique. Two distinct cylindrically symmetric force-free fields are obtained that possess the same normal field component and photospheric shear. The scond field contains a magnetic bubble. As the shear increases beyond a critical value, so the magnetic energy of the first configuration exceeds that of the second. The possibility is therefore suggested of an eruption of the first field outwards towards the second. Such an eruptive instability is proposed as the origin of a two-ribbon solar flare.A new analytic solution to the force-free field equations, of separable form, is discovered and it is pointed out that the existence of shear in a magnetic field does not preclude it from being potential.Now at AWRE, Aldermaston, Reading, Berkshire.     

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.     

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.     

CAN WE EXTRAPOLATE A MAGNETIC FIELD WHEN ITS TOPOLOGY IS COMPLEX?

In order to understand various solar phenomena controlled by the magnetic field, such as X-ray bright points, flares and prominence eruptions, the structure of the coronal magnetic field must be known. This requires a precise extrapolation of the photospheric magnetic field. Presently, only potential or linear force-free field approximations can be used easily. A more realistic modelling of the field is still an active research area because of well-known difficulties related to the nonlinear mixed elliptic-hyperbolic nature of the equations. An additional difficulty arises due to the complexity of the magnetic field structure which is caused by a discrete partition of the photospheric magnetic field. This complexity is not limited to magnetic regions having magnetic nulls (and so separatrices) but also occurs in those containing thin elongated volumes (called Quasi-Separatrix Layers) where the photospheric field-line linkage changes rapidly. There is a wide range for the thickness of such layers, which is determined by the character (bipolar or quadrupolar) of the magnetic region, by the sizes of the photospheric field concentrations and by the intensity of the electric currents. The aim of this paper is to analyse the recent nonlinear force-free field extrapolation techniques for complex coronal magnetic fields.     

Global Solar Free Magnetic Energy and Electric Current Density Distribution of Carrington Rotation 2124

Solar eruptive phenomena, like flares and coronal mass ejections (CMEs), are governed by magnetic fields. To describe the structure of these phenomena one needs information on the magnetic flux density and the electric current density vector components in three dimensions throughout the atmosphere. However, current spectro-polarimetric measurements typically limit the determination of the vector magnetic field to only 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 for global solar atmosphere using nonlinear force-free field (NLFFF) extrapolation codes implemented to a synoptic maps of photospheric vector magnetic field synthesized from the Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) as boundary condition. Using the resulting three-dimensional magnetic field, we calculate the three-dimensional electric current density and magnetic energy throughout the solar atmosphere for Carrington rotation 2124 using our global extrapolation code. We found that spatially, the low-lying, current-carrying core field demonstrates a strong concentration of free energy in the active-region core, from the photosphere to the lower corona (about 70 Mm). The free energy density appears largely co-spatial with the electric current distribution.     

Force-free Magneitic Field Extrapolation for the Complex Sunspot Group of August 1972

Quantitative data on the magnetic field structure at all levels of the solar atmosphere are of basic importance for our understanding of physical processes in active regions in general and in flares in particular. Because photospheric longitudinal magnetograms are the most reliable data available, one has to look for a method of theoretical extrapolation of such data to higher levels. Such a method has been developed for force-free magnetic fields, i.e., ∇ X B = αB , with α = constant, satisfying more realistic boundary conditions as compared with earlier papers; e.g., the net magneitic flux through the magnetogram area is not required to be zero. The method has been used to calculate the magnetic field vector and lines of force in the flare-active region of August 1972. Calculated fields are compared with other observations such as structures in H α . Results of August 3rd show that the loop prominence systems observed during the flares of August 2nd and 4th are represented by a force-free field with positive α rather than by a current-free field ( α = 0) The extractable energy supply of this force-free field is of the order of magnitude of maximum flare demand (10³² erg); the height dependence of the magnetic field strengths agrees with that from radio and X-ray estimates. Similar results are obtained for the August 7th magnetic field structure.     

How To Use Magnetic Field Information For Coronal Loop Identification

The structure of the solar corona is dominated by the magnetic field because the magnetic pressure is about four orders of magnitude higher than the plasma pressure. Due to the high conductivity the emitting coronal plasma (visible, e.g., in SOHO/EIT) outlines the magnetic field lines. The gradient of the emitting plasma structures is significantly lower parallel to the magnetic field lines than in the perpendicular direction. Consequently information regarding the coronal magnetic field can be used for the interpretation of coronal plasma structures. We extrapolate the coronal magnetic field from photospheric magnetic field measurements into the corona. The extrapolation method depends on assumptions regarding coronal currents, e.g., potential fields (current-free) or force-free fields (current parallel to magnetic field). As a next step we project the reconstructed 3D magnetic field lines on an EIT-image and compare with the emitting plasma structures. Coronal loops are identified as closed magnetic field lines with a high emissivity in EIT and a small gradient of the emissivity along the magnetic field.     

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

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

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.     

Magnetic field configurations associated with polarity intrusion in a solar active region

This paper presents a new class of exact solutions describing the non-linear force-free field above a spatially localized photospheric bipolar magnetic region. An essential feature is the variation in all three Cartesian directions and this could not be modelled adequately with previously known symmetric force-free fields. Sequences of force-free fields are constructed and analyzed to simulate the slow growth of a pair of spots on the photosphere. The axis connecting the spots executes rotational motion, distorting the photospheric neutral line separating fluxes of opposite signs. We show directly from the analytic solutions that the resulting reversal of the positions of the spots relative to the background field is associated with (i) the creation of magnetic free energy, (ii) the severe shearing of localized low-lying loops in the vicinity where the photospheric transverse field aligns with the photospheric neutral line, and (iii) the emergence and disappearance of flux from the photosphere at these highly stressed regions. The model relates theoretically for the first time these different magnetic field features that have been suggested by observation and theoretical considerations to be flare precursors. A general formula, based on the virial theorem, is also given for the free energy of a force-free field, strictly in terms of the field value at the photosphere. This formula has obvious practical application.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

A class of force-free magnetic fields for modeling pre-flare coronal magnetic configurations

Three-dimensional, quasi-static evolutions of coronal magnetic fields driven by photospheric flux emergence are modeled by a class of analytic force-free magnetic fields. Our models relate commonly observed photospheric magnetic phenomena, such as the formation and growth of sunspots, the emergence of an X-type separator, and the collision and merging of sunspots, to the three-dimensional magnetic fields in the corona above. By tracking the evolution in terms of a continuous sequence of force-free states, we show that flux emergence and submergence along magnetic neutral lines in the photosphere are essential processes in all these photospheric phenomena. The analytic solutions we present have a parametric regime within which the magnetic energy attained by an evolving force-free field may be of the order of 10³⁰ ergs to several 10³¹ ergs, depending on the magnetic environment into which an emerging flux intrudes. The commonly used indicators of magnetic shear in magnetogram interpretation are discussed in terms of field connectivity in our models. It is demonstrated that the crossing angle of the photospheric transverse magnetic field with the neutral line may not be a reliable indicator of the magnetic shear in the coronal field above, due to the complexity of three-dimensionality. The poorly understood constraint of magnetic-helicity conservation on the availability of magnetic free energy for a flare is briefly discussed.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

The calculation of force-free fields from discrete flux distributions

This paper presents particularly simple mathematical formulas for the calculation of force-free fields of constant α from the distribution of discrete sources on a flat surface. The advantage of these formulas lies in their physical simplicity and the fact that they can be easily used in practice to calculate the fields. The disadvantage is that they are limited to fields of ‘sufficiently small α’. These formulas may be useful in the study of chromospheric magnetic fields by the comparison of high-resolution Hα photographs and photospheric magnetograms.     

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

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

Magnetic Fields and Solar Activities

We discuss the study of solar magnetic fields based on the photospheric vector magnetograms of solar active regions which were obtained at Huairou Solar Observing Station near Beijing in the period of 22nd and 23th solar cycles. The measurements of the chromospheric magnetic field and the spatial configuration of the field at the lower solar atmosphere inferred by the distribution of the solar photospheric and chromospheric magnetic field. After the analysis on the formation process of delta configuration in some super active regions based on the photospheric vector magnetogram observations, some results are obtained: (1) The analysis of magnetic writhe of whole active regions cannot be limited in the strong field of sunspots, because the contribution of the fraction of decayed magnetic field is non-negligible. (2) The magnetic model of kink magnetic ropes, proposed to be generated in the subatmosphere, is not consistent with the evolution of large-scale twisted photospheric transverse magnetic field and the relationship with magnetic shear in some delta active regions completely. (3) The proposition is that the large-scale delta active regions are formed from contribution by highly sheared non-potential magnetic flux bundles generated in the subatmosphere. We present some results of a study of the magnetic helicity. We also compare these results with other data sets obtained by magnetographs (or Stokes polarimeters) at different observatories, and analyze the basic chirality of the magnetic field in the solar atmosphere.     

Magnetic diffusion and flare energy buildup

Photospheric motion shears or twists solar magnetic fields to increase magnetic energy in the corona, because this process may change a current-free state of a coronal field to force-free states which carry electric current. This paper analyzes both linear and nonlinear two-dimensional force-free magnetic field models and derives relations of magnetic energy buildup with photospheric velocity field. When realistic data of solar magnetic field (B ₀ 10³ G) and photospheric velocity field (v _max 1 km s^–1) are used, it is found that 3–4 hours are needed to create an amount of free magnetic energy which is of the order of the current-free field energy. Furthermore, the paper studies situations in which finite magnetic diffusivities in photospheric plasma are introduced. The shearing motion increases coronal magnetic energy, while the photospheric diffusion reduces the energy. The variation of magnetic energy in the coronal region, then, depends on which process dominates.     

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.     

Magnetic Energy Storage and Current Density Distributions for Different Force-Free Models

In the last decades, force-free-field modelling has been used extensively to describe the coronal magnetic field and to better understand the physics of solar eruptions at different scales. Especially the evolution of active regions has been studied by successive equilibria in which each computed magnetic configuration is subject to an evolving photospheric distribution of magnetic field and/or electric-current density. This technique of successive equilibria has been successful in describing the rate of change of the energetics for observed active regions. Nevertheless the change in magnetic configuration due to the increase/decrease of electric current for different force-free models (potential, linear and nonlinear force-free fields) has never been studied in detail before. Here we focus especially on the evolution of the free magnetic energy, the location of the excess of energy, and the distribution of electric currents in the corona. For this purpose, we use an idealised active region characterised by four main polarities and a satellite polarity, allowing us to specify a complex topology and sheared arcades to the coronal magnetic field but no twisted flux bundles. We investigate the changes in the geometry and connectivity of field lines, the magnetic energy and current-density content as well as the evolution of null points. Increasing the photospheric current density in the magnetic configuration does not dramatically change the energy-storage processes within the active region even if the magnetic topology is slightly modified. We conclude that for reasonable values of the photospheric current density (the force-free parameter α<0.25 Mm⁻¹), the magnetic configurations studied do change but not dramatically: i) the original null point stays nearly at the same location, ii) the field-line geometry and connectivity are slightly modified, iii) even if the free magnetic energy is significantly increased, the energy storage happens at the same location. This extensive study of different force-free models for a simple magnetic configuration shows that some topological elements of an observed active region, such as null points, can be reproduced with confidence only by considering the potential-field approximation. This study is a preliminary work aiming at understanding the effects of electric currents generated by characteristic photospheric motions on the structure and evolution of the coronal magnetic field.