Search

In this paper we propose that flux cancellation on small granular scales ($\lesssim 1000$ km) ubiquitously drives reconnection at a multitude of sites in the low solar atmosphere, contributing to chromospheric/coronal heating and the generation of the solar wind. We analyse the energy conversion in these small-scale flux cancellation events using both analytical models and three-dimensional, resistive magnetohydrodynamic (MHD) simulations. The analytical models -- in combination with the latest estimates of flux cancellation rates -- allow us to estimate the energy release rates due to cancellation events, which are found to be of the order $10^6-10^7$ erg cm$^{-2}$sec$^{-1}$, sufficient to heat the chromosphere and corona of the Quiet Sun and active regions, and to power the solar wind. The MHD simulations confirm the conversion of energy in reconnecting current sheets, in a geometry representing a small-scale bipole being advected towards an intergranular lane. A ribbon-like jet of heated plasma is accelerated upwards, that could also escape the Sun as the solar wind in an open-field configuration. We conclude that through two phases of atmospheric energy release -- pre-cancellation and cancellation -- the cancellation of photospheric magnetic flux fragments and the associated magnetic reconnection may provide a substantial energy and mass flux contribution to coronal heating and solar wind generation., Comment: 19 pages, 11 figures, published in ApJ

The slow solar wind is generally believed to result from the interaction of open and closed coronal magnetic flux at streamers and pseudostreamers. We use 3-dimensional magnetohydrodynamic simulations to determine the detailed structure and dynamics of open-closed interactions that are driven by photospheric convective flows. The photospheric magnetic field model includes a global dipole giving rise to a streamer together with a large parasitic polarity region giving rise to a pseudostreamer that separates a satellite coronal hole from the main polar hole. Our numerical domain extends out to 30 solar radii and includes an isothermal solar wind, so that the coupling between the corona and heliosphere can be calculated rigorously. This system is driven by imposing a large set of quasi-random surface flows that capture the driving of coronal flux in the vicinity of streamer and pseudostreamer boundaries by the supergranular motions. We describe the resulting structures and dynamics. Interchange reconnection dominates the evolution at both streamer and pseudostreamer boundaries, but the details of the resulting structures are clearly different from one another. Additionally, we calculate in situ signatures of the reconnection and determine the dynamic mapping from the inner heliosphere back to the Sun for a test spacecraft orbit. We discuss the implications of our results for interpreting observations from inner heliospheric missions, such as Parker Solar Probe and Solar Orbiter, and for space weather modeling of the slow solar wind., Comment: This article was submitted under duress from an unreasonable reviewer for a separate article

Observations of solar flare ribbons show significant fine structure in the form of breaking wave-like perturbations and spirals. The origin of this structure is not well understood, but one possibility is that it is related to the tearing instability in the flare current sheet. Here we study this connection by constructing an analytical three-dimensional magnetic field representative of an erupting flux rope with a flare current sheet below it. We introduce small-scale flux ropes representative of those formed during a tearing instability in the current layer, and use the squashing factor on the solar surface to identify the shape of the presumed flare ribbons and fine structure. Our analysis suggests there is a direct link between flare-ribbon fine structure and flare current sheet tearing, with the majority of the ribbon fine structure related to oblique tearing modes. Depending upon the size, location and twist of the small-scale flux ropes, breaking wave-like and spiral features within the hooks and straight sections of the flare ribbon can be formed that are qualitatively similar to observations. We also show that the handedness of the spirals/waves must be the same as the handedness of the hooks of the main ribbon. We conclude that tearing in the flare current layer is a likely explanation for spirals and wave-like features in flare ribbons., Comment: ApJ, in press, 16 pages, 14 figures

Aims: Emission line profiles from solar coronal loops exhibit properties that are unexplained by current models. We investigate the non-thermal broadening associated with plasma heating in coronal loops that is induced by magnetic field line braiding. Methods: We describe the coronal loop by a 3D magnetohydrodynamic model of the turbulent decay of an initially-braided magnetic field. From this we synthesize the Fe XII line at 193 Angstroms that forms around 1.5 MK. Results: Key features of current observations of extreme UV lines from the corona are reproduced in the synthesised spectra: (i) Typical non-thermal widths range from 15 to 20 km/s. (ii) The widths are approximately independent of the size of the field of view. (iii) There is a correlation between line intensity and the non-thermal broadening. (iv) Spectra are found to be non-Gaussian, with enhanced power in the wings of order 10-20%. Conclusions: Our model provides an explanation that self-consistently connects the heating process to the observed non-thermal line broadening. The non-Gaussian nature of the spectra is a consequence of the non-Gaussian nature of the underlying velocity fluctuations, interpreted as a signature of intermittency in the turbulence., Comment: Published in A&A

(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., Comment: Accepted in A&A, 5 Figures (1 animated)

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., Comment: Published open access in Physics of Plasmas - go to main journal for animations https://doi.org/10.1063/1.5035489

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., Comment: Accepted in ApJ, will be published gold open access, refer to main journal

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, Comment: Accepted in ApJ, will be published gold open access, refer to main journal for animated content

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., Comment: Pre-print version of article accepted for publication in Solar Physics

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., Comment: To appear in Plasma Physics and Controlled Fusion

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., Comment: To appear in ApJ. 20 pages, 10 figures

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., Comment: Submitted to Physics of Plasmas (12 pages, 12 figures)

The manner in which the rate of magnetic reconnection scales with the Lundquist number in realistic three-dimensional (3D) geometries is still an unsolved problem. It has been demonstrated that in 2D rapid non-linear tearing allows the reconnection rate to become almost independent of the Lundquist number (the `plasmoid instability'). Here we present the first study of an analogous instability in a fully 3D geometry, defined by a magnetic null point. The 3D null current layer is found to be susceptible to an analogous instability, but is marginally more stable than an equivalent 2D Sweet-Parker-like layer. Tearing of the sheet creates a thin boundary layer around the separatrix surface, contained within a flux envelope with a hyperbolic structure that mimics a spine-fan topology. Efficient mixing of flux between the two topological domains occurs as the flux rope structures created during the tearing process evolve within this envelope. This leads to a substantial increase in the rate of reconnection between the two domains., Comment: Submitted to Physics of Plasmas (11 pages, 9 figures)

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., Comment: To appear in ApJ

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., Comment: Accepted for publication in the Astrophysical Journal

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., Comment: 12 pages, 14 figures, to appear in Physics of Plasmas

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., Comment: Accepted for publication in ApJ

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., Comment: 16 pages, 13 Figures. To be published in Astronomy and Astrophysics

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.

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., Comment: 9 pages, 11 figures, accepted for publication in A&A

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., Comment: Accepted for publication in Advances in Space Research. (Based on paper presented at 38th COSPAR meeting, Bremen, 2010.)

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., Comment: To be published in A&A

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., Comment: 18 pages, 9 figures

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., Comment: 16 pages, 10 figures. Accepted for publication in Physics of Plasmas

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., Comment: 23 pages, 9 figures

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., Comment: 17 pages, 8 figures

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., Comment: Submitted to ApJ

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., Comment: 19 pages with 7 figures

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., Comment: Accepted for publication in JGR (Space Physics)

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., Comment: to appear in Physics of Plasmas. This version contains updated figures and references, additional discussion, and typos are fixed. This is the second in a series of papers - the first of which (by the same authors) is located at astro-ph/0701462. A version with higher quality figures can be found at http://www.maths.dundee.ac.uk/~dpontin/

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., Comment: to appear in Phys. Plasmas. A version with higher quality figures can be found at http://www.maths.dundee.ac.uk/~dpontin/ In this replacement version, typos have been corrected, and in addition references and some further discussion added