Deriving the strength and direction of the three-dimensional (3D) magnetic field in the solar atmosphere is fundamental for understanding its dynamics. Volume information on the magnetic field mostly relies on coupling 3D reconstruction methods with photospheric and/or chromospheric surface vector magnetic fields. Infrared coronal polarimetry could provide additional information to better constrain magnetic field reconstructions. However, combining such data with reconstruction methods is challenging, e.g., because of the optical-thinness of the solar corona and the lack and limitations of stereoscopic polarimetry. To address these issues, we introduce the Data-Optimized Coronal Field Model (DOCFM) framework, a model-data fitting approach that combines a parametrized 3D generative model, e.g., a magnetic field extrapolation or a magnetohydrodynamic model, with forward modeling of coronal data. We test it with a parametrized flux rope insertion method and infrared coronal polarimetry where synthetic observations are created from a known "ground truth" physical state. We show that this framework allows us to accurately retrieve the ground truth 3D magnetic field of a set of force-free field solutions from the flux rope insertion method. In observational studies, the DOCFM will provide a means to force the solutions derived with different reconstruction methods to satisfy additional, common, coronal constraints. The DOCFM framework therefore opens new perspectives for the exploitation of coronal polarimetry in magnetic field reconstructions and for developing new techniques to more reliably infer the 3D magnetic fields that trigger solar flares and coronal mass ejections.

The Coronal Multichannel Polarimeter (CoMP) routinely performs coronal polarimetric measurements using the Fe XIII 10747 $\AA$ and 10798 $\AA$ lines, which are sensitive to the coronal magnetic field. However, inverting such polarimetric measurements into magnetic field data is a difficult task because the corona is optically thin at these wavelengths and the observed signal is therefore the integrated emission of all the plasma along the line of sight. To overcome this difficulty, we take on a new approach that combines a parameterized 3D magnetic field model with forward modeling of the polarization signal. For that purpose, we develop a new, fast and efficient, optimization method for model-data fitting: the Radial-basis-functions Optimization Approximation Method (ROAM). Model-data fitting is achieved by optimizing a user-specified log-likelihood function that quantifies the differences between the observed polarization signal and its synthetic/predicted analogue. Speed and efficiency are obtained by combining sparse evaluation of the magnetic model with radial-basis-function (RBF) decomposition of the log-likelihood function. The RBF decomposition provides an analytical expression for the log-likelihood function that is used to inexpensively estimate the set of parameter values optimizing it. We test and validate ROAM on a synthetic test bed of a coronal magnetic flux rope and show that it performs well with a significantly sparse sample of the parameter space. We conclude that our optimization method is well-suited for fast and efficient model-data fitting and can be exploited for converting coronal polarimetric measurements, such as the ones provided by CoMP, into coronal magnetic field data.

In this paper we continue our theoretical studies on addressing what are the possible consequences of magnetic helicity accumulation in the solar corona. Our previous studies suggest that coronal mass ejections (CMEs) are natural products of coronal evolution as a consequence of magnetic helicity accumulation and the triggering of CMEs by surface processes such as flux emergence also have their origin in magnetic helicity accumulation. Here we use the same mathematical approach to study the magnetic helicity of axisymmetric power-law force-free fields, but focus on a family whose surface flux distributions are defined by self-similar force-free fields. The semi-analytical solutions of the axisymmetric self-similar force-free fields enable us to discuss the properties of force-free fields possessing a huge amount of accumulated magnetic helicity. Our study suggests that there may be an absolute upper bound on the total magnetic helicity of all bipolar axisymmetric force-free fields. And with the increase of accumulated magnetic helicity, the force-free field approaches being fully opened up, with Parker-spiral-like structures present around a current-sheet layer as evidence of magnetic helicity in the interplanetary space. It is also found that among the axisymmetric force-free fields having the same boundary flux distribution, the one that is self-similar is the one possessing the maximum amount of total magnetic helicity. This possibly gives a physical reason why self-similar fields are often found in astrophysical bodies, where magnetic helicity accumulation is presumably also taking place.

This paper follows up on a previous study showing that in an open atmosphere such as the solar corona the total magnetic helicity of a force-free field must be bounded and the accumulation of magnetic helicity in excess of its upper bound would initiate a non-equilibrium situation resulting in an expulsion such as a coronal mass ejection (CME). In the current paper, we investigate the dependence of the helicity bound on the boundary condition for several families of nonlinear force-free fields. Our calculation shows that the magnitude of the helicity upper bound of force-free fields is non-trivially dependent on the boundary condition. Fields with a multipolar boundary condition can have a helicity upper bound ten times smaller than those with a dipolar boundary condition when helicity values are normalized by the square of their respective surface poloidal fluxes. This suggests that a coronal magnetic field may erupt into a CME when the applicable helicity bound falls below the already accumulated helicity as the result of a slowly changing boundary condition. Our calculation also shows that a monotonic accumulation of magnetic helicity can lead to the formation of a magnetic flux rope applicable to kink instability. This suggests that CME initiations by exceeding helicity bound and by kink instability can both be the consequences of helicity accumulation in the corona. Our study gives insights into the observed associations of CMEs with the magnetic features at their solar surface origins.

A loss of magnetic field confinement is believed to be the cause of coronal mass ejections (CMEs), a major form of solar activity in the corona. The mechanisms for magnetic energy storage are crucial in understanding how a field may possess enough free energy to overcome the Aly limit and open up. Previously, we have pointed out that the accumulation of magnetic helicity in the corona plays a significant role in storing magnetic energy. In this paper, we investigate another hydromagnetic consequence of magnetic-helicity accumulation. We propose a conjecture that there is an upper bound on the total magnetic helicity that a force-free field can contain. This is directly related to the hydromagnetic property that force-free fields in unbounded space have to be self-confining. Although a mathematical proof of this conjecture for any field configuration is formidable, its plausibility can be demonstrated with the properties of several families of power-law, axisymmetric force-free fields. We put forth mathematical evidence, as well as numerical, indicating that an upper bound on the magnetic helicity may exist for such fields. Thus, the accumulation of magnetic helicity in excess of this upper bound would initiate a non-equilibrium situation, resulting in a CME expulsion as a natural product of coronal evolution.