The role of three-dimensional transport in driving enhanced electron acceleration during magnetic reconnection

The role of three-dimensional transport in driving enhanced electron acceleration during magnetic reconnection


Magnetic reconnection is an important driver of energetic particles in many astrophysical phenomena. Using kinetic particle-in-cell (PIC) simulations, we explore the impact of three-dimensional reconnection dynamics on the efficiency of particle acceleration. In two-dimensional systems, Alfvénic outflows expel energetic electrons into flux ropes where they become trapped and disconnected from acceleration regions. However, in three-dimensional systems these flux ropes develop axial structure that enables particles to leak out and return to acceleration regions. This requires a finite guide field so that particles may move quickly along the flux rope axis. We show that greatest energetic electron production occurs when the guide field is of the same order as the reconnecting component: large enough to facilitate strong transport, but not so large as to throttle the dominant Fermi mechanism responsible for efficient electron acceleration. This suggests a natural explanation for the envelope of electron acceleration during the impulsive phase of eruptive flares.


I Introduction

Magnetic reconnection is thought to be an important driver of energetic particles in astrophysical plasmas, releasing stored magnetic energy via efficient acceleration of a non-thermal population. Reconnection-associated energetic particle production has been well-observed in solar flares Lin et al. (2003) and magnetospheric storms (Øieroset et al., 2002). Solar flare observations in particular indicate that reconnection-driven acceleration can be very efficient, driving a large non-thermal electron component with a total energy content comparable to that of the energy in the initial magnetic field Krucker et al. (2010); Oka et al. (2013). This mechanism therefore provides a promising explanation for a variety of astrophysical phenomena characterized by energetic particle production, including stellar flares, gamma-ray bursts Drenkhahn and Spruit (2002); Michel (1994), and gamma-ray flares in pulsar wind nebulae Tavani et al. (2011); Abdo et al. (2011).

Electron acceleration by magnetic reconnection has attracted significant interest, e.g. Hoshino et al. (2001); Zenitani and Hoshino (2001); Drake et al. (2005); Pritchett (2006); Egedal et al. (2009); Oka et al. (2010); Hoshino (2012). Two specific processes have received the most attention. The first is acceleration by electric fields parallel to the local magnetic field () (Litvinenko, 1996; Drake et al., 2005; Uzdensky, Cerutti, and Begelman, 2011; Egedal, Daughton, and Lê, 2012). However, the number of electrons that can be accelerated through this mechanism can be limited because during magnetic reconnection non-zero typically only occur near X-lines and separatrices. Additionally, acceleration by parallel electric fields has a weak energy-scaling Dahlin, Drake, and Swisdak (2015) ( with the particle energy) and characteristically drives bulk electron heating rather than a non-thermal component.

In the second process (Drake et al., 2006), charged particles gain energy as they reflect from the ends of contracting magnetic islands. (An analogous process occurs during the acceleration of cosmic rays by the first-order Fermi mechanism.) In contrast to the localization of , this can occur wherever there are contracting field lines, including the merging of magnetic islands and in the outflows of single X-line reconnection (Drake et al., 2006; Oka et al., 2010; Drake et al., 2010; Drake, Swisdak, and Fermo, 2013; Hoshino, 2012) and in turbulent reconnecting systems where magnetic field lines are stochastic and conventional islands do not exist Dahlin, Drake, and Swisdak (2015). This mechanism is therefore volume-filling and can accelerate a large number of particles. This mechanism scales strongly with the particle energy () and preferentially energizes non-thermal particles Dahlin, Drake, and Swisdak (2015, 2016).

Several recent studies of two-dimensional reconnection Dahlin, Drake, and Swisdak (2014); Wang et al. (2016); Dahlin, Drake, and Swisdak (2016) found that the guide field (the magnetic component parallel to the reconnection axis) controls which mechanisms contribute to electron energy gain. In the antiparallel (small guide field) regime Fermi reflection dominates Dahlin, Drake, and Swisdak (2014, 2016); Guo et al. (2014); Li et al. (2015); Wang et al. (2016), whereas in reconnection where the guide field is much larger than the reconnecting component parallel electric fields drive essentially all of the electron energy gain Numata and Loureiro (2015); Dahlin, Drake, and Swisdak (2016). In the latter (strong guide field) regime energetic electron production is weak, indicating that parallel electric fields are ineffective drivers of energetic electrons in reconnection Dahlin, Drake, and Swisdak (2016).

Studies of particle acceleration in reconnection have primarily been based on 2D simulations, in which accelerated particles are typically localized near the X-line, along magnetic separatrices and within magnetic islands (Drake et al., 2005; Guo et al., 2014). There are some observations with small ambient guide fields (Chen et al., 2008; Retinò et al., 2008; Huang et al., 2012) that support such a picture. A notable exception are Wind observations in which energetic electrons up to 300 keV are seen for more than an hour in an extended region around the reconnection region (Øieroset et al., 2002). These observations correspond to reconnection with a strong guide field.

Two-dimensional simulations impose limitations on the magnetic topology as well as the available spectrum of instabilities. In the presence of an ambient guide field, 3D reconnection can become turbulent as a result of the generation of magnetic islands along separatrices and adjacent surfaces Schreier et al. (2010); Daughton et al. (2011). While test particle trajectories in MHD fields have been used to explore acceleration in such systems(Onofri, Isliker, and Vlahos, 2006; Kowal, de Gouveia Dal Pino, and Lazarian, 2011), the absence of feedback of energetic particles on the reconnection process in such models limits their applicability to physical systems. Recent 3D studies of kinetic reconnection examined particle acceleration in pair plasmas (Sironi and Spitkovsky, 2014; Guo et al., 2014). However, these studies focused on relativistic regimes where the magnetic energy per particle exceeds the rest mass energy and included no ambient guide field.

In a recent kinetic study of nonrelativistic reconnection, we showed that energetic electron production was greatly enhanced in three-dimensional systems Dahlin, Drake, and Swisdak (2015). This occurs because two-dimensional magnetic islands trap particles, limiting energy gain, whereas three-dimensional reconnection generates a stochastic field that enables electrons to access volume-filling acceleration regions. The relative enhancement was found to increase with the size of the simulation domain, suggesting that that the impact of three-dimensional dynamics is robust for astrophysical characterized by spatial and temporal scales that are much larger than kinetic scales.

In this article, we extend this study in several key ways. We begin by reviewing the theory of particle acceleration in reconnection (section II) and describing our kinetic particle-in-cell (PIC) simulations (section III). We then explore the physics of two and three-dimensional reconnection to highlight the remarkable similarity of many of the bulk properties (section IV). In section V, we review the physics of electron acceleration enhancement in 3D reconnection and demonstrate, by varying the spatial length in the third dimension, that transport and enhanced acceleration are intrinsically linked. We show in section VI that a magnetic guide field plays an important role in facilitating three-dimensional transport, and in section VII introduce an ‘injection criterion’ that explains why energetic electrons are enhanced but protons are not. We discuss the astrophysical implications of these results in section VIII.

Ii Particle Acceleration in the Guiding-Center Limit

In order to examine electron acceleration we assume a guiding-center approximation relevant for a strong guide field (Northrop, 1963; Dahlin, Drake, and Swisdak, 2015). In this limit, the evolution of the kinetic energy of a single electron can be written as:


where is the parallel electric field, is the magnetic moment, , and is the magnetic curvature. The velocity components parallel and perpendicular to the magnetic field are and , respectively; is the relativistic Lorentz factor, and is the unit vector in the direction of the local magnetic field.

The first term on the right-hand-side of the equation corresponds to acceleration by parallel electric fields, which are typically localized near the reconnection X-line and along separatrices. The second term corresponds to betatron acceleration associated with conservation in a temporally and spatially varying magnetic field. Because reconnection releases a system’s magnetic energy, this typically causes electron cooling Dahlin, Drake, and Swisdak (2014). The last term corresponds to Fermi reflection of particles from contracting magnetic field lines (Drake et al., 2006, 2010; Hoshino, 2012; Dahlin, Drake, and Swisdak, 2014). Both and Fermi reflection change the parallel energy of the particles, while betatron acceleration changes the perpendicular energy. The term corresponds to local field line contraction: (where is the field line length) and is linked to the conservation of the parallel adiabatic invariant . The guiding-center approximation given in Eq. (1) is accurate when electrons are well-magnetized. In the weak-guide field regime, other terms such as the polarization drift may be significant (compare Li et al.Li et al. (2015)). However, the polarization drift gives the change in the electron bulk flow energy which is typically small for a realistic electron-to-ion mass ratio.

Iii Particle-in-Cell Simulations

We explore particle acceleration in reconnection via simulations using the massively parallel 3D particle-in-cell (PIC) code p3d (Zeiler et al., 2002). Particle trajectories are calculated using the relativistic Newton-Lorentz equation, and the electromagnetic fields are advanced using Maxwell’s equations. The time and space coordinates are normalized, respectively, to the proton cyclotron time and inertial length . The typical grid cell width , where is the electron inertial length. The time step is , where is the electron cyclotron frequency.

All simulations are initialized with a force-free configuration and use periodic boundary conditions. This is chosen as the most generic model for large-scale systems such as the solar corona where the density jump between the current layer and upstream plasma is not expected to be important. The magnetic field is given by and , corresponding to an asymptotic guide field . We include two current sheets at and to produce a periodic system, and choose . This initial configuration is not a kinetic equilibrium, which would require a temperature anisotropy Bobrova et al. (2001), but is in pressure balance. We use at least 50 particles per cell per species. The initial electron and proton temperatures are equal and isotropic with , and the initial density and pressure are constant so that . The speed of light is , where is the Alfvén speed based on the reconnecting component of the magnetic field.

Table 1 lists the simulation configurations discussed in this paper. We focus on configurations SM (‘medium’) and SL (‘large’) with spatial dimensions and , respectively. The larger simulations (SL) have more magnetic flux to reconnect, and can therefore run for a longer time and generate many more energetic electrons. However, the SM configuration is times less expensive and therefore better suited for parameter scans. Except where noted, the simulations are performed with a guide field and the three-dimensional simulations use (as in Dahlin et al., 2015Dahlin, Drake, and Swisdak (2015)). The simulation with in configuration SM uses . The electron-positron (pair plasma) configuration (S1) uses and .

The PIC formulation results in some numerical heating that can differ between 2D and 3D domains. To simplify comparisons of particle energization, we perform a set of simulations with so that the numerical heating is the same as in 3D, yet the reconnection physics remains essentially 2D (see Section V and Fig. 7). To reduce computational expense, configurations SM and SL have an artificial proton-to-electron mass ratio . Simulations with (S100) and (S1) are also presented to explore the impact of the separation between electron and proton scales.

SM 51.2 25.6 2D, 1.6 - 25.6 25 0 - 1.5
SL 102.4 51.2 2D, 1.6, 25.6 25 0, 0.5, 1.0
S1 102.4 51.2 1.6, 51.2 1 1.0
S100 51.2 25.6 0.8, 25.6 100 1.0
Table 1: Simulation Parameters.

Iv Overview of reconnection in 2D and 3D domains

Reconnection develops from particle noise via the tearing instability, which generates interacting flux ropes that grow and merge until they reach the system size. These tearing modes grow wherever (Ref. Furth, Killeen, and Rosenbluth, 1963). In a slab equilibrium with , such modes are characterized by a wavevector and grow on flux surfaces defined by . The pitch of the unstable mode is given by a characteristic angle with respect to the reconnecting field . In a two-dimensional system where , tearing modes are constrained to grow at the center of the current sheet where . However, in three-dimensional systems the nonlinear interaction of modes with different pitches destroys flux surfaces and generates a stochastic, turbulent magnetic field that facilitates enhanced particle transportRechester and Rosenbluth (1978); Daughton et al. (2011).

This stochastic magnetic structure is illustrated in Figs. 1a-b, which show isosurfaces of one component of the electron current density in configuration SM. At , several tearing structures with are visible (compare Fig. 1 in Daughton et al. 2011Daughton et al. (2011)). The filamentary current distribution at showcases the late-time nonlinear development. A different view of the filamentary structure, which emphasizes the stochastic structure of the magnetic field, can be observed in the Poincaré surface-of-section shown in Figure 1c. There is a clear boundary between the stochastic reconnecting region (disordered punctures) and the asymptotic, laminar field.

The stochastic 3D dynamics do not substantially impact the magnetic energy release (see Fig. 2a), as has been noted by Daughton et al.Daughton et al. (2014) Another diagnostic for the energy release is the field-line contraction that drives Fermi acceleration according to Eq. (1). The spatial average , ( calculated over the stochastic reconnection region described below) is shown in Fig. 2b. Although the 2D simulation is relatively bursty, the overall time-evolution is comparable to that in 3D. In both simulations, decreases in time as islands grow and the typical radius of curvature increases. Figure 2c shows the probability distribution function (pdf) of inside the stochastic reconnecting region at . The pdf in the 3D system is symmetric for small values of , and is consistent with a double-exponential distribution where is the Alfvén crossing time. The symmetric component can only produce net acceleration through a second-order Fermi process. The positive mean value is due to the large positive tail visible for . Although the characteristic scales for first and second-order components are comparable, the first-order mechanism is far more efficient, and hence is the dominant driver of particle acceleration in this system. The symmetric (second-order Fermi) component is consistent with Alfvénic fluctuations where the flow and curvature are out of phase corresponding to no net field-line contraction, i.e. .

It has been shown previously that the development of pressure anisotropy with causes the cores of magnetic islands to approach firehose marginal stability, where the tension driving magnetic reconnection ceases, thereby throttling reconnection Drake et al. (2006, 2010); Schoeffler, Drake, and Swisdak (2011). Figure 3a,c show that a significant anisotropy develops, as is the typical case in 2D (not shown). This suggests that the plasma heating and energization occurs in similar ways in 2D and 3D, and that the turbulent magnetic field generated in 3D does not isotropize the plasma. Phase space plots of temperature anisotropy and are shown in Fig. 4 for three values of the guide field in 2D and 3D, along with marginal stability boundaries for the firehose and mirror instabilities (bottom and top, respectively) in a similar format to that used previously in analyzing solar wind data Kasper, Lazarus, and Gary (2002); Bale et al. (2009). For the strong guide field cases (4e,f) there is insufficient free energy in the reconnecting field to drive the system either firehose or mirror unstable. However, for the simulations with (Fig. 4a-d) the phase space brushes up against the stability boundaries. Hence, feedback of heating on reconnection via anisotropy-driven instabilities can occur in 3D systems.

Reconnection primarily drives electron acceleration parallel to the local magnetic field, generating superthermal particles that fill the stochastic reconnecting regionDahlin, Drake, and Swisdak (2015). The presence of reconnection-accelerated electrons is therefore a useful proxy for the reconnection region. This is similar to the “electron mixing” described by Daughton et al. (Ref. Daughton et al., 2014). Indeed, the volume defined by (where is the parallel energy density of electrons exceeding ) corresponds well to the reconnecting region as indicated by the electron pressure and current density (see Fig. 3). Using this marker for designating the reconnection domain allows us to estimate the reconnected volume and can be used to determine a characteristic width in the 3D system: . The 2D analogue is the area inside the separatrices of the primary X-line. A mean inflow velocity can then be determined from , yielding comparable for the 3D and 2D simulations respectively ( is averaged from to ). A calculation of the 2D reconnection rate, determined from the time rate of change of the flux function at the primary X-line, yields a nearly identical inflow velocity . Figure 2d shows that the width of the reconnecting region is comparable to the mean radius of curvature, consistent with . In summary, the total magnetic energy conversion, field line contraction (Fermi drive), and reconnection rate provide strong evidence for the remarkable simlarity of bulk properties of the 2D and 3D reconnection.

V Electron Acceleration

In a previous study Dahlin, Drake, and Swisdak (2015) we found that electron acceleration was enhanced in a 3D reconnecting system. These results are summarized in Fig. 5. Although there is substantial acceleration in both systems, the fraction of electrons with energy exceeding is more than an order of magnitude larger than in the 2D simulation (Fig. 5a). However, as noted in Section IV, the magnetic energy dissipation is comparable in 2D and 3D systems. This suggests that the increased energetic electron production in the 3D system is due to enhanced acceleration efficiency rather than an increase in the total energy imparted to the plasma. According to equation (1) the acceleration mechanisms have different scalings with the particle energy: the Fermi reflection term is second-order in the parallel velocity, whereas the parallel electric field term is only first-order. The instantaneous average acceleration per particle for both and Fermi reflection in configuration SM, () is shown in Fig. 5b. The bulk thermal electrons (low energies) are primarily accelerated by , whereas Fermi reflection is more important at high energies. The energetic electrons are primarily accelerated in the parallel direction so that the momentum distribution exceeds (Fig. 5c), consistent with acceleration via Fermi reflection and . To summarize: electron acceleration, primarily driven by field-line contraction, is enhanced in 3D systems. However, this is not due to greater energy release, so must instead be due to enhancement of the acceleration efficiency.

As was discussed in Section IV, the stochastic structure of the magnetic field in 3D systems allows field-line-following particles to wander throughout the chaotic reconnecting region Jokipii and Parker (1969). However, in 2D systems reconnected field lines form closed loops (islands) that trap particles. The impact of topology on transport is reflected in the spatial distribution of the most energetic particles (shown in Fig. 6a,b). These particles occupy narrow bands well inside the islands in the 2D simulation, but are distributed throughout the reconnecting region in the 3D simulation. The most efficient electron acceleration regions are near the X-lines and at the ends of islands (Fig. 6c,d). In the 2D system, trapped energetic practicles are unable to access these regions, and the overall acceleration efficiency is suppressed with respect to the 3D system where the energetic particles wander the reconnecting region and undergo continuous acceleration. Figure 6e,f shows the resulting rate of Fermi energization for electrons with , revealing an order-of-magnitude difference between the two systems.

In order to examine the transition from 2D to 3D reconnection, we performed a set of simulations with different values of , from a quasi-2D system with to a simulation with . Figure 7a-c shows the spatial distribution of the energetic electrons () for several of these simulations. Surprisingly, there is a sharp transition at : below this threshold, energetic electrons are trapped inside islands, whereas above this threshold the energetic electrons are space-filling. Electron energy spectra exhibit the same transition (Fig. 7d). Simulations with do not show enhancement with respect to the 2D result, whereas simulations are consistent with the 3D result. This reinforces the correlation between enhanced transport and acceleration in 3D systems.

Vi The Role of the Guide Field

In a recent study of two-dimensional kinetic reconnection (Dahlin et al., 2016 Dahlin, Drake, and Swisdak (2016)), we found that the magnetic guide field was a vital parameter that controls the efficiency of electron acceleration. In a system with a guide field much smaller than the reconnecting component, the dominant electron accelerator was a Fermi-type mechanism that preferentially energizes the most energetic particles. In the strong guide field regime, however, the field-line contraction that drives Fermi reflection was negligible. Instead, parallel electric fields () were primarily responsible for driving electron heating. Electron acceleration was suppressed in the systems with a strong guide field. We argued that this was due to the the weaker energy scaling of acceleration.

To probe the role of the guide field in three-dimensional transport and particle acceleration, we performed several three-dimensional simulations in the configuration SM with and compared these results with a set of quasi-2D simulations (). Selected electron energy spectra from these simulations are shown in Figure 8a-c. Figure 8d shows the number of electrons exceeding , revealing that energetic electron production varies strongly with the guide field. The efficiency of the Fermi mechanism that drives energetic electron production weakens with increasing guide fieldDahlin, Drake, and Swisdak (2016). The quasi-2D spectra (Fig. 8c) are consistent with this result, explaining why the energetic electron production diminishes as . The decreasing energetic electron production for must then be due to three-dimensional effects. Indeed, Fig. 8b shows that the relative 3D enhancement () increases with until it saturates above (there is little difference between ).

Representative field lines for several of the 3D simulations are shown in Fig. 10, illustrating the differing field line structures (Fig. 10a shows field lines from a 2D simulation). In Fig. 10b, the field lines are clearly stochastic, and do not show clear flux rope structures that would trap particles; hence, electrons are free to return to acceleration regions. In weak guide-field reconnection, however (Figs. 10c-d), field lines wrap around clearly-defined flux ropes, and in the case with , the field lines approximately close on themselves. This structure inhibits particle escape from islands, similar to the 2D structure shown in Fig. 10a, where flux surfaces are closed and particles become topologically disconnected from acceleration regions. In 3D, the guide field plays a role in breaking the 2D symmetry and allowing particles to escape along the flux rope axis. This explains why the three-dimensional enhancement increases with the guide field. The saturation above can be explained by noting that magnetic structures are typically elongated along the guide field for so that particles must move farther along the axis to escape the island. Hence, further 3D enhancement over what is shown in Fig. 8b should not occur.

The convolution of the Fermi acceleration efficiency and the effectiveness of three-dimensional transport results in a peak energetic electron production at . Results from a set of 3D simulations with (configuration SL) are shown in Fig. 9 (dashed lines indicate an earlier time). These simulations show that the enhancement increases as the spectra extend to higher energies, suggesting that three-dimensional transport will be even more important in physical systems such as the corona where the length scales . The most efficient guide field, in these simulations , will likely depend both on the system size and on other plasma parameters such as the plasma beta, which can impact the relative efficiency of Fermi and -driven acceleration.

Vii An ‘Injection Criterion’ for Enhanced Acceleration

A limitation of the present simulations is the use of an artificial mass ratio, which reduces the separation between proton and electron scales. To examine how the mass ratio impacts particle acceleration, we performed simulations with (configurations S1, SM, and S100) and . Figure 11 shows the relative enhancement of the energy spectra in the three-dimensional simulations (). For the electron-positron case (), there is only a slight enhancement () in the energetic tail for both species. For the electron-proton cases (), the energetic electrons are enhanced whereas the energetic ions are suppressed. The enhancement (suppression) of the energetic electrons (ions) is greater for the more realistic mass ratio. This trend, along with the weak enhancement for the electron-positron case, suggesting that the separation of scales between species plays an important role in 3D particle acceleration and that the impact of 3D transport should be robust for the physical mass ratio.

We propose the following explanation: in order for a charged particle to accelerate multiple times, it must propagate upstream against the Alfvénic outflow that ejects plasma from the energy release regions near the X-line and at the ends of islands. The condition then acts as an ‘injection criterion’ for efficient acceleration (this is analogous to the injection problem in shock-driven particle acceleration). Heavy species (protons for , and both electrons and positrons for ) are responsible for the bulk inertia and hence the reconnection outflow velocity . The characteristic velocity of these particles is therefore of the same order as the Alfvén speed (), so that the bulk particles do not meet the injection criterion and do not experience enhanced 3D acceleration. A few particles in the tail of the distribution satisfy the criterion in the electron-positron case, explaining the small () enhancement at high energies. For the electron-proton simulations () the electron thermal velocity greatly exceeds the Alfvén speed. Recent numerical studies suggest that even in environments with , reconnection heats electrons to an appreciable fraction of the available magnetic energy density Shay et al. (2014); Haggerty et al. (2015): , corresponding to so that essentially all reconnection-heated electrons will satisfy the criterion, independent of the initial temperature. In contrast, ions are typically sub-Alfénic and would require an injection mechanism (e.g. Drake et al. (2009a, b); Knizhnik, Swisdak, and Drake (2011)) to undergo continuous acceleration. However, the suppression of ion acceleration in 3D systems is surprising. While the relative increase in energy going into energetic electrons may play some role, it is not clear that this should preferentially impact the energetic ions. Further treatment is beyond the scope of this paper.

The injection criterion may require modification in the large guide field limit (). In this limit, the electrons dominantly stream in the -direction, so the relevant velocity to compare to the outflow speed is is so the injection criterion becomes . In strongly relativistic systems, all velocities approach , so that the injection criterion cannot be met. This suggests that enhanced 3D acceleration should not occur for relativistic reconnection in either pair Guo et al. (2014); Sironi and Spitkovsky (2014) or electron-proton plasmas.

Viii Discussion

Electron acceleration in three-dimensional systems is a complex problem that intrinsically depends on the transport properties of reconnection-generated magnetic fields. The picture that emerges from this set of simulations is that particle acceleration is efficient in a three-dimensional system when the energetic population can freely access acceleration sites and thereby achieve continuous energy gain. This requires both topological access to energy release regions and a super-Alfvénic particle velocity in order to explore the open topology at a faster time scale than the system evolves (most easily understood as the ejection of flux from the X-line). While electrons satisfy this condition rather easily, heavy species such as protons would require an injection mechanism in order to be able to propagate upstream against the reconnection outflow.

Efficient transport requires a strong guide field. The field structure in antiparallel reconnection is quasi-laminar, so that energetic particles are still well-trapped in islands. Propagation upstream against the Alfvén velocity is not possible in strongly relativistic reconnection, where all characteristic velocities approach the speed of light. This is consistent with studies by Guo et al.Guo et al. (2014, 2015, 2016), that exhibit no substantial difference in energetic particle production between two and three-dimensional simulations in the relativistic regime.

The nonthermal electron spectra in both simulations do not assume a power law form as is frequently observed in nature. The maximum energy gain is limited due to the modest size of the simulations; previous 2D simulations have shown that the total energy gain is greater in larger systems Dahlin, Drake, and Swisdak (2014). An additional issue is that these simulations have periodic boundary conditions that prevent particle loss from the system. Solar observations suggest that electrons are confined in regions of energy release in the corona Krucker et al. (2010); possible mechanisms for this confinement include mirroring and double layersLi, Drake, and Swisdak (2012). This could suggest that particle loss is not an important concern. On the other hand, it has been suggested that the development of a power law requires a loss mechanism in addition to an energy drive Drake, Swisdak, and Fermo (2013). However, several recent simulations suggest that power-law spectra may still develop in the absence of a loss mechanism Guo et al. (2014); Sironi and Spitkovsky (2014); Li et al. (2015); Werner et al. (2016). The set of conditions under which power-law spectra form in kinetic reconnection simulations remains an open issue.

The simulations have a number of numerical limitations. These include comparatively small spatial and temporal scales, and an artificial electron-to-proton mass ratio. However, as was discussed in Section VII, the greater characteristic velocity of a realistically ‘light’ electron facilitates transport in the stochastic topology; the efficiency of Fermi acceleration does not directly dependent on the particle mass. The largest simulations show that three-dimensional dynamics are increasingly important at larger scales; in contrast the diminishing frequency of island mergers leads to suppression of further acceleration in two-dimensional systems.

The role of the guide field in magnetic reconnection has broad implications for reconnection-driven particle acceleration in astrophysical systems. Electron acceleration is most efficient in the regime where both (a) the Fermi mechanism operates and (b) strong three-dimensional transport exists. The former requires , the latter , suggesting that reconnection with a magnetic field of the same order as the reconnecting component () will yield the most efficient energetic electron production. This result is especially relevant for solar flares, where the shear in the magnetic configuration typically diminishes during the impulsive phase (see Fletcher et al. Fletcher et al. (2011) and references therein). This corresponds to a transition from strong to weak guide field reconnection, and could explain why electron acceleration (as inferred from hard X-ray emissions) is typically confined to the impulsive phase.

This work has been supported by NSF Grants AGS1202330 and PHY1102479, and NASA grants NNX11AQ93H, APL-975268, NNX08AV87G, NAS 5- 98033, and NNX08AO83G. J.T.D. acknowledges support from the NASA LWS Jack Eddy Fellowship administered by the University Corporation for Atmospheric Research in Boulder, Colorado. Simulations were carried out at the National Energy Research Scientific Computing Center.
Figure 1: Results for configuration SM. (a,b) Isosurfaces of at respectively, illustrating the nonlinear filamentary current structure. (c) Poincaré surface-of-section at . We trace a set of field lines beginning at , , and plot where they puncture the plane . The surface-of-section shows a clear boundary between the stochastic field lines inside the reconnecting region and the laminar unreconnected fields.
Figure 2: Results for configuration SL: 3D (green) and 2D (red). (a) Magnetic energy release vs. time. (b) Spatially averaged over the reconnecting region, where is the Alfvén crossing time. (c) Probability distribution function of at , with a superimposed double-exponential fit (blue dotted line). (d) Reconnection region half-width and radius of curvature .
Figure 3: Results for configuration SL. (a)(c) slices at of parallel and perpendicular electron pressure. (b) (d) -directed electron current density in planes defined by and , respectively. White and black outlines indicate the contour , a marker for the region of reconnected magnetic field.
Figure 4: Results for configuration SM. Phase space distribution of temperature ratio vs. . Black lines indicate marginal stability conditions for the ideal firehose (bottom) and mirror (top) instabilities.
Figure 5: Results for configuration SL: 3D (green) and 2D (red). (a) Energetic spectra at (dashed) and (solid) (b) Average energy gain due to Fermi acceleration (red) and (blue) in simulation 4a at . (c) Parallel and perpendicular momentum spectra at .
Figure 6: Results for configuration SL: 3D (left) and 2D (right) at in the plane . (a-b) Parallel energy density for electrons with . (c-d) Field-line contraction rate . (e-f) Fermi acceleration for electrons with . The similarity between (c) and (e) reflects the nearly uniform spatial distribution of the energetic electrons in the reconnecting region. Dashed boxes outline one X-line in each panel.
Figure 7: Results for configuration SM, . (a-c) Spatial distribution of electrons with . (d) Electron energy spectra normalized to the 2D spectrum at .
Figure 8: Results for configuration SM at . Three-dimensional electron energy spectra normalized to the initial spectrum (a) and quasi-2D spectra (b). (c) quasi-2D spectra normalized to the initial spectrum . (d) Electrons exceeding vs. guide field. The system with generates the greatest number of energetic electrons.
Figure 9: Results for configuration SL. Three-dimensional electron energy spectra normalized to the initial spectrum (a) and quasi-2D spectra (b) at (solid) and (dashed). (c) quasi-2D spectra normalized to the initial () spectrum. (d) Electrons exceeding at vs. guide field. The system with generates the greatest number of energetic electrons.
Figure 10: Results from configuration SM. Single field lines (blue) for simulations with different guide fields. Two-dimensional () slices of the energetic electron density are shown at the bottom of each panel.
Figure 11: Results from configurations SM, S1, and S100. Three-dimensional enhancement ) for simulations with different mass ratios . For the pair plasma (), both species exhibit a small enhancement , whereas for simulations with , energetic electrons are enhanced whereas energetic ions are suppressed.


  1. preprint: AIP/123-QED


  1. R. P. Lin, S. Krucker, G. J. Hurford, D. M. Smith, H. S. Hudson, G. D. Holman, R. A. Schwartz, B. R. Dennis, G. H. Share, R. J. Murphy, A. G. Emslie, C. Johns-Krull,  and N. Vilmer, Ap. J. 595, L69 (2003).
  2. M. Øieroset, R. P. Lin, T. D. Phan, D. E. Larson,  and S. D. Bale, Phys. Rev. Lett. 89, 195001 (2002).
  3. S. Krucker, H. S. Hudson, L. Glesener, S. M. White, S. Masuda, J.-P. Wuelser,  and R. P. Lin, Ap. J. 714, 1108 (2010).
  4. M. Oka, S. Ishikawa, P. Saint-Hilaire, S. Krucker,  and R. P. Lin, The Astrophysical Journal 764, 6 (2013).
  5. G. Drenkhahn and H. C. Spruit, Astronomy & Astrophysics 391, 1141 (2002).
  6. F. C. Michel, Ap. J. 431, 397 (1994).
  7. M. Tavani, A. Bulgarelli, V. Vittorini, A. Pellizzoni, E. Striani, P. Caraveo, M. C. Weisskopf, A. Tennant, G. Pucella, A. Trois, E. Costa, Y. Evangelista, C. Pittori, F. Verrecchia, E. Del Monte, R. Campana, M. Pilia, A. De Luca, I. Donnarumma, D. Horns, C. Ferrigno, C. O. Heinke, M. Trifoglio, F. Gianotti, S. Vercellone, A. Argan, G. Barbiellini, P. W. Cattaneo, A. W. Chen, T. Contessi, F. D’Ammando, G. DeParis, G. Di Cocco, G. Di Persio, M. Feroci, A. Ferrari, M. Galli, A. Giuliani, M. Giusti, C. Labanti, I. Lapshov, F. Lazzarotto, P. Lipari, F. Longo, F. Fuschino, M. Marisaldi, S. Mereghetti, E. Morelli, E. Moretti, A. Morselli, L. Pacciani, F. Perotti, G. Piano, P. Picozza, M. Prest, M. Rapisarda, A. Rappoldi, A. Rubini, S. Sabatini, P. Soffitta, E. Vallazza, A. Zambra, D. Zanello, F. Lucarelli, P. Santolamazza, P. Giommi, L. Salotti,  and G. F. Bignami, Science 331, 736 (2011).
  8. A. A. Abdo, M. Ackermann, M. Ajello, A. Allafort, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, B. Berenji, R. D. Blandford, E. D. Bloom, E. Bonamente, A. W. Borgland, A. Bouvier, T. J. Brandt, J. Bregeon, A. Brez, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, A. Cannon, P. A. Caraveo, J. M. Casandjian, Ã. Ã‡elik, E. Charles, A. Chekhtman, C. C. Cheung, J. Chiang, S. Ciprini, R. Claus, J. Cohen-Tanugi, L. Costamante, S. Cutini, F. D’Ammando, C. D. Dermer, A. de Angelis, A. de Luca, F. de Palma, S. W. Digel, E. do Couto e Silva, P. S. Drell, A. Drlica-Wagner, R. Dubois, D. Dumora, C. Favuzzi, S. J. Fegan, E. C. Ferrara, W. B. Focke, P. Fortin, M. Frailis, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Gehrels, S. Germani, N. Giglietto, F. Giordano, M. Giroletti, T. Glanzman, G. Godfrey, I. A. Grenier, M.-H. Grondin, J. E. Grove, S. Guiriec, D. Hadasch, Y. Hanabata, A. K. Harding, K. Hayashi, M. Hayashida, E. Hays, D. Horan, R. Itoh, G. Jóhannesson, A. S. Johnson, T. J. Johnson, D. Khangulyan, T. Kamae, H. Katagiri, J. Kataoka, M. Kerr, J. Knödlseder, M. Kuss, J. Lande, L. Latronico, S.-H. Lee, M. Lemoine-Goumard, F. Longo, F. Loparco, P. Lubrano, G. M. Madejski, A. Makeev, M. Marelli, M. N. Mazziotta, J. E. McEnery, P. F. Michelson, W. Mitthumsiri, T. Mizuno, A. A. Moiseev, C. Monte, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, T. Nakamori, M. Naumann-Godo, P. L. Nolan, J. P. Norris, E. Nuss, T. Ohsugi, A. Okumura, N. Omodei, J. F. Ormes, M. Ozaki, D. Paneque, D. Parent, V. Pelassa, M. Pepe, M. Pesce-Rollins, M. Pierbattista, F. Piron, T. A. Porter, S. Rainò, R. Rando, P. S. Ray, M. Razzano, A. Reimer, O. Reimer, T. Reposeur, S. Ritz, R. W. Romani, H. F.-W. Sadrozinski, D. Sanchez, P. M. S. Parkinson, J. D. Scargle, T. L. Schalk, C. Sgrò, E. J. Siskind, P. D. Smith, G. Spandre, P. Spinelli, M. S. Strickman, D. J. Suson, H. Takahashi, T. Takahashi, T. Tanaka, J. B. Thayer, D. J. Thompson, L. Tibaldo, D. F. Torres, G. Tosti, A. Tramacere, E. Troja, Y. Uchiyama, J. Vandenbroucke, V. Vasileiou, G. Vianello, V. Vitale, P. Wang, K. S. Wood, Z. Yang,  and M. Ziegler, Science 331, 739 (2011).
  9. M. Hoshino, T. Mukai, T. Terasawa,  and I. Shinohara, J. Geophys. Res. 106, 25,979 (2001).
  10. S. Zenitani and M. Hoshino, Ap. J. Lett. 562, L63 (2001).
  11. J. F. Drake, M. A. Shay, W. Thongthai,  and M. Swisdak, Phys. Rev. Lett. 94, 095001 (2005).
  12. P. L. Pritchett, J. Geophys. Res. 111, A10212 (2006), 10.1029/2006JA011793.
  13. J. Egedal, W. Daughton, J. F. Drake, N. Katz,  and A. Lê, Phys. Plasmas 16, 050701 (2009), 10.1063/1.3130732.
  14. M. Oka, T.-D. Phan, S. Krucker, M. Fujimoto,  and I. Shinohara, Ap. J. 714, 915 (2010).
  15. M. Hoshino, Phys. Rev. Lett. 108, 135003 (2012).
  16. Y. E. Litvinenko, Ap. J. 462, 997 (1996).
  17. D. A. Uzdensky, B. Cerutti,  and M. C. Begelman, The Astrophysical Journal Letters 737, L40 (2011).
  18. J. Egedal, W. Daughton,  and A. Lê, Nature Phys. 8, 321 (2012).
  19. J. T. Dahlin, J. F. Drake,  and M. Swisdak, Physics of Plasmas 22, 100704 (2015),
  20. J. F. Drake, M. Swisdak, H. Che,  and M. A. Shay, Nature 443, 553 (2006).
  21. J. F. Drake, M. Opher, M. Swisdak,  and J. N. Chamoun, Ap. J. 709, 963 (2010).
  22. J. F. Drake, M. Swisdak,  and R. Fermo, The Astrophysical Journal Letters 763, L5 (2013).
  23. J. T. Dahlin, J. F. Drake,  and M. Swisdak, Phys. Plasmas  (2016), submitted, arXiv:1607.03857.
  24. J. T. Dahlin, J. F. Drake,  and M. Swisdak, Phys. Plasmas 21, 092304 (2014).
  25. H. Wang, Q. Lu, C. Huang,  and S. Wang, The Astrophysical Journal 821, 84 (2016).
  26. F. Guo, H. Li, W. Daughton,  and Y.-H. Liu, Phys. Rev. Lett. 113, 155005 (2014).
  27. X. Li, F. Guo, H. Li,  and G. Li, The Astrophysical Journal Letters 811, L24 (2015).
  28. R. Numata and N. F. Loureiro, Journal of Plasma Physics 81, 305810201 (17 pages) (2015).
  29. L.-J. Chen, A. Bhattacharjee, P. A. Puhl-Quinn, H. Yang, N. Bessho, S. Imada, S. Mühlbachler, P. W. Daly, B. Lefebvre, Y. Khotyaintsev, A. Vaivads, A. Fazakerley,  and E. Georgescu, Nature Physics 4, 19 (2008).
  30. A. Retinò, R. Nakamura, A. Vaivads, Y. Khotyaintsev, T. Hayakawa, K. Tanaka, S. Kasahara, M. Fujimoto, I. Shinohara, J. P. Eastwood, M. André, W. Baumjohann, P. W. Daly, E. A. Kronberg,  and N. Cornilleau-Wehrlin, Journal of Geophysical Research: Space Physics 113 (2008), 10.1029/2008JA013511, a12215.
  31. S. Y. Huang, A. Vaivads, Y. V. Khotyaintsev, M. Zhou, H. S. Fu, A. Retinò, X. H. Deng, M. André, C. M. Cully, J. S. He, F. Sahraoui, Z. G. Yuan,  and Y. Pang, Geophysical Research Letters 39 (2012), 10.1029/2012GL051946, l11103.
  32. R. Schreier, M. Swisdak, J. F. Drake,  and P. A. Cassak, Phys. Plasmas 17, 110704 (2010), 10.1063/1.3494218.
  33. W. Daughton, V. Roytershteyn, H. Karimabadi, L. Yin, B. J. Albright, B. Bergen,  and K. J. Bowers, Nature Phys. 7, 539 (2011).
  34. M. Onofri, H. Isliker,  and L. Vlahos, Phys. Rev. Lett. 96, 151102 (2006), 10.1103/PhysRevLett96.151102.
  35. G. Kowal, E. M. de Gouveia Dal Pino,  and A. Lazarian, The Astrophysical Journal 735, 102 (2011).
  36. L. Sironi and A. Spitkovsky, The Astrophysical Journal Letters 783, L21 (2014).
  37. T. G. Northrop, Reviews of Geophysics 1, 283 (1963).
  38. A. Zeiler, D. Biskamp, J. F. Drake, B. N. Rogers, M. A. Shay,  and M. Scholer, J. Geophys. Res. 107, 1230 (2002).
  39. N. A. Bobrova, S. V. Bulanov, J. I. Sakai,  and D. Sugiyama, Physics of Plasmas 8 (2001).
  40. H. P. Furth, J. Killeen,  and M. N. Rosenbluth, Physics of Fluids 6 (1963).
  41. A. B. Rechester and M. N. Rosenbluth, Phys. Rev. Lett. 40, 38 (1978).
  42. W. Daughton, T. K. M. Nakamura, H. Karimabadi, V. Roytershteyn,  and B. Loring, Physics of Plasmas 21, 052307 (2014),
  43. K. M. Schoeffler, J. F. Drake,  and M. Swisdak, Ap. J. 743, 70 (2011).
  44. J. C. Kasper, A. J. Lazarus,  and S. P. Gary, Geophysical Research Letters 29, 20 (2002), 1839.
  45. S. D. Bale, J. C. Kasper, G. G. Howes, E. Quataert, C. Salem,  and D. Sundkvist, Phys. Rev. Lett. 103, 211101 (2009), 10.1103/PhysRevLett.103.211101.
  46. J. R. Jokipii and E. N. Parker, The Astrophysical Journal 155, 777 (1969).
  47. M. A. Shay, C. C. Haggerty, T. D. Phan, J. F. Drake, P. A. Cassak, P. Wu, M. Oieroset, M. Swisdak,  and K. Malakit, Phys. Plasmas 21, 122902 (2014).
  48. C. C. Haggerty, M. A. Shay, J. F. Drake, T. D. Phan,  and C. T. McHugh, Geophysical Research Letters 42, 9657 (2015), 2015GL065961.
  49. J. F. Drake, M. Swisdak, T. D. Phan, P. A. Cassak, M. A. Shay, S. T. Lepri, R. P. Lon, E. Quataert,  and T. H. Zurbuchen, J. Geophys. Res. 114, A05111 (2009a), 10.1029/2008JA013701.
  50. J. F. Drake, P. A. Cassak, M. A. Shay, M. Swisdak,  and E. Quataert, Ap. J. 700, L16 (2009b).
  51. K. Knizhnik, M. Swisdak,  and J. Drake, Ap. J. Lett. 743 (2011), 10.1088/2041-8205/743/2/L35.
  52. F. Guo, Y.-H. Liu, W. Daughton,  and H. Li, The Astrophysical Journal 806, 167 (2015).
  53. F. Guo, H. Li, W. Daughton, X. Li,  and Y.-H. Liu, Physics of Plasmas 23, 055708 (2016),
  54. T. C. Li, J. F. Drake,  and M. Swisdak, The Astrophysical Journal 757, 20 (2012).
  55. G. R. Werner, D. A. Uzdensky, B. Cerutti, K. Nalewajko,  and M. C. Begelman, The Astrophysical Journal Letters 816, L8 (2016).
  56. L. Fletcher, B. R. Dennis, H. S. Hudson, S. Krucker, K. Phillips, A. Veronig, M. Battaglia, L. Bone, A. Caspi, Q. Chen, P. Gallagher, P. T. Grigis, H. Ji, W. Liu, R. O. Milligan,  and M. Temmer, Space Science Reviews 159, 19 (2011).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description