Fast Magnetic Reconnection: “Ideal” Tearing and the Hall Effect
Abstract
One of the main questions in magnetic reconnection is the origin of triggering behavior with on/off properties that accounts, once it is activated, for the fast magnetic energy conversion to kinetic and thermal energies at the heart of explosive events in astrophysical and laboratory plasmas. Over the past decade progress has been made on the initiation of fast reconnection via the plasmoid instability and what has been called “ideal” tearing, which sets in once current sheets thin to a critical inverse aspect ratio : as shown by Pucci and Velli (2014), at the time scale for the instability to develop becomes of the order of the Alfvén time and independent of the Lundquist number (here defined in terms of current sheet length ). However, given the large values of in natural plasmas, this transition might occur for thicknesses of the inner resistive singular layer which are comparable to the ion inertial length . When this occurs, Hall currents produce a three-dimensional quadrupole structure of magnetic field, and the dispersive waves introduced by the Hall effect accelerate the instability. Here we present a linear study showing how the ”ideal” tearing mode critical aspect ratio is modified when Hall effects are taken into account, including more general scaling laws of the growth rates in terms of sheet inverse aspect ratio: the critical inverse aspect ratio is amended to , at which point the instability growth rate becomes Alfvénic and does not depend on either of the (small) parameters . We discuss the implications of this generalized triggering aspect ratio for recently developed phase diagrams of magnetic reconnection.
1. Introduction
Magnetic reconnection is generally believed to be the mechanism responsible for explosive events in astrophysical, space and laboratory plasmas Yamada et al. (2010) Gonzalez and Parker (2016). Indeed this process allows the magnetic field to access new topologies, dissipating magnetic energy into heat and accelerated particles. However, in order for the energy to be initially stored in the field - a requirement for bursty, intermittent energy release to be possible - reconnection can not be occurring all the time: as measured with a clock based on the characteristic ideal dynamical time-scale, reconnection and resistive instabilities must have an off-on character. A complete understanding of magnetic reconnection therefore requires explaining how energy accumulates in the magnetic field, how the configuration suddenly become unstable, and how magnetic energy release occurs on very fast time scales; this implies a study of dynamics over multiple length and time scales. Since natural plasmas are characterized by magnetic Reynolds numbers which can vary by many orders of magnitude, e.g. from in Tokamaks up to in the solar corona, the dissipation scales can often approach kinetic (ion-inertial or gyroradius) scales, so a lot of work has been devoted to the role of kinetic effects in explaining fast magnetic reconnection Hosseinpur et al. (2009). Recently it has been shown Pucci and Velli (2014) that in high Reynolds number plasmas, and in particular for the ideal limit, , resistive instabilities survive and can become ideal, i.e. grow on the Alfvén time-scale, once an inverse aspect ratio is reached. For such current sheets, which are thicker than the Sweet-Parker sheets, (), the maximum growth rate does not depend on S, which means that the trigger condition is, in some sense, ”independent of the environment” (IT, for “ideal” tearing). This was confirmed by numerical simulations Landi et al. (2015), Tenerani et al. (2015). In particular in Tenerani et al. (2015), where a collapsing current sheet at high Reynolds numbers was simulated, it was shown that, once the critical aspect ratio is reached, the instability takes place on ideal timescales developing multiple plasmoids, from which a hierarchy of secondary IT instabilities takes place.
Even if the inverse aspect ratio of a collapsing current sheet never reaches the Sweet-Parker thickness, kinetic effects may become important at large enough Lundquist numbers if the typical collisional resistive scales are smaller than either the ion inertial length or the thermal ion gyroradius Cassak et al. (2006); Malakit et al. (2009). When characteristic scales approach the ion inertial length, the Hall effect starts to play a role, and a quadrupolar magnetic field emerges, affecting magnetic reconnection dynamics. Indeed there is now copious evidence that Hall reconnection occurs, both from magnetospheric observations Mozer et al. (2002); Cattell et al. (2005); Eastwood et al. (2007); Frank et al. (2016), at the dayside magnetopause Vaivads et al. (2004a) and in the near-Earth magnetotail Vaivads et al. (2004b); Borg et al. (2005); Nakamura et al. (2006), as well as in laboratory experiments Cothran et al. (2005); Ren et al. (2005); Yamada et al. (2006); Tharp et al. (2013); Kaminou et al. (2016) .
It is important to recognize that the Hall effect intervenes well before the overall current sheet thickness approaches , because the Hall term will affect the dynamics already when the inner, tearing mode resistive singular layer thickness approaches ion scales. As we shall see in the following sections, this leads to results which remain consistent with previous works, by reason of the properties of “ideal” tearing. In this paper we will consider the case of a planar configuration without a guide field: the effect of a finite ion inertial length in the framework of the tearing instability was first carried out by Terasawa (1983), where a schematic illustration of the effect of a quadrupole magnetic field on the planar configuration was shown.
Terasawa (1983) demonstrated that the Hall effect produced a growth rate enhancement at high Lundquist numbers, but the time for the instability to develop was still very slow because the current sheet thickness was assumed to be macroscopic. On the other hand, when the ion inertial length becomes of the order of the length scale which characterizes the equilibrium magnetic field gradient , ions and electrons decouple and whistler waves form, so two-fluid effects should be taken into account, and the growth rate should be scaled using typical whistler frequencies. Starting from these considerations, the regime where the inner resistive layer thickness becomes of the same order of the ion inertial length is investigated here, in sheets where the inverse aspect ratio scales as powers both of the inverse Lundquist number and of the normalized ion inertial length . The idea is to recover a specific scaling for which the growth rate becomes independent both of the inverse Lundquist number and the ion inertial length, i.e. intrinsically fast, independently of the dominant small parameter driving the instability.
The plan of the paper is the following: in section two we introduce the Hall effect and summarize the scaling relations of the fastest growing modes for current sheets with fixed aspect ratios. Section 3 then generalizes the “ideal” tearing instability criterion, discussing how the corresponding scaling relations are affected by a finite ion inertial length, and confirming them with numerical solutions of the eigenvalue equations. In the conclusions we discuss this result and place them into context.
2. The tearing instability in the presence of the Hall effect.
In this paragraph we summarize the classical stability problem of an equilibrium magnetic field configuration aligned along the axis ( ) and dependent on the perpendicular coordinate ( ) in the form of a Harris current sheet,
(1) |
where . If there is no out of plane (guide) field, a corresponding equilibrium pressure profile is required to guarantee equilibrium. Quantities are assumed to be uniform in the direction perpendicular to the plane of the equilibrium magnetic field, i.e. but the perturbed fields are 3D, so the reconnection region may develop a spatial 3D structure. That the velocity and magnetic field perturbation along and directions couple, may be seen by writing out the generalized Ohm’s law including the Hall term,
(2) |
where is the electric field, the velocity, the current, the magnetic diffusivity and the electron and ion density. In this form, Ohm’s law includes the Hall effect and collisional resistivity but neglects electron pressure and electron inertia: the effect of these terms on“ideal” reconnection regimes has been discussed in Del Sarto et al. (2016).
Upon linearization, the equations are non-dimensionalized using the magnetic field intensity (the perturbed magnetic field ), the magnetic field gradient scale (essentially the Harris sheet thickness), and the Alfvén time where the Alfvén speed . Wavenumbers along are also scaled with and we introduce the non-dimensional displacement . The growth rate is normalized to ; consequently, the Lundquist number is defined as . Denoting derivatives with respect to the non-dimensional coordinate with , the tearing mode equations become:
(3) | |||
(4) | |||
(5) | |||
(6) |
Here identifies the Hall coefficient which couples the and components of the perturbed magnetic field, and when the classic tearing mode equations (FKR) are recovered. As mentioned above, we expect the Hall term to become important when becomes of the same order as the thickness of the internal, singular layer describing the resistive tearing mode. Since for the classic resistive tearing mode such thickness (also normalized to the shear-length ) , for the fastest growing mode, scales as , we define
(7) |
so that means the Hall effect is no longer negligible.
The system of equations (3-6) is a sixth order two-point eigenvalue problem (for given ) in which the solutions develop large gradients in around the -axis with increasing . As mentioned above, the problem was first studied by Terasawa (1983), who showed that the tearing mode develops three characteristic layers: in addition to the sheet thickness and to the inner, singular layer, familiar from the purely resistive tearing mode Furth et al. (1963), an intermediate layer arises, in which the Hall current effect is also essential. This intermediate layer complicates the asymptotic analysis of the problem, introducing a dependence of the parameter , defined below, on : it is therefore best to resolve the eigenvalue problem numerically using the Lentini-Pereyra method, Lentini and Pereyra (1974), and compare the results to analytic estimations derived from an heuristic generalization of the classic tearing asymptotic matching (in the vein of Terasawa (1983)).
We first illustrate the two-layer development with finite by examining the behavior of the eigenfunctions and . Recalling that is an odd function of , it is easy to see that the ideal, marginal form of eqs. (3-6) yields solutions for , , , that are respectively even, odd and odd, even. Also, at great distances from the -axis, solutions decay exponentially. As proxies for the intermediate, Hall layer, and the inner singular layer around the -axis, we use the distance between two peaks of the displacement and of the Hall generated field , plotted in Fig.(1), left and right panels.
The left panel shows the eigenfunction for fixed Lundquist number and values , in red, blue, green respectively: one sees how the profile widens with increasing . Plotted in the right panel is the eigenfunction for , which displays a different behavior, with a peak to peak central thickness which decreases with increasing . It is not that does not display signatures of the internal, singular layer, it is just that it is less apparent, coming as it does in the form of an abrupt change in its gradient, rather than the more visible maximum/minimum that displays.
The thicknesses of both the intermediate layer and the inner resistive layer are plotted as a function of the parameter in Fig.2.
For very small , one expects to recover the thickness of the resistive layer , while for increasing the intermediate layer thickness increases (almost linearly) while the resistive, singular layer is found to decrease with as . The two curves in the figure do not appear to intersect at small because of the logarithmic scale (convergence would be at on this scale).
To understand the thinning of the inner resistive layer let’s generalize the classical tearing mode asymptotic matching theory heuristically, taking into account the effects of (). In the classic, resistive tearing mode without Hall effect, the maximum growth rate may be obtained by matching the scalings obtained in the so-called large and small regimes. is defined by searching for the solution of the linearized momentum equation neglecting the growth rate, i.e. assuming marginal stability:
Because the tearing mode is centered where the equilibrium magnetic field changes sign, the correct solution vanishes at large . It is easy to see then that the solution has a discontinuity in derivative as , and
which for the particular equilibrium at hand gives . The two regimes of the tearing mode are determined by whether the product (the large regime, also known as the resistive kink regime) or (small or the classical tearing mode regime): depending on which of the two regimes obtains, the second derivative of within the inner singular layer scales either as (large ) or (small ). As a result, the relationships between obtained by matching the inner to outer layer together with may be summarized as follows:
(8) | |||||
(9) | |||||
(10) |
Eq.(8) comes from the ideal terms (outer region) of the induction equation, eq.(9) from the inner, resistive layer of the induction equation, while eq.(10) comes from the momentum equation estimated within the inner, resistive layer. From these equations one immediately finds the growth rate scalings for the large and small regimes respectively, , . With increasing at any fixed modes transition from large to small , and the growth rate of the fastest growing mode at can be found by matching the two dispersion relations, taking into account that in the large regime one has : ,.
Extending this reasoning to include Hall terms requires analyzing the equation for the Hall field, though one must be careful in how to approximate and at the edge of the inner, resistive layer of thickness . Because has the same symmetry as , the dimensional estimate for the second derivative is . A dimensional estimate from eqs.(3-6) then shows, following the strategy of eqs. (8-10) that the Hall term in the induction equation, i.e. the last term in eq.(5) has magnitude
(11) |
In the denominators above, the third term is always negligible, while the term is either , in the small regime, or O(1), in the large regime. In both cases, once , the contribution of this Hall term contribution to eq.(5) dominates compared to the ideal MHD (convective term, first one on the rhs) contribution, leading to dispersion relations of the form
(12) | |||||
in the large and small regimes respectively. Again matching the two to obtain the scaling of for the maximum growth rate, taking into account that in the large regime one has , one finds
showing both that the growth rate of the instability is enhanced and that the resistive internal layer width decreases with increasing . The above scalings have been obtained neglecting the influence of on , which, as already observed by Terasawa (1981) is too strong an approximation. Indeed the direct numerical resolution of the eigenvalue problem shows that the growth rate does not quite follow the derived scalings. To fit the numerical results, we use the approximation for maximum growth rate
(13) |
where the well known scaling in the IT regime when Hall is negligible implies that . Numerical results for the maximum value of the growth rate as a function of , with four different fixed values of the Lundquist number (i.e. variable ) are shown in Fig.3.
Dashed lines join points with the same value of the Hall coefficient : note that the dashed lines in the regimes of small and large , while separately parallel, are not parallel across the transition at . This is because the Hall term influences the scaling with via rather than simply . Fitting the curve in Fig.(3), for we obtain
(14) |
where and . A quick check with the heuristic estimate shows that while the exponent for is the same, is significantly different, i.e. rather than .
We want to also remark that for all of our calculations , i.e. , and , which means energy conversion and outflow from the resistive region occur on Alfvénic timescales. The wavevector for which we have the maximum growth rate is not affected by , i.e. even for , .
3. Fast resistive magnetic reconnection in the presence of the Hall effect.
We come now to the question of how the critical resistive IT aspect ratio is modified once the resistive layer thickness becomes comparable to the ion inertial length. Following our previous work, we will consider current sheets with a macroscopic length , which will therefore be used to define the Lundquist number and Alfvén time: , . Let us first consider the resistive “ideal” tearing mode at the critical current sheet aspect ratio, scaling as . It was shown (Pucci and Velli, 2014) that in this case resistive inner layer scales as . It therefore seems logical to generalize the parameter defining the relevance of the Hall effect to
(15) |
To confirm this, we first consider a sequence of equilibria at various aspect ratios at large and scaling as , but fixing the value of . As increases, we expect the maximum growth rate to be constant until , at which point the Hall effect acceleration should lead to an increase of the growth rate. This is because once is passed, the scaling thins the current sheet too much, leading to the same paradox of the plasmoid instability on SP sheets,Loureiro et al. (2007), namely a growth rate which diverges with increasing . This is confirmed by the numerical solution of the eigenvalue equations shown in Fig.4, where at first (dark blue through dark green lines) one sees that the maximum growth rate is the same, independently of , yet once is surpassed, the maximum growth rate rises again (light green through pink curves), rapidly increasing with increasing . What this means is that once the Hall effect becomes important, the critical aspect ratio should no longer depend only on , but also on the parameter or , as there are now two asymptotic parameters at work.
We therefore generalize to an inverse aspect ratio which varies in the parameter space , scaling as
(16) |
and search for exponents such that the time for the instability to develop becomes independent of the parameters of the system. Starting from Eq.(2), we must now renormalize all quantities to rather than :
(17) |
where was determined in the previous section numerically (. From the definition eq.(15), inserting the explicit aspect ratio dependence from eq.(16) we obtain
(18) |
We now search for the values of that cancel the growth rate dependence on . First, notice that the growth rate is independent of the Lundquist number for , independently of the value of . Then, for we may write:
(19) |
which allows us to calculate the exponents for which :
(20) |
and the error coming from the numerical determination of is then .
4. Results.
The correctness of the prediction for the scaling relation of the inverse aspect ratio eq.(20) was verified numerically: first, in Fig.5 we show solutions for the growth rate as a function of wave-number for different values of and a constant value (). The asymptotic value of the growth rate is , as expected for a growth rate independent of the parameters.
Second, we show the values of the growth rate, again as a function of , at constant values of but different values of (obtained by varying ), in Fig.6. Notice how remains constant,, and there is no shift in wavenumber of the maximum growth rate with changes in .
The fact that the coefficient does not have direct physical significance, stemming as it does from our definitions of aspect ratio scaling: as one may immediately verify, one could change the specific values of these coefficients by redefining the aspect ratio scaling by including an arbitrary multiplicative constant. What is important of course is the value of the scaling exponents which determine the critical aspect ratios beyond which current sheets become so strongly unstable that they will never form.
Our choice for the normalization time, , is valid only if so that (the whistler timescale). We want to verify this hypothesis is satisfied by our critical aspect ratios. Starting from the Eq.(16), we have
(21) |
which means that for our hypothesis is verified, while for we should normalize the growth rate to the whistler time. These quantities have been listed in Tab. 1 for some astrophysical as well as for laboratory plasmas, where we assumed and, if , the inverse aspect ratio , while if , .
Plasma | Solar Cromosphere | Solar Corona | Solar Wind | Magnetotail | MRX |
---|---|---|---|---|---|
n | |||||
L | |||||
B | |||||
T | |||||
S | |||||
/L | |||||
(0.4-1) | 0.2-0.5 | ||||
5. Comments and conclusions.
We have studied the linear resistive tearing instability for current layers whose thickness approaches the ion inertial length, in which the the instability growth rate and parameters are modified by the Hall effect. We have generalized the “ideal” tearing criterion, taking into account a finite ion inertial length, ending up with a trigger relation for the aspect ratio which varies in the parameter space , depending on the specific plasma parameters. The result is a couple of values and defining a critical inverse aspect ratio, scaling as , below which the reconnection process becomes explosive. Recently phase diagrams involving Lundquist number and the macrosopic system size in units of the ion inertial lenght (or ion sound gyroradius, if a guide field is present) have been created, which summarize the essential dynamics of the plasma for a wide range of parameters Ji and Daughton (2011). Fig. 7 summarizes the results obtained using the IT criterion, in the absence of a guide field (in the figure, ). The parameter space has an extension on the left due to the fact that we also take into account the effect of finite electron skin depth in the RMHD case discussed in Del Sarto et al. (2016), where . Indeed the Hall effect on its own can not break the frozen in conditions, i.e. collisionless reconnection must be triggered by other effects, and the inertial terms in Ohm’s law are proportional to . As in Ji and Daughton (2011) we have a region (blu one) where the single X-line reconnection occurs i.e. for Lundquist numbers smaller than the critical one determined by inflows and outflows, see e.g. Tenerani et al. (2016a); Shi et al. (2016). We have investigated the purple region in Pucci and Velli (2014) and the white region within this work. While the orange region has been investigated in Del Sarto et al. (2016) (even if the Hall effect is negligible in the adopted frame), the green one still has to be explored and with that, a possible critical value of for which collisionless reconnection can present single or multiple x-points.
These are useful to understand the ongoing process regime in astrophysical plasmas as well as in laboratory plasmas, and in the latter case to visualize the parameter space of an experimental facility, and in particular can be applied to the multi-point observations of MMS, to compare theoretical predictions the spatial structure of the Hall magnetic and electric fields surrounding the diffusion region. The next step for linear studies is to include the 3D structures which naturally arise in the presence of a mean magnetic field in the direction orthogonal to the plane where magnetic reconnection occurs. In this case other effects may occur due the dependence on the third direction and the system of equations would be more complicated, involving also the complex component of the eigenfunctions. Then, electron pressure terms should be included. Such terms will introduce the dependence of critical aspect ratio also on the thermal ion gyroradius. From the non-linear evolution point of view, it would be of great interest to follow a collapsing current sheet (initially with 0 guide field), into the Hall regime. The resistive internal singular layer, which in 2 Dimensions becomes the thickness of secondary current sheets in nonlinear evolution (Tenerani et al. (2016b)), and where the Hall magnetic field is different from zero, tends to shrink for larger values of the Hall parameter, making it a challenging problem which we hope to address in the near future.
6. Acknowledgements
This work was supported by the NASA Solar Probe Plus Observatory Scientist grant and the NSF - DOE partnership in basic plasma science and engineering.
Footnotes
- affiliationmark:
- affiliationmark:
- affiliationmark:
- affiliationmark:
- affiliationmark:
- affiliationmark:
- affiliationmark:
References
- M. Yamada, R. Kulsrud, and H. Ji, Reviews of Modern Physics 82, 603 (2010).
- W. Gonzalez and E. Parker, eds., Astrophysics and Space Science Library, Astrophysics and Space Science Library, Vol. 427 (2016).
- M. Hosseinpur, N. Bian, and G. Vekstein, Phys. Plasmas 16, 012104 (2009).
- F. Pucci and M. Velli, ApJ 780, L19 (2014).
- S. Landi, L. Del Zanna, E. Papini, F. Pucci, and M. Velli, ApJ 806, 131 (2015), arXiv:1504.07036 [astro-ph.SR] .
- A. Tenerani, A. F. Rappazzo, M. Velli, and F. Pucci, ApJ 801, 145 (2015), arXiv:1412.0047 [physics.plasm-ph] .
- P. A. Cassak, J. F. Drake, and M. A. Shay, ApJ 644, L145 (2006), physics/0604001 .
- K. Malakit, P. A. Cassak, M. A. Shay, and J. F. Drake, Geophys. Res. Lett. 36, L07107 (2009).
- F. S. Mozer, S. D. Bale, and T. D. Phan, Physical Review Letters 89, 015002 (2002).
- C. Cattell, J. Dombeck, J. Wygant, J. F. Drake, M. Swisdak, M. L. Goldstein, W. Keith, A. Fazakerley, M. André, E. Lucek, and A. Balogh, Journal of Geophysical Research (Space Physics) 110, A01211 (2005).
- J. P. Eastwood, T.-D. Phan, F. S. Mozer, M. A. Shay, M. Fujimoto, A. Retinò, M. Hesse, A. Balogh, E. A. Lucek, and I. Dandouras, Journal of Geophysical Research (Space Physics) 112, A06235 (2007).
- A. G. Frank, A. V. Artemyev, and L. M. Zelenyi, Soviet Journal of Experimental and Theoretical Physics 123, 699 (2016).
- A. Vaivads, M. André, S. C. Buchert, J.-E. Wahlund, A. N. Fazakerley, and N. Cornilleau-Wehrlin, Geophys. Res. Lett. 31, L03804 (2004a).
- A. Vaivads, M. André, S. C. Buchert, J.-E. Wahlund, A. N. Fazakerley, and N. Cornilleau-Wehrlin, Geophys. Res. Lett. 31, L03804 (2004b).
- A. L. Borg, M. Øieroset, T. D. Phan, F. S. Mozer, A. Pedersen, C. Mouikis, J. P. McFadden, C. Twitty, A. Balogh, and H. Rème, Geophys. Res. Lett. 32, L19105 (2005).
- R. Nakamura, W. Baumjohann, Y. Asano, A. Runov, A. Balogh, C. J. Owen, A. N. Fazakerley, M. Fujimoto, B. Klecker, and H. RèMe, Journal of Geophysical Research (Space Physics) 111, A11206 (2006).
- C. D. Cothran, M. Landreman, M. R. Brown, and W. H. Matthaeus, Geophys. Res. Lett. 32, L03105 (2005).
- Y. Ren, M. Yamada, S. Gerhardt, H. Ji, R. Kulsrud, and A. Kuritsyn, Physical Review Letters 95, 055003 (2005).
- M. Yamada, Y. Ren, H. Ji, S. Gerhardt, M. Inomoto, R. Kulsrud, S. Dorfman, and Y. Wang, in APS Meeting Abstracts (2006).
- T. D. Tharp, M. Yamada, H. Ji, E. Lawrence, S. Dorfman, C. Myers, J. Yoo, Y.-M. Huang, and A. Bhattacharjee, Physics of Plasmas 20, 055705 (2013).
- Y. Kaminou, M. Inomoto, and Y. Ono, Plasma and Fusion Research 11, 2401052 (2016).
- T. Terasawa, Geophys. Res. Lett. 10, 475 (1983).
- D. Del Sarto, F. Pucci, A. Tenerani, and M. Velli, Journal of Geophysical Research (Space Physics) 121, 1857 (2016), arXiv:1511.00035 [physics.plasm-ph] .
- H. P. Furth, J. Killeen, and M. N. Rosenbluth, Physics of Fluids 6, 459 (1963).
- M. Lentini and V. Pereyra, Math. Comp. 28, XXXXXX (1974).
- N. F. Loureiro, A. A. Schekochihin, and S. C. Cowley, Physics of Plasmas 14, 100703 (2007), astro-ph/0703631 .
- M. J. Aschwanden, Physics of the Solar Corona (Praxis Publishing Ltd, 2004).
- R. Esser and S. R. Habbal, Geophys. Res. Lett. 22, 2661 (1995).
- V. Angelopoulos, A. Runov, X.-Z. Zhou, D. L. Turner, S. A. Kiehas, S.-S. Li, and I. Shinohara, Science 341, 1478 (2013), http://science.sciencemag.org/content/341/6153/1478.full.pdf .
- V. A. Sergeev, D. G. Mitchell, C. T. Russell, and D. J. Williams, J. Geophys. Res. 98, 17345 (1993).
- M. G. Kivelson and C. T. Russell, Introduction to Space Physics, Edited by Margaret G. Kivelson and Christopher T. Russell, pp. 586. ISBN 0521451043. Cambridge, UK: Cambridge University Press, April 1995. (1995) p. 586.
- M. Yamada, J. Yoo, J. Jara-Almonte, H. Ji, R. M. Kulsrud, and C. E. Myers, Nature Communications 5, 4774 (2014).
- H. Ji and W. Daughton, Physics of Plasmas 18, 111207 (2011), arXiv:1109.0756 [astro-ph.IM] .
- A. Tenerani, M. Velli, and C. Shi, Abstract SH51B-2594 presented at 2016 Fall Meeting (2016a).
- C. Shi, A. Tenerani, and Velli, Abstract SH51B-2596 presented at 2016 Fall Meeting (2016).
- A. Tenerani, M. Velli, F. Pucci, S. Landi, and A. F. Rappazzo, Journal of Plasma Physics 82, 535820501 (2016b), arXiv:1608.05066 [astro-ph.SR] .