the magnetic field evolution in plasmas such as the solar corona and the earth s magnetosphere often follows an ideal evolution in which all topological features of the field remain unchanged over time .
this invariance provides a motivation to investigate topological properties of the magnetic field such as magnetic null points , flux surfaces and periodic field lines .
some of these topological features , particularly magnetic null points , are surprisingly stable and turn out to be conserved not only under an ideal plasma dynamics but also for a wide range of non - ideal flows .
mathematically this is reflected in the existence of an index theorem ( greene , 1992 ) for magnetic null points which allows only a small number of bifurcation processes to either generate or destroy magnetic null points .
this underlines the importance of magnetic null points and their properties for the structure and dynamics of astrophysical plasmas .
generic magnetic null points ( i.e.  null points where all eigenvalues of the linearisation of the magnetic field at the null have non - zero real part ) are associated with certain distinguished field lines and flux surfaces , namely spines ( or @xmath0 lines in the terminology of lau & finn , 1990 ) and fan surfaces ( or @xmath1 surfaces ) as well as , sometimes , magnetic separators , the intersection lines of two fan surfaces .
combined , these features are sometimes termed the magnetic skeleton ( e.g.  priest & titov , 1996 ) . in dynamical systems terminology ,
null points are hyperbolic fixed points of a volume preserving ( divergence - free ) vector field and the spine and fan of a null are stable or unstable invariant manifolds .
a null whose fan surface corresponds to a stable manifold and the spine line to an unstable manifold is designated as type  a while the converse is known as type  b. the magnetic separator is a heteroclinic intersection of fan surfaces , lying in the fan plane of both a type  a and type  b null as shown in figure  [ fig : initialstate2 ] .
conservation of topology has important implications for the plasma evolution , limiting the evolutions that may occur and the amount of free magnetic energy bound in the configuration .
however , even for plasmas which are largely ideal , the frozen - in condition can be violated in localised regions ( such as current sheets ) where non - ideal effects become important .
the key requirement for such a process to change the topology is to have a non - zero electric field component in the direction parallel to the magnetic field ( hesse & schindler , 1988 ) .
such a component ( for further details see sections  4 and 5 ) enables a change in connectivity of the magnetic field lines and so a change of topology , the process being known as magnetic reconnection .
reconnection may allow for significant magnetic energy release , both locally at the reconnection site and globally in the previously topologically bound energy .
historically , models of magnetic reconnection were solely two - dimensional .
theoretical considerations show that in two - dimensions reconnection can only occur at a hyperbolic ( x - type ) null - point of the field ( for details see e.g.  hornig , 2007 ) .
such null - points , lying at the intersection of four topologically distinct flux domains , are also likely sites of current sheet formation .
the resultant reconnection is now fairly well understood ( for reviews see , e.g.  biskamp 2000 , priest & forbes 2000 ) , although modelling still required to fully determine the importance of the complete physics including kinetic effects ( birn & priest 2007 ) .    in three - dimensions ,
the typical situation in astrophysical plasmas , there are no fundamental restrictions on where reconnection occurs , the only requirement being the presence of a localised non - ideal term in ohm s law ( hesse & schindler , 2007 ; schindler _ et al . _  1988 ) .
accordingly , reconnection may take place wherever such non - ideal terms become important .
three - dimensional null - points are thought to be one such location ( klapper _ et al . _
various forms of reconnection can occur at 3d nulls ( for a review see priest & pontin 2009 ) and the nature of the reconnection has been analysed in some detail ( e.g.  pontin _ et al . _
2004 , 2005 ; galsgaard & pontin 2011 ) .
current sheets may also form away from null points ( e.g.  titov _ et al . _  2002
; browning _ et al . _  2008
; wilmot - smith _ et al . _
models for the local reconnection process here ( e.g.  hornig & priest 2003 , wilmot - smith _ et al . _  2006 , 2009 ) show that , just as in the 3d null case , there are many distinct 3d characteristics of reconnection , some of which we discuss below .
one topologically distinguished location where current sheets may form away from nulls is at separators which , like the 2d null point , lie at the intersection of four flux domains .
it is this characteristic that has been used to argue why separators should be prone to current sheet formation ( lau & finn 1990 ) .
the minimum energy state for certain field configurations has the current lying along separators ( longcope 1996 , 2001 ) and current sheet formation at separators has been observed in numerical simulations ( e.g.  galsgaard & nordlund , 1997 ; haynes _ et al . _  2007 ) .
reconnection taking place at separators is known as separator reconnection .    a recent series of papers ( including galsgaard & parnell 2005 ; haynes _ et al . _  2007
; parnell _ et al . _  2008
; parnell _ et al . _
2010a ) examine the reconnection taking place in a ` fly - by ' experiment where two magnetic flux patches on the photosphere are moved past each other and the resultant evolution in the corona is followed .
the authors find strong current concentrations and therefore reconnection around each of a number of separators ( up to five in some time frames , haynes _
et al . _
these separators are created in the reconnection process and connect the same two nulls ( see also longcope & cowley , 1996 ) .
magnetic flux is found to evolve in a complex manner , being reconnected multiple times through these distinct reconnection sites ( parnell _ et al . _
2008 ) .
separators were also found to be important in a simulation of magnetic - flux emergence into the corona ( parnell _ et al . _
2010b ) where current sheets between the emerging and pre - existing coronal field were shown to be threaded by numerous separators , a feature of complex magnetic fields suggested by albright ( 1999 ) .
et al . _  ( 2010b ) found the regions of highest integrated parallel electric field to be at or in association with the separators , so arguing that reconnection occurs at the separators themselves .
separator reconnection also appears to be a key process at the magnetopause where magnetic nulls or even clusters of nulls appear in the northern and southern polar cusp regions ( dorelli _ et al . _  2007
these nulls are linked by a separator along the dayside magnetopause which will therefore be involved in reconnection with the interplanetary magnetic field .
indeed reconnection in such a configuration has been implied by observations from the cluster spacecraft ( xiao _ et al . _  2007
these studies partially motivate the present work where we want to investigate how and why multiple separators can be involved in reconnection .
although separator configurations have been analysed for some time ( e.g.  chance _ et al . _  1992
; craig _ et al . _
1999 ) , there remain several unanswered basic questions surrounding the nature of separator reconnection , particularly involving the way magnetic flux evolves , i.e.  exactly how the changes in field line connectivity occur .
there have been suggestions that the reconnection involves a ` cut and paste ' of field lines at the separator itself in a manner similar to the 2d null case ( e.g. lau & finn 1990 ; priest & titov 1996 ; longcope _ et al . _
however , previous detailed investigations into the nature of individual 3d reconnection events both at and away from null points have shown the flux evolution to be quite distinct from the 2d picture ( hornig & priest 2003 ; pontin _ et al . _
2004 , 2005 ) and , in the light of these models , such a simple flux evolution would be surprising .
indeed a number of key features of 3d reconnection are now known which show many of its features are quite distinct to the 2d case .
accordingly it would appear useful to re - examine the fundamental nature of separator reconnection and consider whether it has its own particular distinguishing characteristics .
in doing so we also wish to address a conclusion of parnell et al .
( 2010a ) who suggested that separator reconnection does not appear to involve the nulls that lie at both ends of the separator as well as the recent findings ( parnell _ et al .
_  2010b ) of very large numbers of separators sometimes appearing in numerical mhd experiments .
accordingly the aim of the present work is to describe in detail the nature of an isolated 3d reconnection event in the vicinity of a separator .
we do so using a simple analytical model which is described in detail in section  [ ref : themodel ] .
the model allows us to consider typical magnetic field connectivities resulting from separator reconnection in section  [ sec : bif ] and the nature of the magnetic flux evolution in section  [ sec : evolve ] .
we consider how reconnection rates may be determined in separator configurations in section  [ sec : rates ] before discussing our findings and concluding in sections  [ sec : discussion ] & [ sec : conc ] .
a fully self - consistent model for reconnection must incorporate a dynamic evolution which generates current sheet(s ) as well as the reconnection that takes place at those current sheets and changes the magnetic field topology .
an example in the solar corona is the emergence of a magnetic flux tube from the convection zone and reconnection with the pre - existing coronal magnetic field .
inherent in such events is an enormous separation of scales between the global dynamic process and the local reconnection events .
accordingly , a typical approach to model reconnection itself is to start with a local magnetic field configuration that is considered susceptible to current sheet formation ( such as , in two dimensions , an x - type null point of the field ) . in simulations the magnetic field
is then confined to a finite region and the boundaries driven in such a manner as to initiate a reconnection event which can then be studied in detail .
determining physically realistic boundary conditions is just one of the obstacles in this modelling technique .    here
we also take a simplified approach to consider the local reconnection process , aiming to model the effect of reconnection on the field topology in a three - dimensional magnetic separator configuration .
as discussed in section  1 , such configurations are thought to be likely sites for current sheet formation and associated reconnection in the solar atmosphere . to construct our model
we take advantage of a generic feature of reconnection , namely , the presence of an electric field with a component parallel to the magnetic field in a localised region . such an electric field will , from faraday s law , generate a magnetic flux ring around it . the flux ring will grow in strength until the reconnection ceases .
this situation is illustrated in figure  [ fig : addstwist ] and can be further motivated by considering the equation for the evolution of magnetic helicity . expressing the electric field , @xmath2 , as @xmath3 ( where @xmath4 is the vector potential for the magnetic field @xmath5
, @xmath6 is a potential and @xmath7 represents time ) , we have that the helicity density evolves as @xmath8 in general three - dimensional situations the condition for magnetic reconnection to occur is the existence of isolated regions at which @xmath9 . such regions act as source terms for magnetic helicity ( see the right - hand side of equation  [ eq : helicityevolution ] ) thus imparting a localised twist to the configuration at a reconnection site .
we may investigate the effect that reconnection has on a certain magnetic field topology by determining the effect of additional localised twist within the configuration .
adding any new field component in a manner consistent with maxwell s equations provides a basis for realistic field modifications .
although the question of which states could be accessed in a dynamic evolution is outwith the scope of the model , comparison of results with large - scale numerical simulations may help to determine the plausibility of results ( see section  [ sec : discussion ] for a discussion of this point ) .    moving on to the specific details of our model
, we take as a basic ( pre - reconnection ) state a potential magnetic field configuration in which two magnetic null points at @xmath10 are connected by a separator and is given by @xmath11 here @xmath12 determines the location of the two nulls along the @xmath13-axis while @xmath14 and @xmath15 give the characteristic field strength and length scale , respectively . here
we set @xmath16 , @xmath17 , and @xmath18 .
this configuration is exactly that illustrated in figure  [ fig : initialstate2 ] .
the spine of the upper null point ( itself located at @xmath19 @xmath20 ) is the ( one - dimensional ) unstable manifold of the null and lies along the line @xmath21 .
the fan surface of the upper null , the ( two - dimensional ) stable manifold lies in the plane @xmath22 and is bounded below by the spine of the lower null - point .
the spine of this lower null point ( itself located at @xmath19 , @xmath23 ) is the ( one - dimensional ) stable manifold of the null and lies along the line @xmath24 while the fan surface is the ( two - dimensional ) unstable manifold .
it lies in the plane @xmath25 and is bounded above by the spine of the upper null point .
the two null points are connected by a separator at the intersection of the two fan surfaces , @xmath26 , @xmath27 $ ] . throughout
we consider the magnetic field over the spatial domain @xmath28 $ ] and @xmath29 $ ] .    to simulate the topological effect of reconnection in this configuration we add a magnetic flux ring of the form @xmath30 here the flux ring is centred at @xmath31 , the parameter @xmath32 relates to the radius of the ring , @xmath33 to the height and @xmath34 to the field strength . in the primary model considered here we choose to centre the flux ring along the separator ( at @xmath35 ) and take the parameters @xmath36 , @xmath37 and @xmath38 giving the particular flux ring
@xmath39    we add this flux ring to the potential field in a smooth manner , taking a time evolution satisfying faraday s law , @xmath40 specifically we set @xmath41 with @xmath42 so that the time evolution takes place in @xmath43 .
( note that the gradient of a scalar could also be added to the electric field which could allow for the superimposition of a stationary ideal flow .
we have , for simplicity , neglected this possibility . )
similar evolutions can be obtained for the addition of the more general flux ring to the potential field .
the addition of the magnetic flux ring creates a localised region of twist in the centre of the domain and we now wish to determine whether and how the magnetic field topology is changed . in order to do so first notice that the flux ring is sufficiently localised that the magnetic field at the null points remains unchanged during the time evolution .
accordingly for each magnetic null point we may trace magnetic field lines in the fan surface in the neighbourhood of the null point and out into the volume .
this method allows us to determine how the fan surfaces are deformed by the reconnection and to locate any intersections of these surfaces , i.e.  separators .
we describe our findings on the field topology in the following section .
in order to follow the evolution of magnetic field topology during the reconnection we begin by showing in figure  [ fig : fansz0 ] the intersection of the fan surfaces with the @xmath44 plane .
the strength of the flux ring increases linearly in time up to the final state which we have normalised as @xmath45 .
we show the intersections by tracing field lines from each fan surface in the close neighbourhood of the corresponding null which is possible since the disturbance is localised near the centre of the domain and so the eigenvalues associated with the nulls do not change in time . in the images
the fan surface of the lower null is coloured blue and that of the upper null in orange .
we first note that in the initial phase of the process the angle between the fan surfaces decreases .
if the reconnection is weak the process can stop in this phase without leading to any change in the magnetic skeleton .
however , a stronger reconnection event can lead to a further closing of the angle between the fan surfaces until they intersect ( at about @xmath46 in our model ) . recall that crossings of the fan surfaces give the location of magnetic separators in that particular plane .
hence the process creates two new separators and correspondingly two new magnetic flux domains . in order to properly identify the various flux domains we label each of the distinct topological regions with the numbers i  vi , as shown in the lower - right  hand image of figure  [ fig : fansz0 ] .
we shall examine the nature of these flux domains later in this same section , but at this point it is already possible to make some general statements about this type of bifurcation .
obviously the manner in which the fan planes can fold and intersect leads to the process always creating ( or the reverse process annihilating ) separators in pairs .
there is no way in which we can create a single new separator as long as the reconnection is localised .
( this excludes that the whole domain under consideration is non - ideal and that separators enter or leave the domain across the boundary ) .
also shown in figure  [ fig : fansz0 ] are components of the vector field @xmath47 in the @xmath44 plane which is the vector field perpendicular to the central separator at @xmath48 .
this field structure is initially hyperbolic but becomes elliptic as the reconnection continues .
additionally a separator typically has both elliptic and hyperbolic field regions along its length .
these findings coincide with those of parnell _ et al . _
( 2010 ) who examined the local magnetic field structure along separators in a three - dimensional mhd simulation and found both elliptic and hyperbolic perpendicular field components .    a three - dimensional view of the fan surfaces at @xmath49 ( an instant when three separators are present ) is shown in figure  [ fig : fans07 ] . since each separator begins and ends at the null points , the surfaces must fold in a complex way such that the three intersections merge and coincide at the nulls .
this folding takes place as the fan surface of one null approaches the spine of the other null and creates ` pockets ' of magnetic flux running parallel to the spine , as shown in figure  [ fig : fans07 ] where the flux passing through the @xmath50 boundary can be seen .
the new flux domains alone are shown separately in the right - hand image of figure  [ fig : fans07 ] and correspond to the two new flux domains ( types v and vi ) that appear in the first bifurcation as shown in figure  [ fig : fansz0 ] .
we now examine in more detail the connectivity of the magnetic flux in the topologically distinct flux domains i  vi ( figure  [ fig : fansz0 ] ) , focusing on the flux that passes through the central plane , @xmath44 ( flux above the upper null and below the lower null is not of interest here ) .
the flux through the central plane can be distinguished in two ways .
firstly we can follow field lines in the negative direction to determine on which side of the fan plane of the lower null they end ( @xmath51 or @xmath52 ) and secondly we can follow the field lines in positive direction to determine on which side of the fan plane of the upper null they end ( @xmath53 or @xmath54 ) .
combined this technique shows the nature of the field line connectivity in each of the four different flux domains .
flux in the domain labelled i enters the domain through the boundary @xmath55 and leaves through the boundary @xmath56 . flux in domain ii enters through @xmath57 and leaves through @xmath56 , that in domain iii enters through @xmath57 and leaves through @xmath50 while that in domain iv enters through @xmath55 and leaves through @xmath50 .
note that this distinction of fluxes is based on topological features of the null points and is independent of the particular choice of our boundary .
following the first bifurcation at @xmath58 the two new flux domains v and vi are created .
field lines in flux domain v have the same basic connectivity type as those in flux domain iv ( i.e.  they enter through through the @xmath55 boundary and leave through @xmath50 boundary ) .
field lines in flux domain vi have the same basic connectivity type as those in flux domain ii ( i.e.  they enter through the @xmath57 boundary and leave through @xmath56 boundary ) .
although these two basic connectivity types existed prior to the bifurcation , the flux is enclosed within topologically distinct regions as it threads the new ` pockets ' previously described ( see figure  [ fig : fans07 ] ) .
as the reconnection continues , further distortion of the surfaces gives another bifurcation at @xmath59 generating an additional pair of separators , as shown in figure  [ fig : fansend ] .
thus at the end of the reconnection considered here five separators are present in the configuration , each connecting the same pair of nulls .
figure  [ fig : fansend ] also shows a colour coding of the various flux domains according to connectivity type .
this is shown at the end of the reconnection and gives a summary of the way the flux domains relate to each other ( i.e.  field lines in flux domains shown in the same colour have the same connectivity type ) .
we envisage this type of reconnection as occurring , for example , at the dayside magnetopause due to interaction with the interplanetary magnetic field .
accordingly , a crucial question is exactly how the changes in flux connectivity occur and how flux is transferred between topologically distinct domains . in the magnetopause example
this would have implications for how the solar wind plasma can interact with that in the magnetosphere .
we begin to address this question in the following section .
in order to track the magnetic flux in time and determine how the changes in connectivity occur we examine some particular _ magnetic flux velocities _ and so begin by briefly discussing their motivation . recall that an ideal evolution of the magnetic field is one satisfying @xmath60 ( where @xmath61 is the plasma velocity ) the curl of which gives @xmath62 in such a situation the magnetic field and a line element have same evolution equation and so flux is ` frozen - in ' to the plasma and the magnetic topology is conserved .
a real plasma evolution has @xmath63 where @xmath64 represents some non - ideal term ( such as @xmath65 in a resistive mhd evolution where @xmath66 is the resistivity and @xmath67 the electric current ) which is typically localised to some region of space . in this paper
the effect of the non - ideal term is modelled by the addition of a magnetic flux ring .
even with the inclusion of the non - ideal term we may sometimes still find a velocity @xmath68 with respect to which the magnetic flux is frozen - in if can be written as @xmath69 where @xmath70 is an arbitrary function , since taking the curl again and using faraday s law gives an equation of the form .    in three - dimensional magnetic reconnection
we have localised regions where @xmath71 and , as a result , a unique flux velocity @xmath68 can not be found ( hornig & priest 2003 ) .
instead we may consider field lines as they enter and leave the non - ideal ( @xmath71 ) region , fixing them on one of these sides .
we choose a transversal surface that lies below the non - ideal region and integrate along magnetic field lines into and out of the non - ideal region .
parameterising a magnetic field line by @xmath72 and starting the integration from the point @xmath73 ( with @xmath74 ) on the transversal surface we may determine @xmath70 as @xmath75 and subsequently @xmath68 as @xmath76    carrying out the integration for @xmath70 until field lines meet a boundary and leave the domain and noting that @xmath77 at these boundaries we can find a particular flux velocity @xmath78 and visualise the effect of reconnection on the magnetic flux by showing its component perpendicular to the boundaries .
this procedure can only be carried out in the presence of a single magnetic separator ; multiple separators lead to ambiguities in the potential @xmath70 .
note that in the example considered here the electric field which is responsible for the non - ideal evolution decays exponentially away from its centre . at @xmath79 the electric field has fallen to a value of @xmath80 of the maximum in the domain , sufficiently small to be considered an ideal environment and so we choose this as our transversal surface for the calculation of @xmath70 .    in the cartoon in figure  [ fig : sephft2 ] ( left - hand image ) the direction of the flux velocity @xmath81 is shown for a surface of field lines bounding the non - ideal region .
the situation considered here can be compared with that without null points ( hornig & priest 2003 ) which is indicated in figure  [ fig : sephft2 ] ( right - hand image ) . in the non - null case
the flux velocity has the well - known feature of counter - rotating flows . with the addition of null points ,
the cross - section of the flux tube which bounds the non - ideal region is split into two separate domains connected by a singular line , the spine , for each of the two nulls . along the spine
the flux velocity becomes infinite .
this feature can be also found in figure  [ fig : fluxvel ] , which shows the component of @xmath82 perpendicular to the @xmath50 boundary at @xmath83 .
the left - hand image shows streamlines of @xmath82 ( grey lines ) and the intersection of the fan surface of the lower null with the boundary ( black line ) .
the flow direction is counter - clockwise .
the magnitude of @xmath84 is shown in the right - hand image .
a logarithmic scale is taken and the flow becomes infinite at the point of intersection of the upper null s spine ( @xmath85 ) . the flux velocity shown in figure  [ fig : fluxvel ] is typical of that in the early evolution when only one separator is present .
a rotational flux velocity is found within the flux tube passing through the non - ideal region .
additionally , the fan surface behaves as if advected by this flux velocity .    with these considerations in mind we may consider the nature of the reconnection to be as follows .
for exactness in the discussion we assume the field lines are fixed below the non - ideal region although a symmetric situation occurs in the alternative case with field lines fixed above the non - ideal region .
a magnetic field line threading the non - ideal region leaves the domain through a side boundary , @xmath50 , say ( that for which the flux velocity is illustrated in figure  [ fig : fluxvel ] ) .
there are two distinct ways the field line topology changes , illustrative phases of which are shown in figure  [ fig : recntypes ] :    1 .   the first situation is shown in the upper panel of figure  [ fig : recntypes ] .
the field line initially lies to the right of the fan surface of the lower null ( flux domain iv ) .
the motion of this unanchored end is towards the spine of the upper null .
it reaches the spine in a singular moment at which the flux velocity is infinite . here
the field line is connected to the upper null , lying in the fan surface of that null ( this is a dynamic situation with the fan surfaces moving in time ; the fan sweeps across the anchored end of the field line ) .
the field line then moves through the spine and is connected to the opposite side boundary as it leaves the domain ( flux domain i ) . in this process reconnection is occurring continuously while the field line is connected to the non - ideal region .
the ` flipping ' of of the magnetic field line and change in its topology occurs at a singular moment during the reconnection .
the second situation is shown in the lower panel of figure  [ fig : recntypes ] .
the field line initially lies to the left of the fan surface of the lower null ( flux domain iii ) .
the motion of the unanchored end is again towards the spine of the upper null .
the fan surface of the upper null is moving away from the anchored end of the field line while order to change its global connectivity ( pass through the upper null ) the entire field line must be lying in the fan plane of that null . accordingly
a ` flipping ' of the field line can only occur at or after the first bifurcation of separators . in this bifurcation a fold ( pocket ) in the fan surface is created which sweeps up and over the anchored end of the field line . in the moment at which the global field line connectivity changes the field line is connected to the ( rising ) pocket of the fan and the upper null . in the following instant the unanchored end of the field line moves to the opposite side boundary ( @xmath56 ) and the field connectivity is the new type vi .    while the separator reconnection configuration is sometimes viewed as a three - dimensional analogue of two - dimensional reconnection case at an x - type null point
, this example illustrates that the behaviour of the two is , in many ways , quite distinct .
in the two - dimensional case an analysis of the magnetic flux velocities ( see , for example , hornig 2007 ) shows that the reconnection of field lines occurs _ at _ the magnetic null point where the magnetic flux velocity is infinite and the field lines are ` cut ' and rejoined .
our analysis demonstrates that in the separator case a magnetic null is also key to the reconnection with field lines passing through the null at the moment their global connectivity changes , again when the flux velocity is infinite .
however , this is a non - local process with the null itself far removed from the reconnection site ( indeed the null may lie in an ideal environment ) .
the magnetic separator itself has only an indirect role in the process .
we suggest that , physically , the locations of the singularity in the flux velocity may be associated with regions of strong particle acceleration in real separator reconnection events .    having considered the way in which field line reconnection occurs in the separator configuration we proceed next to a quantitative analysis where we ask how to measure and interpret the rate of reconnection in the configuration with one or more separators . the question here is whether it is necessary to know the global field topology ( including the location and number of magnetic separators ) in order to determine the reconnection rate .
in a two - dimensional configuration where reconnection takes place at an x - type null point of the field the reconnection rate is given by the value of the electric field at the null point and measures the rate at which magnetic flux is transferred between the four topologically distinct flux domains . in order to express the rate as a dimensionless quantity that electric field
is normalised to a characteristic convective electric field and so the reconnection rate measured in terms ( fractions ) of the alfvn mach number .    in three - dimensions
we also have a measure for the rate of reconnection .
this is given by the maximum integrated parallel electric field over all field lines that thread the non - ideal region ( schindler _ et al . _
1988 ) : @xmath86 however , the interpretation of the ( maximum ) integral as a unique reconnection rate relies on the assumption that the topology of the magnetic field in the reconnection region is simple .
this assumption is justified in our case only up to the time of the bifurcation of separators .
the formulation is consistent with the two - dimensional measure of the electric field at the null with the two - dimensional reconnection rate being the three - dimensional reconnected flux per unit length in the invariant direction .
note that the question of how or whether to normalise the three - dimensional reconnection rate has not been properly addressed .
in contrast to the two - dimensional case a ` high ' reconnection rate can be obtained by having a strong value of the electric field or a long non - ideal region ( long path of the integral ) and the question of what constitutes ` typical ' convective electric field , or field strength to normalise it to , is not easily answered .    to examine the reconnection rate in the model presented here ,
we first consider only the early stages of the evolution , @xmath87 $ ] , when a single separator is present ( along the @xmath19 line and with @xmath88 $ ] ) . during this time
the field line with the field line with the maximum integrated parallel electric field along it is the separator field line itself .
the linear increase in time of the flux ring ( simulating the reconnection ) implies that the electric field at this line is constant in time and so the reconnection rate of the configuration for @xmath87 $ ] is given by @xmath89 in this case the coincidence of the separator and the line of maximum parallel electric field gives the clear and intuitive interpretation of the reconnection rate as the rate at which flux is transferred between the topologically distinct flux domains ( from regions ii and iv and into regions i and iii ) .
following the first bifurcation three separators are present in the domain ( and five following the second bifurcation ) ; we aim to determine how the reconnection rate should be measured and interpreted in this situation .
we label the central separator @xmath90 and the two separators lying off the central axis following the first bifurcation @xmath91 and @xmath92 ( by symmetry the order is not important ) .
the values of the quantities @xmath93 over time are shown in figure  [ fig : recnratedata ] along with the associated cumulative reconnected fluxes ( time integrals of the reconnection rates ) . for
@xmath94 $ ] a total of @xmath95 units of flux are reconnected through each of the separators @xmath91 and @xmath92 .    the physical interpretation for the fluxes @xmath96 and @xmath97 comes from considering the difference in reconnected flux between the separators @xmath90 and @xmath91 ( or , equivalently , @xmath92 ) in the interval @xmath98 $ ] .
this is given by @xmath99 this value is that of the magnetic flux passing through the surface bounded by the separators @xmath91 and @xmath90 at the end of the reconnection and so the flux contained in each of the new flux domains ( regions v and vi ) at the end of the reconnection ( @xmath45 ) .
overall , the situation is illustrated in figure  [ fig : recnrateinterpret ] where the flux transport between domains ( across fan surfaces ) is shown . in the figure
the black arrows show flux transport in the situation where field lines are considered as as fixed from below the non - ideal region .
the alternative , the flux transport obtained by considering field lines as fixed above the non - ideal region , is shown in grey .
when one separator is present ( left - hand image ) the integrated parallel electric field along that separator gives information on flux transport between the four domains in the relatively simple way already discussed .
when multiple separators are present ( right - hand image ) the separator with the highest integrated parallel electric field provides the primary transport of flux while the secondary separators provide additional information on the amount of flux in each domain . considering flux in domain vi for example ( the lower closed domain ) , the transport of flux into the domain ( thick grey and black arrows ) associated with the primary , central separator is faster than transport out of the domain ( thinner grey and black arrows ) associated with the secondary bounding separator .
this results in a net transport of flux into this region .    when three separators are present the reconnection rate with respect to the four basic flux connectivity types ( corresponding to flux domains i , ii / vi , iii , iv / v , as described in section  [ sec : bif ] ) is given by @xmath100 for @xmath101 , the situation considered here , this rate is lower than @xmath102.= the reduction in the rate of change of flux between domains compared with the reconnection rate as given by the maximum integrated parallel electric field across the region ( @xmath102 ) is a result of recursive reconnection between the multiple separators .
the recursive nature of the reconnection can be seen in figure  [ fig : recnrateinterpret ] . as an example considering flux as fixed from below the non - ideal region ( i.e.  with the grey arrows ) , the amount of flux in domain i is being reduced at a rate @xmath102 by reconnection into region vi but also increased at a rate @xmath103 from domain vi and rate @xmath104 from domain ii .    note that the situation @xmath105 is also conceivable . in this case
the flux in the regions @xmath106 and @xmath107 decreases in time and the additional flux coming from these regions can push the reconnection rate [ eq : threeseprate ] above the maximum of @xmath108 .
both situations show that due to the non - trivial topology of the magnetic field the reconnection rate with respect to the flux domains can differ significantly from the reconnection rate [ eq:3drec ] .
this is of particular relevance for the measurement and interpretation of reconnection since separators can be hard to detect in real systems ( and numerical simulations ) .
we remark further that if one is , for example , interested in particle acceleration rather than flux evolution , it may not be important that the reconnection processes can cancel each other in a topological sense . in this case the sum of all the reconnection rates , @xmath109 would give a better measure for the efficiency of the reconnection process .
the preceding analysis examined a rather particular situation in which a reconnection process was centred exactly on a magnetic separator that was identified as the reconnection line .
a natural question arises : is a real reconnection event likely to be centred in this way or will it perhaps encompass a separator but with a different magnetic field line being the reconnection line ?
it seems likely that both cases will occur in real reconnection events and so we now consider how our findings are altered by centring the reconnection away from the initial separator .
such a modification to the model can easily be made by taking @xmath110 , @xmath111 , @xmath112 as not all zero in equation   and making the appropriate adjustments to equation  .
a first question surrounds our finding that reconnection at a separator can create new separators . by choosing different locations ( @xmath113 )
at which to centre the reconnection region we can further address the question of whether reconnection in the vicinity of a separator tends to create new separators .
we find that this is indeed the case so long as the reconnection region overlaps the separator to some degree .
such an example is shown in figure  [ fig : offaxis ] ( left ) . if the reconnection region lies on or about a fan surface but not including a separator then the effect of the reconnection is to twist or wrap up the fan surface ( figure  [ fig : offaxis ] , right ) .
then only a very particular combination of reconnection events that distort two separate fan surfaces could lead to an intersection of the surfaces and the creation of a new separator .
we conclude that the creation of new separators by a reconnection event that includes a separator is a generic process .
the finding may help to explain the large number of separators found in mhd reconnection simulations such as those of parnell _ et al . _
( 2010b ) .
next we consider the nature of the magnetic flux evolution where a single separator is present but the separator itself does not give the maximum integrated parallel electric field .
figure  [ fig : sephft2 ] shows the magnetic flux evolution in the case where the separator and the maximum integrated parallel electric field do coincide ( the reconnection region is centred at the separator ) . here a rotational flux evolution in the magnetic flux tube threading the non - ideal region splits into four ` wings ' from the presence of the nulls .
a similar evolution occurs when the reconnection region is offset from the separator . in that case
the flux tube is split into four wings of unequal sizes but , nevertheless , the primary features of a continuous change in field line connectivity , a flipping of magnetic field lines as their topology is changed and the non - local involvement of the null all remain .
that is , it is not necessary for a separator to be the reconnection line in order for the the reconnection to have these distinguishing features of separator reconnection . in this way we can think of separator reconnection as a process that takes place when a magnetic separator passes through a reconnection region .
the separator need not itself be associated with the maximum parallel electric field .
separator reconnection is a ubiquitous process in astrophysical plasmas ( e.g.  longcope _ et al . _  2001
; haynes _ et al . _  2007
; dorelli _ et al . _
2007 ) but several fundamental details of how it takes place including how flux evolves in the process are not yet well understood .
here we have introduced a simple model for reconnection at an isolated non - ideal region threaded by a magnetic separator in an effort to better understand the basics of the process .
our first new finding is that reconnection events in the neighbourhood of a separator can ( and if sufficiently strong do ) create new separators .
this occurs even though the null points themselves remain in an ideal region .
such a bifurcation of separators has to occur in pairs ( from 1 to 3 , 5 , 7 , ... separators ) and the reverse process is , of course , also possible .
the fan planes of the two null points divide the space into distinct regions which can be distinguished due to the connectivity of their field lines with respect to any surface enclosing the configuration .
each bifurcation introduces a new pair of flux domains with a connectivity different from their neighbouring fluxes .
the rate of change of flux between four neighbouring flux domains is usually found by the integral of the parallel electric field along their dividing separator . however , in cases where multiple separators thread the non - ideal region the change in flux between domains is more complicated , as expressed here by equation  , and it will typically be less than the maximum of rates along the separators due to recursive reconnection between separators .
furthermore , considerations of asymmetric cases show that the maximum of @xmath114 does not have to occur along a separator and hence even for a single separator and a single isolated reconnection region the maximum integrated parallel electric field may not correspond to the rate of change of flux between domains . in three dimensions
the topology of the magnetic field must be known before a meaningful reconnection rate can be determined .    finally , the nature of the magnetic flux evolution in a separator reconnection process shows several distinguishing characteristics including rotational flux velocities that are split along the separator .
these characteristics are present regardless of whether the separator itself possesses the maximum integrated parallel electric field .
this analytical model is necessarily limited , one such limitation being in not addressing the way in which reconnection could be initiated at a separator .
simulations do show current layers forming in the neighbourhood of separators ( e.g. parnell _ et al . _  2010b ) but the extent of the layers may be variable .
for example , we have assumed here that the non - ideal region is localised to the central portion of the separator and the nulls remain in an ideal environment .
if a non - ideal region were to extend further along the domain and to include the nulls then further bifurcation processes can occur due to bifurcations of null points and a considerably more complex topology can evolve .      the authors would like to thank an anonymous referee for pointing out the possibility of having a reconnection rate with respect to the flux domains which is greater than any of the rates along separators .
hornig , g. 2007 , ` fundamental concepts ' in ` reconnection of magnetic fields : magnetohydrodynamics and collisionless theory and observations ' , eds .
birn , j. and priest , e.r .
, cambridge university press .
longcope , d. 2007 , ` separator reconnection ' in ` reconnection of magnetic fields : magnetohydrodynamics and collisionless theory and observations ' , eds .
birn , j. and priest , e.r . , cambridge university press .                                       plane at times indicated . superimposed is the vector field @xmath47 of horizontal magnetic field components at the relevant times .
a folding of the fan surfaces creates new separators in pairs as the fan surfaces intersect .
the lower - right  hand image indicates our labelling of the topologically distinct regions.,title="fig:",scaledwidth=26.0% ]   plane at times indicated .
superimposed is the vector field @xmath47 of horizontal magnetic field components at the relevant times .
a folding of the fan surfaces creates new separators in pairs as the fan surfaces intersect .
the lower - right  hand image indicates our labelling of the topologically distinct regions.,title="fig:",scaledwidth=26.0% ]   plane at times indicated .
superimposed is the vector field @xmath47 of horizontal magnetic field components at the relevant times .
a folding of the fan surfaces creates new separators in pairs as the fan surfaces intersect .
the lower - right  hand image indicates our labelling of the topologically distinct regions.,title="fig:",scaledwidth=26.0% ]   plane at times indicated .
superimposed is the vector field @xmath47 of horizontal magnetic field components at the relevant times .
a folding of the fan surfaces creates new separators in pairs as the fan surfaces intersect .
the lower - right  hand image indicates our labelling of the topologically distinct regions.,title="fig:",scaledwidth=26.0% ]   plane at times indicated .
superimposed is the vector field @xmath47 of horizontal magnetic field components at the relevant times .
a folding of the fan surfaces creates new separators in pairs as the fan surfaces intersect .
the lower - right  hand image indicates our labelling of the topologically distinct regions.,title="fig:",scaledwidth=24.8% ]    .
three magnetic separators are present where the fan surfaces intersect .
the location of the nulls is marked with black spheres .
the folding and intersection of the fan surfaces has created two new flux domains .
the corresponding flux tubes are also shown separately ( _ right _ ) where the three separators are marked with black lines ( the aspect ratio of the image has been further distorted to allow the main features to be identified).,title="fig:",scaledwidth=35.0% ] .
three magnetic separators are present where the fan surfaces intersect .
the location of the nulls is marked with black spheres .
the folding and intersection of the fan surfaces has created two new flux domains .
the corresponding flux tubes are also shown separately ( _ right _ ) where the three separators are marked with black lines ( the aspect ratio of the image has been further distorted to allow the main features to be identified).,title="fig:",scaledwidth=32.0%,height=377 ]     plane at @xmath45 . superimposed is the vector field @xmath47 .
five separators are present at this time . were
the reconnection process to continue , further folding of the fan surfaces would continue to generate separators in pairs in this manner .
( _ right _ ) at @xmath45 the various topologically flux domains are coloured according to their boundary connectivity types ( for the main regions compare with figure  3 ) showing the relation between connectivities of the various domains.,title="fig:",scaledwidth=42.5% ]   plane at @xmath45 . superimposed is the vector field @xmath47 .
five separators are present at this time . were
the reconnection process to continue , further folding of the fan surfaces would continue to generate separators in pairs in this manner .
( _ right _ ) at @xmath45 the various topologically flux domains are coloured according to their boundary connectivity types ( for the main regions compare with figure  3 ) showing the relation between connectivities of the various domains.,title="fig:",scaledwidth=40.0% ]        boundary at @xmath83 .
the left - hand image shows streamlines of the flux velocity ( grey lines ) and the location of the fan surface of the lower null ( black line ) . the right - hand image indicates the magnitude of the flow which becomes singular at the spine ( @xmath85 ) .
the direction of the flow is anti - clockwise .
, scaledwidth=70.0% ]     ( upper panel ) and @xmath115 ( lower panel ) . ,
title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) .
, title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) .
, title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) .
, title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) . ,
title="fig:",scaledwidth=19.0% ]     ( upper panel ) and @xmath115 ( lower panel ) . ,
title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) . ,
title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) . , title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) .
, title="fig:",scaledwidth=19.0% ]   ( upper panel ) and @xmath115 ( lower panel ) . , title="fig:",scaledwidth=19.0% ]     ( solid black line ) and the separators @xmath91 , @xmath92 ( solid grey line ) that form in the first bifurcation .
the cumulative reconnected flux across the separators is shown in the dashed lines of corresponding colour .
the dotted vertical line indicates the moment of bifurcation of separators.,scaledwidth=99.0% ]         . in the left - hand image
the reconnection region is centred at @xmath116 and overlaps the initial separator .
two more separators are created as the fan surfaces intersect . in the right - hand image
the reconnection region is centred at @xmath117 on a fan surface but not including the separator .
no new separator is created in this process.,title="fig:",scaledwidth=42.0% ] . in the left - hand image
the reconnection region is centred at @xmath116 and overlaps the initial separator .
two more separators are created as the fan surfaces intersect . in the right - hand image
the reconnection region is centred at @xmath117 on a fan surface but not including the separator .
no new separator is created in this process.,title="fig:",scaledwidth=45.1% ]