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.     

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.     

Magnetic stars with an external non-linear force-free field

The possible existence of strong magnetic fields in stars is discussed and a method of constructing highly distorted models of magnetic, rotating stars developed. For stars with both poloidal and toroidal fields at the surface a force-free outer boundary condition is necessary. Non-linear solutions of the force-free equations must be used. The force-free equations and the structure equations for a white dwarf are solved simultaneously by a finite difference method.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

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.     

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.     

Force-free magnetic fields on curvilinear coordinate surfaces

To seek nonlinear solutions of force-free magnetic fields, some symmetries or approximations are usually invoked. We consider magnetic fields lying on coordinate surfaces of an orthogonal curvilinear coordinate system. We conclude that only fields on parallel planes or spherical shells can be expressed in the form provided by Low in 1980s. These force-free fields are stable against small perturbations with rigid boundaries. Fields on cylindrical shells are also considered.     

Force-free magnetic fields with singular current density surfaces and their stability

Starting from Bernstein's principle of magnetohydrodynamic energy, a general analysis is presented for the stability of a kind of 1-D force-free magnetic fields with singular current density surfaces and a single parameter in cylindrical coordinates. It is found that in the parameter space of this kind of force-free magnetic fields there simultaneously exist stable and unstable regions. Their stability is solely determined by the radial distribution of the magnetic pitch in the neighborhood of the cylinder axis, and is independent of the presence of singular current density surface at the boundary of the field.     

Equilibrium and stability of force-free magnetic field

Force-free magnetic fields (f.f.f) are considered as the first approximation of magnetic hydrodynamic equations in the case when the energy of the field exceeds the thermal energy of the medium. Such a relation of energies takes place in the upper atmosphere of the Sun in active regions.The consequence of the virial theorem obtained shows that for any solution of the corresponding non-linear system of equations only two cases are possible: either the total energy of the field is given by a divergent integral, or in some regions the force-free character of field is destroyed. This permits the conclusion that it is impossible to build f.f. current systems everywhere, and therefore boundary problems for this type of fields are of the same importance as for harmonic fields.Integral relations are obtained which are the necessary conditions for the solution of boundary problems. According to the classical principle of Thompson the harmonic fields are always stable, while f.f.f. may be stable or unstable.It is shown that: (1) arbitrary f.f.f. are stable to small changes of boundary conditions; (2) among f.f.f. the hydrodynamically stable configurations exist.The hydrodynamic stability condition restricts the size of force-free currents in such configurations.     

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.     

The topological association of Hα structures and magnetic fields

The topological associations between H structures and magnetic fields are examined for an active region observed on two different dates. The structures seen in the on and off band of H filtergrams are compared with the contour maps of magnetic fields at the level of magnetogram observations. Similar comparisons are made also with the configurations of force-free magnetic lines of force at various heights evaluated with the use of formulations developed previously by Nakagawa and Raadu (1972).Among the results of significance, we may note that (1) H plages could be identified with regions of magnetic field larger than ±80 G, (2) the network of bright dots seen in H -1 Å filtergrams follow closely ±80 G contours. (3) stable prominences lie along either neutral lines or valleys of magnetic fields, (4) the configuration of magnetic lines of force shows discrete domain structures suggesting bipolar nature of local magnetic fields, and (5) within a domain the configuration is governed apparently by evalutional consequences. Details of analyses are described with discussions on the limitations and possible future improvements.The National Center for Atmospheric Research is sponsored by the National Science Foundation.Operated by Association of Universities for Research in Astronomy, Inc., under contract with 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.     

Topological differences and similarities between force-free and potential models of coronal magnetic fields

Potential fields and linear force-free fields are often used as models for the magnetic field of the Sun's corona. They can be written as analytical expressions in terms of boundary values at the photosphere. Because of their relative simplicity compared with nonlinear force-free fields, these two models are of particular importance in topological analysis of solar phenomena. However, it has been suggested by Hudson and Wheatland (1999) that the topologies of potential and force-free models are in general not even qualitatively equivalent. In this paper, their example is re-examined and it is found that the opposite conclusions hold. In general, potential and force-free fields are topologically similar sufficiently close to localized sources. The exception to this are structurally unstable states, such as bifurcation states, where a small change of current can produce a significant change of topology.     

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.     

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.     

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.     

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.     

Magnetodynamic stability of a fluid cylinder under the lundquist force-free magnetic field

The magnetodynamic (in)stability of a conducting fluid cylinder subject to the capillarity and electromagnetic forces has been developed. The cylinder is pervaded by a uniform magnetic field but embedded in the Lundquist force-free varying field that allows for flowing a current surrounding the fluid. A general eigenvalue relation is derived based on a study of the equilibrium and perturbed states. The stability criterion is discussed analytically in general terms. The surface tension is destabilizing for small axisymmetric mode and stable for all others. The principle of the exchange of stability is allowed for the present problem due to the non-uniform behaviour of the force-free field. Each of the axial and transverse force-free fields separately exerts a stabilizing influence in the most dangerous mode but the combined contribution of them is strongly destabilizing. Whether the model is acted upon the electromagnetic force (with the Lundquist field) the stability restrictions or/and the capillarity force are identified.Several reported works can be recovered as limiting cases with appropriate simplifications.     

A Test on Two Codes for Extrapolating Solar Linear Force-free Magnetic Fields

We evaluate the performances of two newly-implemented codes for extrapolating the solar linear force-free magnetic fields, by measuring their quantified responses to the lower boundary vector field data on a finite region using analytical models. The codes are based on two boundary integral formulas with different mechanisms in utilizing the transverse boundary field: the first one only utilizes the transverse boundary field to derive the value of the force-free parameter, while the other one explicitly utilizes the whole transverse boundary field, in addition to the vertical field component. Studies on the test cases show that both of the codes could reproduce the analytical model fields with reasonable accuracy within the valid domain, provided a sufficient amount of data is available at the lower boundary. The code utilizing explicitly all three components of the boundary field shows generally better performances and requires a smaller boundary-data area in order to achieve the same degree of accuracy in the metric quantities such as the normalized vector error, vector correlation, vector angular difference, and magnetic energy; however, the accuracy in the divergence-free metric and especially the force-free metric are less ideal.