Effective Resistivity in Collisionless Magnetic Reconnection

An effective resistivity relevant to collisionless magnetic reconnection (MR) in plasma is presented. It is based on the argument that pitch angle scattering of electrons in the small electron diffusion region around the X line can lead to an effective, resistivity in collisionless plasma. The effective resistivity so obtained is in the form of a power law of the local plasma and magnetic field parameters. Its validity is confirmed by direct collisionless particle-in-cell (PIC) simulation. The result agrees very well with the resistivity (obtained from available data) of a large number of environments susceptible to MR: from the intergalactic and interstellar to solar and terrestrial to laboratory fusion plasmas. The scaling law can readily be incorporated into existing collisional magnetohydrodynamic simulation codes to investigate collisionless MR, as well as serve as a guide to ab initio theoretical investigations of the collisionless MR process.

Magnetic reconnection (MR) often occurs in plasmas containing sheared magnetic fields and can efficiently convert magnetic energy into the kinetic and thermal energies of the charged particles [1][2][3] . The process plays important roles in the evolution of the solar corona 4,5 , the geomagnetic tail 6,7 , the magnetosphere 8,9 , the intergalactic and interstellar, as well as laboratory plasmas 10,11 . In particular, collisionless or fast MR (FMR) on time scales much less than the inter-particle collision time can occur. FMR has often been attributed to anomalous resistivity arising from local current-instability driven turbulence in the small electron diffusion region of the MR 12,13 . However, the evolution and effect of the turbulence during the MR are difficult to follow and remain unclear. Existing studies 14,15 have noted that the lifetime of a particle in the diffusion region can be considered as an effective collision time, since in terms of its momentum and energy changes, the electron dynamics in the electron diffusion region resembles that of electrons being scattered by collisions, except that here the scattering partners are the local magnetic and electric fields. In particular, Speiser 14 introduced an ad hoc resistivity (more precisely, conductivity) based on the behavior of the electric current flow in the diffusion region. However, the problem remains unclear and no general conclusion can be drawn [14][15][16] .
It is well known that pitch-angle scattering of electrons in highly bending magnetic fields such as that in the diffusion region around the X point of MR can lead to particle momentum transfer from the parallel to the perpendicular initial current direction. In this paper, we reconsider the dynamics of electrons in this small region. The transit times of typical electrons in the region near the X line are determined by following their motion as the FMR process evolves. An effective resistivity in the form of power-law scaling of the most relevant local plasma parameters is obtained by replacing the mean-free-time in the expression for the collisional resistivity by an ensemble averaged electron transit time that depends on the local plasma and field parameters in the diffusion region. Validity of our approach is confirmed by full particle-in-cell (PIC) simulation of the FMR. Moreover, when compared with a large number of plasmas susceptible to MR: from the intergalactic and interstellar space to solar and terrestrial, as well as fusion, plasmas, it is found that the effective resistivity agrees very well with that estimated from the known parameters of these plasmas. The scaling law can readily be incorporated into the existing macroscopic MHD simulation codes 17 for investigating FMR in complex space and fusion plasmas, as well as serve as guide for detailed theoretical investigation of the FMR physics.

Analytical formulation of effective resistivity
Accordingly, collisionless MR can be investigated by replacing the mean free time τ mft between collisions in the collisional resistivity η coll = m e /ne 2 τ mft , where n, e, and m e are the electron number density, charge, and mass, respectively, by the mean transit time τ transit of electrons in the diffusion region around the X line (see Fig. 1 and MHD simulation codes. In the following, we shall obtain τ transit by concentrating only on, in our opinion, the most relevant physics involved. The local magnetic and electric fields in the electron diffusion region of the MR ( Fig. 1) can be approximated by where L x and L z are the characteristic lengths of the electron diffusion region, respectively, and the coefficients B x , B z (<<B x ), and E y are constants. When L z << L x , the current sheet becomes elongated. This MR geometry is ofter referred to as of Y type. The electron trajectory near the X line is then governed by where for simplicity we shall assume (consistent with the PIC simulation results below) that the electron velocity v y in the current sheet varies only a little from the average value 〈v y 〉 ~ −J y /ne.
Equations (3) and (4) yield x , and (x 0 , z 0 ) and (v x0 , v z0 ) are the initial position and velocity of the electron.
Equations (5) and (6) show that the electron is accelerated in the x direction but it only oscillates in the z direction. We can thus consider the transit time τ transit as the time for the electron to traverse the diffusion region in the x direction. Accordingly, from equation (5) we get can be chosen to be at the edge of the simulation box.
Near the X line, we can also reasonably assume that thermal effects can be neglected and the initial in-plane electron velocity is nearly zero. Considering that the electron motion in the z direction is oscillatory, for the transit time we only need to follow its motion in the x direction. Accordingly, equation (7) becomes Since the initial position x 0 of an electron can be anywhere inside the diffusion region, the mean transit time is x and induction electric field =Ê E y y in the diffusion region. The current sheet is in the y direction. The magnetic field increases (from null) with the distance from the X line, which is in the y direction and appears in the x, z plane here as the X point at (0, 0). Several idealized electron trajectories are shown as dashed curves.
So that the effective collisionless resistivity is 3 2 In deriving the electron transit time, we have made the reasonable but unsubstantiated assumption on the existence of a collisionless, or effective, resistivity that imitates the function of the classical collisional resistivity in the fluid description of the plasma. For verification, we next carry out full PIC simulations of the FMR in collisionless plasma.

Particle-in-Cell simulation
We have performed 2.5D full PIC simulations for plasma particle motion in the diffusion region by assuming ∂ y = 0. For simplicity, we use the charge-conservation scheme (CCS) instead of solving the Poisson equation, and the finite difference time domain (FDTD) method to solve the other Maxwell's equations. The equations used in the PIC simulations are where c is the light speed, v j and p j = m j v j are the particle velocity and momentum, respectively. The variables are normalized as follows: Ai0 , the ion-to-electron mass ratio μ = m i /m e is from 25 to 400, and the initial ion-to-electron temperature ratio is T i /T e = 5. Our simulation domain is −D x /2≤ x ≤D x /2, −D z /2 ≤ z ≤D z /2, where D x = 12.8d i0 , D z = 6.4d i0 , dx = dz = 0.01d i0 and the time step is ω ci0 Δt = 0.0002. Periodic and closed boundary conditions are adopted in the x and z directions, respectively. Nearly 82 million simulation particles for each species are used.
We use the Harris equilibrium as the initial configuration. The initial magnetic field is given by In the simulation, the reconnection process is initiated by a small perturbation of the magnetic field. Pressure balance yields where P and B are the local thermal pressure and magnetic field, β = P/(B 2 /2), and P is normalized by μ B /2 0 2 0 . Here we set β = 0.2.

Comparison of the analytical and simulated resistivities
We first examine the y component of the velocity v y of electrons entering and leaving the electron diffusion region. For ion-to-electron mass ratio μ = 400, during the peak reconnection period (from t = 16 to 17) we found that the average change of v y is less than 10%. That is, v y is indeed roughly constant, as assumed in the evaluation of (5) and (6).
Next we compare the resistivities from our analytical model and the PIC simulation. Figure 2 shows the evolution of the resistivities in the electron diffusion region. In the analytical formula for η eff , the electron number density n, velocity v y , B z , and L x have been replaced by, as calculated from the PIC simulation results, the average electron number density n and electron velocity v y in the electron diffusion region (of size of d e ), the maximum B z , and the characteristic length L x of B z in the X-point region, respectively. On the other hand, the resistivity from the PIC simulation is calculated directly from the relation η s = E y /J y by substituting the measured values of the average electric field intensity E y and sheet-current density J . We see that our effective resistivity agrees quite well with the collisionless resistivity obtained from the PIC simulations. In particular, in the fast reconnection phase both resistivities increase rapidly and in a similar manner, thereby verifying the scaling of the field parameters in our model effective resistivity. We can also see that the peak value of the resistivity decreases with increase of the mass ratio μ, and the analytical and simulation results approach each other.

Discussion and Summary
In collisional plasma, the characteristic thickness of the current sheet taking into account magnetic field diffusion is 18,19 η where L is the plasma size, v A is the Alfven speed, and η spz is the Spitzer resistivity (based on Coulomb collisions). For the parameters of the Earth's magnetopause and magnetotail 20 , the half-thickness of the thinnest current sheet during the nonlinear stage of magnetic reconnection can be of the order 10 m that is four to five orders of magnitude smaller than that obtained from the satellite data during magnetic reconnection: namely, ~100 km in the magnetopause 21,22 and ~1000 km in the magnetotail 23,24 . Such a huge descrepancy suggests that Coulomb collisions are not the dominant dissipation mechanism for MR in the magnetopause and magnetail.
The effective resistivity in our model can be rewritten in the form of a power law: x z z x between the local magnetic field components near the X line. The parameter α is thus the ratio of the characteristic lengths of B x and B z . Simulations have shown that in the MR configuration α increases as reconnected magnetic field B z increases with time in the diffusion region and α is of order 0.1 when MR gets into the nonlinear stage. In the magnetotail 20 on Earth's nightside, the electron density is n e (n i ) ≈ 0.3cm −3 , the electron temperature is T e ≈ 600eV, the magnetic field is B 0 ≈ 2 × 10 −8 T, and plasma size L ≈ 100R E , the current sheet thickness is about Δ SP ≈ 0.15R E ≈ 953km in the nonlinear phase of magnetic reconection with α = 0.1. On the other hand, for the magnetopause 20 on the Earth's dayside, with electron density n e (n i ) ≈ 10cm −3 , electron temperature T e ≈ 300eV, magnetic field B 0 ≈5 × 10 −8 T, and plasmas size L ≈ 10R E , the current sheet thickness is Δ SP ≈ 0.021R E ≈ 137km. Thus, for both the magnetopause and magnetotail, the current sheet thicknesses as predicted by our model are in good agreement with that from the satellite observations. Moreover, one can easily show that for the parameters of the experimental device MRX 20 , our model yields Δ SP ≈ 3.38 cm, which is in good agreement with that from the direct laboratory measurement.
It is of interest to make a broader comparison of the results from our model with that from existing data on space and laboratory plasmas where MR is observed or expected to exist. Since classical collisions can be important or relevant in some of the environments 20 , it is useful to introduce the total resistivity η tot =η spz +η eff . Figure 3 shows the plot of η tot normalized by the Spitzer resistivity η spz versus the mean-free-path λ mfp normalized by the Sweet-Parker current sheet thickness Δ SP = 2L z . From Fig. 3, it is clearly shown that all data points are distributed near the best-fitting line.  collisionless (λ mfp /Δ SP > 1) plasmas. The slope of the best-fit (red dashed) line in Fig. 3 is C = 0.1, which is in the middle of its range. We note that the data points only slightly deviate from the best-fit line, and can be attributed to uncertainties in the obervational data. The somewhat larger deviation for the data points from the magnetic-confinement-fusion devices ITER and TFTR can be attributed to the strong guiding magnetic field. It is well known that a guide field can suppress magnetic reconnection or reduce the effective resistivity, which is consistent with the fact that the ITER and TFTR data points are located below the best-fit line. The results here suggest that pitch angle scattering of electrons due to bending magnetic field lines, which is the basic assumption of our model, may be responsible for fast MR in collisionless plasmas.
The effective resistivity in our 2D model is mainly a result of diversion of electrons in the electron diffusion region, which has very small spatial scale. Three-dimensional (3D) effects are ignorable if the spatial scale of the magnetic field in the third direction is larger than the electron inertia length. This condition is usually valid for space and laboratory plasmas, which explains why the effective resistivity from our model agrees very well with that of a large number of environments, namely from intergalactic and interstellar to solar and terrestrial to laboratory fusion plasmas. The reconnection rate in general depends on the ratio between the thickness and length of the diffusion region, or the current sheet. The thickness is usually determined by dissipation, or resistivity, of the system. Therefore, for given resistivity, the reconnection rate can increase with decrease of the current sheet length, say by an external driving force. For example, in tokamak plasmas, the slower tearing modes correspond to spontaneous MR and the sawtooth oscillations correspond to FMR driven by internal kink instabilities.
In summary, for understanding FMR we have introduced an effective resistivity that contains no free parameters. The effective resistivity is based on self-consistent scattering or acceleration of electrons by bending of magnetic field lines, and it agrees in magnitude with that of a large number of environments where MR is observed or suspected. It can also be readily adapted in existing collisional-fluid simulation codes for investigating collisionless FMR. The present work can also serve as a guide for a formal first-principles derivation of such an effective resistivity. Finally, it may be of interest to point out that our results are clearly also applicable to very small scale and very fast MR in the absence of ion dynamics. In fact, such novel ultrafast (45 millisecond) MR phenomena have been recently reported to be occuring within the entangled magnetic fields in the Earth's turbulent magnetosheath. 25 Figure 3. Plot of the normalized total resistivity η tot /η spz for space and laboratory plasmas exhibiting MR versus the normalized mean-free-path λ mfp /Δ SP . Here, the thickness of the current sheet according to the Sweet-Parker model is Δ =L L S 2 2 / , SP z L 1/2 where S L = μ 0 Lv A /η tot is the Lundquist number, or the ratio of the resistive to Alfvèn times, and L is the overall plasma size. Note that the normalization parameters η spz and Δ SP are different for different data points. The solid circles and triangles are for collisional (λ mfp /Δ SP <1) and collisionless (λ mfp / Δ SP >1) plasmas, respectively. The corresponding current sheet thicknesses Δ SP based on our effective resistivity are given in the insets. The red dashed line is the best fit for the data points in the collisionless plasma regime (λ mfp /Δ SP > 1).