Why are flare ribbons associated with the spines of magnetic null points generically elongated?
Abstract
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 in the corona. 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, . In this paper we make a detailed analysis of the distribution of in the presence of magnetic nulls. While is formally infinite on both the spine and fan of the null, the decay of away from these structures is shown in general to depend strongly on the null point structure. For the generic case of a nonradiallysymmetric null, decays most slowly away from the spine/fan in the direction in which increases most slowly. In particular, this demonstrates that the extended, elliptical high halo around the spine footpoints observed by \citeauthoryearMasson et al.2009 (Astrophys. J., 700, 559, \citeauthoryearMasson et al.2009) is a generic feature. The asymmetry of the halo of contours around the spine/fan is shown to be strongest for the highest contours, and increases as the null point asymmetry increases. This extension of the halos around the spine/fan footpoints is in general 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 the extent to which 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. In summary, we conclude that both the physics in the vicinity of the null and how this is related to the extension of away from the spine/fan can be used in tandem to understand observational signatures of reconnection at coronal null points.
addressref=aff1,corref,email=d.i.pontin@dundee.ac.uk]\initsD.I.\fnmDavid \lnmPontin\orcid0000000210899270addressref=aff2,email=kg@nbi.ku.dk]\initsK.\fnmKlaus \lnmGalsgaard\orcid0000000188821708addressref=aff3,email=Pascal.Demoulin@obspm.fr]\initsP.\fnmPascal \lnmDémoulin\orcid0000000182156532
1 Introduction
As new generations of solar telescopes allow ever more detailed views of the Sun’s atmosphere, the link between magnetic topological structures and observed sites of energy release becomes increasingly apparent. The magnetic structure of the corona is highly complex over a broad range of scales, as a result of the complex array of magnetic polarities that appear in a continually evolving pattern on the photosphere. The magnetic flux from each polarity region on the photosphere generically connects to many other flux patches of opposite polarity. The structure of the associated coronal magnetic field can appear bewilderingly complex, but advances in theory, modelling and observations have allowed a characterisation of the key features of the 3D structure – such as likely sites for dynamic events to take place. One particular tool for analysing the coronal field structure is the magnetic field line mapping between positive and negative polarity regions of the photosphere. In particular, field lines along which this mapping is discontinuous – usually separatrix surfaces associated with magnetic null points – or has strong gradients – at quasiseparatrix layers (QSLs) – are now known to be likely sites for current accumulation and energy dissipation. These structures are defined and discussed in the following section.
In the last 20 years or so, a wealth of observational evidence has accumulated for energy release at both magnetic null points and QSLs in the form of flares, jets, and bright points \citep[e.g.][]demoulin1997,fletcher2001,mandrini2006,luoni2007,masson2009,liu2011,zhang2012. In each of these cases, the magnetic field structure in the corona must be inferred by employing some extrapolation method that uses the observed photospheric field as a boundary condition. One particular recent focus has been to understand in flaring regions how the observed flare ribbons map to the coronal magnetic field structure – and what one can subsequently deduce about the flare energy release process. In configurations containing QSLs, the footprints of these QSLs have been shown to be colocated with observed ribbons \citepdemoulin1993,dalmasse2015,savcheva2015,savcheva2016. In specific cases a magnetic null point is also present. It has a fan separatrix and spines which define the topology of the magnetic configuration, and since the vicinity of the null is a preferential site for current accumulation and reconnection, so the footpoints of the spine and fan structures are often where the flare ribbons are located. While the fan surface footprint naturally defines elongated ribbons, the spine lines define locally compact regions, so that one would expect compact ribbons – however they are also observed to be elongated. A link is observed between this elongation and the squashing degree surrounding the fan and spines \citep[e.g.][]mandrini2006 – the reason for this link is explored herein. [\citeauthoryearMasson et al.2009] observed a socalled circular flare ribbon associated with the footprint of the separatrix surface of a coronal null point – and they noted again the elongation of the spine footpoint ribbons. Since this observation, a number of further studies have confirmed that these findings are generic \citepwang2012b,yang2015,liu2015.
In this paper we make more concrete the link between the null point magnetic field structure, the geometries of associated features in the field line mapping, and the expected locations of flare ribbons. This is done by analysing the field line mapping in the vicinity of generic 3D null point structures, and relating this to the known properties of current sheet formation and magnetic reconnection around these nulls. The paper is organised as follows. We start in Section 2 by discussing some necessary background on 3D magnetic topology and reconnection. In Sections 3 and 5 the field line mapping in a linear null configuration and coronal separatrix dome configuration are analysed, respectively. In Section 4 the field line mapping is studied in the context of magnetic reconnection around the null point in MHD simulations. We end in Section 6 with a discussion.
2 Background
2.1 Magnetic null points
In this paper we deal with magnetic null points in the solar corona. Such coronal null points have been demonstrated to exist in abundance by various surveys of coronal magnetic field extrapolations, both potential and forcefree \citepregnier2008,edwards2015,freed2015. A magnetic null is simply a location in space at which the magnetic field strength is exactly zero, , and in three dimensions (3D) this condition is met generically only at isolated points. The magnetic field in the vicinity of the null is characterised by a pair of spine field lines that asymptotically approach (or recede from) the null and a fan surface within which field lines radiate away from (or approach) the null point – see [\citeauthoryearLau and Finn1990, \citeauthoryearParnell et al.1996] for a full description, and Figure 1 for a visualisation. The fan surface forms a separatrix surface in the field, distinguishing two volumes of magnetic flux within which the field line connectivity is topologically distinct.
The simplest generic coronal structure involving a magnetic null point occurs when a ‘parasitic’ polarity region on the photosphere is surrounded by polarity region(s) of the opposite sign (and greater total flux). In this configuration a magnetic null is located where the field contributions from the two polarities cancel. The fan separatrix surface then forms a dome structure and separates flux connecting the dominant polarity to the parasitic polarity (beneath the dome) from that which connects from the dominant polarity to locations further away on the photosphere. The effect of reconnection in such a separatrix dome configuration has been considered by [\citeauthoryearEdmondson et al.2010, \citeauthoryearPontin, Priest, and Galsgaard2013]. Null points and their associated separatrix surfaces also occur in many more complicated topological configurations involving multiple null points, separatrix surfaces, and separators (separatrix surface intersections) – see for example [\citeauthoryearPlatten et al.2014]. However, here we restrict our analysis to a single null point, considering a separatrix dome configuration in Section 5.
3D null points are one of the preferential sites for reconnection in the corona. This is because in the perfectlyconducting limit, singular current layers are known to form at the null when rather general perturbations are applied \citeppontincraig2005. Therefore no matter how small the dissipation, nonideal processes will eventually become important as the field around the null point collapses. Note that in an equilibrium there should be zero current at the null, since in general a pressure gradient cannot balance the Lorentz force in the vicinity of the null \citepparnell1997. Any nonzero current at the null will in general lead to a Lorentz force that drives the null point to collapse to form a current sheet \citepklapper1996,pontincraig2005,fuentes2012,craig2014. The resulting reconnection may take different forms, the most general mode of reconnection being spinefan reconnection that is associated with transfer of magnetic flux across the separatrix surface – this process permitting in principle the release of significant stored magnetic energy \citepantiochos1996,pontinbhat2007a,pariat2009,pontin2013. There are also two other modes of reconnection at 3D nulls – torsionalspine and torsionalfan reconnection, that involve a rotational slippage of field lines around the spine but involve no flux transfer across the separatrix \citeppontin2011b.
2.2 The squashing factor and quasiseparatrix layers
The principal reason why magnetic null points were first proposed as sites of current accumulation and therefore magnetic reconnection in 3D is that at the null the field line mapping is discontinuous. It is now well established that 3D reconnection may also occur in the absence of null points or separatrices, and in particular natural sites for the formation of intense current layers are regions in which the field line mapping exhibits strong gradients. Analysis of these gradients is typically performed by evaluating the (covariant) squashing degree, defined for planar boundaries by
(1) 
where and are field line footpoints on the ‘launch’ boundary, and and are the footpoint locations on the ‘target’ boundary, see [\citeauthoryearTitov, Hornig, and D\a’emoulin2002]. The general expression for nonplanar boundaries can be found in Equations (1114) of [\citeauthoryearTitov2007]. Note that the denominator in Eq. (1) can also be represented by where and are the field components normal to the boundaries at the launch and target footpoints, respectively. Numerically it is usually more stable to use this expression in practice. One potential weakness of as defined in Eq. (1) is that the values obtained depend on the orientation at which field lines intersect the launch and target boundaries. An alternative formulation is the perpendicular covariant squashing factor, , as defined in equations (3036) of [\citeauthoryearTitov2007], which removes such projection effects by evaluating the mapping deformation for infinitesimal perpendicular planes at the locations of the launch and target surfaces.
Bundles of magnetic field lines along which or have large values (, the minimum value, obtained for a uniform field) are known as quasiseparatrix layers, or QSLs. It has been demonstrated that these are natural locations for accumulation of intense currents, using both modelling approaches \citeptitov2003,galsgaard2003,aulanier2005,effenberger2011 and solar observations \citep[e.g.][]demoulin1997,mandrini2006. The term QSL comes from the fact that true separatrices can be thought of as a limiting case of a QSL [?, see]for a detailed exposition]demoulin2006. In particular, due to the discontinuity in the field line mapping at a separatrix, is by definition infinite there. In addition, must also be large but finite in the region adjacent to the separatrix – and on this we focus in Section 3.
2.3 The nature of 3D reconnection
In order to understand energy release mediated by magnetic reconnection in 3D, it is important to understand a key property of 3D reconnection. Specifically, in contrast to the 2D case reconnection in 3D always occurs in a finite volume. That is, rather than field lines breaking and rejoining at a single point (the Xpoint) as in 2D, field lines change connectivity continuously throughout the (finitesized) nonideal region \citeppriest2003a. This nonideal region is in general any region within which the electric field component parallel to the magnetic field () is nonzero, and for which , the integral being evaluated along field lines \citepschindler1988. A consequence of the breaking and rejoining of field lines throughout the nonideal region is as follows. For reconnection in the absence of separatrices, for example at QSLs, there is an everywhere continuous ‘flipping’ or ‘slipping’ of reconnecting field lines \citeppriest1995. One can also distinguish this flipping motion further, to slipping or sliprunning, depending on whether the velocity of the apparent field line motion is sub or superAlfvénic, respectively \citepaulanier2006. Such slipping motions are now observed during energy release in the corona \citepaulanier2007,sun2013,li2014,dudik2014.
All of the above statements hold true for 3D null point reconnection. In particular, the reconnection happens not only at the null point itself, but throughout the nonideal region (current sheet) surrounding the null \citeppontin2004,pontinhornig2005. Therefore what we call ‘null point reconnection’ is more precisely reconnection that occurs within a finite region (the current layer) surrounding a null – the importance of the null being that it is a favourable site for intense currents to develop (as described above). There is still a continuous change of field line connectivity within the current sheet, and an apparent flipping motion of field lines. However, in the presence of a null there is also one discontinuous jump of connectivity, for every field line that is reconnected through either the spine of the null or the separatrix (fan) surface. For a given reconnection event, if there is a null point in the nonideal region the field line velocity must by necessity be ‘sliprunning’ near the null, since it is infinite at the spine/fan. For the case of a single null point in a separatrix dome configuration, that we examine below, the expected patterns of the field line slippage motions were described in detail by [\citeauthoryearPontin, Priest, and Galsgaard2013].
3 Linear null point
3.1 Preliminaries
In the following section we examine the distribution of the squashing degree, , in the vicinity of null points. First we consider the simplest case of a linear null point. Note that any generic 3D null point can be represented locally (i.e. sufficiently close to the null) by this linearisation – and conversely for the null point to be topologically stable, the linear term in the Taylor expansion of about the null point must be nonzero \citephornig1996. (Topological stability implies that an arbitrary perturbation does not destroy the topology, in contrast to the case where the first nonzero term is the quadratic term – those higher order nulls are topologically unstable since an arbitrary perturbation of changes the topology.) Note that there are two factors that influence the variation of in the vicinity of a null point and associated separatrix: one is that is formally infinite for spine and fan field lines due to the discontinuity in the field line mapping. There is then a characteristic decay of away from these field lines. The other factor is that field lines in the fan surface typically become oriented parallel to one of the fan eigenvectors (corresponding to the largest fan eigenvalue) at larger distances. This naturally means there is a stronger divergence of field lines away from the fan eigenvector direction corresponding to the smaller fan eigenvalue. This leads to a rotational asymmetry of the decay of away from the spine and fan – as demonstrated below.
Evaluation of the squashing factor requires that we select two surfaces that each field line intersects once and only once. We consider two cases: in the first case we take both boundaries to be planar (Figure 1a), and in the second case we take a plane of constant and a circular cylinder surface (see Figure 1b). When both surfaces are planar we can obtain exact expressions for and , whereas for the cylindrical boundary we must evaluate them numerically.
3.2 Squashing factor between two planar boundaries
Consider an equilibrium magnetic null point (zero current). The field can be represented by
(2) 
with , where we have chosen to orient the coordinate system such that the spine lies along the axis, the fan surface is coincident with the plane, and the two eigenvectors of are parallel to the and axes. The corresponding eigenvalues are unequal when , and in the plane the field strength increases most quickly away from the origin along the direction of the eigenvector associated with the largest eigenvalue – the direction for . We refer below to this direction as the strong field direction in the fan, and correspondingly to the orthogonal direction as the weak field direction. In the corona, one would in general expect field lines in the strong field direction to connect to the strongest nearby photospheric flux concentrations (though this need not necessarily be the case). Taking two planar boundaries as shown in Figure 1(a), we can obtain a mathematical expression for how decays away from (say) the spine along the strong and weak field directions in the fan.
Here we take the ‘launch plane’ for field lines to be and the target plane to be ( and constant). That is, we consider the mapping . For simplicity we only examine how decays from the spine along the axis (setting ). This corresponds to the strong field direction in the fan if , and the weak field direction for . We examine the decay in arbitrary directions in the next section. Setting , we obtain the following expression for as a function of , the distance of the launch footpoint from the spine
(3) 
see Appendix A. We see that, as expected, as (the spine footpoint). We are interested in the behaviour for small , and since , the relevant term close to the spine is the first one, so that
(4) 
This implies that in the weak field region () should be larger (for fixed ) since we have a larger negative exponent. This is to be expected given the strong field line divergence in the plane in this region. Hence a level curve of would be expected to be elongated along the weak field direction. Note that for the rotationally symmetric case () we have that decays like .
Consider now the distribution of (see Section 2.2). Evaluating on the axis as before we obtain the following expression
(5) 
(Appendix A). For , we have the same dominant scaling in as before, specifically
(6) 
That the behaviour of and is identical in this plane is expected since in this plane as we get close to the spine and fan the field lines intersect the boundaries approximately perpendicular.
We now examine the dimensions of contours of on the spine boundary (). Specifically, we rearrange Equation (3) to find the radius at which . Assuming that we can directly invert Equation (6) to obtain
(7) 
One can obtain from this the same information along the axis by making the replacement to give .
As shown in Figure 2(a), for the spacing of contours drops off more slowly in the weak field direction, consistent with above.
Now let us analyse the asymmetry in the contour produced by the asymmetry in null eigenvalues in more detail, by examining the relative decay of from the spine footpoint along the two axes. Using Eq. (7) we have that
(8) 
In Figure 2(b) we plot the ratio for and three particular values of . We see that, as expected, this aspect ratio increases as the magnetic field asymmetry increases. Furthermore, we observe that for high values the distribution is highly ‘noncircular’ even for moderate values of , while lower level contours of only show high asymmetry when is far from 0.5 (note that Figure 2(b) is a loglog plot). Note also that a full range of degrees of null point eigenvalue ratios (corresponding to ) is obtained in solar extrapolations, and even for a moderate asymmetry (values of close to 0.5) the eccentricity is significant. For example [\citeauthoryearDémoulin, Hénoux, and Mandrini1994] analysed 6 nulls with – , the configuration analysed by [\citeauthoryearMasson et al.2009] contained a null point with , and [\citeauthoryearFreed, Longcope, and McKenzie2015] carried out an extensive survey of potential field extrapolations from three years’ worth of magnetogram data, identifying 1924 coronal null points – choosing the orientation such that their data shows a relatively uniform distribution of values between 0.1 and 0.5, with a mean value of of 0.26 and a standard deviation of 0.11.
3.3 Squashing factor for a surface encircling the fan
The calculations of the previous section allow us to visualise the distribution along the coordinate directions. However, since the vertical planar boundary intersects only a subset of the fan field lines, they do not allow us to visualise the distribution all around the spine or fan footpoints. In order to do this we must choose a target plane that intersects all fan field lines (such that all field lines passing close to the null intersect both the launch and target surfaces).
As such, we now calculate between a planar ‘launch’ surface at intersecting the spine, and a cylindrical ‘target’ surface at encircling the fan, constant (Figure 1b). That is, we study the mapping generated by the field lines . In this case we are unable to obtain a full analytical expression for the field line mapping and its inverse, and thus we evaluate numerically. We integrate between and field lines from a rectangular grid of starting points at to obtain their intersections with , then perform derivatives of the mapping using a fourthorderaccurate centred difference over that grid. In this case we evaluate , since especially in the regions of strongly diverging field lines (the weak field region, around the axis for ) the field lines intersect the circular cylinder surface far from perpendicular.
The maps at for different are presented in Figure 3. As predicted by the planar boundary analysis, we see a stretching of the contours along the weak field direction in the fan ( for ). We also see that the contours do not form simple ellipses, but are pinched in the middle so that their maximum extension is at finite values. Therefore in order to measure the asymmetry of the distribution around the spine, it is arguably most useful to determine the largest extent of these contours along the and directions, rather than simply examining the profile along the coordinate axes. In Figure 4 we plot (with circles and solid lines) the contour aspect ratio defined as the maximum contour extent along divided by the maximum extent along (over all ).
The plot demonstrates the same trends as observed in the previous section, specifically that the aspect ratio is greater for both increased null point asymmetry (smaller for ) and for higher levels. Note that the qualitative features discussed above are present if one considers instead of .
We can now perform the same calculation procedure but with the launch and target boundaries reversed, in order to find the pattern of in the vicinity of the fan surface – shown in Figure 5. We observe that decreases in some locations and increases in others as we increase the asymmetry of the null. The weak field region is along the axis, which corresponds to , by the usual convention. Again the widest contours are located in the vicinity of these weak field directions. That is, for a given distance from the separatrix (height ), is largest along the weak field direction.
The results above are entirely consistent with those of e.g. [\citeauthoryearMasson et al.2009]. In particular, we see that one obtains extended (or ) contours around the spine for even a moderate degree of null point asymmetry. Of course other global features of the field could well also contribute, but the figures show that the high region in which the null is embedded is not a special additional feature of the particular field studied by [\citeauthoryearMasson et al.2009], but is rather a natural consequence of having a coronal null that is not rotationally symmetric.
4 Current layer formation in MHD simulations, and flare ribbon locations
4.1 Simulation setup and results
In this section we ask the question: what is the relation between current layers that form at 3D nulls and the distribution of the squashing factor identified in the previous section, and to what extent can we expect the (or ) profile to predict where flare ribbons might be observed? We consider resistive MHD simulations similar to those of [\citeauthoryearGalsgaard and Pontin2011]. Specifically, at in our simulations we have a linear magnetic null point of the form
(9)  
which reduces to Eq. (2) for . For , the fan plane eigenvectors are rotated by an angle with respect to the coordinate axes. The simulation domain is and we set . We apply a driving velocity on the domain boundaries that advects the spine footpoints in opposite directions on opposite boundaries, in the direction. Specifically, we let
(10)  
with . Outside these regions on all boundaries and is linetied. The plasma density and pressure are initially uniform, , , so that the driving velocity is highly subsonic and subAlfvénic. In addition, the resistivity , also spatially uniform, throughout.
As the simulation proceeds the stress injected by the boundary driving focusses around the null point – the null point collapses and a current sheet forms around it \citep[see][]pontinbhat2007a. We examine three different simulations; in the first the null point is rotationally symmetric, while in the second and third we take an asymmetric null with corresponding to frame (c) in Figures 3 and 5. Specifically, we use parameters in Eq. (9) as follows: simulation 1: ; simulation 2: ; simulation 3: . Hereafter we analyse the state reached in the simulations at time , this constituting a representative time by which the current layer has formed and reconnection is underway.
The current layer formed at is shown for each simulation in Fig. 6. In each case the current is maximum at the null point. Note that the current layer extends from the high region into the low region in all simulations. In cases (a) and (b) the current extends perpendicular to the direction of the driving motion, which is also the weak field direction for case (b), while in case (c) there is a competition between the driving and weak field directions. However for this moderate value of , the extension is not very pronounced. It is already known that when the null point is asymmetric, the current tends to spread preferentially along the weak field direction in the fan plane, at least when the driving has a nonzero component perpendicular to this direction \citepalhachami2010,galsgaard2011b, since the weak field in this region is less able to withstand the field collapse. Note that one can in general expect currents to extend along this weak field direction also due to the strong values further from the separatrix there, which reflect the variety of field line connectivities present nearby. Such field lines are anchored in distant locations at the boundary, so they typically experience different magnetic stress from the boundary motions, so different perturbed . However, at this time in our simulations there has not been time for any communication with the and boundaries, and so the current accumulation is associated only with stresses being applied from the spine boundaries and the local collapse dynamics.
4.2 Predicting flare ribbon locations
A proper diagnosis of expected locations of flare ribbons would require a selfconsistent modelling of particle acceleration in a null point current sheet. This is yet to be done – most existing studies use simplified analytical models that do not properly represent the structure of the magnetic field and current layer. We emphasise then that the following analysis based on the MHD approximation is a crude first step toward predicting the expected location of energetic particles.
First, we should understand the mechanisms by which the acceleration could occur. Null points have been proposed as efficient particle acceleration sites first because the geometry of the field around the null naturally allows for magnetic mirroring, and particle acceleration by gradient and curvature drifts \citepvekstein1997,petkaki2007,guo2010,stanier2012. The particle dynamics in the vicinity of the null may indeed be inherently chaotic \citepmartin1986. In addition, when a current layer is present at the null during reconnection, there can be direct acceleration by the associated electric field. This was observed to be the dominant acceleration mechanism in the PIC simulations of [\citeauthoryearBaumann, Haugbølle, and Nordlund2013], who studied null point reconnection in the corona using a configuration similar to that of [\citeauthoryearMasson et al.2009]. To understand the resulting particle deposition patterns, one must first understand the structure of the electric current layer at the null. This electric current distribution is determined in general by a combination of factors. During spinefan reconnection, the spine and fan of the null point locally collapse towards one another (see Figure 7) as in our simulations. This collapse of the null occurs in general when a shear perturbation of either the spine or the fan occurs \citeppontinbhat2007a. The plane in which this collapse occurs (plane that contains the deformed spine line – see Figure 7) is determined both by the perturbation that drives the collapse and the null point structure. The associated current sheet that forms has a current vector that at the null is oriented perpendicular to the plane of collapse, see Figure 7. Thus (in resistive MHD) the parallel electric field is oriented along the fan surface, perpendicular to the spine and the plane of the null point collapse \citeppontinbhat2007a. Hence, we expect a strong acceleration layer near the null in the fan plane, and thus particle deposition in the vicinity of the fan surface footprint (of oppositely charged particles on opposite sides). It is also possible that particles accelerated towards the null in this layer may follow the field lines out along the spines. However, when they reach the vicinity of the null point they become effectively demagnetised, and most particles are simply accelerated ‘across’ the current layer, and out along the fan (rather than being deflected up the spine). Exactly this effect was observed by [\citeauthoryearBaumann, Haugbølle, and Nordlund2013], who noticed very few particles accelerated out along the spines.
We should note that there are at least three factors that could cause enhanced acceleration along the spine structures as well, to create the spine footpoint ribbons observed by [\citeauthoryearMasson et al.2009] and others. First, at solar parameters the reconnection process around the site of the original null is likely to be significantly more complex than in the simple models where a single laminar current layer is present. Indeed this current sheet is susceptible to a tearingtype instability that leads to a fragmented current layer containing many nulls, as described by [\citeauthoryearWyper and Pontin2014]. In such a configuration, the vicinity of the original null becomes highly turbulent, and one would expect efficient particle scattering along both the largescale spine and fan directions. Second, if there is some largescale rotational external motion, this can drive torsional spine reconnection, associated with a component of current parallel to the spine \citeppontingalsgaard2007 which can accelerate particles along the spine \citephosseinpour2014. Finally, one could expect strong mirroring of particles close to the fan footpoints to lead to a distribution of particles also around the spine footpoints (note that the PIC simulations of [\citeauthoryearBaumann, Haugbølle, and Nordlund2013] did not cover the domain all the way to the photosphere).
Based on the above considerations, there exists no unequivocal way of diagnosing general expected particle deposition footprints. However, independent of the details of the acceleration mechanism, if the acceleration happens during reconnection at the null point, then as particles move away from the null they will become magnetised, and would be expected to be observed in the vicinity of either the spine or fan footpoints. For our simulations we make the following basic assumption: Particles will be accelerated in some manner in the vicinity of the current sheet around the null. Thus, we estimate expected deposition patterns by tracing field lines from the current sheet to find their intersections with the boundaries. To predict particle deposition locations on the spine boundaries (), we therefore perform the following procedure: We select an array of points that lie on a given current contour level within the domain, and trace field lines from each of these points to the spine boundaries. We then compare the intersections of these field lines with the distribution on the boundary to determine whether they match, i.e. whether the distribution can be expected to give a good prediction of the geometry of particle deposition signatures.
We first evaluate in our simulations using launch boundary and target boundary a circular cylinder of radius 2.8, as in Section 3.3. These contours are plotted, together with footpoint locations for field lines threading the current layer, in Fig. 8. Consider first simulation 1 where the field is initially rotationally symmetric (Figs. 8a, 6a). We see that there is some  asymmetry in the 3D current distribution (Fig. 6a) that is due to the orientation of the boundary driver. Here the null point field was initially rotationally symmetric, and when we map field lines from the current layer to the boundary we see a similar degree of asymmetry in the 3D current layer and its 2D projection on the boundary (Fig. 8a).
As a contrast, now compare Figs. 6(b), 8(b). There is little difference in the asymmetry of the 3D current density distribution from Fig. 6(a). However, the field geometry means that field lines approach the spine more slowly along the weak field direction, so that the projection of is elongated along this direction. This reinforces the fact that the current preferentially spreads along that direction – but it is the 3D field line geometry rather than the current layer geometry that has the major effect on the projected map. This projection of along the field lines now has a comparable geometry to the contours, although the projected map does not exhibit such high eccentricity/asymmetry. We observe a similar behaviour when we examine Figure 8(c), except that in this case (where the boundary driving is at a finite angle to the null point eigenvectors) the two distributions appear to be rotated with respect to one another – more pronounced for low current contour levels. These results indicate that a complicated combination of the driving geometry, the field geometry, and the current intensity in the current sheet (itself dependent on plasma parameters and the driving of the system) will influence the expected particle precipitation locations (even using this simple estimate for these locations).
Now consider the boundaries intersected by fan field lines. The same method as applied before using field line mapping is not useful in determining the angular distribution of expected particle locations, since by definition every field line of the fan connects back to the null and therefore the maximum current region. What is more important for this angular distribution is the orientation of the electric field in the acceleration region together with the global field structure. Let us make the following simple considerations. Particle acceleration along fan field lines will occur through direct acceleration by the DC electric field (or other mechanisms as mentioned above). It is expected that at the null this electric field is directed predominantly towards the weak field direction (perpendicular to the plane of null collapse, see above). However, the field lines diverge away from this direction, and converge towards the strong field region. Thus the particles – which becomes remagnetised as they are accelerated away from the null – are naturally channelled along the field lines into the neighbourhood of the strong field direction. This is the region in which (or ) falls off more quickly away from the separatrix – and thus we expect that particles accumulate around the strong field regions of the fan footprint, corresponding to the locations where the contours are narrowest about the fan.
5 Effect of the global field
The above sections showed that the asymmetry of the field in the local vicinity of the null can have a profound effect on , and on the mapping of field lines from the current layer. However, other features of the global field can clearly distort this picture. In this section we return to an equilibrium field and examine the effect of the global coronal geometry, to determine whether the above results regarding the distribution asymmetry carry through beyond the linear null point field. We consider a null point in a separatrix dome configuration, as shown in Figure 9(a).
The magnetic field is potential, and on the photosphere corresponds to a magnetic dipole, one polarity of which contains an embedded ‘parasitic polarity’. This field is constructed by placing four magnetic point charges at locations outwith our domain of interest. Specifically, we restrict our studies to the halfspace , where represents the photosphere, and place all point charges at . The magnetic field is given by
(11) 
where are the locations and are the strengths of the point charges. Here we take and , , , . The charges located at and are associated with the parasitic polarity around . The parameter controls the separation of these two charges. When they are coincident, and the parasitic polarity is approximately circular. However, as is increased, the parasitic polarity becomes increasingly stretched (see Figure 9). This elongation of the parasitic polarity means that the field strength becomes less homogeneous in the vicinity of the fan footprint, and also that the field asymmetry at the null – as measured by the fan plane eigenvalues – increases. Here we analyse the magnetic topology for three different values of : leads to an approximately symmetric field in the local vicinity of the null that corresponds closely to the symmetric case for the linear field, that gives fan eigenvalues that correspond to for the linear null, and that corresponds to .
As shown in Figure 9(b), when the null point is asymmetric the level curves of on the photosphere form extended structures. As is increased, so too the asymmetry of these level curves increases – both around the inner and outer spine footpoints – see Figure 10. To determine whether the same scaling with the null asymmetry as above is observed we measure the aspect ratio of these ribbons, defined as follows. The ‘length’ of a given contour is defined as the maximum extent along any line passing through the spine line footpoint (location at which ). The width is then defined as the maximum extent of the contour along any line perpendicular to this, analogous to section 3.3. The aspect ratio, being the ratio of the length over the width, is calculated for different contours for each value of , and the results are plotted in Figure 4. What we see is that there are clear differences in the values of the calculated aspect ratio, both between the inner and outer spine footpoints, and to the results for the linear field (crosses, squares, and circles in Figure 4, respectively). However, as shown in Eq. (8) the exact value of the aspect ratio depends on the locations of the launch and target footpoints (for the linear null and ), thus we would not expect an exact agreement. Moreover, we note that the overall scaling of the aspect ratio with the level is rather well reproduced between the linear field and coronal null point field, and that for the higher contour levels considered the aspect ratios are in rather good agreement.
6 Discussion and conclusions
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 are preferential locations for reconnection due to the discontinuity of the field line mapping. This field line mapping also exhibits strong gradients adjacent to the separatrix and spine field lines, that we have analysed here using the squashing factor (and ). Understanding the distribution of in the presence of separatrices is of timely importance due to the increasing use of calculation of maps in analysing the coronal field topology. While a map of the distribution shows the presence of both true separatrices and finite QSLs, one should note that the physics of current layer formation / energy storage is critically different between a high region containing a separatrix and one that does not. In particular, current singularities are known to form in the ideal limit in the presence of separatrices \citeppontincraig2005. Thus reconnection onset is inevitable irrespective of the dissipation (though may be ‘slow’ in an energy storage phase). By contrast, the current layers that form at QSLs are probably finite \citepcraig2014,effenberger2016, with the onset of reconnection at coronal parameters then requiring a thinning of the QSL and current layer during the energy storage phase \citepaulanier2005,demoulin2006. What is clear is that in the case of both null points (separatrices) and QSLs, the current layer formation and eventual dynamics are crucially dependent on the driving of the system, for example from the photosphere.
In this paper, we have made a detailed analysis of the distribution of in the presence of magnetic nulls and their associated separatrices. The main results can be summarised as follows.

It is generically the case that is not uniformly distributed around the spine and fan footpoints. Specifically, a generic null point is not rotationally symmetric, and while is infinite formally on both the spine and fan of the null, it decays most rapidly away from the spine/fan in the direction in which increases most rapidly. When a linearisation of the null is performed (this linearisation characterising the local topology of the field for any topologically stable null \citephornig1996), this direction corresponds to the eigenvector of the largest fan eigenvalue.

The result of the above is that contours of are broadest along the direction of the eigenvalue with smallest fan eigenvalue – denoted herein as the ‘weak field direction’. In particular, this demonstrates that the extended, ellipticallike high halo around the spine footpoints observed by, e.g., [\citeauthoryearMasson et al.2009, \citeauthoryearSun et al.2013] is not a special feature of the particular observations, but is a generic feature when a coronal null is present whose fan eigenvalues are not equal (i.e. when the field strength is not homogeneous around the fan footprint).

The asymmetry of the halo of contours around the spine/fan increases as the null point asymmetry (measured by the ratio of the eigenvalues) increases. Furthermore, for a given null point asymmetry, the stretching of the contours is most extreme for the highest contour levels.

When the global field geometry (beyond the linear field region) is considered, the exact aspect ratios of the contours are modified from the simple linear null case, but the core of the distribution of still reflects the conditions around the null. This is especially true for high contour levels.

As a first approximation for understanding why the geometry of flare ribbons is observed to agree well with the geometry of the halo in circular ribbon flares \citep[e.g.][]masson2009, we analysed MHD simulations of null point reconnection. We traced field lines through the current layer, and analysed the relationship between their intersections with the boundary and the contours on the boundary. While no ontoone relation was found, we showed that field lines traced from the core of the current layer match rather well with the highest contours. Thus, particularly for the kernels of the flare ribbons, the distribution should in general be expected to predict well the location and orientation of the ribbons.
It is well established that an understanding of the null point structure and its relation to the driving of the system is crucial for determining the current layer formation at the null and associated dynamics. We have shown here that this null point structure, defined by its local eigenvectors and eigenvalues, is intrinsically linked to the distribution of away from the spine/fan. Furthermore, this extension of the halos around the spine/fan footpoints is in general important for diagnosing the regions of the photosphere that are magnetically connected to any current layer that forms at the null. If we hypothesise this current layer to be a primary site of particle acceleration, this provides predictive properties for e.g. flare ribbon formation. We conclude that the physics in the vicinity of the null and how this is related to the extension of away from the spine/fan can be used in tandem to understand observational signatures of reconnection at coronal null points.
The authors acknowledge fruitful discussions with G. Valori, E. Pariat, P. Wyper, E. Priest and M. Janvier. D.P. is grateful for financial support from the Leverhulme Trust and the UK’s STFC (grant number ST/K000993). The work by K.G. was supported by a research grant (VKR023406) from VILLUM FONDEN.
Appendix A Calculation of and for the linear null
For the linear null point magnetic field of Equation (2), the field line equations may be solved to obtain parametric equations for the field lines;
(12) 
where at . Now, set on the plane , intersecting the spine, so that . Then we can eliminate in the above equations, to obtain
(13) 
We now choose the ‘target plane’ to be , constant. Finally, identifying in Eq. (1) with , the required derivatives may be obtained. A little algebra leads to the following expression for :
(14) 
which reduces to Equation (3) for .
Evaluation of for the same planar boundaries as above requires that we calculate
(15) 
where is the Kronecker delta, and summation over repeated indices is assumed \citeptitov2007. In addition, , , are evaluated at the target boundary , and are evaluated at the launch plane . The resulting expression is too lengthy to reproduce here, but reduces to Equation (5) for .
References
 AlHachami, A.K., Pontin, D.I.: 2010, Magnetic reconnection at 3D null points: effect of magnetic field asymmetry. Astron. Astrophys. 512, A84. DOI.
 Antiochos, S.K.: 1996, Solar Drivers of Space Weather. In: Balasubramaniam, K.S., Keil, S.L., Smartt, R.N. (eds.) Solar Drivers of the Interplanetary and Terrestrial Disturbances, Astronomical Society of the Pacific Conference Series 95, 1.
 Aulanier, G., Pariat, E., Démoulin, P.: 2005, Current sheet formation in quasiseparatrix layers and hyperbolic flux tubes. Astron. Astrophys. 444, 961. DOI.
 Aulanier, G., Pariat, E., Démoulin, P., Devore, C.R.: 2006, Sliprunning reconnection in quasiseparatrix layers. Solar Phys. 238, 347. DOI.
 Aulanier, G., Golub, L., DeLuca, E.E., Cirtain, J.W., Kano, R., Lundquist, L.L., Narukage, N., Sakao, T., Weber, M.A.: 2007, Slipping magnetic reconnection in coronal loops. Science 318, 1588. DOI.
 Baumann, G., Haugbølle, T., Nordlund, Å.: 2013, Kinetic Modeling of Particle Acceleration in a Solar Nullpoint Reconnection Region. Astrophys. J. 771, 93. DOI.
 Craig, I.J.D., Pontin, D.I.: 2014, Current Singularities in Linetied Threedimensional Magnetic Fields. Astrophys. J. 788, 177. DOI.
 Dalmasse, K., Chandra, R., Schmieder, B., Aulanier, G.: 2015, Can we explain atypical solar flares? Astron. Astrophys. 574, A37. DOI.
 Démoulin, P.: 2006, Extending the concept of separatrices to QSLs for magnetic reconnection. Adv. Space Res. 37, 1269. DOI.
 Démoulin, P., Hénoux, J.C., Mandrini, C.H.: 1994, Are null magnetic points important in solar flares? Astron. Astrophys. 285, 1023.
 Démoulin, P., van DrielGesztelyi, L., Schmieder, B., Henoux, J.C., Csepura, G., Hagyard, M.J.: 1993, Evidence for magnetic reconnection in solar flares. Astron. Astrophys. 271, 292.
 Démoulin, P., Bagala, L.G., Mandrini, C.H., Hénoux, J.C., Rovira, M.G.: 1997, Quasiseparatrix layers in solar flares. II. Observed magnetic configurations. Astron. Astrophys. 325, 305.
 Dudík, J., Janvier, M., Aulanier, G., Del Zanna, G., Karlický, M., Mason, H.E., Schmieder, B.: 2014, Slipping Magnetic Reconnection during an Xclass Solar Flare Observed by SDO/AIA. Astrophys. J. 784, 144. DOI.
 Edmondson, J.K., Antiochos, S.K., DeVore, C.R., Lynch, B.J., Zurbuchen, T.H.: 2010, Interchange Reconnection and Coronal Hole Dynamics. Astrophys. J. 714, 517. DOI.
 Edwards, S.J., Parnell, C.E.: 2015, Null Point Distribution in Global Coronal Potential Field Extrapolations. Solar Phys. 290, 2055. DOI.
 Effenberger, F., Craig, I.J.D.: 2016, Simulations of 3D Magnetic Merging: Resistive Scalings for Null Point and QSL Reconnection. Solar Phys. 291, 143. DOI.
 Effenberger, F., Thust, K., Arnold, L., Grauer, R., Dreher, J.: 2011, Numerical simulation of current sheet formation in a quasiseparatrix layer using adaptive mesh refinement. Physics of Plasmas 18(3), 032902. DOI.
 Fletcher, L., Metcalf, T.R., Alexander, D., Brown, D.S., Ryder, L.A.: 2001, Evidence for the flare trigger site and threedimensional reconnection in multiwavelength observations of a solar flare. Astrophys. J. 554, 451. DOI.
 Freed, M.S., Longcope, D.W., McKenzie, D.E.: 2015, ThreeYear Global Survey of Coronal Null Points from PotentialFieldSourceSurface (PFSS) Modeling and Solar Dynamics Observatory (SDO) Observations. Solar Phys. 290, 467. DOI.
 FuentesFernández, J., Parnell, C.E.: 2012, Magnetohydrodynamics dynamical relaxation of coronal magnetic fields. iii. 3d spiral nulls. Astron. Astrophys. 544, A77. DOI.
 Galsgaard, K., Pontin, D.I.: 2011, Current accumulation in an asymmetric 3D null point: Dependence on the angle between the spine driver and the strongest fan eigenvector. Astron. Astrophys. 534, A2. DOI.
 Galsgaard, K., Titov, V.S., Neukirch, T.: 2003, Magnetic pinching of hyperbolic flux tubes. II. Dynamic numerical model. Astrophys. J. 595, 506. DOI.
 Guo, J.N., Büchner, J., Otto, A., Santos, J., Marsch, E., Gan, W.Q.: 2010, Is the 3D magnetic null point with a convective electric field an efficient particle accelerator? Astron. Astrophys. 513, A73. DOI.
 Hornig, G., Schindler, K.: 1996, Magnetic topology and the problem of its invariant definition. Phys. Plasmas 3, 781.
 Hosseinpour, M., Mehdizade, M., Mohammadi, M.A.: 2014, Comparison of test particle acceleration in torsional spine and fan reconnection regimes. Phys. Plasmas 21(10), 102904. DOI.
 Klapper, I., Rado, A., Tabor, M.: 1996, A lagrangian study of dynamics and singularity formation at magnetic null points in ideal threedimensional magnetohydrodynamics. Phys. Plasmas 3(11), 4281.
 Lau, Y.T., Finn, J.M.: 1990, Three dimensional kinematic reconnection in the presence of field nulls and closed field lines. Astrophys. J. 350, 672.
 Li, T., Zhang, J.: 2014, Slipping Magnetic Reconnection Triggering a Solar Eruption of a Triangleshaped Flag Flux Rope. Astrophys. J. Lett. 791, L13. DOI.
 Liu, C., Deng, N., Liu, R., Lee, J., Pariat, É., Wiegelmann, T., Liu, Y., Kleint, L., Wang, H.: 2015, A Circularribbon Solar Flare Following an Asymmetric Filament Eruption. Astrophys. J. Lett. 812, L19. DOI.
 Liu, W., Berger, T.E., Title, A.M., Tarbell, T.D., Low, B.C.: 2011, Chromospheric Jet and Growing “Loop” Observed by Hinode: New Evidence of Fanspine Magnetic Topology Resulting from Flux Emergence. Astrophys. J. 728, 103. DOI.
 Luoni, M.L., Mandrini, C.H., Cristiani, G.D., Démoulin, P.: 2007, The magnetic field topology associated with two M flares. Adv. Space Res. 39, 1382. DOI.
 Mandrini, C.H., Démoulin, P., Schmieder, B., Deluca, E.E., Pariat, E., Uddin, W.: 2006, Companion event and precursor of the X17 flare on 28 October 2003. Solar Phys. 238, 293. DOI.
 Martin, R.F.J.: 1986, Chaotic particle dynamics near a twodimensional magnetic neutral point with application to the geomagnetic tail. J. Geophys. Res. 91, 11985.
 Masson, S., Pariat, E., Aulanier, G., Schrijver, C.J.: 2009, The nature of flare ribbons in coronal nullpoint topology. Astrophys. J. 700, 559. DOI.
 Pariat, E., Antiochos, S.K., DeVore, C.R.: 2009, A model for solar polar jets. Astrophys. J. 691, 61. DOI.
 Parnell, C.E., Smith, J.M., Neukirch, T., Priest, E.R.: 1996, The structure of threedimensional magnetic neutral points. Phys. Plasmas 3(3), 759.
 Parnell, C.E., Neukirch, T., Smith, J.M., Priest, E.R.: 1997, Structure and collapse of threedimensional magnetic neutral points. Geophys. Astrophys. Fluid Dynamics 84, 245.
 Petkaki, P., MacKinnon, A.L.: 2007, Particle acceleration by fluctuating electric fields at a magnetic field null point. Astron. Astrophys. 472(2), 623. DOI.
 Platten, S.J., Parnell, C.E., Haynes, A.L., Priest, E.R., Mackay, D.H.: 2014, The solar cycle variation of topological structures in the global solar corona. Astron. Astrophys. 565, A44. DOI.
 Pontin, D.I.: 2011, Threedimensional magnetic reconnection regimes: A review. Adv. Space Res. 47, 1508. DOI.
 Pontin, D.I., Craig, I.J.D.: 2005, Current singularities at finitely compressible threedimensional magnetic null points. Phys. Plasmas 12, 072112. DOI.
 Pontin, D.I., Galsgaard, K.: 2007, Current amplification and magnetic reconnection at a 3D null point. Physical characteristics. J. Geophys. Res. 112, A03103. DOI.
 Pontin, D.I., Bhattacharjee, A., Galsgaard, K.: 2007, Current sheet formation and nonideal behaviour at threedimensional magnetic null points. Phys. Plasmas 14, 052106. DOI.
 Pontin, D.I., Hornig, G., Priest, E.R.: 2004, Kinematic reconnection at a magnetic null point: Spinealigned current. Geophys. Astrophys. Fluid Dynamics 98, 407. DOI.
 Pontin, D.I., Hornig, G., Priest, E.R.: 2005, Kinematic reconnection at a magnetic null point: Fanaligned current. Geophys. Astrophys. Fluid Dynamics 99, 77. DOI.
 Pontin, D.I., Priest, E.R., Galsgaard, K.: 2013, On the Nature of Reconnection at a Solar Coronal Null Point above a Separatrix Dome. Astrophys. J. 774, 154. DOI.
 Priest, E.R., D\a’emoulin, P.: 1995, Threedimensional magnetic reconnection without null points. 1. Basic theory of magnetic flipping. J. Geophys. Res. 100, 23443.
 Priest, E.R., Hornig, G., Pontin, D.I.: 2003, On the nature of threedimensional magnetic reconnection. J. Geophys. Res. 108(A7), 1285. DOI.
 Régnier, S., Parnell, C.E., Haynes, A.L.: 2008, A new view of quietSun topology from Hinode/SOT. Astron. Astrophys. 484, L47. DOI.
 Savcheva, A., Pariat, E., McKillop, S., McCauley, P., Hanson, E., Su, Y., Werner, E., DeLuca, E.E.: 2015, The Relation between Solar Eruption Topologies and Observed Flare Features. I. Flare Ribbons. Astrophys. J. 810, 96. DOI.
 Savcheva, A., Pariat, E., McKillop, S., McCauley, P., Hanson, E., Su, Y., DeLuca, E.E.: 2016, The Relation between Solar Eruption Topologies and Observed Flare Features. II. Dynamical Evolution. Astrophys. J. 817, 43. DOI.
 Schindler, K., Hesse, M., Birn, J.: 1988, General magnetic reconnection, parallel electric fields, and helicity. J. Geophys. Res. 93(A6), 5547.
 Stanier, A., Browning, P., Dalla, S.: 2012, Solar particle acceleration at reconnecting 3D null points. Astron. Astrophys. 542, A47. DOI.
 Sun, X., Hoeksema, J.T., Liu, Y., Aulanier, G., Su, Y., Hannah, I.G., Hock, R.A.: 2013, Hot Spine Loops and the Nature of a Latephase Solar Flare. Astrophys. J. 778, 139. DOI.
 Titov, V.S.: 2007, Generalized squashing factors for covariant description of magnetic connectivity in the solar corona. Astrophys. J. 660, 863. DOI.
 Titov, V.S., Galsgaard, K., Neukirch, T.: 2003, Magnetic pinching of hyperbolic flux tubes. i. basic estimations. Astrophys. J. 582, 1172. DOI.
 Titov, V.S., Hornig, G., D\a’emoulin, P.: 2002, The theory of magnetic connectivity in the corona. J. Geophys. Res. 107, 1164. DOI.
 Vekstein, G.E., Browning, P.K.: 1997, Electricdrift generated trajectories and particle acceleration in collisionless magnetic reconnection. Phys. Plasmas 4(6), 2261.
 Wang, H., Liu, C.: 2012, Circular Ribbon Flares and Homologous Jets. Astrophys. J. 760, 101. DOI.
 Wyper, P.F., Pontin, D.I.: 2014, Nonlinear tearing of 3D null point current sheets. Phys. Plasmas 21(8), 082114. DOI.
 Yang, K., Guo, Y., Ding, M.D.: 2015, On the 2012 October 23 Circular Ribbon Flare: Emission Features and Magnetic Topology. Astrophys. J. 806, 171. DOI.
 Zhang, Q.M., Chen, P.F., Guo, Y., Fang, C., Ding, M.D.: 2012, Two types of magnetic reconnection in coronal bright points and the corresponding magnetic configuration. Astrophys. J. 746, 19. DOI.