Abstract
In near-collisionless plasmas, which are ubiquitous in astrophysics, entropy production relies on fully-nonlinear processes such as turbulence and reconnection, which lead to particle acceleration. Mechanisms for turbulent reconnection include multiple magnetic flux ropes interacting to generate thin current sheets which undergo reconnection, leading to mixing and magnetic merging and growth of coherent structures, unstable reconnection current layers that fragment and turbulent reconnection outflows. All of these processes act across, and encompass, multiple reconnection sites. We use Magnetospheric Multi Scale four-point satellite observations to characterize the magnetic field line topology within a single reconnection current layer. We examine magnetopause reconnection where the spacecraft encounter the Electron Diffusion Region (EDR). We find fluctuating magnetic field with topology identical to that found for dynamically evolving vortices in hydrodynamic turbulence. The turbulence is supported by an electron-magnetohydrodynamic (EMHD) flow in which the magnetic field is effectively frozen into the electron fluid. Accelerated electrons are found in the EDR edge where we identify a departure from this turbulent topology, towards two-dimensional sheet-like structures. This is consistent with a scenario in which sub-ion scale turbulence can suppress electron acceleration within the EDR which would otherwise be possible in the electric field at the X-line.
Similar content being viewed by others
Introduction
Collisionless magnetic reconnection1 and plasma turbulence2,3 are fundamental mechanisms that transfer energy across scales and between electromagnetic fields and particles4,5. Magnetic field topology is central to both these processes. In two dimensions, reconnection is mediated by structures and interactions that are topologically constrained. In a turbulent flow, multiple reconnection X-lines can co-exist6. In such a two-dimensional flow, magnetic flux islands interact to generate thin current sheets7 which undergo reconnection, leading to magnetic merging and growth of coherent structures8,9,10. A given reconnection current layer is also unstable and tends to fragment into a population of secondary magnetic islands and current filaments11,12,13, even when embedded in a laminar flow. These processes are seen both in direct numerical simulations11,14 and in in-situ observations15,16. Three-dimensional reconnection differs significantly from its two-dimensional counterpart due to turbulence6,17,18 and asymmetries of the reconnecting current sheet19,20. Fully kinetic three dimensional simulations of reconnection reveal a highly tangled flux rope structure21 and outflow regions that can become turbulent22. All of these processes are either turbulence in an accelerated outflow, or are the generation of, or interaction between, multiple X-lines or magnetic null points.
Here, we classify the magnetic field topology observed as the four MMS spacecraft fly through a well resolved reconnection site. These observations allow us to quantify the topology of turbulence within a single reconnection layer. In-situ field and plasma observations from a tetrahedron of four spacecraft can be combined to directly classify the magnetic field line topology23,24,25,26 which is fully characterised by the eigenvalues of the magnetic field gradient tensor near a magnetic null27. The MMS spacecraft separation defines a spatial ‘yardstick’, which is of order of the ion inertial range \(d_i\), for sampling magnetic field topology. Crucially, whilst the magnetic field gradients, and the corresponding magnetic topology, are directly spatially resolved on the scale of the spacecraft separation, how the topology varies spatially is indirectly captured on a much finer spatial scale with the fast sampling in time of the magnetic field simultaneously at all four spacecraft. Within a given plasma region this gives many such samples since the cadence of the observations is fast (magnetic field measurements at 8192 Hz). We will use this high cadence to track how the magnetic field topology varies as the spacecraft move through the reconnection region. We obtain the time evolution of the magnetic field gradient tensor (MFGT) invariants, including, for the first time, those related to the magnetic field deformation (strain rate) tensor. When the spacecraft encounter the Electron Diffusion Region (EDR), we identify a particular evolution of the topology of the disordered magnetic field, for which the MFGT invariants explore a trajectory that parallels that seen in systems governed by the Euler equations such as experimental realizations of hydrodynamic turbulence28,29. In the sub-MHD scales of the MMS reconnection encounter, plasma turbulence is supported by an electron-magnetohydrodynamic (EMHD) flow in which the magnetic field is effectively frozen into the electron fluid7.
In this MMS reconnection encounter, accelerated electrons are found at the edge of the EDR, where the magnetic topology becomes sheet-like rather than turbulent. The absence of turbulence would permit direct acceleration of suprathermal electrons at the EDR edge. Heliospheric observations do not unambiguously confirm strong reconnection-driven electron acceleration; whilst it is seen at solar flares30,31 and in Earth’s magnetotail32, magnetosheath reconnection produces only mildly energised electron jets at a few Alfvén velocities which are often observed near, but not at, the reconnection X-line33,34. Our results are consistent with a scenario in which sub-ion scale turbulence can suppress electron acceleration within the EDR which would otherwise be possible in the electric field at the X-line.
Results
Figure 1 is an overview of the MMS fly-through of a magnetosheath reconnection site35,36,37,38. The left panel sketches a simplified 2D schematic showing regions of plasma inflow and outflow centered on a magnetic field null or X-line. The overall physical picture of this isolated reconnection site is of a magnetohydrodynamic (MHD) flow which sweeps magnetic flux into the inflow region, where MHD breaks down and the flux reconnects at the X-line within the small EDR marked by blue shaded rectangle. The process generates stretched field lines in the outflow region which relax, accelerating a plasma exhaust. Apart from energization in the outflow region, particle acceleration relies on physics on sub-MHD scales within the central current sheet. On scales between the ion and electron inertial lengths, \(d_i\!>\! d \!> \!d_e\), marked by the green shaded rectangle, the ions are demagnetized and the magnetic flux is “frozen” into the EMHD flow. Magnetic field observations from the full spacecraft tetrahedron, with spacing \(D_{sc}\! \sim \!70\) km or \(\sim \! d_i\), then characterize the topology of the magnetic field embedded in this EMHD flow. For this event, the spacecraft formation had high Tetrahedron Quality Factor at 0.95. Coloured circles with numbers indicate the positions of MMS spacecraft 1–4 at one time snapshot. MMS 4 is the closest to the Earth and may have encountered the EDR prior to the other spacecraft, however, not all data signatures are conclusive in supporting this encounter. A previous study found evidence of secondary flux ropes produced by the Kelvin–Helmholtz instability during the ion diffusion region crossing by MMS 438, which could imply the existence of multiple EDRs within the interval considered here. While results presented below are consistent with the presence of such flux ropes, the accuracy of the magnetic field measurements does not allow us to draw strong conclusions as to the presence of the multiple X-lines at these times. Among MMS spacecraft 1, 2 and 3, MMS 3 encounters the reconnection layer first. The right-hand panels of Fig. 1 show the MMS 3 orbit and the spacecraft formation on the day of interest. The red triangle in panels (a) and (b) shows the MMS3 location at the time of the reconnection event. The spacecraft were located towards the dusk flank of the magnetopause, at the distance of approximately \(8.8R_E\) from the Earth.
The time series from a single spacecraft is at sufficiently high cadence to capture features at the scale of the electron inertial length. The time series from spacecraft MMS 3 is shown in Fig. 2. This MMS 3 encounter with the EDR on the 19 of September 2015 at around 07:43:30 has been previously documented35,36,37,38, see Methods for details. The relevant features of this encounter with the reconnection region are as follows. The magnetic field null is encountered by MMS 3 at 07:43:30.488, and this time is used throughout to define \(t=0\). The MMS encounter with the reconnection site, that is, the time interval in which at least one MMS spacecraft is within the EDR region, is indicated by green shading across all panels. Blue shading across all panels indicate times when MMS 3 is within the EDR.
In Fig. 2, panel (a) plots the magnetic field magnitude (black), and the electric field component \(E_N\) (blue) in the event LMN coordinates defined as described in35, a more general description can be found in18. In the vicinity of the reconnection region, there is a coherent negative excursion of the electric field which defines the EDR35, this is offset from the magnetic field null (yellow triangle). There are strong fluctuations in both the electric and magnetic field in the vicinity of the EDR. The fluctuations in the magnetic field components (blue) and magnitude (black) of the combined FGM/SCM data, captured by a band pass filter with frequency range 64–256 Hz, are plotted in panel (b). These fluctuations are a combination of turbulent electron motion as well as whistler waves37, as evident from the power spectral density of magnetic field components (see figure 4 of the SI document). The magnetic components fluctuate strongly in the vicinity of the EDR (time interval indicated by red dashed lines) and there is a suppression of magnetic field magnitude variation within the EDR. Panel (c) plots the electric field fluctuations obtained with the same method as those in panel (b). The electric field fluctuates strongly in the vicinity of the EDR and these fluctuations are suppressed within the EDR. The bulk electron flow speed is enhanced in the EDR whereas the bulk proton flow remains approximately constant and sub-Alfvénic throughout, as shown in panel (d). There is evidence of electron heating at the edge of the EDR seen in \(\textbf{j} \cdot \textbf{E}_e\) calculated in the electron frame of reference (see Methods), shown in panel (e). Suprathermal electrons coincide with the excursion of the electric field within the EDR, shown in panel (e) which plots differential electron flux integrated over all angles and over the energy range 800–\(1\times 10^4\) eV. Finally, the bulk electron anisotropy is enhanced where the magnetic field most strongly fluctuates, just outside the EDR. Within the EDR, the bulk electrons are essentially isotropic. The observed change in the electron temperature anisotropy is rather modest, with a maximum at about \(60\%\) of the mean. Finally, Panel (h) plots time series of the MFGT invariants Q and R which we discuss next.
This set of time series from MMS 3 provides a clear chronology of the distinct physical regions of a well resolved single reconnection layer. However, if we combine magnetic field observations from all four MMS spacecraft, we can classify the magnetic topology sampled as the spacecraft fly through the distinct physical regions, using timings from MMS 3. The magnetic field observed at four MMS spacecraft is used to construct the MFGT \(\varvec{M}(t)\). We have applied two different methods, one based on a direct matrix inversion39, 40 and another based on nonlinear least squares optimisation41,42 (the lsqnonlin function in Matlab). We have found that the least squares method gives an MFGT which satisfies \(\nabla \cdot \textbf{B} = 0\) to much higher accuracy. All results presented below are obtained using the least squares optimisation. We give a more detailed description of these two approaches in the Magnetic Field Gradient Tensor section. The topological invariants, R and Q, of this tensor characterise the magnetic field line topology23,25,26.
The symmetric part of \(\varvec{M}\), \(\varvec{S}\), with invariants \(R_s\) and \(Q_s\), relates to a rate of strain of magnetic field, that is, all curl-free deformations of the magnetic field lines which are not related to a current. All invariants are normalised by the power in the current \((\mu _0 j)^2\) (see Methods). As the four spacecraft travel through the reconnection region, they provide multiple samples of the magnetic field topological invariants and these are plotted in Fig. 3. Figure 3 is structured as follows: Panel (a) plots the magnetic field magnitude measured by MMS 3 over a slightly extended time interval compared with Fig. 2a, with the field null (yellow triangle) at \(t=0\). Green shading again indicates the time interval in which at least one MMS spacecraft samples the EDR. The EDR region is indicated by blue shading. The magnetic field magnitude time trace is color coded to indicate sub-intervals of the encounter: dark green, blue and red traces respectively indicate the leading edge of the EDR, the EDR, and the trailing edge of the EDR, respectively. The regions before and after the EDR is encountered are indicated by grey and black time traces, respectively. These colors are used throughout Fig. 3 to refer to MMS’s encounter with different regions of the reconnection site.
This high resolution time series of topological invariants can then be plotted in the invariant space, for the field gradient tensor \(\varvec{M}\) in R, Q space (panel (b)) and for the deformation tensor \(\varvec{S}\), in \(Q_s(t),R_s(t)\) (panel (c)). The overall magnetic field topology in R, Q space is discriminated by the contour \(\Delta \!=\!(27/4)R^2+Q^3\) (see Methods). Elliptic magnetic field lines, that is, flux ropes, satisfy \(\Delta >0\), whereas for \(\Delta <0\) the field lines are hyperbolic, forming a multi-pole separatrix structures consistent with 3D reconnection (3D X-line). In terms of the eigenvalues \(\lambda ^s_1\), \(\lambda ^s_2\), \(\lambda ^s_3\) of the deformation tensor S, \(R_s \!= \!-\lambda ^s_1 \lambda ^s_2 \lambda ^s_3\) and \(Q_s \!= \! \lambda ^s_1 \lambda ^s_2 + \lambda ^s_2 \lambda ^s_3 + \lambda ^s_1 \lambda ^s_3\). In \(R_s\), \(Q_s\) space, the contour \(\Delta _s=0\) (where \(\Delta _s \!=\!(27/4)R_s^2+Q_s^3\)) is the boundary of curl-free deformations of the magnetic field, corresponding to eigenvalue ratios 2:\(-1\):\(-1\) (\(R_s<0\)) and 1:1:\(-2\) (\(R_s>0\)). The topology tends to two-dimensional, stretched field lines in the vicinity of \(R=0\) and \(R_s=0\). The sign of R and \(R_s\) is determined by the overall magnetic field polarity through the reversal.
We now use Fig. 3 to track how the magnetic field topology evolves as the MMS spacecraft tetrahedron transits the reconnection site. Panel (b) tracks the overall topology and at times well outside the EDR and the encounter (grey and black points), the topology of structures in the magnetic field have \(\Delta\) just above zero and \(R\simeq 0\), that is, they are quasi 2D flux ropes. This is consistent with the previous findings of secondary flux ropes at the time of MMS 4 encounter with the ion diffusion region38 (grey points). The brief incursion of black points below the \(\Delta =0\) separatrix corresponds to a second minimum of the magnetic field strength at \(t\! \approx \! 2\) seconds after the magnetic null is encountered by MMS 3. Within the encounter, at the leading edge of the EDR (green trace), the topology tracks the \(\Delta =0\) separatrix and then transitions into hyperbolic field topology region \(\Delta <0\). Within the EDR region (blue traces) the topology is consistently X-line type (\(\Delta <0\)) and is closer to \(R\simeq 0\), indicating reduced dimensionality of magnetic field in one of the directions. The field null (yellow triangle) is also well within (\(\Delta <0\)). The overall topology shown in panel (b) is then just that expected for an encounter with a reconnection region, outside of the encounter, it is flux-rope like, and during the encounter, separatrix-like. The main MFGT invariants are not expected to vary in time at magnetic null point25, as can be seen to be the case in the time trace of these invariants plotted in panel (h) of Fig 2.
The invariants \(R_s,Q_s\) of curl-free deformations are well known to be informative in the study of hydrodynamic turbulence28,29. Here, we obtain these for the first time for the magnetic field; these are plotted in panel (c). During the MMS encounter with the reconnection site (green shading), we then find two distinct field deformation topologies. Turbulent-type topology is seen during the MMS 3 encounter with the leading (green) and trailing (red) edges of the EDR. For this topology, the invariants track the magenta contour of triaxial deformations with eigenvalue ratios − 3:\(-\,\)1:4 (\(R_s<0\)) and 3:1:\(-4\) (\(R_s>0\)), including the points near the magnetic null indicated by a yellow triangle. There is a departure from this contour, to dimensionally reduced topology when MMS 3 is within the EDR (blue), during which time it is confined to a line of constant \(R_s\).
The turbulent type topology that is found in the magnetic field during the reconnection encounter coincides with that seen in the deformations of the velocity stream lines of neutral fluid turbulence on dissipative scales, obtained from the observed velocity gradient tensor28,29. As detailed in the Discussion below, this suggests a turbulent Euler-EMHD flow supported by the bulk, thermal electron population. Dimensionally reduced topology is seen in the departure from the magenta contour when MMS 3 is within the EDR (blue). During this time the magnetic field deformations plotted in panel (b) are instead confined to a line of constant \(R=0\). This indicates that the variation of the magnetic field in one of the directions is fixed by the eigenvalues of the remaining two directions. The third eigenvalue is then fixed by \(\lambda _3 = K(\lambda _1 \lambda _2)^{-1}=0\). It is a field line stretching which will render magnetic structures more sheet-like and will suppress the 3D EMHD supported turbulence. This is a path by which the curl-free, and perforce divergence free, magnetic field deformations can evolve toward a field null.
Anisotropy in the bulk electron thermal populations coincides with electric and magnetic field fluctuations at the edge of the EDR, whereas suprathermal electrons are seen within the EDR which is a region of persistently non-zero electric field (Fig. 2f,g). This is consistent with both simulations43 and observations33,34. Within the reconnection encounter, we have identified the topology of EMHD turbulent-type fluctuations. These turbulent fluctuations could quench the direct electric field acceleration of electrons to suprathermal energies. The overall effect would be to suppress electron acceleration throughout the reconnection encounter in regions where the topology tracks the turbulent type contour in Fig. 3c. Enhanced anisotropy of the bulk thermal electron population, and both electric and magnetic field fluctuations coincide with this turbulent-type topology of the magnetic field. Within the EDR, the topology is dimensionally reduced and departs from the signature of turbulence-type topology; this coincides with enhanced suprathermal electron flux (Fig. 2f) and isotropic bulk electrons.
We verify this scenario in Fig. 4 where we directly map the observed electron suprathermal energy, \(\mathcal {E}_e\), (integrated differential flux in the range of 800 eV to \(10^4\) eV) and electron temperature anisotropy \(\alpha _e\) onto the observed magnetic deformations in the same \(R_s,Q_s\) space as Fig. 3c. The red contour again plots \(\Delta _s=0\), the boundary of curl-free deformations of the magnetic field and the magenta contour indicates turbulent type topology identified in Fig. 3(c). The details of how these plots are constructed are given in the Methods. We consider three time intervals: (i) prior to the EDR encounter by any of the spacecraft, (ii) during the encounter (as in Fig. 3) and (iii) after the encounter with the EDR. These intervals are indicated on the MMS 3 time series in panels (a) and (b) of Fig. 4 where we show differential energy flux of suprathermal electrons, and the electron temperature anisotropy. Suprathermal electrons appear intermittently, and at the leading edge of the EDR rather than at the magnetic field null, as reported previously (eg44).
Panels (ci), (ciii) and (di), (diii) map quantities \(\mathcal {E}_e\) and \(\alpha _e\) onto \(Q_s,R_s\) invariants space when all MMS spacecraft are outside of the EDR. The colour scale is common to all three panels in each row. We see that for time intervals (i) and (iii) the suprathermal energy is low and the temperature anisotropy is approximately 1, and there is no clear organization of energization with magnetic field line topology. During the reconnection encounter, time interval (ii), we again see the points track the turbulent like topology contour except during the EDR encounter where the topology is dimensionally reduced. The peak integrated energy flux \(\mathcal {E}_e\) (cii) of suprathermal electrons coincides with dimensionally reduced field topology, whereas the peak of the thermal (bulk) electron temperature anisotropy \(\alpha _e\) (dii) coincides with turbulent type magnetic field topology. The time at which the magnetic null point is encountered, indicated by the yellow triangle, does not coincide with a flux of energetic electrons. The field null is located at the transition between dimensionally reduced, and turbulent type topology. Figure 4 required temporal averaging of the topological invariants, they are low-pass filtered to match the timescale of the electron plasma observations (see Methods).
Particle distribution functions corroborate these findings and these are plotted in the SI. These plots show non-Maxwellian distributions for both the parallel and perpendicular velocity VDF. The parallel velocity VDF show anisotropic heating at the time when the peak in the electron temperature anisotropy is found in Fig. 2 panel (g). The perpendicular velocity VDF does not change significantly until the MMS 3 is in the EDR, when it shows both an shift in its mean and elevated tails.
Discussion and conclusions
High cadence, four spacecraft observations have allowed us to directly determine the topology of the magnetic field in-situ within a reconnection site. The magnetic field gradients, and the corresponding magnetic topology, are directly spatially resolved on the scale of the spacecraft separation. However, how the topology varies spatially is indirectly captured on a much finer spatial scale with the fast sampling in time of the magnetic field simultaneously at all four spacecraft. We resolve the trajectory in the invariant space that identifies this topology on the fast sampling timescale. Thus we can directly compare how the magnetic field topology is changing with signatures of electron heating and acceleration which are also available at high time cadence.
Previously, mechanisms proposed for turbulent reconnection6,7,8,10,11,12,13, are processes that act across, and encompass, multiple reconnection sites. We have identified the topological signature of fluid-like turbulence within a single reconnection site. We find that the topology of turbulent deformation of magnetic field lines directly orders electron acceleration. Fluid-like turbulent topology is found where suprathermal electrons are absent; when the topology transitions to reduced dimensionality and is closer to sheet-like in the electron diffusion region (EDR), suprathermal electrons are encountered. This suggest that the magnetic topological features that we have identified are thus an integral part of how electrons are accelerated by reconnection.
The topology of turbulent deformation of magnetic field lines that we have identified coincides with that seen in an ideal turbulent hydrodynamic flow. The correspondence between magnetic field rotation-free deformations and the behaviour of the dissipative structures in HD turbulence can be understood as follows. The spacecraft tetrahedron spans scales on which the ions are effectively demagnetised, so EMHD will apply. For a near-collisionless plasma, considered here, the magnetic field is effectively frozen into the electron fluid45, and if we can neglect electron inertia, the induction equation then becomes:
where \(\varvec{v}_e\) is the velocity field of the electrons. We then expect the topology of fluctuations in the magnetic field, and its scaling properties, to track that of the electron fluid. When the scales of the fluctuations become very short, on scales \(d \ll d_e\), where the electron vorticity \(w_e\gg j/d_e\), EMHD reduces to the Euler equation for the electron fluid velocity field7. The topology of this turbulent electron flow will then be identical to ideal hydrodynamic flow, as they are both constrained by the same Euler equation. It is the clear time-ordered flow of the invariants in their phase space near the EDR (Fig. 3c) that supports this hypothesis. Away from the encounter, we do not detect such ordered flow and the invariants show more random behaviour, indicative of sampling uncorrelated structures with different topology. Taking the maximal electron fluid velocity of \(\sim \!800\) km s\(^{-1}\) as measured by the MMS 3 at the EDR, and the sampling frequency of the magnetic field at 8192 Hz, we obtain the corresponding spatial scale of the temporal measurement to be \(d \approx \! 0.1\) km, as compared to the electron inertial length \(d_e = 1.2\) km.
An alternative scenario is fast self-similar current sheet fragmentation process via finite length plasmoid production. The lack of lasting elliptic magnetic field line features, indicative of plasmoids, within the encounter suggests that this hypothesis is less likely.
Our results provide direct observational constraints on the topology that can discriminate these mechanisms. In this first study we have analysed one well-sampled reconnection event. Routine computation of field topological invariants in the studies of reconnection, both in observations, and simulations, would establish whether this behaviour is ubiquitous, or indeed, universal, or under what conditions it occurs.
Methods
Reconnection event overview
There are several quantities which indicate an encounter of the spacecraft with the EDR and these are shown in panels (a)–(h) of Fig. 2 in the main text. The extent of the EDR is marked by the blue shaded region near \(t=0\) for each panel. Panel (a) of Fig. 2 shows the decrease in magnetic field magnitude near the X-point crossed at 07:43:30.488 and this crossing defines the epoch time \(t=0\). We use merged FGM46 and SCM47 instrument magnetic field data with the sampling frequency 8192 Hz. The blue trace shows one electric field component, \(E_y\) in GSE coordinates, obtained from the EDP instrument48,49 with sampling frequency of 8192 Hz. A secular negative excursion of this field component in the vicinity of \(t=0\) defines the EDR35. The magnetic and electric field magnitude and component fluctuations, shown in panels (b) and (c), are calculated using a band pass filter with frequency range 64–256 Hz. Panel (d) shows the expected behaviour of the ion and electron velocities, normalised to the downstream Alfvén speed, \(v_{A,ds} = B_{ds} / \sqrt{\mu _0 n_{p,ds} m_p}\), with electrons accelerated by a factor of few Alfvén speeds within a narrow region within the EDR. The electron velocity distributions and their moments are from FPI-DEM instrument50 with sampling frequency of 33 Hz. Panel (e) shows the power density transferred between particles and fields, \(P^m_e = \varvec{j} \cdot \varvec{E}_e\), where \(\varvec{E}_e\) is the electric field in the electron frame of reference51, \(\varvec{E}_e = \varvec{E} + \varvec{v}_e \times \varvec{B}\) and the current density obtained from particle data, \(\varvec{j}=q n_e(\varvec{v}_i - \varvec{v}_e)\). A positive \(P^m_e\) indicates a loss of electromagnetic energy, which is converted to macroscopic electron flows. Where particles give their energy back to the electromagnetic fluctuations \(P^m_e\) is negative. Overall we observe a net transfer of fields’ energy to particles in the shaded region. The MMS mission measures the differential electron flux as a function of energy and pitch angle. Panel (f) shows the omni-directional electron flux integrated over a range of energies, \(\mathcal {E}_e\), defined as
where \(E_0=800\) eV and \(E_1=10^4\) eV are the chosen minimum and the maximum energy values and \(\Phi _e\) is the differential electron flux in the units of [(cm\(^2\) s sr)\(^{-1}\)]. We observe a modest increase in the electron flux for higher energies, and critically the maximum is not at the suspected X-line position \(t=0\). The reversal of the magnetic field near the X-point, shown in panel (a), must generate an electric field52 which, at least in principle, can accelerate electrons to extremely large energies. Three-dimensional PIC simulation shows that such direct acceleration is not efficient when the guide field is weak, as it is in our case53. On larger scales, the fast release of the local magnetic field tension re-configures the field and these changes propagate at approximately Alfvén speed. Away from the X-point particles are still magnetized and are accelerated by this motion of the newly generated magnetic field lines. This should lead to increase in the electron temperature anisotropy, \(\alpha _e\), which is shown in the panel (g). These quantities are also plotted in Fig. 4a,b. In order to project the quantities \(\mathcal {E}_e\) and \(\alpha _e\) onto the phase space of the magnetic field topological invariants (\(Q_s\), \(R_s\)), we average the invariants between each two consecutive electron measurements (over approximately 248 magnetic data points) obtaining integrated electron flux and the electron temperature anisotropy associated with an averaged value of invariants (\({\overline{R}}_s\) and \({\overline{Q}}_s)\)), \(\mathcal {E}_e({\overline{R}}_s,{\overline{Q}}_s)\) and \(\alpha _e({\overline{R}}_s, {\overline{Q}}_s)\).
Magnetic field gradient tensor
We have obtained the MFGT using two approaches, based on the Taylor expansion of the magnetic field at a point which retains only linear terms in the spacecraft position. A magnetic field vector is then written as
where \(\varvec{r}\) represents spacecraft position, in our case, in GSE coordinates.
The first approach uses differences between spacecraft pairs and expresses the relation between magnetic field increment and spacecraft separation vectors as a matrix equation39,40. Defined with respect to the centroid, \(\varvec{\rho }_{1}\!=\!(\varvec{r}_1-\varvec{r}_2)/\sqrt{2}\), \(\varvec{\rho }_{2}\!=\!(\varvec{r}_1+\varvec{r}_2-2\varvec{r}_3)/\sqrt{6}\), and \(\varvec{\rho }_{3}\!=\!(\varvec{r}_1+\varvec{r}_2+\varvec{r}_3-3\varvec{r}_4)/\sqrt{12}\). Similarly, \(\varvec{b}_{1}\!=\!(\varvec{B}_1-\varvec{B}_2)/\sqrt{2}\), \(\varvec{b}_{2}\!=\!(\varvec{B}_1+\varvec{B}_2-2\varvec{B}_3)/\sqrt{6}\), \(\varvec{b}_{3}\!=\!(\varvec{B}_1+\varvec{B}_2+\varvec{B}_3-3\varvec{B}_4)/\sqrt{12}\), and \(\varvec{B}_s\) is the magnetic field vector at spacecraft s. Inverting the relation \(\textbf{b}\!=\!\varvec{A}\varvec{\rho }\) with the constraint \(\nabla \cdot \varvec{b}=0\) gives the tensor elements \(M_{ij}\!=\!(\varvec{b}\varvec{\rho }^{-1})_{ij}-\delta _{ij}{\text {tr}} \left( \varvec{b} \varvec{\rho }^{-1} \right) /3\). For the interval considered here, this method produced MFGTs with the average trace values of \(\sim \! 1 \times 10^{-4}\), while appropriately normalised tensor invariants of interest reached values of orders \(0.1-1\).
The second approach is based on nonlinear least squares fits of the polynomial (3) to the measured magnetic field41. Note that this method has 12 parameters, 3 components of \(\varvec{B}(\varvec{r}_0)\) and 9 elements of the MFGT. Given 4 spacecraft, each measuring 3 components of the magnetic field vector, we have 12 equations to be solved, thus the problem is well-posed. We have used Matlab’s lsqnonlin function to perform the least squares fit with the constraint of the divergence-free solution. This method produced an MFGT with trace no larger then \(\sim \! 1 \times 10^{-14}\), while appropriately normalised tensor invariants were still of order 0.1–1. Recently, second order methods have been proposed, which use particle currents at each spacecraft in addition to magnetic field data41,42. These methods may give more accurate results when the current is strong, which is not the case in our interval.
Figure 5 shows differences between the phase spaces (R, Q) and (R\(_s\), Q\(_s\)) obtained by these two different approaches. It is clear that the overall flow of invariants in their phase space is the same, but there are small differences in how the invariants approach 2D line of \(R=0\) in panels (a) and (c). The nonlinear least squares method shows that the magnetic field topology has reduced dimensionality at the EDR, that is, the magnetic field variation in one of the direction is negligible. The dynamics of the points within the EDR, those shown in blue, is also different. The matrix inversion method shows almost no variation in the values of R, while the least squares method does show some variations.
The topological invariants, R and Q, of this tensor characterise the magnetic field line topology23,25,26. Given three eigenvalues of the MFGT, \(\lambda _1\), \(\lambda _2\) and \(\lambda _3\), the main invariants can be computed from the following expressions54:
The symmetric part of \(\varvec{M}\), \(\varvec{S}=(\varvec{M}+\varvec{M}^T)/2\), relates to a rate of strain of magnetic field lines, that is, it accounts for all deformations of the magnetic field topology which are not related to a current. The invariants of \(\varvec{S}\) will be indicated by \(R_s\), \(Q_s\). The antisymmetric current density rate tensor \(\varvec{J}=(\varvec{M}-\varvec{M}^T)/2\) relates to current density, \(j_k\!=\!\epsilon _{ijk} J_{ij}\) via the Faraday law \(\nabla \times \varvec{B}=\mu _0 \varvec{j}\). All invariants are normalised by the power in \(j^2\), for example, we infer the topology from \(Q(t) \!=\! \tilde{Q}(t)/\langle j^2 \rangle\) and \(R(t) \!=\! \tilde{R}(t)/\langle j^2 \rangle ^{3/2}\), where \(\langle j^2 \rangle\) is the time average.
The sign of the quantity \(\Delta = (27/4)R^2+Q^3\) defines two regions in the (R, Q) plane. Invariants that give \(\Delta >0\) indicate elliptic magnetic field lines, while for \(\Delta <0\) the field lines are hyperbolic, forming a multi-pole separatrix structure consistent with 3D reconnection (3D X-line). When \(\Delta >0\), two eigenvalues of \(\varvec{M}\) are complex and one is real, indicative of the magnetic flux rope topology that generalises a two-dimensional O-point in three dimensions. We note that quasi-2D flux ropes with a negligible magnetic field variation along the structure on the scale of the spacecraft tetrahedron have \(R\! \sim \! 0\), while finite length, self-contained structures, such as three-dimensional plasmoids, have \(|R| \gg 0\). For \(\Delta <0\) all three eigenvalues are real and such hyperbolic magnetic field line topology is the 3D generalisation of a 2D X-point. The sign of R discriminates between two different directions of the magnetic field vector for each topology.
The deformation tensor invariants are defined in a similar way. If \(\lambda ^s_1\), \(\lambda ^s_2\), and \(\lambda ^s_3\) are the eigenvalues of tensor \(\varvec{S}\), then
Applying the same normalisation as was used for \(\tilde{R}\) and \(\tilde{Q}^3\), we obtain invariants \(R_s\), \(Q_s\). Recall that the eigenvalues of a symmetric tensor tensor must be real. This reduces the available phase space of the invariants \(R_s\), \(Q_s\), to a region \(\Delta _s = (27/4)R_s^2+Q_s^3 \le 0\). The red contour in Fig. 3c is the boundary of this available phase space, and corresponds to the axis-symmetric contraction (\(R_s<0\)) with the eigenvalue ratios of 2:\(-1\):\(-1\), or axis-symmetric expansion (\(R_s>0\)) with the eigenvalue ratios 1:1:\(-2\). The magenta contour in this panel coincides with the triaxial deformations given by eigenvalue ratios \(-3\):\(-1\):4 (\(R_s<0\)) and 3:1:\(-4\) (\(R_s>0\)). The eigenvalue ratio 3:1:\(-4\) (\(R_s>0\)) has been found to characterise the deformations of the velocity field gradient tensor in a three dimensional hydrodynamic turbulent flow28,29. The opposite sign of this ratio found in our work reflects a particular magnetic field polarity in the reconnection region.
Data availibility
The datasets generated and/or analysed during the current study are available in the MMS Science Data Center repository, https://lasp.colorado.edu/mms/sdc/public/.
References
Birn, J. & Priest, E. R. Reconnection of Magnetic Fields: Magnetohydrodynamics and Collisionless Theory and Observations (Cambridge University Press, New York, 2007).
Krommes, J. A. Fundamental statistical descriptions of plasma turbulence in magnetic fields. Phys. Rep. 360(1–4), 1–352 (2002).
Matthaeus, W. H. & Velli, M. Who needs turbulence?. Space Sci. Rev. 160(1), 145–168 (2011).
Lemoine, M. Particle acceleration in strong MHD turbulence. Phys. Rev. D 104(6), 063020 (2021).
Egedal, J., Daughton, W. & Le, A. Large-scale electron acceleration by parallel electric fields during magnetic reconnection. Nat. Phys. 8(4), 321–324 (2012).
Lazarian, A., Eyink, G., Vishniac, E. & Kowal, G. Turbulent reconnection and its implications. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 373(2014), 20140144 (2015).
Biskamp, D., Schwarz, E. & Drake, J. F. Two-fluid theory of collisionless magnetic reconnection. Phys. Plasmas 4(4), 1002–1009 (1997).
Linton, M. G., Dahlburg, R. B. & Antiochos, S. K. Reconnection of twisted flux tubes as a function of contact angle. Astrophys. J. 553(2), 905 (2001).
Karimabadi, H. et al. Coherent structures, intermittent turbulence and dissipation in high-temperature plasmas. Phys. Plasmas 20, 012303. https://doi.org/10.1063/1.4773205 (2013).
Zhou, M., Wu, D. H., Loureiro, N. F. & Uzdensky, D. A. Statistical description of coalescing magnetic islands via magnetic reconnection. J. Plasma Phys. 87(6), 905870620 (2021).
Drake, J. F., Swisdak, M., Schoeffler, K. M., Rogers, B. N. & Kobayashi, S. Formation of secondary islands during magnetic reconnection. Geophys. Res. Lett. 33, L13105. https://doi.org/10.1029/2006GL025957 (2006).
Daughton, W., Scudder, J. & Karimabadi, H. Fully kinetic simulations of undriven magnetic reconnection with open boundary conditions. Phys. Plasmas 13, 072101 (2006).
Lu, S. et al. Particle-in-cell simulations of secondary magnetic islands: Ion-scale flux ropes and plasmoids. Astrophys. J. 900(2), 145 (2020).
Del Sarto, D., Califano, F. & Pegoraro, F. Secondary instabilities and vortex formation in collisionless-fluid magnetic reconnection. Phys. Rev. Lett. 91(23), 235001 (2003).
Eriksson, S., Newman, D. L., Lapenta, G. & Angelopoulos, V. On the signatures of magnetic islands and multiple X-lines in the solar wind as observed by ARTEMIS and WIND. Plasma Phys. Control. Fusion 56(6), 064008 (2014).
Zhong, Z. H. et al. Stacked electron diffusion regions and electron Kelvin–Helmholtz vortices within the ion diffusion region of collisionless magnetic reconnection. Astrophys. J. Lett. 926(2), L27 (2022).
Strauss, H. R. Turbulent reconnection. Astrophys. J. 326, 412–417 (1988).
Eastwood, J. P., Phan, T. D., Bale, S. D. & Tjulin, A. Observations of turbulence generated by magnetic reconnection. Phys. Rev. Lett. 102(3), 035001 (2009).
Eastwood, J. P., Shay, M. A., Phan, T. D. & Oieroset, M. Asymmetry of the ion diffusion region Hall electric and magnetic fields during guide field reconnection: Observations and comparison with simulations. Phys. Rev. Lett. 104(20), 205001 (2010).
Fuselier, S. A. et al. Large-scale characteristics of reconnection diffusion regions and associated magnetopause crossings observed by MMS. J. Geophys. Res. Space Phys. 122(5), 5466–5486 (2017).
Daughton, W. et al. Role of electron physics in the development of turbulent magnetic reconnection in collisionless plasmas. Nat. Phys. 7(7), 539–542 (2011).
Leonardis, E., Chapman, S. C., Daughton, W., Roytershteyn, V. & Karimabadi, H. Identification of intermittent multi-fractal turbulence in fully kinetic simulations of magnetic reconnection. Phys. Rev. Lett. 110, 205002. https://doi.org/10.1103/PhysRevLett.110.205002 (2013).
Dallas, V. & Alexakis, A. Structures and dynamics of small scales in decaying magnetohydrodynamic turbulence. Phys. Fluids 25(10), 105106 (2013).
Fu, H. S. et al. How to find magnetic nulls and reconstruct field topology with MMS data?. J. Geophys. Res. Space Phys. 120(5), 3758–3782 (2015).
Quattrociocchi, V., Consolini, G., Marcucci, M. F. & Materassi, M. On geometrical invariants of the magnetic field gradient tensor in turbulent space plasmas: Scale variability in the inertial range. Astrophys. J. 878(2), 124 (2019).
Hnat, B., Chapman, S. C. & Watkins, N. W. Magnetic topology of actively evolving and passively convecting structures in the turbulent solar wind. Phys. Rev. Lett. 126(12), 125101 (2021).
Greene, J. M. Geometrical properties of three-dimensional reconnecting magnetic fields with nulls. J. Geophys. Res. Space Phys. 93(A8), 8583–8590 (1988).
Ashurst, W. T., Kerstein, A. R., Kerr, R. M. & Gibson, C. H. Alignment of vorticity and scalar gradient with strain rate in simulated Navier–Stokes turbulence. Phys. Fluids 30, 2343–2353 (1987).
Tsinober, A., Kit, E. & Dracos, T. Experimental investigation of the field of velocity gradients in turbulent flows. J. Fluid Mech. 242, 169–192 (1992).
Lin, R. P. Energy release and particle acceleration in flares: Summary and future prospects. Space Sci. Rev. 159(1), 421–445 (2011).
Zhang, Q., Guo, F., Daughton, W., Li, H. & Li, X. Efficient nonthermal ion and electron acceleration enabled by the flux-rope kink instability in 3D nonrelativistic magnetic reconnection. Phys. Rev. Lett. 127(18), 185101 (2021).
Øieroset, M., Lin, R. P., Phan, T. D., Larson, D. E. & Bale, S. D. Evidence for electron acceleration up to 300 keV in the magnetic reconnection diffusion. Phys. Rev. Lett. 89, 195001 (2002).
Imada, S., Hoshino, M. & Mukai, T. Average profiles of energetic and thermal electrons in the magnetotail reconnection regions. Geophys. Res. Lett.https://doi.org/10.1029/2005GL022594 (2005).
Phan, T. D. et al. Electron magnetic reconnection without ion coupling in Earth’s turbulent magnetosheath. Nature 557(7704), 202–206 (2018).
Chen, L. J. et al. Electron energization and mixing observed by MMS in the vicinity of an electron diffusion region during magnetopause reconnection. Geophys. Res. Lett. 43(12), 6036–6043 (2016).
Fuselier, S. A. et al. Magnetospheric multiscale science mission profile and operations. Space Sci. Rev. 199(1), 77–103 (2016).
Cao, D. et al. MMS observations of whistler waves in electron diffusion region. Geophys. Res. Lett. 44(9), 3954–3962 (2017).
Zhong, Z. H. et al. Evidence for secondary flux rope generated by the electron Kelvin–Helmholtz instability in a magnetic reconnection diffusion region. Phys. Rev. Lett. 120(7), 075101 (2018).
Paschmann, G. & Daly, P.W. (Eds). Analysis methods for multi-spacecraft data. ISSI Scientific Report SR-001 (1998)
Chertkov, M., Pumir, A. & Shraiman, B. I. Lagrangian tetrad dynamics and the phenomenology of turbulence. Phys. Fluids 11(8), 2394–2410 (1999).
Denton, R. E. et al. Polynomial reconstruction of the reconnection magnetic field observed by multiple spacecraft. J. Geophys. Res. Space Phys. 125(2), e2019JA027481 (2020).
Torbert, R. B. et al. A new method of 3-D magnetic field reconstruction. Geophys. Res. Lett. 47, e2019GL085542. https://doi.org/10.1029/2019GL085542 (2020).
Lapenta, G., Berchem, J., El Alaoui, M. & Walker, R. Turbulent energization of electron power law tails during magnetic reconnection. Phys. Rev. Lett. 125(22), 225101 (2020).
Haynes, C. T., Burgess, D. & Camporeale, E. Reconnection and electron temperature anisotropy in sub-proton scale plasma turbulence. Astrophys. J. 783(1), 38 (2014).
Lighthill, M. J. Studies on magneto-hydrodynamic waves and other anisotropic wave motions. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 252(1014), 397–430 (1960).
Russell, C. T. et al. The magnetospheric multiscale magnetometers. Space Sci. Rev. 199, 189–256 (2016).
Le Contel, O. et al. The search-coil magnetometer for MMS. Space Sci. Rev. 199(1), 257–282 (2016).
Ergun, R. E. et al. The axial double probe and fields signal processing for the MMS mission. Space Sci. Rev. 199, 167–188 (2016).
Lindqvist, P. A. et al. The spin-plane double probe electric field instrument for MMS. Space Sci. Rev. 199, 137–165 (2016).
Pollock, C. et al. Fast plasma investigation for magnetospheric multiscale. Space Sci. Rev. 199, 331–406 (2016).
Zenitani, S., Hesse, M., Klimas, A. & Kuznetsova, M. New measure of the dissipation region in collisionless magnetic reconnection. Phys. Rev. Lett. 106, 195003 (2011).
Pritchett, P. L. & Mozer, F. S. The magnetic field reconnection site and dissipation region. Phys. Plasmas 16(8), 080702 (2009).
Dahlin, J. T., Drake, J. F. & Swisdak, M. The role of three-dimensional transport in driving enhanced electron acceleration during magnetic reconnection. Phys. Plasmas 24(9), 092110 (2017).
Meneveau, C. Lagrangian dynamics and models of the velocity gradient tensor in turbulent flows. Annu. Rev. Fluid Mech. 43(1), 219–245 (2011).
Fu, H. S., Khotyaintsev, Y. V., Vaivads, A., Retinò, A. & André, M. Energetic electron acceleration by unsteady magnetic reconnection. Nat. Phys. 9(7), 426–430 (2013).
Acknowledgements
We thank the entire MMS team and instrument leads for data access and support. The data presented in this paper can be accessed from MMS Science Data Center (https://lasp.colorado.edu/mms/sdc/public/). The authors acknowledge AFOSR grant FA8655-22-1-7056. SCC acknowledges support from ISSI via the J. Geiss fellowship.
Funding
AFOSR grant FA8655-22-1-7056.
Author information
Authors and Affiliations
Contributions
B.H. performed the analysis. B.H. and S.C.C. designed the analysis. All authors co-wrote the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hnat, B., Chapman, S. & Watkins, N. Topology of turbulence within collisionless plasma reconnection. Sci Rep 13, 18665 (2023). https://doi.org/10.1038/s41598-023-45650-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-45650-x
- Springer Nature Limited