Three-dimensional solar active region magnetohydrostatic models and their stability using Euler potentials

Active regions (ARs) are magnetic structures typically found in the solar atmosphere. We calculated several magnetohydrostatic (MHS) equilibrium models that include the effect of a finite plasma-beta and gravity and that are representative of AR structures in three dimensions. The construction of the models is based on the use of two Euler potentials, alpha and beta, that represent the magnetic field as B = del alpha x del beta. The ideal MHS nonlinear partial differential equations were solved numerically using finite elements in a fixed 3D rectangular domain. The boundary conditions were initially chosen to correspond to a potential magnetic field (current-free) with known analytical expressions for the corresponding Euler potentials. The distinctive feature of our model is that we incorporated the effect of shear by progressively deforming the initial potential magnetic field. This procedure is quite generic and allowed us to generate a vast variety of MHS models. The thermal structure of the ARs was incorporated through the dependence of gas pressure and temperature on the Euler potentials. Using this method, we achieved the characteristic hot and overdense plasma found in ARs, but we demonstrate that the method can also be applied to study configurations with open magnetic field lines. Furthermore, we investigated basic topologies that include neutral lines. Our focus is on the force balance of the structures, and we do not consider the energy balance in the constructed models. In addition, we addressed the difficult question of the stability of the calculated 3D models. We find that if the plasma is convectively stable, then the system is not prone, in general, to develop magnetic Rayleigh-Taylor instabilities. However, when the plasma-beta is increased or the density at the core of the AR is high, then the magnetic configuration becomes unstable due to magnetic buoyancy.

Automated summary

In this study, we calculated three-dimensional equilibrium models of active region structures that take into account finite plasma-beta and gravity effects. By gradually deforming the initial potential magnetic field, we successfully incorporated the effect of shear and generated a variety of magnetohydrostatic models. The thermal structure of the active regions was achieved through the dependence of gas pressure and temperature on the Euler potentials. Additionally, we investigated the basic topologies that include neutral lines and analyzed the stability of the calculated 3D models.