Onset of fast “ideal” tearing in thin current sheets: dependence on the equilibrium current profile

Onset of fast “ideal” tearing in thin current sheets: dependence on the equilibrium current profile

F. Pucci fpucci@nifs.ac.jp, fpucci@princeton.edu National Institute for Fusion Science, National Institutes of Natural Sciences, Toki 509-5292, Japan Princeton Plasma Physics Laboratory, Princeton University, Princeton, New Jersey, USA    M. Velli mvelli@ucla.edu University of California Los Angeles, USA    A. Tenerani annatenerani@epss.ucla.edu University of California Los Angeles, USA    D. Del Sarto daniele.del-sarto@univ-lorraine.fr Institut Jean Lamour UMR 7198 CNRS, Université de Lorraine, Nancy, France
July 30, 2019

In this paper we study the scaling relations for the triggering of the fast, or “ideal”, tearing instability starting from equilibrium configurations relevant to astrophysical as well as laboratory plasmas that differ from the simple Harris current sheet configuration. We present the linear tearing instability analysis for equilibrium magnetic fields which a) go to zero at the boundary of the domain and b) contain a double current sheet system (the latter previously studied as a cartesian proxy for the m=1 kink mode in cylindrical plasmas). More generally, we discuss the critical aspect ratio scalings at which the growth rates become independent of the Lundquist number , in terms of the dependence of the parameter on the wavenumber of unstable modes. The scaling with at small is found to categorize different equilibria broadly: the critical aspect ratios may be even smaller than with originally found for the Harris current sheet, but there exists a general lower bound .

Suggested keywords
Valid PACS appear here
preprint: AIP/123-QED

I Introduction

It is widely accepted that the mechanism responsible for explosive events, both in space and laboratory plasmas, is magnetic reconnection, a complex phenomenon that involves several scales and physical effects in which the dynamical evolution is affected by the initial configuration. One way in which reconnection may be triggered is via the tearing instability of thin current sheets. The latter is characterized by a fundamental parameter , which depends uniquely on the equilibrium magnetic field through the current density and its gradient, as well as the specific boundary conditions to which the plasma is subject. The parameter is also important for the subsequent non-linear evolution of magnetic islands in different geometrical configurations Furth, Rutherford, and Selberg (1973); Connor et al. (1988); Ottaviani et al. (2004), though on its own does not suffice to fully determine the dynamics of the system Militello, Grasso, and Borgogno (2014). Studying the instability of a Harris current sheet, Pucci and Velli (2014) (heareafter PV14) showed that for thin current sheets of thickness and length (assumed to be of the same order as the macroscopic length scale of the system), the specific dependence of the inverse aspect ratio on the Lundquist number determines when the instability becomes fast or “ideal”, where by the latter term it is implied that the growth rate no longer depends explicitly on the Lundquist number. In other words, PV14 showed that the limit of resistive MHD is not ideal MHD, but rather a magnetohydrodynamics where the frozen in condition breaks down in thin current sheets on ideal timescales , with the characteristic Alfvén speed. This suggests that reconnection survives in the limit but that no divergence of growth rates occurs, the latter result being implicated by the studies of the plasmoid instability (Loureiro, Schekochihin, and Cowley (2007); Huang and Bhattacharjee (2016)), based on the Sweet-Parker current sheet and corresponding island number scaling laws. More specifically, by assuming the form , PV14 showed that the critical inverse aspect ratio for current sheets was reached at , at which point the fastest growing tearing mode has a time-scale comparable to the Alfvén crossing time and, more fundamentally, is independent of . Therefore, the value separates slowly unstable current sheets (i.e. with growth rate scaling as a negative, fractional exponent of the Lundquist number) from those so violently unstable (with a growth rate scaling as a positive exponent of the Lundquist number, including the famed Sweet-Parker configuration) that they will never form in the first place. Simulations of a collapsing current sheet by Tenerani et al. (2015) demonstrated the fundamental correctness of the above scenario for a configuration comprised of an isolated current sheet of the Harris type.

PV14 also discussed a number of other consequences of “ideal” tearing: namely that the singular layer of the tearing mode has a thickness which scales as and a number of islands born out of the linear instability that scales as . The corollary of the result that at the fastest growing mode no longer scales with is that such modes, close to the maximum growth rates, remain the only growing modes at asymptotically large , as modes at any other values of the wave number end up decaying in time. Such result is now corroborated by other simulations (Huang et al., 2017), which include the effects of flows and the dynamics of current sheet thinning, even though flows are known to play a stabilizing effect at low Lundquist numbers (Chen et al. 2018, Bulanov et al. 1978).

The scaling in analysis in PV14 was fundamentally dependent on the details of the chosen equilibrium profile, the Harris current sheet with characteristic gradient scale . On the other hand isolated Harris-type current sheet are probably a rare occurrence in natural plasmas, where magnetic configurations more complex than the planar sheets are continuously perturbed by background turbulence and/or by macroscopic forcings, which may drive the field to ideal instability before reconnection takes place. Multiple current sheet configurations commonly occur in astrophysical plasmas: rapid rotators with an inclined magnetic axis, such as the case of pulsar magnetospheres Contopoulos (2007), lead to a waving current sheet which piles up into a series of alternating current layers almost orthogonal to the rotation axis. In the solar corona, the helmet streamers where the closed field is pinched off into the solar wind, are not always completely dipolar, but may contain multiple polarities, leading to the formation of multiple current sheets as the streamer is stretched into the wind Klein and Burlaga (1980); Woo (1994); Crooker et al. (1993); Dahlburg and Karpen (1995). Multiple current sheet configurations have been observed in laboratory plasmas as well: as pointed out e.g. by Ono et al. (1997), both for Reversed Field Pinches and Tokamaks, the periodic and/or continuous self-organization of their magnetic configurations produces multiple reconnection points Longcope and Klapper (2002); Bierwage et al. (2005).

In addition, recent attempts incorporate the “ideal” tearing concept in models of fully developed MHD turbulence Loureiro and Boldyrev (2017); Mallet, Schekochihin, and Chandran (2017): in these calculations the scaling relations for current sheet stability modify the energy spectral index of the cascade. But the current sheets arising in a turbulent cascade most likely have a whole spectrum of current profiles crossing them, so that modifications of the critical aspect sheet ratio dependencies arising from the detailed current profiles will also end up affecting the energy cascade. We will further comment on this latter aspect in the conclusion, but it is clearly of great interest to extend the stability study to more general configurations with differing magnetic field profiles.

Understanding the dependence of the classic tearing mode instability on the structure of any given equilibrium depends (at least assuming free or distant boundaries) on one fundamental parameter, the so called delta prime or . The latter is determined by resolving the first order perturbed momentum equation across the current sheet in the ideal limit. The perturbed, reconnecting magnetic field component (which will be called ) is found to be continuous across the sheet, but its first derivative retains a finite jump, (that will be ultimately resolved by including resistive diffusion in the calculation):


here depends only on the coordinate across the sheet and represents the typical equilibrium scale-length (or current sheet half-thickness).

The parameter depends on the normalized wave-number (along the sheet) of the reconnecting velocity and magnetic field perturbations. In this paper we investigate how, in the limit , different dependencies determine the corresponding “ideal tearing” scaling thresholds, i.e. the aspect ratio scaling index in , at which the maximum growth rate no longer depends on . This affects not only the critical aspect ratio, but also the number of islands that emerge from the linear phase once the critical aspect ratio is reached and the tearing instability becomes fast.

We will present detailed solutions for two cases: a configuration where the magnetic field far from the sheet vanishes and a double current sheet configuration. The plan of the paper is the following. In the next section we recall some terminology and properties related with the concept of the ”ideal” tearing (IT hereafter) instability and we describe it in terms of the parameters that define the equilibrium configurations. In section three we discuss the scaling relations for an equilibrium magnetic field which goes to zero at the boundary of the domain, recovering the trigger condition for the ”ideal” tearing mode to be developed. In section four we will discuss the transition to different regimes of a double current sheet configuration. Finally in the last section we discuss the results and the more general relevance to nonlinear plasma dynamics.

Ii Instability regimes in slab tearing reconnection.

Consider a general slab magnetic field configuration described as:


where is a function of the coordinate , defining the scale of the gradient of magnetic field and is a field which may depend on to guarantee equilibirum. We study the linearized resistive MHD equations, where the frozen in condition for magnetic field lines is locally violated due to the effect of a finite magnetic diffusivity ; we neglect kinetic scale effects in Ohm’s law (see Del Sarto et al. (2016); Pucci, Velli, and Tenerani (2017) for the extension of the “ideal” tearing scenario including such terms in Ohm’s law). Quantities are assumed to be uniform in the direction perpendicular to the plane of the equilibrium magnetic field, i.e. . The equations are non-dimensionalized by normalizing the magnetic field and its perturbation in terms of the mean field , the wavenumber along with and we introduce the non-dimensional displacement and . The growth rate is normalized to the macroscopic Alfvén time as previously defined, where the Alfvén speed ( is the plasma density); the Lundquist number is defined as . We therefore define a new independent variable , and denote derivatives with respect to with . We will consider an equilibrium current sheet aspect ratio scaling as . The tearing mode equations become:


where is the non dimensional growth rate of the instability. To recover the classic tearing results, in which , while the Lundquist number and the Alfvén speed must be replaced with and respectively. The asymptotic, large Lundquist number analysis of the solutions to the above equations, subject to the boundary conditions that the velocity and magnetic field perturbations vanish far from the sheet, show that the solution consists of two regions: a boundary layer of thickness around the center of the current sheet (which we will center around ), and an outer regions where diffusion may be neglected. Asymptotic matching then requires the jump in derivative of as calculated from the outer solution only, the parameter defined above, to equal the jump in derivative of calculated moving outwards from the inner solution.

The continuous spectrum of tearing unstable modes is then described in the two limits of large and small regimes (also known as the constant- ordering), in terms of whether or . The maximum growth rate can be found by matching the two regimes, i.e. where , Bhattacharjee et al. (2009); Loureiro, Schekochihin, and Uzdensky (2013); Del Sarto et al. (2016); Tenerani et al. (2016). We recall the results (where we have defined )


where is a coefficient given in terms of the values of Euler’s -function. If we introduce a dependence of on the wavenumber at small as a power law, and assume this to hold also at the maximum growth wave-number ,


where and are specified by the magnetic field equilibrium profile, we can introduce this into the small dispersion relation regime, Del Sarto et al. (2016). Assuming that the two regimes are correct at maximum growth, using Eq.(6-5) and we find


and for the singular resistive layer . For the classical Harris-pinch equilibrium, plotted in Fig.1, the estimation obtained in the instability range of the corresponding reveals to be quite reasonable, providing the known results and Furth, Killeen, and Rosenbluth (1963), verified numerically in Pucci and Velli (2014). As discussed in Landi et al. (2015), even in the presence of the general equilibrium in Eq.2 it is easy to show that the linear analysis of the instability is unchanged with respect to Pucci and Velli (2014).

The relationships in Eq.(7) can be generalized to a current sheet characterized by a thickness and a length and in particular, since and , we have and Pucci and Velli (2014) Tenerani et al. (2016). While the IT condition separates slowly unstable current sheets from supertearing instabilities Loureiro, Schekochihin, and Cowley (2007); Huang and Bhattacharjee (2016), the condition to have plasmoids, i.e. more than one island along the current sheet length associated to the fastest mode, since , is having an integer . This means , i.e. for an inverse aspect ratio


i.e. for instance if , .

Iii Boundary Magnetic Field

Consider an equilibrium magnetic field of the form


showed in Fig.1, where the equilibrium magnetic field is localized in the direction.

Figure 1: Typical profiles of the different magnetic field equilibria considered here: Harris equilibrium (red), magnetic field which goes to zero at the boundary of the domain (blue) and a double current sheet configuration (green).

For this equilibrium magnetic field profile, the parameter is


where and the instability regime occurs for , Porcelli et al. (2002). For Eq.(10) can be express in the form Eq.(6), for . From Eq.(7), for the maximum growth rate we find


which is verified numerically as shown in Fig.2.

Figure 2: Scaling of the maximum growth rate as a function of the Lundquist number and wave vector (in the inset) corresponding to the maximum growth rate for the classic tearing mode, i.e. inverse aspect ratio , for an equilibrium profile of the form Eq.(9).

Applying the renormalization criterion, as described in Sect.II, neglecting the dependence on the coefficients and , one gets


and assuming a scaling for the inverse aspect ratio , the condition


As shown in Fig.3 for this value of for the growth rate reaches an asymptotic value of , independent of the resistivity. For arbitrary , Eq. (14) generalizes to


so . Notice that , i.e. the critical inverse aspect ratio is larger in this case with respect to an initial Harris type current sheet, so still much thicker than the Sweet Parker current sheet configuration.

Figure 3: Growth rate normalized to the Alfvén time () as a function of the wave vector for a current sheet with an inverse aspect ratio , for an equilibrium profile of the form Eq.(9). Notice the maximum growth rate does not depend on the Lundquist number and is of the order of an Alfvén time.

Iv The double tearing mode.

Consider now a double current sheet configuration, described by


where is the semi-distance between the two rational surfaces and the thickness of the current sheet, see Fig.1. The mutual interaction between the rational surfaces may lead to a coupled tearing mode so that, depending on the distance which separates the two current sheets, the dispersion relation presents two asymptotic regions Pritchett, Lee, and Drake (1980). For the classic double tearing mode instability the latter paper shows that, setting , when


the growth rate scales as


The lower limit in Eq.(17) guarantees that , i.e. the inner resistive layer is much smaller than the distance between the two rational surfaces, which allows the matching between the inner resistive region and the outer ideal region. The upper limit is provided by where , where is the double kink mode growth rate, described by Eq.(17) of Pritchett, Lee, and Drake (1980). The latter depends on the equilibrium configuration, in particular for Eq.(16):


i.e. , so that provides the right inequality in Eq.(17). In the opposite case, for


the growth rate scales as


where, is the Euler gamma function. The upper limit is necessary to provide the solution in the ideal region for the asymptotic matching. The transition between the two regimes, for , occurs for . The regime identified by Eq.(18) appears to be the fastest one, in particular for a fixed distance between the two resonant surfaces, we have the maximum growth rate at the transition between the two regimes, i.e. for


We verified numerically the scaling for the maximum growth rate as a function of the Lundquist number, as well as the wavevector corresponding to the maximum growth rate, as shown in Fig.4

Figure 4: Scaling of the maximum growth rate as a function of the Lundquist number and wave vector corresponding to the maximum growth rate for the configuration described by Eq.(16).

If we want to generalize the IT condition for the double current sheet configuration we must take into account the condition , so


i.e. the distance between the two resonant surfaces is fixed and of the order of . Assuming a scaling for the inverse aspect ratio , the condition


verified numerically, as shown in Fig.5.

Figure 5: Growth rate for the configuration described by Eq.(16), as a function of the wave vector for . Different colors code different values of the Lundquist number.

In terms of wavevector the coupled tearing mode regime occurs for


Using (24):


where in our case. Since the minimum available wavevector is , for the coupled tearing regime is realized for a large number of modes (for instance the condition holds ). The Eq.(26) can be interpreted as a condition on the distance between the two rational surfaces in order to have a coupled tearing mode regime, i.e . For , for instance in the case of , the analysis in Pritchett, Lee, and Drake (1980) is invalid yet a coupled configuration should still exist. Indeed we estimated numerically the value of in the case finding that .

V Conclusions

In this paper we have discussed the critical aspect ratios required for the fastest tearing mode growth times to reach values comparable to ideal time-scales for current sheet equilibria differing from the Harris current sheet. Under fairly general conditions, it was shown that the crucial parameter affecting the scaling of the critical aspect ratio with the Lundquist number was the power-law scaling of the parameter in the large regime, . For the case of a single current sheet profile, such as the one presented in Sect.III, deviations from the limit of Eq.(6), modify the scaling relations of the classic tearing mode instability resulting ins modification of the trigger coefficient , constraining it to be in the range , i.e. with a thickness of the limiting current sheet greater than the Sweet Parker sheet type.

In the case of a double layer configuration, the distance between the two resonant surfaces, determines wether the tearing mode instability develops in a coupled or decoupled regime, while the maximum growth rate occurs at the transition between the two. The rescaling argument provides a value for the trigger aspect ratio very similar to the value estimated in the linear numerical solutions presented in Baty (2017).

As discussed in Tenerani et al. (2015, 2016) the critical aspect ratio criterion for fast reconnection is also important because it affects the subsequent nonlinear evolution of a current sheet or current sheet system in the nonlinear phase. More specifically, the criterion can be applied to understand the trigger of nonlinear secondary reconnection and the subsequent self-similar collapse, at least in 2D, giving insight into turbulence initiated by tearing in an initially laminar current sheet. With application to the secondary reconnecting instabilities in the nonlinear regime of the m=1 resistive kink mode in tokamaks, Del Sarto and Ottaviani (2017) have discussed the nonlinear transition to fast reconnection stressing the re-scaling (or embedding) effects due to the decreased magnetic field on the reconnecting current sheet compared to the overall global field intensity. They found subsequent non-linear secondary instability to evolve in this case on a very weakly dependent time-scale.

Recently, attempts have also been made to incorporate the effects of the fast tearing of thin current sheets into models of fully developed MHD turbulence (Loureiro and Boldyrev (2017), Mallet, Schekochihin, and Chandran (2017)). In these models, the authors assume - as has been often conjectured and been previously shown in many numerical simulations (see, e.g., the simulations of Rappazzo et al. (2008) and Yang et al. (2017)) - that the fundamental, dissipative coherent structure evolving out of the anisotropic nonlinear cascade is a current sheet (with aspect ratios depending on the scale and precise nature of the anisotropy). The idea is that the subsequent destabilization of the current sheet, once its aspect ratio has become too large, affects the subsequent nonlinear cascade and therefore also the inertial range spectrum. In this case, the trigger condition is given by the fastest tearing mode growth time becoming comparable to the nonlinear transfer time at the given scale (both time-scales independent of the Lundquist number). This gives a condition on the scale at which the inertial range spectrum is modified: the actual aspect ratio at which the tearing mode loses its scaling with the Lundquist number therefore also affects the subsequent inertial range spectrum. The results of the present paper are relevant in that they bracket the ranges of critical aspect ratios for a specific set of current profiles. It is important to note that in the turbulent cascade additional details will almost certainly become relevant, from the combined flow-current profile of the individual current sheet (coherent structure), to the precise topological nature of the structure itself (neutral vs non-neutral current sheets, for example). Finally, though the trigger aspect-ratios are fundamental in the 2D nonlinear evolution, a 3D system may evolve directly through ideal secondary instabilities permitted by the initial reconnection-determined topology change: for example, kink-type instabilities may take over. This would modify the nonlinear recursive resistive evolution discussed by Tenerani et al. (2015). We plan to address some of these issues in upcoming papers.

Vi Acknowledgements

This work was supported by the NASA Parker Solar Probe Observatory Scientist grant NNX15AF34G and the NSF-DOE Partnership in Basic Plasma Science and Engineering Award 1619611. The research presented also benefited greatly from discussions at the ISSI team meetings on “Explosive Processes in the Magnetotail: Reconnection Onset and Associated Plasma Instabilities” held in Bern, Switzerland from 17–21 October 2016 and 23-27 October 2017.


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