The Creation of Twist by Reconnection of Flux Tubes

A fundamental process in a plasma is the magnetic reconnection of one pair of flux tubes (such as solar coronal loops) to produce a new pair. During this process magnetic helicity is conserved, but mutual helicity can be transformed to self-helicity, so that the new tubes acquire twist. However, until recently, when Wright (Astrophys. J.878, 102, 2019) supplied a solution, the partition of self-helicity between the two tubes was an outstanding puzzle. Here we examine Wright’s result in detail and apply it to a variety of cases. The simplest case, which Wright himself used to illustrate the result, is that of thin ribbons or flux sheets. We first explicitly apply his method to the usually expected standard case (when the tubes approach one another without twisting before reconnection) and confirm his result is valid for flux sheaths and tubes as well as sheets. For the reconnection of sheets, it is shown that the orientation of the sheets needs to be chosen carefully. For flux sheaths and tubes, Wright’s results are demonstrated to hold for the standard case. There is both a local and a global aspect to the effect of reconnection. The local effect of reconnection is to produce an equipartition of the added self-helicity (and therefore of twist), but the extra global effect of the location and orientation of the feet of the sheet, shell or tube in general adds different amounts of magnetic helicity to the two structures. It is important, as Wright realized, to account for any twist or writhe already existing in the fluxes prior to reconnection. Here we show explicitly that, if a section of a flux sheet is twisted by a multiple of π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\pi $\end{document} with its ends held fixed and is then reconnected with another sheet, then the effect of the reconnection is to add that multiple of π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\pi $\end{document} to one sheet and subtract it from the other, while conserving the total helicity. If, on the other hand, the central part of a flux sheath is twisted before reconnection by any angle, then the effect of reconnection is to add that amount of twist to one sheath and subtract it from the other, while conserving the total helicity. Thus, for the local part of the process in both sheets and sheaths, there is no longer helicity equipartition. Finally, we apply Wright’s results explicitly to flux tubes, and in particular to determine the twist that is acquired by an erupting flux rope due to reconnection during an eruptive solar flare or coronal mass ejection.


Introduction
The reconnection of neighbouring magnetic flux tubes in the solar atmosphere is a very common event. It is central on a large scale to the creation of twist in eruptive flux tubes, which show up as erupting filaments in eruptive solar flares and coronal mass ejections. It also has been proposed to take place during so-called zipper reconnection as an explanation for the creation of flare ribbons (Priest and Longcope, 2017). Furthermore, it is highly likely to be at the core of coronal heating events, both in X-ray bright points (Priest, Parnell, and Martin, 1994;Parnell, Priest, and Titov, 1994;Parnell and Priest, 1995) and also in nanoflares due to multiple current sheets produced by braiding of a uniform field (Parker, 1972;Galsgaard and Nordlund, 1996) or of a field produced by the magnetic carpet (Priest, Heyvaerts, and Title, 2002) or by relaxation of a braided field (Wilmot-Smith, Hornig, and Pontin, 2009). Indeed, Close et al. (2004) used quiet-Sun magnetograms from the Michelson Doppler Imager (MDI) instrument on the Solar and Heliospheric Observatory (SOHO) to track the motion of individual magnetic fragments and deduced that the time for all the field lines in the quiet Sun to change their connections is only 1.4 hours.
Magnetic helicity is a topological invariant (Woltjer, 1958;Moffatt, 1969;Arnol'd, 1974), which is crucial to three-dimensional reconnection, since it is conserved to a very high degree during the reconnection process and can be used to determine the lowest-energy equilibrium magnetic field configuration (Taylor, 1974). Fundamental advances in understanding magnetic helicity were made by Berger and Field (1984), who introduced a gaugeindependent measure of magnetic helicity in a volume not bounded by a magnetic surface. They also realized that magnetic helicity arises partly from the internal structure within a flux tube (such as twist and kinking) and partly from external relations between tubes (i.e. linking and knotting), and that reconnection can transfer internal helicity into external helicity and vice versa. These two types of magnetic helicity were given the names mutual helicity (H m ) and self-helicity (H s ) by Wright and Berger (1989), who also applied the ideas to flux transfer events, showing how reconnection can generate twist in a flux tube.
Self-helicity therefore refers to the twisting and kinking of a flux tube, while mutual helicity represents the linkage between different flux tubes. In turn, the self-helicity of an isolated flux tube of magnetic flux F , say, consists of the sum of its writhe helicity and twist helicity (Berger, 1999;Berger and Prior, 2006), namely, where the writhe (W) measures the winding and kinking of the axis of the tube, while the twist ( ) represents the mean twist of field lines in the flux tube about that axis. W is measured in units such that it equals unity for one complete writhe, while = 2π for one complete twist. Whereas the self-helicity of a flux tube is a topological invariant under ideal motions, neither the writhe nor the twist is invariant. Thus, a tube with high twist but very little writhe may be transformed by an ideal motion into a tube with small twist and large writhe, as indicated in Figure 1. Suppose a plane surface is threaded by N thin flux tubes of magnetic flux F i and that each footpoint possesses a uniform rotation at rate ω i together with a translation. Then a useful expression for the rate at which the magnetic helicity changes is where dθ ij /dt represents the rate of change of the angle between the footpoints i and j (Berger and Field, 1984). The first term represents the rate at which twist is injected into each footpoint (the self-helicity) and the second term represents the braiding of the thin flux tubes (their mutual helicity). Woltjer (1958) suggested the lowest-energy state that conserves magnetic helicity is a linear force-free field, and its role has been studied in laboratory plasmas (Taylor, 1974) and dynamo theory (Moffatt, 1978). Later, Heyvaerts and Priest (1984) suggested the importance of magnetic helicity in the solar corona, both in coronal heating and coronal mass ejections, a role that has been further developed by others, such as Pevtsov, Canfield, and Metcalf (1995), Pevtsov (1998), andPevtsov et al. (2008). An important development was the proposal of a gauge-invariant relative magnetic helicity where B 0 = ∇ × A 0 is a reference potential field with the same normal component on the boundary as the magnetic field B (Berger and Field, 1984), A is the vector potential. Our aim here is to apply Wright (2019)'s recent innovative suggestion for calculating the changes in magnetic helicity in two initially untwisted coronal loops when they reconnect.
First, we set up the model and calculate the effect of magnetic helicity conservation on the twists that are produced in the reconnected loops as mutual helicity is converted into self-helicity (Section 2).
Then we introduce Wright's proposal for a second relation between the twists in the reconnected loops (Section 3).
Next, we consider the reconnection of coronal loops modelled as magnetic flux sheets (Section 4) or flux shells (Section 5) or flux tubes (Section 6). By explicit calculation, we confirm that Wright's technique is valid in each case and gives the simplest and most likely partition of helicities between the reconnected tubes, but that other partitions are possible when pre-reconnection twist or writhe is taken into account.
Finally, we use Wright's technique to estimate the twist produced in an eruptive flare (Section 7). connections of an initial pair of untwisted flux tubes AB and CD, viewed from above and lying above a polarity inversion line (PIL). (b) The effect of reconnection is to create an overlying twisted flux rope AD (denoted by R) together with underlying flux rope BC, denoted by U. (c) and (d) show perspective views, while (e) and (f) show the same situation as (c) and (d) for sheets that lie in a direction parallel to AB.

Effect of Magnetic Helicity Conservation on the Reconnection of Two Coronal Loops
Consider a pair of flux tubes side by side, straddling a polarity inversion line indicated by a dotted line in the Figure 2. The footpoints of one of the flux tubes are located at points A and B, while those of the other tube are at C and D. A perspective view is shown in Figure 2c. We suppose that reconnection creates a long twisted flux rope denoted by R and joining A to D, together with an underlying flux tube denoted by U and joining B to C, as shown in Figures 2b and d. The distance BD is denoted by L, while the perpendicular distance between AC and BD (namely, AE) is w and the shear distance EB is s (see Figure 2a). We assume that initially both flux tubes have flux F and are untwisted, so that they possess zero self-helicity. The initial mutual helicity (H i m ) is, however, non-zero when there is shear -i.e. when s = 0. It has the form where the angles θ 1 and θ 2 are indicated in Figure 2a.
where the angle θ 3 is indicated in Figure 2b The final self-helicity has contributions from the overlying and underlying flux ropes and may be written where R is the twist in the overlying flux rope, while U is the twist in the underlying flux. We need two conditions to determine the twists R and U produced by reconnection, the first of which is magnetic helicity conservation, so that the initial mutual helicity will equal the sum of the final mutual and self-helicities, namely, For the second condition, Priest, Longcope, and Janvier (2016) and Priest and Longcope (2017) suggested helicity equipartition as a hypothesis when applying the theory to the creation of twist in an eruptive flare. This hypothesis was motivated by considering the region in which the reconnection occurs. In the local neighbourhood of the reconnection region, there is no obvious difference between the two reconnection by-products, so it seemed reasonable to assume they each receive the same self-helicity. The flaw in this reasoning is that helicity is not a local quantity, so even a local process like reconnection can produce unbalanced results. Indeed, Wright (2019) derived an alternative condition instead, which for our configuration may be written as proved in the next section.

The Twist of Sheets with Different Orientation
The method of Wright (2019) for calculating the self-helicity of the reconnected fluxes is as follows: (i) start with equal fluxes that are untwisted and have no writhe (Figure 2c), and suppose that they lean over to reconnect, as shown in Figure 2d; (ii) consider the local region around the reconnection site, and allow a contribution to the self-helicity equivalent to a half twist for each of the reconnected fluxes; (iii) the total self-helicity of a given flux can be determined by identifying the footpoint rotations that are required to leave the flux untwisted and with no writhe via the self-helicity removed by these footpoint rotations; (iv) if there is twist or writhe in the fluxes prior to reconnection, this must be accounted for in addition to the calculation in steps (i)-(iii). Wright (2019) demonstrates the method by applying it explicitly to flux in the form of sheets like those in Figure 2e that are aligned parallel to AB and CD. Suppose that the sheets lean over to reconnect, as shown in Figure 2f. Then the problem is to determine the twist and therefore the self-helicity of the flux rope AD in Figure 2f. Here we first consider flux sheets and reproduce Wright (2019)'s analysis for sheets of different orientations, before discussing more details of flux sheet reconnection (Section 4) and extending it to sheaths (Section 5) and tubes (Section 6).

Sheets Oriented Parallel to AB
First of all, suppose the feet of the flux tube AD are rotated to lie along lines perpendicular to AD, as shown in Figure 3a, in which case its twist would simply be π . However, in order to transform from the configuration of Figure 2f into that of Figure 3a, we would need to rotate footpoint A through an angle of −π/2 + θ A and footpoint D through an angle of π/2 + θ D , as indicated in Figure 3a. In this figure the angles θ A = θ D = θ 3 − θ 1 , so that the extra twist that is added in going from Figure 3a to Figure 2f is θ A + θ D = 2(θ 3 − θ 1 ). Thus, in total, the required twist in the flux rope R stretching between A and D in Figures 2d and 2f is and its appearance is sketched in Figure 3b. A similar argument gives the twist in the underlying flux tube U as which may be proved as follows. We need the feet of the flux tube to be oriented parallel to AB, but suppose to start with that they were rotated to lie along lines perpendicular to BC (Figure 4a), in which case its twist would simply be 0. Now, transform this into the configuration of Figure 2f by rotating both footpoints B and C clockwise through an angle of (π/2 − θ B ). However, from the geometry of the parallelogram ABDC, it can be seen that and so Thus, in total, the twist in the underlying flux tube U stretching between B and C is 2(θ 3 + θ 2 ) − π , as required. One consequence of Equations 10 and 11 for the twists is that magnetic helicity is conserved during the reconnection process, since Equation 8 is indeed satisfied. Another is that R lies between π and 2π , while U lies between 0 and π , with R + U lying between 0 and 2π . The flux rope BC created by reconnection. Assume to start with that the feet B and C are oriented in a direction normal to BC, so that the twist vanishes. Then rotate both footpoints B and C through angles (π/2 − θ B ) as shown, to align them parallel to AB. (b) The flux rope AD created by reconnection. Assume the feet A and D are oriented in a direction normal to AD, and then rotate both footpoints through angles (θ 3 − θ 1 ) to align them perpendicular to AB.

Sheets Oriented in a Direction Perpendicular to AB
Is the method of Wright robust when applied to sheets? Does the orientation of the sheets matter? To address this question we consider alternative scenarios where the sheets are oriented in different directions. We start with the flux tubes in Figures 2c and d and suppose they are flattened out to form sheets whose feet lie perpendicular to AB, as shown in Fig In this case, when the rope AD has its feet rotated to be normal to AD (Figure 4b), the twist is π . In order to recover the tube AD with their original positions, we have to rotate both A and D in the same direction by (θ 3 − θ 1 ), and so the twist becomes R = π + 2(θ 3 − θ 1 ), as before.
The reason for this result being unchanged is that, in going from a tube that is flattened in the AB direction to one flattened in a direction perpendicular to AB, you rotate A by π/2 in Figure 6 An untwisted magnetic flux tube whose axis lies in a vertical plane and which consists of four fibrils joining the footpoints that are shown. a clockwise direction and B by π/2 in a counterclockwise direction, and so do not change the net twist.

Sheets Oriented in Other Directions
So what happens when the flux sheets are oriented in other directions? Consider an untwisted flux tube with no writhe whose axis lies in a vertical plane. Suppose it consists of four fibrils (each of flux F f ), as shown in Figure 6 joining footpoints A 1 to B 1 , A 2 to B 2 , A 3 to B 3 , and A 4 to B 4 . Consider what happens when it is flattened in the direction A 2 A 4 and B 2 B 4 to produce a sheet perpendicular the vertical plane in which the axis of the tube lies. Then the change in magnetic helicity produced by motions of A 1 , A 3 , B 1 , and B 3 towards the centres O A and O B (see Figure 7) are given from Equation 2 by where θ i,j represents the changes in orientation of each of the 8 fibril footpoints with respect to all the others. However, if the separation L of the two feet of the flux tube is much greater than the radius of the flux tube, then the only important contributions arise from the footpoints A 1 , A 2 , A 3 , and A 4 with respect to each other and from B 1 , B 2 , B 3 , and B 4 with respect to each other. Thus, it can be seen from Figure 7a that and so the flattening of the tube in this way does not change the magnetic helicity. Next suppose instead ( Figure 7b) that the flattening occurs in the direction A 1 A 3 to produce a sheet lying in the same plane as the axis of the flux tube ( Figure 6). Then the change in magnetic helicity produced by motions of A 2 , A 4 , B 2 , and B 4 to the centres O A and O B is given from Equation 2 by and so the flattening of the tube in this way again does not change the magnetic helicity. Finally, suppose the flattening is at an arbitrary angle θ to the first case, and consider four fibrils with feet at A 5 , A 6 , A 7 , and A 8 that map to points B 5 , B 6 , B 7 , and B 8 , as shown in Figure 7c. Then the change in helicity is given by where α B 5 B 8 , α B 6 B 7 , α B 7 B 8 , α B 5 B 6 are the angles through which B 5 B 8 , B 6 B 7 , B 7 B 8 , and B 5 B 6 , respectively, rotate. Thus, which in general is non-zero and so flattening in an arbitrary direction adds magnetic helicity to the tube. Note that this amount vanishes when θ = 0 and it is a multiple of 2π when θ = π/2. We therefore conclude that flattening in the way we have suggested can add self-helicity prior to reconnection which would need to be taken into account in addition to the self-helicity associated with reconnection. The exception is flattening parallel or perpendicular to AB when no self-helicity is added. This section has been concerned with describing the self-helicity of a sheet with general orientations of its two feet. We have seen that, when the orientations of a sheet at both ends are θ relative to the line joining the feet, and θ is not equal to 0 or π/2, then the sheet is twisted, and its self-helicity needs to be accounted for. On the other hand, Wright (private correspondence) has noticed that, if the orientations of the two feet are instead θ and −θ , then the sheet is untwisted and so it has no self-helicity.

Reconnection of Flux Sheets
We have seen in the previous section that the way in which magnetic helicity is distributed between two sheets by reconnection depends on two factors, namely, the reconnection process itself and the relative location of the footpoints. We regard the former as a local process and the latter as a global process. In this section, let us isolate the local reconnection effect, by neglecting the effects of the locations of the feet and the orientations of the flux sheet that we considered in Section 3 -i.e. supposing the sheets are of unlimited length, so that there is no reflection at the feet of the torsional Alfvén waves that communicate the effect of the reconnection away from the reconnection site. We shall find that the local process produces an equipartition of helicity by distributing the helicity that is released equally between the two flux sheets. The simplest way to see this is to suppose two sheets are initially crossing and then reconnect.

Local Effect of Reconnection -Equipartition of Helicity
Consider two untwisted flux sheets, each of flux F , intersecting at an arbitrary angle and having their footpoints at arbitrary locations ( Figure 8a). Suppose the footpoints are flattened in directions normal to the field lines and are passive in the sense of having no effect on the helicity.
Suppose each of the flux sheets consists of N fibrils of flux F/N (as shown in Figure 8a for N = 4). Then it may be proved that, in the limit as N → ∞, there is equipartition of selfhelicity, so that each tube acquires a twist of π . The proof of this is given in the appendix of Wright (2019) and we summarise his argument as follows.
The initial mutual helicity for such a tube of flux F crossing another tube of flux F is simply Alternatively, this may be shown by summing the effect of the multiple fibrils. Thus, the mutual helicity of fibril a 1 with b 1 , both having flux F/N, is (F /N ) 2 . But a 1 crosses N other fibrils and so it contributes N(F /N) 2 to the total mutual helicity. The same is true for each of the N other fibrils a 1 , a 2 , . . . , a N , and so in total they all give a mutual helicity of which is the same as before. Now suppose each of a 1 b 1 , a 2 b 2 , . . . , a N b N reconnects, as shown in Figure 8b for the case of N = 4, so that the upper and lower parts of the initial tubes are no longer linked. Consider the upper flux bundle. Its initial magnetic helicity from the a 1 b 1 crossover is (F /N ) 2 , and so after reconnection assume a fraction α 1 , say, of this is given as self-helicity to the upper flux rope and 1 − α 1 is given to the lower flux rope. The upper flux rope consists of 1 2 (N 2 − N) crossings, and so its mutual helicity will be 1 2 (N 2 − N )(F /N ) 2 , making a total helicity of In the limit as N → ∞, this reduces to as required. Note that this is just a local argument for the twist added by reconnection, since it does not take account of the orientations of the feet, which also play a role, as in Figure 2. The reason that Figure 2 gives a different answer is that the flux tube may possess extra writhe or twist (and therefore helicity) between the feet and the reconnection location. In other words, the resulting twist depends both on the local effect of the reconnection and the global geometry of the feet in accordance with Section 2.

Local Reconnection of Flux Sheets Consisting of Two Fibrils
Consider two untwisted sheets each of flux F and each consisting of just two fibrils of flux 1 2 F intersecting at right angles. Is the situation just considered in Figure 8 the only way to reconnect the two sheets? This question has been partially addressed by Wright and Berger (1989). Their Figure 2 considers configurations similar our Figure 9a-d, since the focus of their study was on reconnecting the fluxes to produce completely unlinked fluxes. For context, we reproduce their analysis here before extending it to the situation in Figure 9e and f.
In the initial state (Figure 9a), the mutual magnetic helicity may be regarded as that of two sheet crossings of flux F (namely, F 2 ) or as four fibril crossings of flux 1 2 F , namely, First, suppose that a 1 reconnects with b 1 and a 2 with b 2 , as indicated by the large dots in Figure 9a, to give the situation in Figure 9b. When each of the crossing fibrils of mutual helicity F 2 /4 reconnects, the two new fibrils each gain self-helicity F 2 /8, so that the fibrils in the resulting upper flux rope have self-helicity F 2 /4 and mutual helicity F 2 /4, making a total of 1 2 F 2 , while the same is true of the lower flux rope. Next, suppose the two upper fibrils in Figure 9b reconnect to give Figure 9c, in which the self-helicity of the fibrils has been increased by F 2 /4 to 1 2 F 2 , while their mutual helicity has decreased to zero. Then, in Figure 9d, the same thing happens to the two lower fibrils, increasing their self-helicity also to 1 2 F 2 . Finally, suppose the upper and lower crossings (Figure 9e) are reconnected to give Figure 9f. In this case, the upper and lower fibrils are freed, but the two middle fibrils are still interlinked, so we have not succeeded in reconnecting the whole flux sheet. We conclude that there are several ways to reconnect the initial sheets (Figure 9a) to give the configurations shown in Figures 9b, c and d, but in each case, although the internal arrangement of fibrils is different, there is local self-helicity equipartition with the total self-helicities in the reconnected sheets both being the same ( 1 2 F 2 ). We may regard this as the standard and natural situation.
However, it is possible to obtain an infinite number of different non-standard situations without local self-helicity equipartition, but differing by discrete values of self-helicity. Wright (2019) showed how flipping a flux sheet over at the reconnection site prior to reconnection will add ±1/2 turn of twist to the sections of sheet adjacent to the reconnection site. It is important to recognise that the helicity associated with this twist must be added to the 1/2 turns that result from reconnection. The effect on the flux sheets, therefore, is to leave one reconnected sheet untwisted, whilst the other has a twist of one full turn. Clearly preconditioning the reconnected fluxes affects their final self-helicities.

Figure 10
The same situation as in Figure 9, except that the central part of b 1 b 2 is rotated before reconnection, with the result that a 1 now connects with b 2 instead of b 1 .
We now demonstrate this effect for two fibrils. Suppose that, in place of Figure 9a, the central part of the underlying tube is rotated through half a turn, say, while keeping the footpoints fixed, so that no helicity is added to the tube as a whole -i.e. half a turn is subtracted from the section of the underlying tube that is in the upper right of Figure 10a, while half a turn is added to the section in the lower left. After reconnection, the result is that now fibril a 1 is joined to b 2 instead of b 1 , while a 2 is joined to b 1 . The result is that the upper flux sheet is untwisted and has zero helicity, while the lower sheet is twisted by one full turn (2π ) and so its helicity is F 2 . We still have helicity conservation but no longer equipartition. Thus, as well as Wright's standard result, a range of other scenarios is possible, in which the central portion of a flux sheet is twisted by multiples of π before reconnection.

Local Reconnection of Flux Sheets Consisting of Three Fibrils
Next consider the case of sheets consisting of three fibrils (Figure 11a), each fibril having no twist and containing a flux 1 3 F . The initial mutual helicity is Suppose first that there are reconnections at the crossings marked by dots in Figure 11a to give the configuration shown in Figure 11b. A reconnection at each crossing converts mutual to self-helicity in the same way as in the previous subsection and the whole process is linear, so the total helicity is conserved as before.
But what happens if different crossings are reconnected, such as in Figure 12c. Here A 1 B 1 reconnects with C 2 D 2 instead of C 1 D 1 , while A 2 B 2 reconnects with C 3 D 3 instead of C 2 D 2 , and A 3 B 3 reconnects with C 1 D 1 instead of C 3 D 3 . However, as can be seen, the problem is that now the two new flux bundles are not separated, since A 1 C 2 and A 2 C 3 still link with D 1 B 3 . Of course, we could then allow extra reconnections at the joins of D 1 B 3 with the other two fibrils, but the result is to give two flux ropes with exactly the same total helicity as before, even though the distribution between mutual and self-helicity may be different.
Is there any way to obtain a different amount of helicity in the resulting flux sheets? Yes indeed. Just as in the case described in Figure 10a, pairs of fibrils may have their central parts rotated before reconnection without changing the initial helicity, with the result that helicity equipartition during reconnection no longer applies, since the twists of the resulting reconnected sheets can differ by multiples of a turn.

Reconnection of Flux Shells
Do our results change when we model the flux tubes as shells rather than sheets? The way to address this is to regard the shells as consisting of a discrete number of fibrils (just as in our Figure 11 Two of the different ways in which a pair of flux sheets intersecting at right angles and each consisting of 3 fibrils may reconnect. (a) The locations of the reconnections of fibrils A 1 B 1 with C 1 D 1 , A 2 B 2 with C 2 D 2 , and A 3 B 3 with C 3 D 3 , and (b) the result of the reconnections. (c) The locations of the reconnections of fibrils A 1 B 1 with C 2 D 2 , A 2 B 2 with C 3 D 3 , and A 3 B 3 with C 1 D 1 , and (d) the result of the reconnections. The effect of reconnecting A 1 B 1 to C 3 D 3 , A 2 B 2 to C 4 D 4 , A 3 B 3 to C 1 D 1 , and A 4 B 4 to C 2 D 2 . (c) The resulting two flux tubes. discussion of sheets) and then to carefully follow the effects of reconnecting each fibril. We first of all consider the local reconnection process and discover some differences between sheets and shells. Then we discuss how to recover Wright's result when regarding the structures as flux shells rather than sheets. Although Wright's result is the standard, most natural solution, it was derived for fluxes that have no twisting or writhing prior to reconnection. Naturally, there are, therefore, other scenarios that can produce different distributions of magnetic helicity in the reconnected sheaths by preconditioning the fluxes to have portions of non-zero twist. The self-helicity associated with sections of flux that have twist or writhe prior to reconnection needs to be accounted for in addition to the self-helicity associated with reconnection and footpoint location (see Wright, 2019).

Local Effect of Reconnecting Flux Shells
Consider two intersecting flux tubes, each consisting of, for example, four fibrils that are located in a flux shell, as shown in Figure 12a. The first shell has fibrils A 1 B 1 , A 2 B 2 , A 3 B 3 , and A 4 B 4 , while the second shell has fibrils C 1 D 1 , C 2 D 2 , C 3 D 3 , and C 4 D 4 . Now reconnect A 1 B 1 with C 3 D 3 , A 2 B 2 with C 4 D 4 , A 3 B 3 with C 1 D 1 and A 4 B 4 with C 2 D 2 , in the manner shown in Figure 12b. This produces two flux tubes modelled as flux shells, each of twist π , as shown in Figure 12c. In the first tube, A 1 joins to C 3 , A 2 to C 4 , A 3 to C 1 , and A 4 to C 2 , while in the second tube B 1 joins to D 3 , B 2 to D 4 , B 3 to D 1 , and B 4 to D 2 . We therefore have equipartition of self-helicity, as in the case of flux sheets, so let us next include also the global effect of the orientation of the feet.

Deriving Wright's Result by Modelling the Flux Tubes as Flux Shells instead of Sheets
Consider a pair of flux tubes AB and CD that each consist of four fibrils located in a shell. The first possesses fibrils A 1 B 1 , A 2 B 2 , A 3 B 3 , and A 4 B 4 , while the second has fibrils C 1 D 1 , C 2 D 2 , C 3 D 3 , and C 4 D 4 . Suppose they reconnect to give two flux tubes AD and BC (Figure 13a), in which footpoint A 1 joins to D 3 , A 2 to D 4 , A 3 to D 1 and A 4 to D 2 , while B 1 joins to C 3 , B 2 to C 4 , B 3 to C 1 , and B 4 to C 2 . Then, at first sight, it looks a simple matter to calculate the helicity, since the twist in AD appears to be π , but the problem is that it is not immediately obvious what is the effect of the orientation of the feet. Furthermore, it should be borne in mind that, the writhe and twist are not conserved, since one can be converted to another by an ideal motion while their total, namely, the magnetic self-helicity, remains the same. This means that the extra effect of the footpoint orientation may show up in a given example either as writhe or as extra twist or a combination of both.
One way to calculate the helicity is as follows. First of all, imagine the tube AD has its axis located in a vertical plane through AB, so that the centre of the footpoint D is located at D 0 on the line AB instead of at D (as shown in Figure 14a or, viewed from above in Figure 14b). Since this tube has no writhe, its helicity will come entirely from its twist, namely, in this case π . Next, suppose AD 0 is rigidly rotated through an angle θ 3 − θ 1 so that the footpoint D 0 moves to D. This means that the orientations of the fibril footpoints A 1 , A 2 , A 3 , A 4 and D 1 , D 2 , D 3 , D 4 have all rotated by an angle θ 3 − θ 1 , as shown in Figure 14c. Thus, in order to bring the fibril footpoints to their correct locations as seen in Figure 13a, all of the fibril footpoints need to rotate by (θ 3 − θ 1 ) when viewed from above, to give Figure 14d.
However, when the flux tube is straightened out and viewed in perspective (Figure 14e), rather than being curved and viewed from above (Figures 14b-d), it can be seen that the two ends of the flux tube have been rotated in opposite directions when going from Figure 14c to Figure 14d. This provides an extra twist of 2(θ 3 − θ 1 ) to the ends of the tube in addition to the twist of π provided by the local effect of reconnection, and it may show up either as extra flux rope twist or as writhe. Thus, in total, the self-helicity of the reconnected flux rope AD is and, if there is no writhe, its twist is therefore in agreement with Wright's result (Equation 9).
Another way to calculate the extra helicity above the local reconnection value ( 1 2 F 2 ) due to the footpoint orientations is by supposing this extra helicity is all in the form of writhe. Suppose, as in Figure 13b, that the axis of the flux tube starts out near A pointing upwards Figure 15 (a) The writhed axis AMD of the flux tube AD, whose projection onto the xy-plane makes a curve that starts out at A inclined an angle θ 3 − θ 1 to the x-axis, passes through P and then ends up at D also inclined at θ 3 − θ 1 . (b) The tantrix sphere containing points that map from the ends of the unit vectors that are parallel to the flux tube axis. and in a vertical plane through AB. Then it rotates to be horizontal over the midpoint between A and D, and, during the second half of its path, it ends up approaching D along a curve lying in a vertical plane through D.
Writhe can be calculated in terms of the so-called "tantrix curve" (Berger and Prior, 2006). The unit vector that is tangent to the axis of the writhing flux tube AD maps points on the axis of the tube to points on a unit sphere, called the tantrix sphere. The tantrix curve is then the path that the tip of the tangent vector traces out on the surface of the sphere. The writhe of a closed curve is given in terms of the spherical area A enclosed by the tantrix curve as For a flux tube such as AD, the tantrix curve is completed in two steps: first the tube is extended by adding sections below the xy-plane that change the tangent direction to the vertical and so add sections to the tantrix curve that extend it to the north and south poles along geodesics; secondly, the curve is closed by adding a semicircle. Thus, in Figure 15a the points A, M and D, map to A , M , and D on the tantrix curve ( Figure 15a), with M lying on the curve φ = 0 of longitude. Then the footpoints A and D are continued down to N and S below the xy-plane, which map to N and S on the tantrix. Finally, the curve is completed by a semicircle which maps onto the curve φ = π of longitude on the tantrix. Now the area between the curves φ = 0 and φ = π is 2π and the area we require is that area minus the area between the curves φ = 0 and the curve N A D S (namely, φ = θ 3 − θ 1 ), which gives an area of 4π(θ 3 − θ 1 ). Thus, and so as required.

Generalisations of Wright's Standard Result
We have verified Wright's result for the standard natural reconnection of flux sheaths when there is initial twisting or writhing present. However, it is possible to reconnect the field in other ways, as follows, depending on the direction in which flux sheaths approach one Figure 16 Reconnection of a pair of untwisted flux tubes AB and CD, viewed in perspective (top) or from above (middle) or from the side, in which the first fibrils to reconnect are (a) A 2 B 2 with C 4 D 4 , or (b) A 3 B 3 with C 4 D 4 or (c) A 3 B 3 with C 1 D 1 .
another before reconnecting and also on a possible twisting of part of the sheaths. As mentioned in Wright's (2019) discussion of sheets, any initial twisting or writhing needs to be accounted for in addition to the self-helicities expected from reconnecting fluxes with no twist or writhe. First, consider an arcade of coronal loops, modelled as a pair of untwisted flux sheaths AB and CD, which reconnect to create an overlying twisted flux rope AD together with underlying flux BC ( Figure 16). Suppose the loop AB consists of four fibrils A 1 B 1 , A 2 B 2 , A 3 B 3 , A 4 B 4 , while CD consists of fibrils C 1 D 1 , C 2 D 2 , C 3 D 3 , C 4 D 4 . The simplest and most obvious way for the flux tubes AB and CD to touch and reconnect is in a symmetric manner, which first leads to A 2 B 2 reconnecting with C 4 D 4 , as shown in Figure 16a. However, an alternative is for them to approach asymmetrically, with, for instance, tube CD lying closer to the solar surface than AB and approaching at differing angles. For instance, in Figure 16b the first field lines to reconnect are A 3 B 3 and C 4 D 4 , whereas in Figure 16c the first ones are A 3 B 3 and C 1 D 1 .
Secondly, consider the nature of the local reconnection regardless of the orientation of the feet. We found for flux sheets that rotation of the central part of a sheet before reconnection could lead to a different distribution of twist, and the same is true for flux sheaths, but with an important difference. Consider, for example, the case shown in Figure 17 of two flux sheaths, each consisting of four fibrils, approaching one another and reconnecting. Suppose sheath AB is rotated by π/4, so that A 1 B 1 reconnects with C 2 D 2 instead of C 3 D 3 , while A 2 B 2 reconnects with C 3 D 3 , A 3 B 3 reconnects with C 4 D 4 and A 4 B 4 reconnects with C 1 D 1 . Thus, after reconnection A 1 , A 2 , A 3 , and A 4 connect to C 2 , C 3 , C 4 , and C 1 , respectively, while B 1 , B 2 , B 3 , and B 4 connect to D 2 , D 3 , D 4 , and D 1 , respectively. The result is that Figure 17 Local reconnection of a pair of flux tubes AB and CD in which flux tube AB is first rotated about its centre by π/4 before reconnection, leading to a twist of π/4 in BD and 3π/4 in AC.
instead of both sheaths being twisted by π/2, sheath AC is twisted by 3π/4, while sheath BD is twisted by π/4. This may be generalised to allow the central part of sheath AB to twist by any multiple of π/4 before reconnection, so that after reconnection AC will be twisted by any corresponding multiple of π/4, while BD is twisted by an amount such that magnetic helicity is conserved. A further generalisation is to allow the sheaths to consist of, say, n fibrils, in which case pre-twisting by any multiple of π/n is possible, and so the final twist of AC could be any multiple of π/n. In the limit as n → ∞, rotation by any angle before reconnection is possible and so AC could be twisted by any corresponding amount after reconnection, with again BD twisted by the amount that conserves magnetic helicity. This differs from our previous consideration of flux sheets, when only rotation by discrete multiples of π was allowed.

Reconnection of Flux Tubes
Flux tubes, which consist of a continuum of field lines spread through a volume, rather than along a line or in a shell, are potentially much more complex than flux sheets or flux shells. For example, they may be initially twisted ropes, when Linton, Dahlburg, and Antiochos (2001) and Linton and Priest (2003) showed that one can obtain tunnelling with double reconnection rather than slingshot reconnection. Another complex example is the case of double separator reconnection between two flux tubes that was considered by Priest and Longcope (2017).
However, here we focus on much simpler scenarios, by generalising from the above results for flux sheets and sheaths. So, the question that arises is: do the results for simple two-dimensional sheets or sheaths apply to three-dimensional flux tubes? We would like to suggest that yes they do, on the following grounds.
First of all, let us regard the flux tube as being composed of a series of flux sheaths ( Figure 18). Then reconnection of the flux tube cannot simply be modelled as the sequential reconnection, as in Figure 8, of the set of flux sheaths for the following reason. Suppose we start with two flux sheaths, an inner one and an outer one, as in Figure 18a, with the outer sheaths each consisting of four fibrils, namely, A 1 B 1 , A 2 B 2 , A 3 B 3 , A 4 B 4 in the first tube and C 1 D 1 , C 2 D 2 , C 3 D 3 , A 4 B 4 in the second tube. Then, at first sight, there seems to be no problem in reconnecting these fibrils in the outer sheaths to give the situation shown in Figure 18b. However, a difficulty arises since the fibrils in Figure 18b do not form smooth  sheaths, because, if the number of fibrils is increased, there will eventually be two fibrils between C 1 D 1 and C 4 D 4 in the one sheath and between A 2 B 2 and A 3 B 3 in the other sheath which cannot reconnect, since their paths are blocked by the inner sheaths. That is to say, the first fibril lies above the inner sheath of the upper flux tube and the second fibril lies below the lower flux tube. Reconnection of these two fibrils can, however, still take place after the inner sheaths have reconnected.
As a better alternative, a flux tube may be sliced into a series of flux sheets (Figure 19), so that reconnection of the flux tubes may be modelled as the sequential reconnection, as in Figure 8, of the set of flux sheets. Figure 19 illustrates the case with seven sheets, such that field line A 3 B 3 reconnects with field line C 1 D 1 , and then sheet a 1 reconnects with sheet c 1 , followed by a 2 with c 2 , a 3 with c 3 , a 4 with c 4 , a 5 with c 5 , a 6 with c 6 , and finally field line A 1 B 1 with field line C 3 D 3 . There seems therefore to be no conceptual problem to regarding the reconnection of a three-dimensional flux tube as the reconnection of an infinite sequence of flux sheets.

Figure 20
The variation of (a) R and (b) U with s/w for various values of L/w.

Application of Wright's Method to the Eruption of a Simple Arcade in an Eruptive Solar Flare
In previous publications, we have considered the influence of magnetic helicity conservation on the creation of twist in an eruptive solar flare (Priest, Longcope, and Janvier, 2016;Priest and Longcope, 2017). We included as a hypothesis magnetic helicity equipartition during reconnection, but Wright (2019) has suggested that there will not generally be such equipartition, so here we reconsider the role of helicity conservation and its role in creating twist in eruptive flares and coronal mass ejections by adopting Wright's method.

R and U as functions of s/w and L/w
Following Priest, Longcope, and Janvier (2016), suppose that a preflare sheared arcade is modelled simply as a pair of flux tubes AB and CD with significant mutual helicity, and suppose that the tubes reconnect to give an eruptive flux rope AD together with underlying flux BC (Figure 2). The aim is to calculate the self-helicity in the eruptive flux rope. Wright (2019) has proposed that, for this situation the twist in the eruptive flux tube AD if it has no writhe is given by while the twist in the underlying flux tube U is as proved in Section 3, and so we may calculate its variation with various parameters as follows. R and U are functions of θ 1 , θ 2 and θ 3 (indicated in Figure 2), and therefore of the arcade shear (s/w) and length (L/w) normalised with respect to its width (w). First of all, consider R , where tan θ 1 = w/(L + s) and tan θ 3 = w/s, so that tan(θ 3 − θ 1 ) = tan θ 3 − tan θ 1 1 + tan θ 3 tan θ 1 = w/s − w/(L + s) As s increases from 0 to ∞, the expression 2(θ 3 − θ 1 ) in Equation 29 decreases from 2 arctan(L/w) to 0, while, as L/w decreases from ∞ to 0, 2 arctan(L/w) decreases from π to 0. The resulting behaviour of R with s/w and L/w is plotted in Figure 20a.
Finally, consider for which Thus, as s → 0 or ∞, R + U → 2π , whereas, when L w, R + U → 0. In general, as found by Priest, Longcope, and Janvier (2016), R + U lies between 0 and 2π . The results in Figure 20 can also be compared with some of the properties given in Wright (2019), who found the overlying tube would always have a twist greater than a half turn (i.e., R /π > 1) and the underlying tube would always have a twist less than a half turn (i.e. U /π < 1).
For the specific situation in which the lines AB and CD in Figure 2 are parallel and the fluxes lie adjacent to one another (so are touching and ready to reconnect), we recover the situation in Figure 8b of Wright (2019). In terms of our notation, this corresponds to L/w → ∞, s/w → ∞ and s < L. Figure 20 shows that the overlying and underlying tubes both have a self-helicity of π , equivalent to a twist of a half turn, in agreement with Equation 23 of Wright (2019) when his θ A → 0 and θ C → 0.

Conclusion
A huge puzzle in the basic theory of reconnection has been to understand the nature of the conversion of, say, mutual magnetic helicity into self-helicity during the reconnection of two magnetic flux tubes, which in turn is important in determining the twist that can be added to a flux tube (such as an erupting solar flux rope on a large scale or, say, spicules on a small scale). The main unanswered question is: when two flux tubes reconnect, how is the resulting self-helicity (which can be either twist or writhe) partitioned between the two flux tubes? One indisputable constraint is total magnetic helicity conservation, so that the fall, say, in mutual helicity during reconnection must equal the rise in the sum of the self-helicities of the two flux tubes. But there must be a second condition that determines the distribution of the total self-helicity between the two reconnected flux tubes. Priest, Longcope, and Janvier (2016) and Priest and Longcope (2017), in their consideration of eruptive solar flares, suggested helicity equipartition, so that the same amount of self-helicity is added to the two tubes. More recently, Wright (2019) has proposed a clever technique for determining the second constraint and so deducing the resulting partition of self-helicity. The purpose of this paper has been to investigate and develop his proposal in detail and, in particular, to consider explicit calculations for flux sheets, flux sheaths, and flux tubes.
Our analysis and that of Wright (2019) imply that there is both a local and a global aspect to the effect of reconnection, and so our treatment is split into three parts. The first is to analyse the local effect of reconnection, which is found to produce an equipartition of the self-helicity (and therefore the twist) added to the structures, regardless of whether they are sheets, sheaths or tubes. The second is to consider the additional global effect of the orientation of the feet of the sheet, shell or tube. In general, it is found to add different amounts of magnetic self-helicity to the two structures. The third part is to allow for twisting of the central region of the sheet, sheath or tube before reconnection.
The local effect of reconnection is demonstrated both for sheets, sheaths, and tubes to be one of equipartition. This is achieved by splitting the structure into a series of fibrils and then letting the fibrils reconnect, one after the other.
However, for the global effect, namely, the effect of the orientation of the feet of the structures after reconnection, the extra amounts of self-helicity added to the tubes differ and depend on their orientations. This is because what matters is the orientation of the feet relative to the new orientation of the flux tube after reconnection, which differ from the initial orientation before reconnection. For flux sheets, the orientation before reconnection needs to be chosen carefully. For flux sheaths, Wright's results are also verified, by either regarding the self-helicity solely in terms of twist and considering carefully the locations of the feet of the fibrils or by regarding it as writhe and using the so-called tantrix method.
If the central part of a flux sheet, sheath or tube is twisted before reconnection while the feet are held fixed, then there is no longer local helicity equipartition and the untwisted calculation needs to be modified to include additional contributions. However, sheets behave very differently from sheaths or tubes, since sheets can have their central parts rotated only by a multiple of π , whereas sheaths and tubes can be rotated by arbitrary angles. The result is that after reconnection an extra amount of twist has been added to one flux tube and subtracted from the other.
In the future, it will be interesting to conduct numerical experiments to test the applicability of the result, and in particular to produce much more realistic analysis of the twist acquired by an erupting flux rope during an eruptive solar flare or coronal mass ejection. In addition, the energetics of the configurations needs to be investigated, so as to decide which transitions are favourable and therefore which are likely to lead to reconnection, as discussed in Priest, Longcope, and Janvier (2016) and Priest and Longcope (2017). Furthermore, the results can be applied to build an improved model for zipper reconnection (Priest and Longcope, 2017) in order to explain the early stages of a flare or coronal mass ejection, where reconnection starts at one location and spreads along the polarity inversion line to create flare ribbons and a rising arcade of reconnected coronal loops.