Starting from an exact, steady-state, force-free solution of the magnetohydrodynamic (MHD) equations, we investigate how resistive current layers are induced by perturbing line-tied three-dimensional magnetic equilibria. This is achieved by the superposition of a weak perturbation field in the domain, in contrast to studies where the boundary is driven by slow motions, like those present in photospheric active regions. Our aim is to quantify how the current structures are altered by the contribution of so called quasi-separatrix layers (QSLs) as the null point is shifted outside the computational domain. Previous studies based on magneto-frictional relaxation have indicated that, despite the severe field line gradients of the QSL, the presence of a null is vital in maintaining fast reconnection. Here, we explore this notion using highly resolved simulations of the full MHD evolution. We show that for the null-point configuration, the resistive scaling of the peak current density is close to $J\sim\eta^{-1}$, while the scaling is much weaker, i.e. $J\sim\eta^{-0.4}$, when only the QSL connectivity gradients provide a site for the current accumulation.

The open problem of how singular current structures form in line-tied, three-dimensional magnetic fields is addressed. A Lagrangian magneto-frictional relaxation method is employed to model the field evolution towards the final near-singular state. Our starting point is an exact force-free solution of the governing magnetohydrodynamic (MHD) equations which is sufficiently general to allow for topological features like magnetic nulls to be inside or outside the computational domain, depending on a simple set of parameters. Quasi-separatrix layers (QSLs) are present in these structures and together with the magnetic nulls, they significantly influence the accumulation of current. It is shown that perturbations affecting the lateral boundaries of the configuration lead not only to collapse around the magnetic null, but also to significant QSL currents. Our results show that once a magnetic null is present, the developing currents are always attracted to that specific location and show a much stronger scaling with resolution than the currents which form along the QSL. In particular, the null point scalings can be consistent with models of "fast" reconnection. The QSL currents also appear to be unbounded but give rise to weaker singularities, independent of the perturbation amplitude.

This paper considers the current distributions that derive from finite amplitude perturbations of line-tied magnetic fields comprising hyperbolic field structures. The initial equilibrium on which we principally focus is a planar magnetic X-point threaded by a uniform axial field. This field is line-tied on all surfaces but subject to three-dimensional (3D) disturbances that alter the initial topology. Results of ideal relaxation simulations are presented which illustrate how intense current structures form that can be related, through the influence of line-tying, to the quasi-separatrix layers (QSLs) of the initial configuration. It is demonstrated that the location within the QSL that attracts the current, and its scaling properties, are strongly dependent on the relative dimensions of the QSL with respect to the line-tied boundaries. These results are contrasted with the behavior of a line-tied 3D field containing an isolated null point. In this case, it is found that the dominant current always forms at the null, but that the collapse is inhibited when the null is closer to a line-tied boundary.

Force-free magnetic fields are important in many astrophysical settings. Determining the properties of such force-free fields -- especially smoothness and stability properties -- is crucial to understanding many key phenomena in astrophysical plasmas, for example energy release processes that heat the plasma and lead to dynamic or explosive events. Here we report on a serious limitation on the computation of force-free fields that has the potential to invalidate the results produced by numerical force-free field solvers even for cases in which they appear to converge (at fixed grid resolution) to an equilibrium magnetic field. In the present work we discuss this problem within the context of a Lagrangian relaxation scheme that conserves magnetic flux and div(B) identically. Error estimates are introduced to assess the quality of the calculated equilibrium. We go on to present an algorithm, based on re-writing the curl operation via Stokes' theorem, for calculating the current which holds great promise for improving dramatically the accuracy of the Lagrangian relaxation procedure.

We use a numerical nonlinear multigrid magnetic relaxation technique to investigate the generation of current sheets in three-dimensional magnetic flux braiding experiments. We are able to catalogue the relaxed nonlinear force-free equilibria resulting from the application of deformations to an initially undisturbed region of plasma containing a uniform, vertical magnetic field. The deformations are manifested by imposing motions on the bounding planes to which the magnetic field is anchored. Once imposed the new distribution of magnetic footpoints are then taken to be fixed, so that the rest of the plasma must then relax to a new equilibrium configuration. For the class of footpoint motions we have examined, we find that singular and nonsingular equilibria can be generated. By singular we mean that within the limits imposed by numerical resolution we find that there is no convergence to a well-defined equilibrium as the number of grid points in the numerical domain is increased. These singular equilibria contain current "sheets" of ever-increasing current intensity and decreasing width; they occur when the footpoint motions exceed a certain threshold, and must include both twist and shear to be effective. On the basis of these results we contend that flux braiding will indeed result in significant current generation. We discuss the implications of our results for coronal heating.