Magnetic helicity is conserved under ideal magnetohydrodynamics (MHD) and quasi-conserved even under a resistive process. The standard definition for magnetic helicity cannot be applied directly to an open magnetic field in a volume, because it is gauge-dependent. Instead, the relative magnetic helicity is widely used. We find that the energy of a potential magnetic field in a rectangular domain with periodic lateral boundary conditions is less than that of the field with a fixed normal component on all six boundaries. To make use of this lower energy potential field in the analysis of relative magnetic helicity, we introducing a new definition for magnetic helicity for the magnetic field, which involves the periodic potential field. We apply this definition to a sequence of analytic solutions and a numerical simulation. The results show that our new gauge-invariant helicity is very close to the current-carrying part of the relative magnetic helicity of the original magnetic field. We find also that the ratio between the current-carrying helicity and the relative magnetic helicity for the original and our defined relative helicity show different behavior. It seems that the new helicity is more sensitive to the component of the field due to the electric current in the volume, which is the source for instabilities and solar eruptive phenomena.

Eruptive flares are sudden releases of magnetic energy that involve many phenomena, several of which can be explained by the standard 2D flare model and its realizations in three-dimensions. We analyze a three-dimensional magnetohydrodynamics simulation in the framework of this model that naturally explains the contraction of coronal loops in the proximity of the flare sites, as well as the inflow towards the region above the cusp-shaped loops. We find that two vorticity arcs located along the flanks of the erupting magnetic flux rope are generated as soon as the eruption begins. The magnetic arcades above the flux-rope legs are then subjected to expansion, rotation or contraction depending on which part of the vortex-flow advects them. In addition to the vortices, an inward-directed magnetic pressure gradient exists in the current sheet below the magnetic flux rope. It results in the formation of a sink that is maintained by reconnection. We conclude that coronal loop apparent implosions observed during eruptive flares are the result of hydro-magnetic effects related to the generation of vortex- and sink-flows when a flux rope moves in a magnetized environment.

Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.

We studied the magnetic topology of active region 12158 on 2014 September 10 and compared it with the observations before and early in the flare which begins at 17:21 UT (SOL2014-09-10T17:45:00). Our results show that the sigmoidal structure and flare ribbons of this active region observed by SDO/AIA can be well reproduced from a Grad-Rubin non linear force free field extrapolation method. Various inverse-S and -J shaped magnetic field lines, that surround a coronal flux rope, coincide with the sigmoid as observed in different extreme ultraviolet wavelengths, including its multi-threaded curved ends. Also, the observed distribution of surface currents in the magnetic polarity where it was not prescribed is well reproduced. This validates our numerical implementation and set-up of the Grad-Rubin method. The modeled double inverse-J shaped Quasi-Separatrix Layer (QSL) footprints match the observed flare ribbons during the rising phase of the flare, including their hooked parts. The spiral-like shape of the latter may be related to a complex pre-eruptive flux rope with more than one turn of twist, as obtained in the model. These ribbon-associated flux-rope QSL-footprints are consistent with the new standard flare model in 3D, with the presence of a hyperbolic flux tube located below an inverse tear drop shaped coronal QSL. This is a new step forward forecasting the locations of reconnection and ribbons in solar flares, and the geometrical properties of eruptive flux ropes.

Magnetic flux ropes are topological structures consisting of twisted magnetic field lines that globally wrap around an axis. The torus instability model predicts that a magnetic flux rope of major radius $R$ undergoes an eruption when its axis reaches a location where the decay index $-d(\ln B_{ex})/d(\ln R)$ of the ambient magnetic field $B_{ex}$ is larger than a critical value. In the current-wire model, the critical value depends on the thickness and time-evolution of the current channel. We use magneto-hydrodynamic (MHD) simulations to investigate if the critical value of the decay index at the onset of the eruption is affected by the magnetic flux rope's internal current profile and/or by the particular pre-eruptive photospheric dynamics. The evolution of an asymmetric, bipolar active region is driven by applying different classes of photospheric motions. We find that the critical value of the decay index at the onset of the eruption is not significantly affected by either the pre-eruptive photospheric evolution of the active region or by the resulting different magnetic flux ropes. As in the case of the current-wire model, we find that there is a `critical range' $ [1.3-1.5]$, rather than a `critical value' for the onset of the torus instability. This range is in good agreement with the predictions of the current-wire model, despite the inclusion of line-tying effects and the occurrence of tether-cutting magnetic reconnection.

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

We present a code for solving the nonlinear force-free equations in spherical polar geometry, with the motivation of modeling the magnetic field in the corona. The code is an implementation of the Grad-Rubin method. Our method is applicable to a spherical domain of arbitrary angular size. The implementation is based on a global spectral representation for the magnetic field which makes no explicit assumptions about the form of the magnetic field at the transverse boundaries of the domain. We apply the code to a bipolar test case with analytic boundary conditions, and we demonstrate the convergence of the Grad-Rubin method, and the self-consistency of the resulting numerical solution.

The coronal magnetic field cannot be directly observed, but in principle it can be reconstructed from the comparatively well observed photospheric magnetic field. A popular approach uses a nonlinear force-free model. Non-magnetic forces at the photosphere are significant meaning the photospheric data are inconsistent with the force-free model, and this causes problems with the modeling (De Rosa et al., Astrophys. J. 696, 1780, 2009). In this paper we present a numerical implementation of the Grad-Rubin method for reconstructing the coronal magnetic field using a magnetostatic model. This model includes a pressure force and a non-zero magnetic Lorentz force. We demonstrate our implementation on a simple analytic test case and obtain the speed and numerical error scaling as a function of the grid size.

Presently, many models of the coronal magnetic field rely on photospheric vector magnetograms but these data have been shown to be problematic as the sole boundary information for nonlinear force-free field (NLFFF) extrapolations. Magnetic fields in the corona manifest themselves in high-energy images (X-rays and EUV) in the shapes of coronal loops, providing an additional constraint that at present is not used due to the mathematical complications of incorporating such input into numerical models. Projection effects and the limited number of usable loops further complicate the use of coronal information. We develop and test an algorithm to use images showing coronal loops in the modeling of the solar coronal magnetic field. We first fit projected field lines with field lines of constant-\als force-free fields to approximate the three-dimensional distribution of currents in the corona along a sparse set of trajectories. We then apply a Grad-Rubin-like iterative technique to obtain a volume-filling nonlinear force-free model of the magnetic field, modifying method presented in \citet{Wheatland2007}. We thoroughly test the technique on known analytical and solar-like model magnetic fields previously used for comparing different extrapolation techniques \citep{Schrijver2006, Schrijver2008} and compare the results with those obtained by presently available methods that rely only on the photospheric data. We conclude that we have developed a functioning method of modeling the coronal magnetic field by combining the line-of-sight component of photospheric magnetic field with information from coronal images. Vector magnetograms over the full or partial photospheric boundary of the numerical domain could optionally be used.

The NOAA active region AR 11029 was a small but highly active sunspot region which produced 73 GOES soft X-ray flares. The flares appear to show a departure from the well known power-law frequency-size distribution. Specifically, too few GOES C-class and no M-class flares were observed by comparison with a power-law distribution (Wheatland in Astrophys. J. 710, 1324, 2010). This was conjectured to be due to the region having insufficient magnetic energy to power large events. We construct nonlinear force-free extrapolations of the coronal magnetic field of active region AR 11029 using data taken on 24 October by the SOLIS Vector-SpectroMagnetograph (SOLIS/VSM), and data taken on 27 October by the Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). Force-free modeling with photospheric magnetogram data encounters problems because the magnetogram data are inconsistent with a force-free model, and we employ a recently developed `self-consistency' procedure which addresses this and accommodates uncertainties in the boundary data (Wheatland and Regnier in Astrophys. J. 700, L88, 2009). We calculate the total energy and free energy of the self-consistent solution and find that the free energy was 4x10^29 erg on 24 October, and 7x10^31 erg on 27 October. An order of magnitude scaling between RHESSI non-thermal energy and GOES peak X-ray flux is established from a sample of flares from the literature and is used to estimate flare energies from observed GOES peak X-ray flux. Based on the scaling, we conclude that the estimated free energy of AR 11029 on 27 October when the flaring rate peaked is sufficient to power M-class or X-class flares, and hence the modeling does not appear to support the hypothesis that the absence of large flares is due to the region having limited energy.