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

The theoretical force-free magnetic fields in the first paper of this series, modeling magnetic configurations associated with polarity intrusion in active regions, are established to be all stable to linear ideal hydromagnetic perturbations under the boundary condition that anchors the lines of force rigidly to the photosphere. It is shown first that these force-free fields belong to an even larger class found by Chang and Carovillano (1981). A proof by the energy principle is then given to establish that all force-free magnetic fields in the larger class are absolutely stable. The physical implications of this result are discussed.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Influence of fluctuating field on the configuration of force-free field

The consistent equations for random force-free magnetic fields are obtained from the fluctuation and average equations of both momentum and induction ones based on the statistical approach. Some features of average force-free fields are discussed. This idea may be applied to explain many configurations of magnetic field in the astronomical environment.     

On axisymmetric hydromagnetic equilibria

The physical characteristics of possible axisymmetric equilibria are examined on the basis of the integrals of hydromagnetic equations. It is shown for nearly spherical configurations that a surface differential rotation is possible only in the absence of a meridional circulation with either purely toroidal or purely poloidal magnetic field. In the presence of a meridional circulation, it is shown that no surface rotation or constant rotation is possible if the magnetic field is purely toroidal, and that no rotation is possible if the magnetic field is purely poloidal. A brief discussion is given on the possible solutions including the case of stellar winds with force-free magnetic fields.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Non-ideal evolution of non-axisymmetric, force-free magnetic fields in a magnetar

Recent numerical magnetohydrodynamic calculations by Braithwaite and collaborators support the 'fossil field' hypothesis regarding the origin of magnetic fields in compact stars and suggest that the resistive evolution of the fossil field can explain the reorganization and decay of magnetar magnetic fields. Here, these findings are modelled analytically by allowing the stellar magnetic field to relax through a quasi-static sequence of non-axisymmetric, force-free states, by analogy with spheromak relaxation experiments, starting from a random field. Under the hypothesis that the force-free modes approach energy equipartition in the absence of resistivity, the output of the numerical calculations is semiquantitatively recovered: the field settles down to a linked poloidal–toroidal configuration, which inflates and becomes more toroidal as time passes. A qualitatively similar (but not identical) end state is reached if the magnetic field evolves by exchanging helicity between small and large scales according to an α-dynamo-like, mean-field mechanism, arising from the fluctuating electromotive force produced by the initial random field. The impossibility of matching a force-free internal field to a potential exterior field is discussed in the magnetar context.     

Force-free magnetic fields and flares of August 1972

On the basis of observations (Zirin and Tanaka, 1973) inferring the presence of shear in magnetic fields, the amount of extractable energy stored in a class of force-free magnetic fields is evaluated for the flares of August 1972, using the formulations developed by Nakagawa and Raadu (1972). It is shown that the evaluated energy storage could be built up by the proper motions of sunspots in the active region McMath 11976 during July 31 and August 7. Then for the flare of August 7, a detailed analysis is made of the manner of energy release in the post maximum phases deduced from the configuration of flare loops. It is shown that the observed flare loops could be represented closely by the force-free magnetic fields and that the evaluated rate of energy release is consistent with observed rate given by the soft X-ray emission. The results of analysis suggest that the flare of August 1972 could be identified with the relaxation of an energetic force-free magnetic field towards lower energy states. The limitations and possible future extension of this type of analysis are discussed.Visiting scientist from the Tokyo Astronomical Observatory, Mitaka, Tokyo, Japan.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Differential rotation and the structure and energy content of coronal magnetic fields

It is argued that differential rotation of the photospheric magnetic fields will induce currents in the corona. The work done against surface magnetic stresses will increase the energy content of the coronal magnetic field. The electrical conductivities are high and the foot points of field lines move with the differential rotation. The force-free field equations are solved with this constraint to obtain a minimum estimate of the energy increase for a quadrupole field. During a solar rotation the magnetic energy increases by 25%. Local release of this energy in the corona would have a significant effect. The expansion of field lines as a result of the differential rotation should increase the amount of flux and the field strength in the solar wind region.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.     

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.     

Computation of solar magnetic fields from photospheric observations

The observational difficulties of obtaining the magnetic field distribution in the chromosphere and corona of the Sun has led to methods of extending photospheric magnetic measurements into the solar atmosphere by mathematical procedures. A new approach to this problem presented here is that a constant alpha force-free field can be uniquely determined from the tangential components of the measured photospheric flux alone. The vector magnetographs now provide measurements of both the solar photospheric tangential and the longitudinal magnetic field. This paper presents derivations for the computation of the solar magnetic field from these type of measurements. The fields considered are assumed to be a constant alpha force-free fields or equivalent, producing vanishing Lorentz forces. Consequently, magnetic field lines and currents are related by a constant and hence show an identical distribution. The magnetic field above simple solar regions are described from the solution of the field equations.     

Topology of Magnetic Field and Coronal Heating in Solar Active Regions

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

Calculating and Testing Nonlinear Force-Free Fields

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

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.     

Computing Nonlinear Force-Free Coronal Magnetic Fields in Spherical Geometry

We describe a newly developed code for the extrapolation of nonlinear force-free coronal magnetic fields in spherical coordinates. The program uses measured vector magnetograms on the solar photosphere as input and solves the force-free equations in the solar corona. The method is based on an optimization principle and the heritage of the newly developed code is a corresponding method in Cartesian geometry. We test the newly developed code with the help of a semi-analytic solution and rate the quality of our reconstruction qualitatively by magnetic field line plots and quantitatively with a number of comparison metrics. We find that we can reconstruct the original test field with high accuracy. The method is fast if the computation is limited to low co-latitudes (say 30°≤θ≤150°), but it becomes significantly slower if the polar regions are included.     

On practical representation of magnetic field

Various manners of determination of a magnetic field are reviewed briefly from the standpoint of practicality and uniqueness. Then a practical representation of magnetic fields in terms of a class of force-free magnetic field is described. The proposed scheme is based on the physical consideration that in the chromosphere and lower corona a quasistatic magnetic field must be nearly force-free and that for the class of force-free magnetic field, i.e., ×B=B with = constant, the magnetic field can be determined uniquely from the observed distribution of the vertical component of a magnetic field. The applicability of the representation is demonstrated by examples and the limitations are discussed.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

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.     

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.     

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.     

Suppression of the kink instability for magnetic flux ropes in the chromosphere

Energy storage in chromospheric flux ropes is discussed, in the context of solar flares. The structure is represented by a cylindrically symmetric magnetic field of finite length. The field is assumed to be approximately force-free. The stability of the field to a kink perturbation is investigated. Flux ropes are rooted in dense photospheric material. So the ends of the field lines are taken to be fixed on rigid boundaries for all perturbations. An energy perturbation method is used and the boundary conditions give a stabilizing effect. It is shown that for a moderate degree of twisting the fields are stable to a kink perturbation. Thus energy can be stored in cylindrical fields prior to release in a solar flare.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Comparing Values of the Relative Magnetic Helicity in Finite Volumes

Relative magnetic helicity, as a conserved quantity of ideal magnetohydrodynamics, has been highlighted as an important quantity to study in plasma physics. Due to its nonlocal nature, its estimation is not straightforward in both observational and numerical data. In this study we derive expressions for the practical computation of the gauge-independent relative magnetic helicity in three-dimensional finite domains. The derived expressions are easy to implement and rapid to compute. They are derived in Cartesian coordinates, but can be easily written in other coordinate systems. We apply our method to a numerical model of a force-free equilibrium containing a flux rope, and compare the results with those obtained employing known half-space equations. We find that our method requires a much smaller volume than half-space expressions to derive the full helicity content. We also prove that values of relative magnetic helicity of different magnetic fields can be compared with each other in the same sense as free-energy values can. Therefore, relative magnetic helicity can be meaningfully and directly compared between different datasets, such as those from different active regions, but also within the same dataset at different times. Typical applications of our formulae include the helicity computation in three-dimensional models of the solar atmosphere, e.g., coronal-field reconstructions by force-free extrapolation and discretized magnetic fields of numerical simulations.     

Including stereoscopic information in the reconstruction of coronal magnetic fields

We present a method to include stereoscopic information about the three-dimensional structure of flux tubes into the reconstruction of the coronal magnetic field. Due to the low plasma beta in the corona we can assume a force-free magnetic field, with the current density parallel to the magnetic field lines. Here we use linear force-free fields for simplicity. The method uses the line-of-sight magnetic field on the photosphere as observational input. The value of is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force-free solution constrained by both the photospheric magnetogram and the observed plasma structures. As an example we apply our method to SOHO MDI/EIT data of an active region. In the future it is planned to apply the method to analyse data from the SECCHI instrument aboard the STEREO mission.