In this paper we seek to understand the timescale on which the photospheric motions on the Sun braid coronal magnetic field lines. This is a crucial ingredient for determining the viability of the braiding mechanism for explaining the high temperatures observed in the corona. We study the topological complexity induced in the coronal magnetic field, primarily using plasma motions extracted from magneto-convection simulations. This topological complexity is quantified using the field line winding, finite time topological entropy and passive scalar mixing. With these measures we contrast mixing efficiencies of the magneto-convection simulation, a benchmark flow known as a ``blinking vortex', and finally photospheric flows inferred from sequences of observed magnetograms using local correlation tracking. While the highly resolved magneto-convection simulations induce a strong degree of field line winding and finite time topological entropy, the values obtained from the observations from the plage region are around an order of magnitude smaller. This behavior is carried over to the finite time topological entropy. Nevertheless, the results suggest that the photospheric motions induce complex tangling of the coronal field on a timescale of hours.

We study the effect of photospheric footpoint motions on magnetic field structures containing magnetic nulls. The footpoint motions are prescribed on the photospheric boundary as a velocity field which entangles the magnetic field. We investigate the propagation of the injected energy, the conversion of energy, emergence of current layers and other consequences of the non-trivial magnetic field topology in this situation. These boundary motions lead initially to an increase in magnetic and kinetic energy. Following this, the energy input from the photosphere is partially dissipated and partially transported out of the domain through the Poynting flux. The presence of separatrix layers and magnetic null-points fundamentally alters the propagation behavior of disturbances from the photosphere into the corona. Depending on the field line topology close to the photosphere, the energy is either trapped or free to propagate into the corona.

We investigate the existence of magnetohydrostatic equilibria for topologically complex magnetic fields. The approach employed is to perform ideal numerical relaxation experiments. We use a newly-developed Lagrangian relaxation scheme that exactly preserves the magnetic field topology during the relaxation. Our configurations include both twisted and sheared fields, of which some fall into the category for which Parker (1972) predicted no force-free equilibrium. The first class of field considered contains no magnetic null points, and field lines connect between two perfectly conducting plates. In these cases we observe only resolved current layers of finite thickness. In further numerical experiments we confirm that magnetic null points are loci of singular currents.

By defining an appropriate field line helicity, we apply the powerful concept of magnetic helicity to the problem of global magnetic field evolution in the Sun's corona. As an ideal-magnetohydrodynamic invariant, the field line helicity is a meaningful measure of how magnetic helicity is distributed within the coronal volume. It may be interpreted, for each magnetic field line, as a magnetic flux linking with that field line. Using magneto-frictional simulations, we investigate how field line helicity evolves in the non-potential corona as a result of shearing by large-scale motions on the solar surface. On open magnetic field lines, the helicity injected by the Sun is largely output to the solar wind, provided that the coronal relaxation is sufficiently fast. But on closed magnetic field lines, helicity is able to build up. We find that the field line helicity is non-uniformly distributed, and is highly concentrated in twisted magnetic flux ropes. Eruption of these flux ropes is shown to lead to sudden bursts of helicity output, in contrast to the steady flux along the open magnetic field lines.

We examine the dynamics of magnetic flux tubes containing non-trivial field line braiding (or linkage), using mathematical and computational modelling, in the context of testable predictions for the laboratory and their significance for solar coronal heating. We investigate the existence of braided force-free equilibria, and demonstrate that for a field anchored at perfectly-conducting plates, these equilibria exist and contain current sheets whose thickness scales inversely with the braid complexity - as measured for example by the topological entropy. By contrast, for a periodic domain braided exact equilibria typically do not exist, while approximate equilibria contain thin current sheets. In the presence of resistivity, reconnection is triggered at the current sheets and a turbulent relaxation ensues. We finish by discussing the properties of the turbulent relaxation and the existence of constraints that may mean that the final state is not the linear force-free field predicted by Taylor's hypothesis.

The Parker hypothesis (Parker (1972)) assumes that heating of coronal loops occurs due to reconnection, induced when photospheric motions braid field lines to the point of current sheet formation. In this contribution we address the question of how the nature of photospheric motions affects heating of braided coronal loops. We design a series of boundary drivers and quantify their properties in terms of complexity and helicity injection. We examine a series of long-duration full resistive MHD simulations in which a simulated coronal loop, consisting of initially uniform field lines, is subject to these photospheric flows. Braiding of the loop is continually driven until differences in behaviour induced by the drivers can be characterised. It is shown that heating is crucially dependent on the nature of the photospheric driver - coherent motions typically lead to fewer large energy release events, while more complex motions result in more frequent but less energetic heating events.

Predicting the final state of turbulent plasma relaxation is an important challenge, both in astrophysical plasmas such as the Sun's corona and in controlled thermonuclear fusion. Recent numerical simulations of plasma relaxation with braided magnetic fields identified the possibility of a novel constraint, arising from the topological degree of the magnetic field-line mapping. This constraint implies that the final relaxed state is drastically different for an initial configuration with topological degree 1 (which allows a Taylor relaxation) and one with degree 2 (which does not reach a Taylor state). Here we test this transition in numerical resistive-magnetohydrodynamic simulations, by embedding a braided magnetic field in a linear force-free background. Varying the background force-free field parameter generates a sequence of initial conditions with a transition between topological degree 1 and 2. For degree 1, the relaxation produces a single twisted flux tube, while for degree 2 we obtain two flux tubes. For predicting the exact point of transition, it is not the topological degree of the whole domain that is relevant, but only that of the turbulent region.

One proposed resolution to the long-standing problem of solar coronal heating involves the buildup of magnetic energy in the corona due to turbulent motions at the photosphere that braid the coronal field, and the subsequent release of this energy via magnetic reconnection. In this paper the ideal relaxation of braided magnetic fields modelling solar coronal loops is followed. A sequence of loops with increasing braid complexity is considered, with the aim of understanding how this complexity influences the development of small scales in the magnetic field, and thus the energy available for heating. It is demonstrated that the ideally accessible force-free equilibrium for these braided fields contains current layers of finite thickness. It is further shown that for any such braided field, if a force-free equilibrium exists then it should contain current layers whose thickness is determined by length scales in the field line mapping. The thickness and intensity of the current layers follow scaling laws, and this allows us to extrapolate beyond the numerically accessible parameter regime and to place an upper bound on the braid complexity possible at coronal plasma parameters. At this threshold level the braided loop contains $10^{26}$--$10^{28}{\rm ergs}$ of available free magnetic energy, more than sufficient for a large nanoflare.

We investigate the evolution of field line helicity for magnetic fields that connect two boundaries without null points, with emphasis on localized finite-B magnetic reconnection. Total (relative) magnetic helicity is already recognized as an important topological constraint on magnetohydrodynamic processes. Field line helicity offers further advantages because it preserves all topological information and can distinguish between different magnetic fields with the same total helicity. Magnetic reconnection changes field connectivity and field line helicity reflects these changes; the goal of this paper is to characterize that evolution. We start by deriving the evolution equation for field line helicity and examining its terms, also obtaining a simplified form for cases where dynamics are localized within the domain. The main result, which we support using kinematic examples, is that during localized reconnection in a complex magnetic field, the evolution of field line helicity is dominated by a work-like term that is evaluated at the field line endpoints, namely the scalar product of the generalized field line velocity and the vector potential. Furthermore, the flux integral of this term over certain areas is very small compared to the integral of the unsigned quantity, which indicates that changes of field line helicity happen in a well-organized pairwise manner. It follows that reconnection is very efficient at redistributing helicity in complex magnetic fields despite having little effect on the total helicity.

We present a new code that performs a relaxation of a magnetic field towards a force-free state (Beltrami field) using a Lagrangian numerical scheme. Beltrami fields are of interest for the dynamics of many technical and astrophysical plasmas as they are the lowest energy states that the magnetic field can reach. The numerical method strictly preserves the magnetic flux and the topology of magnetic field lines. In contrast to other implementations we use mimetic operators for the spatial derivatives in order to improve accuracy for high distortions of the grid. Compared with schemes using direct derivatives we find that the final state of the simulation approximates a force-free state with a significantly higher accuracy. We implement the scheme in a code which runs on graphical processing units (GPU), which leads to an enhanced computing speed compared to previous relaxation codes.

A topological flux function is introduced to quantify the topology of magnetic braids: non-zero line-tied magnetic fields whose field lines all connect between two boundaries. This scalar function is an ideal invariant defined on a cross-section of the magnetic field, whose integral over the cross-section yields the relative magnetic helicity. Recognising that the topological flux function is an action in the Hamiltonian formulation of the field line equations, a simple formula for its differential is obtained. We use this to prove that the topological flux function uniquely characterises the field line mapping and hence the magnetic topology. A simple example is presented.

We introduce a topological flux function to quantify the topology of magnetic braids: non-zero, line-tied magnetic fields whose field lines all connect between two boundaries. This scalar function is an ideal invariant defined on a cross-section of the magnetic field, and measures the average poloidal magnetic flux around any given field line, or the average pairwise crossing number between a given field line and all others. Moreover, its integral over the cross-section yields the relative magnetic helicity. Using the fact that the flux function is also an action in the Hamiltonian formulation of the field line equations, we prove that it uniquely characterizes the field line mapping and hence the magnetic topology.

Magnetic flux tubes reaching from the solar convective zone into the chromosphere have to pass through the relatively cool, and therefore non-ideal (i.e. resistive) photospheric region enclosed between the highly ideal sub-photospheric and chromospheric plasma. It is shown that stationary MHD equilibria of magnetic flux tubes which pass through this region require an inflow of photospheric material into the flux tube and a deviation from iso-rotation along the tube axis. This means that there is a difference in angular velocity of the plasma flow inside the tube below and above the non-ideal region. Both effects increase with decreasing cross section of the tube. Although for characteristic parameters of thick flux tubes the effect is negligible, a scaling law indicates its importance for small-scale structures. The relevance of this "inflow effect" for the expansion of flux tubes above the photosphere is discussed.

Aims. We show how the build-up of magnetic gradients in the Sun's corona may be inferred directly from photospheric velocity data. This enables computation of magnetic connectivity measures such as the squashing factor without recourse to magnetic field extrapolation. Methods.Assuming an ideal evolution in the corona, and an initially uniform magnetic field, the subsequent field line mapping is computed by integrating trajectories of the (time-dependent) horizontal photospheric velocity field. The method is applied to a 12 hour high-resolution sequence of photospheric flows derived from Hinode/SOT magnetograms. Results. We find the generation of a network of quasi-separatrix layers in the magnetic field, which correspond to Lagrangian coherent structures in the photospheric velocity. The visual pattern of these structures arises primarily from the diverging part of the photospheric flow, hiding the effect of the rotational flow component: this is demonstrated by a simple analytical model of photospheric convection. We separate the diverging and rotational components from the observed flow and show qualitative agreement with purely diverging and rotational models respectively. Increasing the flow speeds in the model suggests that our observational results are likely to give a lower bound for the rate at which magnetic gradients are built up by real photospheric flows. Finally, we construct a hypothetical magnetic field with the inferred topology, that can be used for future investigations of reconnection and energy release.

We investigate the relaxation of braided magnetic loops in order to find out how the type of braiding via footpoint motions affects resultant heating of the loop. Two magnetic loops, braided in different ways, are used as initial conditions in resistive MHD simulations and their subsequent evolution is studied. The fields both undergo a resistive relaxation in which current sheets form and fragment and the system evolves towards a state of lower energy. In one case this relaxation is very efficient with current sheets filling the volume and homogeneous heating of the loop occurring. In the other case fewer current sheets develop, less magnetic energy is released in the process and a patchy heating of the loop results. The two cases, although very similar in their setup, can be distinguished by the mixing properties of the photospheric driver. The mixing can be measured by the topological entropy of the plasma flow, an observable quantity.

We present a model for separator reconnection due to an isolated reconnection process. Separator reconnection is a process which occurs in the neighbourhood of a distinguished field line (the separator) connecting two null points of a magnetic field. It is, for example, important for the dynamics of magnetic flux at the dayside magnetopause and in the solar corona. We find that, above a certain threshold, such a reconnection process generates new separators which leads to a complex system of magnetic flux tubes connecting regions of previously separated flux. Our findings are consistent with the findings of large numbers of separators in numerical simulations. We discuss how to measure and interpret the reconnection rate in a configuration with multiple separators.

The definition and measurement of magnetic reconnection in three-dimensional magnetic fields with multiple reconnection sites is a challenging problem, particularly in fields lacking null points. We propose a generalization of the familiar two-dimensional concept of a magnetic flux function to the case of a three-dimensional field connecting two planar boundaries. Using hyperbolic fixed points of the field line mapping, and their global stable and unstable manifolds, we define a unique flux partition of the magnetic field. This partition is more complicated than the corresponding (well-known) construction in a two-dimensional field, owing to the possibility of heteroclinic points and chaotic magnetic regions. Nevertheless, we show how the partition reconnection rate is readily measured with the generalized flux function. We relate our partition reconnection rate to the common definition of three-dimensional reconnection in terms of integrated parallel electric field. An analytical example demonstrates the theory, and shows how the flux partition responds to an isolated reconnection event.

A topological constraint on the dynamics of a magnetic field in a flux tube arises from the fixed point indices of its field line mapping. This can explain unexpected behaviour in recent resistive-magnetohydrodynamic simulations of magnetic relaxation. Here we present the theory for a general periodic flux tube, representing, for example, a toroidal confinement device or a solar coronal loop. We show how an ideal dynamics on the side boundary of the tube implies that the sum of indices over all interior fixed points is invariant. This constraint applies to any continuous evolution inside the tube, which may be turbulent and/or dissipative. We also consider the analogous invariants obtained from periodic points (fixed points of the iterated mapping). Although there is a countably infinite family of invariants, we show that they lead to at most two independent dynamical constraints. The second constraint applies only in certain magnetic configurations. Several examples illustrate the theory.

Aims: Our aim is to investigate the resistive relaxation of a magnetic loop that contains braided magnetic flux but no net current or helicity. The loop is subject to line-tied boundary conditions. We investigate the dynamical processes that occur during this relaxation, in particular the magnetic reconnection that occurs, and discuss the nature of the final equilibrium. Methods: The three-dimensional evolution of a braided magnetic field is followed in a series of resistive MHD simulations. Results: It is found that, following an instability within the loop, a myriad of thin current layers forms, via a cascade-like process. This cascade becomes more developed and continues for a longer period of time for higher magnetic Reynolds number. During the cascade, magnetic flux is reconnected multiple times, with the level of this `multiple reconnection' positively correlated with the magnetic Reynolds number. Eventually the system evolves into a state with no more small-scale current layers. This final state is found to approximate a non-linear force-free field consisting of two flux tubes of oppositely-signed twist embedded in a uniform background field.

The final state of turbulent magnetic relaxation in a reversed field pinch is well explained by Taylor's hypothesis. However, recent resistive-magnetohydrodynamic simulations of the relaxation of braided solar coronal loops have led to relaxed fields far from the Taylor state, despite the conservation of helicity. We point out the existence of an additional topological invariant in any flux tube with non-zero field: the topological degree of the field line mapping. We conjecture that this constrains the relaxation, explaining why only one of three example simulations reaches the Taylor state.

The response of the solar coronal magnetic field to small-scale photospheric boundary motions including the possible formation of current sheets via the Parker scenario is one of open questions of solar physics. Here we address the problem via a numerical simulation. The three-dimensional evolution of a braided magnetic field which is initially close to a force-free state is followed using a resistive MHD code. A long-wavelength instability takes place and leads to the formation of two thin current layers. Magnetic reconnection occurs across the current sheets with three-dimensional features shown, including an elliptic magnetic field structure about the reconnection site, and results in an untwisting of the global field structure.

The squashing factor Q (Titov et al. 2002), a property of the magnetic field line mapping, has been suggested as an indicator for the formation of current sheets, and subsequently magnetic reconnection, in astrophysical plasmas. Here we test this hypothesis for a particular class of braided magnetic fields which serve as a model for solar coronal loops. We explore the relationship between Quasi-Separatrix Layers (QSLs), that is, layer-like structures with high Q value, electric currents and integrated parallel currents, the latter being a quantity closely related to the reconnection rate. It is found that as the degree of braiding of the magnetic field is increased the maximum values of Q increase exponentially. At the same time, the distribution of Q becomes increasingly filamentary, with the width of the high-Q layers exponentially decreasing. This is accompanied by an increase in the number of layers so that as the field is increasingly braided the volume becomes occupied by a myriad of thin QSLs. QSLs are not found to be good predictors of current features in this class of braided fields. Indeed, despite the presence of multiple QSLs, the current associated with the field remains smooth and large-scale under ideal relaxation; the field dynamically adjusts to a smooth equilibrium. Regions of high Q are found to be better related to regions of high integrated parallel current than to actual current sheets.

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.

In this series of papers we examine magnetic reconnection in a domain where the magnetic field does not vanish and the non-ideal region is localised in space. In a previous paper we presented a technique for obtaining analytical solutions to the stationary resistive MHD equations in such a situation and examined specific examples of non-ideal reconnective solutions. Here we further develop the model, noting that certain ideal solutions may be superimposed onto the fundamental non-ideal solutions and examining the effect of imposing various such flows. Significant implications are found for the evolution of magnetic flux in the reconnection process. It is shown that, in contrast to the two-dimensional case, in three-dimensions there is a very wide variety of physically different steady reconnection solutions.

The braiding of the solar coronal magnetic field via photospheric motions - with subsequent relaxation and magnetic reconnection -- is one of the most widely debated ideas of solar physics. We readdress the theory in the light of developments in three-dimensional magnetic reconnection theory. It is known that the integrated parallel electric field along field lines is the key quantity determining the rate of reconnection, in contrast with the two-dimensional case where the electric field itself is the important quantity. We demonstrate that this difference becomes crucial for sufficiently complex magnetic field structures. A numerical method is used to relax a braided magnetic field to an ideal force-free equilibrium; that equilibrium is found to be smooth, with only large- scale current structures. However, the equilibrium is shown to have a highly filamentary integrated parallel current structure with extremely short length- scales. An analytical model is developed to show that, in a coronal situation, the length scales associated with the integrated parallel current structures will rapidly decrease with increasing complexity, or degree of braiding, of the magnetic field. Analysis shows the decrease in these length scales will, for any finite resistivity, eventually become inconsistent with the stability of a force- free field. Thus the inevitable consequence of the magnetic braiding process is shown to be a loss of equilibrium of the coronal field, probably via magnetic reconnection events.

A magnetic helicity integral is proposed which can be applied to domains which are not magnetically closed, i.e. have a non-vanishing normal component of the magnetic field on the boundary. In contrast to the relative helicity integral, which was previously suggested for magnetically open domains, it does not rely on a reference field and thus avoids all problems related to the choice of a particular reference field. Instead it uses a gauge condition on the vector potential, which corresponds to a particular topologically unique closure of the magnetic field in the external space. The integral has additional elegant properties and is easy to compute numerically in practice. For magnetically closed domains it reduces to the classical helicity integral.

Newly emerging magnetic flux can show a complicated linked or interwoven topology of the magnetic field. The complexity of this linkage or knottedness of magnetic flux is related to the free energy stored in the magnetic field. Magnetic reconnection provides a process to release this energy on the time scale of the dynamics. At the same time it approximately conserves the total magnetic helicity. Therefore the conservation of total magnetic helicity is a crucial constraint for the relaxation of complex magnetic fields. However, the total magnetic helicity is only the first, most elementary, quantity of an infinite series of topological invariants of the magnetic field. All these invariants are strictly conserved in ideal magnetohydrodynamics. As an example a preliminary set of these invariants is derived. The relevance of these higher order invariants for the final state of relaxation under magnetic reconnection and their implications for the release of magnetic energy are discussed.

Although the analysis of observational data indicates that quasi-separatrix layers (QSLs) have to play an important role in magnetic configurations of solar flares, the corresponding theory is only at an initial stage so far. In particular, there is still a need in a proper definition of QSL. This problem is analyzed here on the basis of geometrical properties of the mapping produced by the field lines which connect photospheric areas of positive and negative magnetic polarities of active regions. In general, one can find on the photosphere a unique pair of locally perpendicular line elements at one footpoint of a given field line, so that this pair is mapped into a similar pair of elements at the other footpoint. Along the directions of these elements the field line mapping only stretches and compresses the corresponding displacements of the neighboring field-line footpoints. This fact enables us to get a remarkably clear and concise definition of QSL as a volume filled by coronal field lines for which the ratio of the corresponding stretching and compressing coefficients is anomalously large. The new definition is also compared with the ones previously introduced by other authors. The theory is applied to flare events of the so-called sigmoid-type, using an analytical model of a twisted force-free configuration. It is shown that such a configuration may contain a QSL even if genuine separatrix surfaces are absent. Identifying QSLs is an essential prerequisite for understanding the mechanism of magnetic energy release in this type of flares. It is also demonstrated that the magnetic field under study has a Hamiltonian structure, which makes it possible to reveal a geometrical reason for the appearance of QSL in this case.