How to Implement Boundary Conditions for Force-Free Magnetic Field Computations Using Vector Potentials

A force-free field (FFF) is determined solely by the normal components of magnetic field and current density on the entire boundary of the domain. Methods employing three components of magnetic field suffer from overspecification of boundary conditions and/or a nonzero divergence-B problem. A vector potential formulation eliminates the latter issue, but introduces difficulties in imposing the normal component of current density at the boundary. This paper proposes four different boundary treatment methods within the vector potential formulation. We conduct a comparative analysis of the vector potential FFF solvers that we have developed incorporating these methods against other FFF codes in different magnetic field representations. Although the vector potential solvers with the new boundary treatments do not outperform our poloidal-toroidal formulation code, they demonstrate comparable or superior performance compared to the optimization code in SolarSoftWare. The methods developed here are expected to be readily applied not only to force-free field computations but also to time-dependent data-driven simulations.