An optimization principle for the computation of MHD equilibria in the solar corona

AIMS: We develop an optimization principle for computing stationary MHD equilibria. METHODS: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma uses non-force-free MHD equilibria. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements. The program uses photospheric vector magnetograms and coronal EUV images as input. We tested our reconstruction code with the help of a semi-analytic MHD-equilibrium. The quality of the reconstruction was judged by comparing the exact and reconstructed solution qualitatively by magnetic field-line plots and EUV-images and quantitatively by several different numerical criteria. RESULTS: Our code is able to reconstruct the semi-analytic test equilibrium with high accuracy. The stationary MHD optimization code developed here has about the same accuracy as its predecessor, a non-linear force-free optimization code. The computing time for MHD-equilibria is, however, longer than for force-free magnetic fields. We also extended a well-known class of nonlinear force-free equilibria to the non-force-free regime for purposes of testing the code. CONCLUSIONS: We demonstrate that the code works in principle using tests with analytical equilibria, but it still needs to be applied to real data.Comment: 6 pages, 3 figure

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.Comment: 23 pages, 5 figure

Can We Improve the Preprocessing of Photospheric Vector Magnetograms by the Inclusion of Chromospheric Observations?

The solar magnetic field is key to understanding the physical processes in the solar atmosphere. Nonlinear force-free codes have been shown to be useful in extrapolating the coronal field upward from underlying vector boundary data. However, we can only measure the magnetic field vector routinely with high accuracy in the photosphere, and unfortunately these data do not fulfill the force-free condition. We must therefore apply some transformations to these data before nonlinear force-free extrapolation codes can be self-consistently applied. To this end, we have developed a minimization procedure that yields a more chromosphere-like field, using the measured photospheric field vectors as input. The procedure includes force-free consistency integrals, spatial smoothing, and -- newly included in the version presented here -- an improved match to the field direction as inferred from fibrils as can be observed in, e.g., chromospheric H$\alpha$ images. We test the procedure using a model active-region field that included buoyancy forces at the photospheric level. The proposed preprocessing method allows us to approximate the chromospheric vector field to within a few degrees and the free energy in the coronal field to within one percent.Comment: 22 pages, 6 Figur

Quasi-Static 3D-Magnetic Field Evolution in Solar Active Region NOAA 11166 Associated with X1.5 Flare

We study the quasi-static evolution of coronal magnetic fields constructed from the Non Linear Force Free Field (NLFFF) approximation aiming to understand the relation between the magnetic field topology and ribbon emission during an X1.5 flare in active region (AR) NOAA 11166. The flare with a quasi-elliptical, and two remote ribbons occurred on March 9, 2011 at 23:13UT over a positive flux region surrounded by negative flux at the center of the bipolar AR. Our analysis of the coronal magnetic structure with potential and NLFFF solutions unveiled the existence of a single magnetic null point associated with a fan-spine topology and is co-spatial with the hard X-ray source. The footpoints of the fan separatrix surface agree with the inner edge of the quasi-elliptical ribbon and the outer spine is linked to one of the remote ribbons. During the evolution, the slow footpoint motions stressed the fieldlines along the polarity inversion line and caused electric current layers in the corona around the fan separatrix surface. These current layers trigger magnetic reconnection as a consequence of dissipating currents, which are visible as cusped shape structures at lower heights. The reconnection process reorganised the magnetic field topology whose signatures are observed at the separatrices/QSL structure both in the photosphere and corona during the pre-to-post flare evolution. In agreement with previous numerical studies, our results suggest that the line-tied footpoint motions perturb the fan-spine system and cause null point reconnection, which eventually causes the flare emission at the footpoints of the fieldlines.Comment: Accepted, 2014, The Astrophysical Journa

Testing non-linear force-free coronal magnetic field extrapolations with the Titov-Demoulin equilibrium

CONTEXT: As the coronal magnetic field can usually not be measured directly, it has to be extrapolated from photospheric measurements into the corona. AIMS: We test the quality of a non-linear force-free coronal magnetic field extrapolation code with the help of a known analytical solution. METHODS: The non-linear force-free equations are numerically solved with the help of an optimization principle. The method minimizes an integral over the force-free and solenoidal condition. As boundary condition we use either the magnetic field components on all six sides of the computational box in Case I or only on the bottom boundary in Case II. We check the quality of the reconstruction by computing how well force-freeness and divergence-freeness are fulfilled and by comparing the numerical solution with the analytical solution. The comparison is done with magnetic field line plots and several quantitative measures, like the vector correlation, Cauchy Schwarz, normalized vector error, mean vector error and magnetic energy. RESULTS: For Case I the reconstructed magnetic field shows good agreement with the original magnetic field topology, whereas in Case II there are considerable deviations from the exact solution. This is corroborated by the quantitative measures, which are significantly better for Case I. CONCLUSIONS: Despite the strong nonlinearity of the considered force-free equilibrium, the optimization method of extrapolation is able to reconstruct it; however, the quality of reconstruction depends significantly on the consistency of the input data, which is given only if the known solution is provided also at the lateral and top boundaries, and on the presence or absence of flux concentrations near the boundaries of the magnetogram.Comment: 6 pages, 2 figures, Research Not

Nonlinear force-free modelling: influence of inaccuracies in the measured magnetic vector

Context: Solar magnetic fields are regularly extrapolated into the corona starting from photospheric magnetic measurements that can suffer from significant uncertainties. Aims: Here we study how inaccuracies introduced into the maps of the photospheric magnetic vector from the inversion of ideal and noisy Stokes parameters influence the extrapolation of nonlinear force-free magnetic fields. Methods: We compute nonlinear force-free magnetic fields based on simulated vector magnetograms, which have been produced by the inversion of Stokes profiles, computed froma 3-D radiation MHD simulation snapshot. These extrapolations are compared with extrapolations starting directly from the field in the MHD simulations, which is our reference. We investigate how line formation and instrumental effects such as noise, limited spatial resolution and the effect of employing a filter instrument influence the resulting magnetic field structure. The comparison is done qualitatively by visual inspection of the magnetic field distribution and quantitatively by different metrics. Results: The reconstructed field is most accurate if ideal Stokes data are inverted and becomes less accurate if instrumental effects and noise are included. The results demonstrate that the non-linear force-free field extrapolation method tested here is relatively insensitive to the effects of noise in measured polarization spectra at levels consistent with present-day instruments. Conclusions heading: Our results show that we can reconstruct the coronal magnetic field as a nonlinear force-free field from realistic photospheric measurements with an accuracy of a few percent, at least in the absence of sunspots.Comment: A&A, accepted, 9 Pages, 4 Figure

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, \textit{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 magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/Vector Spectromagnetograph(VSM) 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(ARs). We have studied the time evolution of 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.Comment: Solar Physic

Coronal hole boundaries evolution at small scales: I. EIT 195 A and TRACE 171 A view

We aim at studying the small-scale evolution at the boundaries of an equatorial coronal hole connected with a channel of open magnetic flux with the polar region and an `isolated' one in the extreme-ultraviolet spectral range. We intend to determine the spatial and temporal scale of these changes. Imager data from TRACE in the Fe IX/X 171 A passband and EIT on-board Solar and Heliospheric Observatory in the Fe XII 195 A passband were analysed. We found that small-scale loops known as bright points play an essential role in coronal holes boundaries evolution at small scales. Their emergence and disappearance continuously expand or contract coronal holes. The changes appear to be random on a time scale comparable with the lifetime of the loops seen at these temperatures. No signature was found for a major energy release during the evolution of the loops. Although coronal holes seem to maintain their general shape during a few solar rotations, a closer look at their day-by-day and even hour-by-hour evolution demonstrates a significant dynamics. The small-scale loops (10" - 40" and smaller) which are abundant along coronal hole boundaries have a contribution to the small-scale evolution of coronal holes. Continuous magnetic reconnection of the open magnetic field lines of the coronal hole and the closed field lines of the loops in the quiet Sun is more likely to take place.Comment: 7 pages, 6 figures. in press in A&