(Abridged for ArXiv) An injection of energy towards a magnetic null point can drive reversals of current sheet polarity leading to time-dependent Oscillatory Reconnection, which is a possible explanation of how periodic phenomena can be generated when reconnection occurs in the solar atmosphere. However, the details of what controls the period of these oscillations is poorly understood, despite being crucial in assessing whether OR can account for observed periodic behaviour. This paper aims to highlight that different types of reconnection reversal are supported about null points, and that these are distinct from the oscillation on the closed-boundary, linear systems considered in the 1990s. In particular, we explore the features of a nonlinear oscillation local to the null point, and examine the effect of resistivity and perturbation energy on the period, contrasting it to the linear case. It is found that in the linear systems, the inverse Lundquist number dictates the period, provided the perturbation energy is small relative to the inverse Lundquist number defined on the boundary, regardless of the broadband structure of the initial perturbation. However, when the perturbation energy exceeds the threshold required for 'nonlinear' null collapse to occur, a complex oscillation of the magnetic field is produced which is, at best, only weakly-dependent on the resistivity. The resultant periodicity is strongly influenced by the amount of free energy, with more energetic perturbations producing higher-frequency oscillations. Crucially, with regards to typical solar-based and astrophysical-based input energies, we demonstrate that the majority far exceed the threshold for nonlinearity to develop. This substantially alters the properties and periodicity of both null collapse and subsequent OR. Therefore, nonlinear regimes of OR should be considered in solar and astrophysical contexts.

Implosive formation of current sheets is a fundamental plasma process. Previous studies focused on the early time evolution, while here our primary aim is to explore the longer-term evolution, which may be critical for determining the efficiency of energy release. To address this problem we investigate two closely-related problems, namely: (i) 1D, pinched anti-parallel magnetic fields and (ii) 2D, null point containing fields which are locally imbalanced ('null-collapse' or 'X-point collapse'). Within the framework of resistive MHD, we simulate the full nonlinear evolution through three distinct phases: the initial implosion, its eventual halting mechanism, and subsequent evolution post-halting. In a parameter study, we find the scaling with resistivity of current sheet properties at the halting time is in good agreement - in both geometries - with that inferred from a known 1D similarity solution. We find that the halting of the implosions occurs rapidly after reaching the diffusion scale by sudden Ohmic heating of the dense plasma within the current sheet, which provides a pressure gradient sufficient to oppose further collapse and decelerate the converging flow. This back-pressure grows to exceed that required for force balance and so the post-implosion evolution is characterised by the consequences of the current sheet `bouncing' outwards. These are: (i) the launching of propagating fast MHD waves (shocks) outwards and (ii) the width-wise expansion of the current sheet itself. The expansion is only observed to stall in the 2D case, where the pressurisation is relieved by outflow in the reconnection jets. In the 2D case, we quantify the maximum amount of current sheet expansion as it scales with resistivity, and analyse the structure of the reconnection region which forms post-expansion, replete with Petschek-type slow shocks and fast termination shocks.

The topology of magnetic fields near the open-closed flux boundary in the Sun's corona is an important influencing factor in the process of interchange reconnection, whereby plasma is exchanged between open and closed flux domains. Maps of the magnetic squashing factor at the radial outer boundary in coronal field models reveal the presence of the so-called `S-web', and suggest that interchange reconnection could potentially deposit closed coronal material into high-latitude regions far from the heliospheric current sheet. Here we demonstrate that certain features of the S-web reveal the underlying topological structure of the magnetic field. Specifically, in order for the arcing bands of highly squashed magnetic flux of the S-web to terminate or intersect away from the helmet streamer apex, there must be a null spine line that maps a finite segment of the photospheric open-closed boundary up to a singular point in the open flux domain. We propose that this association between null spine lines and arc termination points may be used to identify locations in the heliosphere that are preferential for the appearance of solar energetic particles or slow solar wind plasma with certain characteristics.

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.

Null collapse is an implosive process whereby MHD waves focus their energy in the vicinity of a null point, forming a current sheet and initiating magnetic reconnection. We consider, for the first time, the case of collapsing 3D magnetic null points in nonlinear, resistive MHD using numerical simulation, exploring key physical aspects of the system as well as performing a detailed parameter study. We find that within a particular plane containing the 3D null, the plasma and current density enhancements resulting from the collapse are quantitatively and qualitatively as per the 2D case in both the linear and nonlinear collapse regimes. However, the scaling with resistivity of the 3D reconnection rate - which is a global quantity - is found to be less favourable when the magnetic null point is more rotationally symmetric, due to the action of increased magnetic back-pressure. Furthermore, we find that with increasing ambient plasma pressure the collapse can be throttled, as is the case for 2D nulls. We discuss this pressure-limiting in the context of fast reconnection in the solar atmosphere and suggest mechanisms by which it may be overcome. We also discuss the implications of the results in the context of null collapse as a trigger mechanism of Oscillatory Reconnection, a time-dependent reconnection mechanism, and also within the wider subject of wave-null point interactions. We conclude that, in general, increasingly rotationally-asymmetric nulls will be more favourable in terms of magnetic energy release via null collapse than their more symmetric counterparts.

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.

Here we detail the dynamic evolution of localised reconnection regions about three-dimensional (3D) magnetic null points by using numerical simulation. We demonstrate for the first time that reconnection triggered by the localised collapse of a 3D null point due to an external MHD wave involves a self-generated oscillation, whereby the current sheet and outflow jets undergo a reconnection reversal process during which back-pressure formation at the jet heads acts to prise open the collapsed field before overshooting the equilibrium into an opposite-polarity configuration. The discovery that reconnection at fully 3D nulls can proceed naturally in a time-dependent and periodic fashion is suggestive that oscillatory reconnection mechanisms may play a role in explaining periodicity in astrophysical phenomena associated with magnetic reconnection, such as the observed quasi-periodicity of solar and stellar flare emission. Furthermore, we find a consequence of oscillatory reconnection is the generation of a plethora of freely-propagating MHD waves which escape the vicinity of the reconnection region

Coronal magnetic null points exist in abundance as demonstrated by extrapolations of the coronal field, and have been inferred to be important for a broad range of energetic events. These null points and their associated separatrix and spine field lines represent discontinuities of the field line mapping, making them preferential locations for reconnection. This field line mapping also exhibits strong gradients adjacent to the separatrix (fan) and spine field lines, that can be analysed using the `squashing factor', $Q$. In this paper we make a detailed analysis of the distribution of $Q$ in the presence of magnetic nulls. While $Q$ is formally infinite on both the spine and fan of the null, the decay of $Q$ away from these structures is shown in general to depend strongly on the null-point structure. For the generic case of a non-radially-symmetric null, $Q$ decays most slowly away from the spine/fan in the direction in which $|{\bf B}|$ increases most slowly. In particular, this demonstrates that the extended, elliptical high-$Q$ halo around the spine footpoints observed by Masson et al. (Astrophys. J., 700, 559, 2009) is a generic feature. This extension of the $Q$ halos around the spine/fan footpoints is important for diagnosing the regions of the photosphere that are magnetically connected to any current layer that forms at the null. In light of this, we discuss how our results can be used to interpret the geometry of observed flare ribbons in `circular ribbon flares', in which typically a coronal null is implicated. We conclude that both the physics in the vicinity of the null and how this is related to the extension of $Q$ away from the spine/fan can be used in tandem to understand observational signatures of reconnection at coronal null points.

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.

Global magnetic field extrapolations are now revealing the huge complexity of the Sun's corona, and in particular the structure of the boundary between open and closed magnetic flux. Moreover, recent developments indicate that magnetic reconnection in the corona likely occurs in highly fragmented current layers, and that this typically leads to a dramatic increase in the topological complexity beyond that of the equilibrium field. In this paper we investigate the consequences of reconnection at the open-closed flux boundary ("interchange reconnection") in a fragmented current layer. We demonstrate that it leads to a situation in which magnetic flux (and therefore plasma) from open and closed field regions is efficiently mixed together. This corresponds to an increase in the length and complexity of the open-closed boundary. Thus, whenever reconnection occurs at a null point or separator of the open-closed boundary, the associated separatrix arc of the so-called "S-web" in the high corona becomes not a single line but a band of finite thickness within which the open-closed flux boundary is highly structured. This has significant implications for the acceleration of the slow solar wind, for which the interaction of open and closed field is thought to be important, and may also explain the coronal origins of certain solar energetic particles. The topological structures examined contain magnetic null points, separatrices and separators, and include a model for a pseudo-streamer. The potential for understanding both the large scale morphology and fine structure observed in flare ribbons associated with coronal nulls is also discussed.

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 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.

In this work the dynamic magnetic field within a tearing-unstable three-dimensional (3D) current sheet about a magnetic null point is described in detail. We focus on the evolution of the magnetic null points and flux ropes that are formed during the tearing process. Generally, we find that both magnetic structures are created prolifically within the layer and are non-trivially related. We examine how nulls are created and annihilated during bifurcation processes, and describe how they evolve within the current layer. The type of null bifurcation first observed is associated with the formation of pairs of flux ropes within the current layer. We also find that new nulls form within these flux ropes, both following internal reconnection and as adjacent flux ropes interact. The flux ropes exhibit a complex evolution, driven by a combination of ideal kinking and their interaction with the outflow jets from the main layer. The finite size of the unstable layer also allows us to consider the wider effects of flux rope generation. We find that the unstable current layer acts as a source of torsional MHD waves and dynamic braiding of magnetic fields. The implications of these results to several areas of heliophysics are discussed.

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.

Three-dimensional magnetic null points are ubiquitous in the solar corona, and in any generic mixed-polarity magnetic field. We consider magnetic reconnection at an isolated coronal null point, whose fan field lines form a dome structure. We demonstrate using analytical and computational models several features of spine-fan reconnection at such a null, including the fact that substantial magnetic flux transfer from one region of field line connectivity to another can occur. The flux transfer occurs across the current sheet that forms around the null point during spine-fan reconnection, and there is no separator present. Also, flipping of magnetic field lines takes place in a manner similar to that observed in quasi-separatrix layer or slip-running reconnection.

We report here, for the first time, an observed instability of a Kelvin-Helmholtz (KH) nature occurring in a fully three-dimensional (3D) current-vortex sheet at the fan plane of a 3D magnetic null point. The current-vortex layer forms self-consistently in response to foot point driving around the spine lines of the null. The layer first becomes unstable at an intermediate distance from the null point, with the instability being characterized by a rippling of the fan surface and a filamentation of the current density and vorticity in the shear layer. Owing to the 3D geometry of the shear layer, a branching of the current filaments and vortices is observed. The instability results in a mixing of plasma between the two topologically distinct regions of magnetic flux on either side of the fan separatrix surface, as flux is reconnected across this surface. We make a preliminary investigation of the scaling of the system with the dissipation parameters. Our results indicate that the fan plane separatrix surface is an ideal candidate for the formation of current-vortex sheets in complex magnetic fields and, therefore, the enhanced heating and connectivity change associated with the instabilities of such layers.

The compression or expansion of a magnetic field that is initially potential is considered. It was recently suggested by Janse & Low [2009, ApJ, 690, 1089] that, following the volumetric deformation, the relevant lowest energy state for the magnetic field is another potential magnetic field that in general contains tangential discontinuities (current sheets). Here we examine this scenario directly using a numerical relaxation method that exactly preserves the topology of the magnetic field. It is found that of the magnetic fields discussed by Janse & Low, only those containing magnetic null points develop current singularities during an ideal relaxation, while the magnetic fields without null points relax toward smooth force-free equilibria with finite non-zero current.

Magnetic reconnection is a fundamental process in a plasma that facilitates the release of energy stored in the magnetic field by permitting a change in the magnetic topology. In this article we present a review of the current state of understanding of magnetic reconnection. We discuss theoretical results regarding the formation of current sheets in complex 3D magnetic fields, and describe the fundamental differences between reconnection in two and three dimensions. We go on to outline recent developments in modelling of reconnection with kinetic theory, as well as in the MHD framework where a number of new 3D reconnection regimes have been identified. We discuss evidence from observations and simulations of solar system plasmas that support this theory, and summarise some prominent locations in which this new reconnection theory is relevant in astrophysical plasmas.

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.

Three-dimensional null points are present in abundance in the solar corona, and the same is likely to be true in other astrophysical environments. Recent studies suggest that reconnection at such 3D nulls may play an important role in the coronal dynamics. In this paper the properties of the torsional spine and torsional fan modes of magnetic reconnection at 3D nulls are investigated. New analytical models are developed, which for the first time include a current layer that is spatially localised around the null, extending along either the spine or the fan of the null. These are complemented with numerical simulations. The principal aim is to investigate the effect of varying the degree of asymmetry of the null point magnetic field on the resulting reconnection process - where previous studies always considered a non-generic radially symmetric null. The geometry of the current layers within which torsional spine and torsional fan reconnection occur is found to be strongly dependent on the symmetry of the magnetic field. Torsional spine reconnection still occurs in a narrow tube around the spine, but with elliptical cross-section when the fan eigenvalues are different, and with the short axis of the ellipse being along the strong field direction. The spatiotemporal peak current, and the peak reconnection rate attained, are found not to depend strongly on the degree of asymmetry. For torsional fan reconnection, the reconnection occurs in a planar disk in the fan surface, which is again elliptical when the symmetry of the magnetic field is broken. The short axis of the ellipse is along the weak field direction, with the current being peaked in these weak field regions. The peak current and peak reconnection rate in this case are clearly dependent on the asymmetry, with the peak current increasing but the reconnection rate decreasing as the degree of asymmetry is increased.

Context. We investigate the dynamical evolution of the reconnection process at an initially linear 3D null point that is stressed by a localised shear motion across the spine axis. The fan plane is not rotationally symmetric and this allows for different behaviours depending on the alignment of the fan plane relative to the imposed driver direction. Aims. The aim is to show how the current accumulation and the associated reconnection process depends on the relative orientation between the driver imposed stress across the spine axis of the null and the main eigenvector direction in the fan plane. Methods. The time evolution of the 3D null point is investigated solving the 3D non-ideal MHD equations numerically in a Cartesian box. The magnetic field is frozen to the boundaries and the boundary velocity is only non-zero where the imposed driving is for stressing the system is applied. Results. The current accumulation is found to be along the direction of the fan eigenvector associated with the smallest eigenvalue until the direction of the driver is almost parallel to this eigenvector. When the driving velocity is parallel to the weak eigenvector and has an impulsive temporal profile the null only forms a weak current layer. However, when the null point is stressed continuously boundary effects dominates the current accumulation. Conclusions. There is a clear relation between the orientation of the current concentration and the direction of the fan eigenvector corresponding to the small eigenvalue. This shows that the structure of the magnetic field is the most important in determining where current is going to accumulate when a single 3D null point is perturbed by a simple shear motion across the spine axis. As the angle between the driving direction and the strong eigenvector direction increases, the current that accumulates at the null becomes progressively weaker.

To systematically stress a rotationally symmetric 3D magnetic null point by advecting the opposite footpoints of the spine axis in opposite directions. This stress eventually concentrates in the vicinity of the null point forming a local current sheet through which magnetic reconnection takes place. The aim is to look for a steady state evolution of the current sheet dynamics which may provide scaling relations for various characteristic parameters of the system. The evolution is followed by solving numerically the non-ideal MHD equations in a Cartesian domain. The null point is embedded in an initially constant density and temperature plasma. It is shown that a quasi-steady reconnection process can be set up at a 3D null by continuous shear driving. It appears that a true steady state in unlikely to be realised as the current layer tends to grow until restricted by the geometry of the computational domain and imposed driving profile. However, ratios between characteristic quantities clearly settle after some time to stable values -- so that the evolution is quasi-steady. The experiments show a number of scaling relations, but they do not provide a clear consensus for extending to lower magnetic resistivity or faster driving velocities. More investigations are needed to fully clarify the properties of current sheets at magnetic null points.

The magnetic field in many astrophysical plasmas -- such as the Solar corona and Earth's magnetosphere -- has been shown to have a highly complex, three-dimensional structure. Recent advances in theory and computational simulations have shown that reconnection in these fields also has a three-dimensional nature, in contrast to the widely used two-dimensional (or 2.5-dimensional) models. Here we discuss the underlying theory of three-dimensional magnetic reconnection. We also review a selection of new models that illustrate the current state of the art, as well as highlighting the complexity of energy release processes mediated by reconnection in complicated three-dimensional magnetic fields.

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 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 aim of this paper is to investigate the properties of magnetic reconnection at a 3D null point, with respect to their dependence on the symmetry of the magnetic field around the null. In particular we examine the rate of flux transport across the null point with symmetric/asymmetric diffusion regions, as well as how the current sheet forms in time, and its properties. Mathematical modelling and finite difference resistive MHD simulations are used. It is found that the basic structure of the mode of magnetic reconnection considered is unaffected by varying the magnetic field symmetry, that is, the plasma flow is found cross both the spine and fan of the null. However, the peak intensity and dimensions of the current sheet are dependent on the symmetry/ asymmetry of the field lines. As a result, the reconnection rate is also found to be strongly dependent on the field asymmetry.

Recent advances in theory and computational experiments have shown the need to refine the previous categorisation of magnetic reconnection at three-dimensional null points -- points at which the magnetic field vanishes. We propose here a division into three different types, depending on the nature of the flow near the spine and fan of the null. The spine is an isolated field line which approaches the null (or recedes from it), while the fan is a surface of field lines which recede from it (or approach it). So-called "torsional spine reconnection" occurs when field lines in the vicinity of the fan rotate, with current becoming concentrated along the spine, so that nearby field lines undergo rotational slippage. In "torsional fan reconnection" field lines near the spine rotate and create a current that is concentrated in the fan with a rotational flux mismatch and rotational slippage. In both of these regimes, the spine and fan are perpendicular and there is no flux transfer across spine or fan. The third regime, called "spine-fan reconnection", is the most common in practice and combines elements of the previous spine and fan models. In this case, in response to a generic shearing motion, the null point collapses to form a current sheet that is focused at the null itself, in a sheet that locally spans both the spine and fan. In this regime the spine and fan are no longer perpendicular and there is flux transfer across both of them.

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.

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.

The nature of current sheet formation in the vicinity of three-dimensional (3D) magnetic null points is investigated. The particular focus is upon the effect of the compressibility of the plasma on the qualitative and quantitative properties of the current sheet. An initially potential 3D null is subjected to shearing perturbations, as in a previous paper [Pontin et al., Phys. Plasmas, in press (2007)]. It is found that as the incompressible limit is approached, the collapse of the null point is suppressed, and an approximately planar current sheet aligned to the fan plane is present instead. This is the case regardless of whether the spine or fan of the null is sheared. Both the peak current and peak reconnection rate are reduced. The results have a bearing on previous analytical solutions for steady-state reconnection in incompressible plasmas, implying that fan current sheet solutions are dynamically accessible, while spine current sheet solutions are not.

The nature of the evolution of the magnetic field, and of current sheet formation, at three-dimensional (3D) magnetic null points is investigated. A kinematic example is presented which demonstrates that for certain evolutions of a 3D null (specifically those for which the ratios of the null point eigenvalues are time-dependent) there is no possible choice of boundary conditions which renders the evolution of the field at the null ideal. Resistive MHD simulations are described which demonstrate that such evolutions are generic. A 3D null is subjected to boundary driving by shearing motions, and it is shown that a current sheet localised at the null is formed. The qualitative and quantitative properties of the current sheet are discussed. Accompanying the sheet development is the growth of a localised parallel electric field, one of the signatures of magnetic reconnection. Finally, the relevance of the results to a recent theory of turbulent reconnection is discussed.

The behaviour of magnetic perturbations of an initially potential three-dimensional equilibrium magnetic null point are investigated. The basic components which constitute a typical disturbance are taken to be rotations and shears, in line with previous work. The spine and fan of the null point (the field lines which asymptotically approach or recede from the null) are subjected to such rotational and shear perturbations, using three-dimensional magnetohydrodynamic simulations. It is found that rotations of the fan plane and about the spine lead to current sheets which are spatially diffuse in at least one direction, and form in the locations of the spine and fan. However, shearing perturbations lead to 3D-localised current sheets focused at the null point itself. In addition, rotations are associated with a growth of current parallel to the spine, driving rotational flows and a type of rotational reconnection. Shears, on the other hand, cause a current through the null which is parallel to the fan plane, and are associated with stagnation-type flows and field line reconnection across both the spine and fan. The importance of the parallel electric field, and its meaning as a reconnection rate, are discussed.