A Hard Energy Spectrum in 3D Guide-Field Magnetic Reconnection

Magnetic reconnection has long been known to be the most important mechanism not only for mixing the plasmas by changing the magnetic field topology but also for releasing the magnetic field energy into the plasma kinetic energy. During magnetic energy release, it is possible for some of the heated plasma to be accelerated to energies much higher than the thermal energy. Recently, the energy partitioning of the thermal and the nonthermal energy has been studied by using particle-in-cell (PIC) simulations, and it has been shown that the acceleration efficiency of nonthermal particles increases with increasing the plasma temperature, and the nonthermal energy density occupies more than 90% in the total heated plasma when the Alfven velocity is close to the speed of light c. However, the acceleration efficiency decreases as the guide magnetic field increases. So far the acceleration efficiency has been mainly studied in two-dimensional systems, but it is interesting to study three-dimensional effects where the patchy and turbulent reconnection can dynamically occur. This study explores the effects of three-dimensional relativistic reconnection on a pair plasma with the guide magnetic field, utilizing three-dimensional (3D) PIC simulations. The results indicate that the decrease in nonthermal particle production is smaller in 3D guide-field reconnection compared to 2D. More importantly, contrary to general expectation, 3D reconnection is capable of maintaining a hard nonthermal energy spectrum even in the presence of a strong guide magnetic field.


I. INTRODUCTION
Magnetic reconnection is widely recognized as a crucial mechanism for plasma heating and particle acceleration in various explosive plasma phenomena in the Universe [1][2][3][4][5] .The magnetic field energy stored in the plasma sheet is released, which results in the thermalization of the plasma up to the equivalent temperature of the Alfvén velocity.Furthermore, a portion of the plasma is accelerated to higher energies above the thermal temperature, forming a nonthermal component that can be approximated by a power-law spectrum.
Since there is much observation evidence that the nonthermal particles associated with the hot thermalized plasma formed during the solar flares 6,7 and the magnetospheric substroms in the Earth's magnetosphere [8][9][10] , it is believed that the generation of nonthermal particles in reconnection is a ubiquitous process in the collisionless plasma environment.Then the high-energy phenomena in a relativistic plasma medium observed around black holes, pulsar nebulae, and magnetars, etc., are also thought to be related to magnetic reconnection [11][12][13][14] .Indeed, the numerical simulations have shown that the efficiency of the nonthermal particle acceleration is extremely high in such a relativistic plasma regime with V A ∼ c, and the formation of a hard energy spectrum has been discussed [15][16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31] .
While many observational and theoretical/computational studies have been conducted to study nonthermal particle acceleration in the plasma sheet, it was not necessarily investigated by focusing on the energy partitioning between the thermal energy density and the nonthermal energy density 32,33 .Recently, we discussed the energy partitioning during magnetic reconnection with a pair plasma by using a two-dimensional particle-in-cell (PIC) simulation 34,35 .It is found that the production of nonthermal particles increases with increasing Alfvén speed for magnetic reconnection 34 .In the relativistic reconnection, where the Alfvén speed is close to the speed of light, it has been shown that the nonthermal energy density can occupy more than 90% of the total kinetic plasma energy density for an anti-parallel magnetic field geometry without a guide magnetic field.In addition,we have discussed that the production of nonthermal particles decreases with the increase of a guide magnetic field 35 .
So far our study on the energy partitioning of thermal and nonthermal plasma has been limited to two-dimensional (2D) system 35 (Paper I), but it is important to investigate threedimensional (3D) reconnection.In 2D reconnection, the magnetic field "reconnection" can only occur in the neutral sheet, but 3D reconnection can occur not only in the neutral sheet but also outside the neutral sheet [36][37][38][39] .It has been shown that the three-dimensional reconnection has more dynamic time evolution, leading to patchy and turbulent reconnection [37][38][39] .In addition, it has been discussed that the efficiency of electron acceleration in a 3D and non-relativistic reconnection is enhanced when a guide magnetic field is comparable to the reconnecting magnetic field component 38 .After these previous studies on the importance of 3D turbulent reconnection and plasma transport with a guide magnetic field, it is intriguing to investigate the nonthermal particle acceleration in a relativistic plasma with a more stronger guide magnetic field, because the patchy and turbulent reconnection is expected to be present in a wide range of plasma parameters in the plasma universe.Here we study the efficiency of the nonthermal particle acceleration of the 3D relativistic reconnection as a function of the guide magnetic field by using 3D PIC simulations.We show that the efficiency of the magnetic energy release in 3D is almost the same as that in 2D for a weak guide magnetic field, and that the efficiency of both 2D and 3D reconnection decreases with increasing the guide magnetic field.The decrease in efficiency of 3D, however, is less than that of 2D.More importantly, while the kappa index (the power law index) increases with increasing the guide magnetic field for 2D reconnection, the nonthermal spectra in 3D reconnection with a strong guide magnetic field regime retain a hard spectrum.

II. TURBULENT 3D GUIDE-FIELD RECONNECTION
Our objective is to understand particle acceleration in three-dimensional reconnection with a guide magnetic field, but before discussing the particle acceleration, let us study the difference and similarity between 2D and 3D reconnection.We use basically the same 3D simulation code as that used in the study of the magnetorotational instability 40 .The numerical scheme with a semi-implicit time integration is the same as that used in our previous study of 2D reconnection 35 .We assume a periodic boundary in the x and z directions with the system size of 0 < x/λ < 322.5 and −80.7 < z/λ < 80.7.We use the so-called conducting wall in the y direction with the upper and lower boundaries situated at y = ±80.7λ,where λ is the thickness of the initial plasma sheet.The computational grid size is (N x , N y , N z ) = 4032 × 2016 × 2016.The total number of particles is set as N p = 1.6 × 10 11 , and the average number of particles per grid size at the neutral sheet is 81.For the two-dimensional simulation, we use the same system size in the x and y directions as for the three-dimensional simulation, and the average number of particles per grid size at the neutral sheet is 972.These values are used to accurately calculate the magnetic energy dissipation process in the plasma sheet.
We adopt the Harris solution 41 for the pair plasma with the same mass m and uniform temperature T 0 in space and focused on the energy partitioning of an idealized magnetic reconnection in a collisionless plasma system.The magnetic field B(y) and plasma density N(y) are expressed as and respectively, where B G is the uniform guide magnetic field parallel to the z axis, and N b is the background uniform plasma density to demonstrate continuous plasma injection from the outside plasma sheet to the reconnection downstream.The background plasma density N b is set as 5% of the maximum plasma density N(0) at y = 0.The relationship between the pressure balance B 2 0 /8π = 2N 0 T 0 and force balance 2T 0 /λ = |e|u d B 0 /c is satisfied, where u d is the drift velocity, and the initial electric field E is zero.
In this paper, we fix the plasma temperature T 0 /mc 2 = 1, and investigate the energy partitioning between the thermal and nonthermal populations by changing only the magnitude of the guide magnetic field B G /B 0 from 0 to 5/2.The magnetization parameter σ = B 2 0 /(8πN 0 mc 2 ) = 2 if we use the number density for the value of the Harris plasma sheet N 0 and the magnetic field for the value of the anti-parallel magnetic field component B 0 .The Alfvén speed v A is expressed as v A = c σ/(1 + σ) in the cold plasma limit, where the density and magnetic field are used for the values of the central Harris plasma sheet and outside plasma sheets, respectively.The inertia length c/ω p in the plasma sheet is 6.5 grid cell, and the ratio of the gyro-radius r g for the Harris magnetic field B 0 and the thickness of the plasma sheet λ is set r g /λ = 0.45.The gyro-radius r g decreases with an increase in the guide magnetic field.The drift speed u d is set u d /c = 0.52.The red and blue isosurfaces are for the high and low plasma densities in the plasma sheet, respectively, and the green lines are the magnetic field lines.Three colored contour plots outside the box show the cross section of the plasma density at x = 0, z/λ = 80.6 and y = 0, respectively.One can clearly see that the plasma density structure becomes turbulent for the guide-field reconnection in panel (a), but not for the anti-parallel magnetic field reconnection in panel (b).The xy plane cross section above the box in panel (b) shows a large magnetic island around x/λ = 150 and three small magnetic islands around x/λ = 20, 220, and 310, while the xy plane cross section in panel (a) shows a large magnetic island around x/λ = 200 surrounded by many other turbulent density structures.A similar structure can also be seen in the yz plane cross section.The dynamically turbulent structure seen in the 3D guide-field reconnection is basically the same as that discussed by Daughton et al. 37 .
The bottom two panels (c) and (d) show the comparison of the plasma density between 3D and 2D guide-field reconnection with B G /B 0 = 1.Panel (c) is obtained by integrating the plasma density along the z-axis in panel (a), and panel (d) is the result under the same plasma parameters as panel (a) by using 2D PIC code.Comparing the density structure between the cross section (a) and the integrated structure (c), it can be seen that the patchy and turbulent density structure seen in panel (a) is more or less smeared out in the integrated density in panel (c), and the density in panel (c) is similar to the result of the 2D simulation in panel (d).However, we will show later that the patchy density structure plays an important role in particle acceleration.
The guide magnetic field plays an important role in the patch and turbulent magnetic reconnection, as shown by comparing panels (a) and (b) obtained in 3D simulations.Comparing panels (a) and (c) for the guide-field reconnection, the 3D effect is another important condition for the patch and turbulent reconnection.The guide-field reconnection in 2D system can occur only at the neural sheet, while the guide-field reconnection in 3D system can occur not only at the neutral sheet but also outside the neutral sheet 36,37,39 .This is a reason why the patch and turbulent reconnection can occur under both the effects of 3D geometry and finite guide magnetic field.

III. DENSITY PROFILE AND TURBULENT WAVES
It is important to compare the density profile between 3D and 2D guide-field reconnection.Shown in Figures 2 (a  For reference, the density profile of the 3D anti-parallel magnetic field reconnection at the same time stage t(V A /λ) = 365, which is obtained from panel (b) in Figure 1, is also shown in panel (a) by the red dashed line in Figure 2. Comparing the red solid line with the red dashed line, we have to keep in mind that the growth rate of the anti-parallel magnetic field reconnection is larger than that of the guide magnetic field reconnection .In fact, we find that the plasma density of the 3D anti-parallel reconnection outside plasma sheet is slightly lower than the initial plasma density N b /N(0) = 0.05, because the reconnection proceeds faster and the outer propagating slow mode rarefaction wave pulls the outside plasma into the central plasma sheet.Nevertheless, the density profile with the anti-parallel reconnection (red dashed line) is narrower than that with the guide magnetic field (red solid line).In addition, the density profile around the neutral sheet shows a flat shape for the guide-field reconnection, suggesting strong diffusion in the central plasma sheet for the guide-field reconnection.
The turbulent and broad density profile in the 3D guide-field reconnection may be related to the excitation of the oblique tearing mode [42][43][44]47 . It s interesting to examine the wave spectrum associated with the turbulent reconnection.While the tearing mode propagates parallel to the x-axis for the anti-parallel magnetic field reconnection, the oblique tearing mode can be excited outside the neutral sheet for the guide-field reconnection to satisfy where k is the wave vector of the tearing mode.Since the ambient magnetic field B is a function of y, the propagation of the tearing mode also changes.Furthermore, when the sizes of these magnetic islands due to the oblique tearing mode/reconnection become large enough to overlap, the islands are destroyed and the particles begin to wander stochastically.This process can lead to the diffusion of the plasma density.In addition to the wave spectra of B y (k x , y, k z ) showing the oblique tearing mode excitation, we have also studied the wave spectra of the plasma density and the electric field E z , which can capture other instabilities such as the drift-kink mode 48,49 .The profile of these wave spectra in the nonlinear phase was basically the same as that shown in Figure 3.
However, in the linear phase for the weak guide magnetic field cases with B G /B 0 < 1/4, we clearly observed two unstable modes: one is the tearing mode with the maximum growth rate along the x-axis and the other is the drift-kink mode with the maximum growth rate along the z-axis (not shown here).
The intensity of these oblique waves is another key parameter controlling the formation of the broad density profile.Figure 4 shows the comparison of the waves spectra of B y (k x ) between the 3D guide-field reconnection in Figures 1 (a reconnection is stronger than that of the 3D reconnection.At the edge of the plasma sheet, however, the wave intensity of 2D reconnection becomes weaker, while the wave intensity of 3D reconnection at y = −20λ remains almost at the same level of the wave intensity at y = 0. Therefore, we think that the broad density profile is generated by the particle scattering by the oblique tearing mode existing in a wide range of the plasma sheet for the 3D guide-field reconnection.

IV. TIME EVOLUTION OF PLASMA HEATING
In the previous section, we discussed that the 3D guide-field reconnection generates the dynamic and turbulent structure in the plasma sheet.It is interesting to study whether or not the dimensionality of 3D/2D makes a difference to the plasma heating and particle acceleration.Figures 5 (a

Lorentz factor γ
The total kinetic energy is normalized to the initial total kinetic energy at t(V A /λ) = 0.In our simulation study, the initial plasma parameters including the initial kinetic energy are all the same except for the magnitude of the guide magnetic field.In Figure 5, five different cases for the guide-field reconnection are plotted: the solid blue, purple, red, orange, and green lines show the cases with the guide magnetic field of B G /B 0 = 0, 1/4, 1/2, 1 and 2, respectively.
The kinetic energy at t(V A /λ) = 1000 becomes almost four times higher than that at the initial state for the case with no guide magnetic field, but the energy gain decreases with increasing the guide magnetic field for both 2D and 3D cases.For the guide magnetic field B G /B 0 = 2, the energy gain remains at a level of 50%.The linear growth rate of the collisionless tearing mode is known to decrease with increasing the guide magnetic field [42][43][44][45][46][47] , because the resonant particles around the magnetic diffusion region can be easily accelerated along the finite guide magnetic field.In fact, the resonant particles can be trapped by the guide magnetic field, and then the effective conductivity becomes high [42][43][44][45] .The tendency of the kinetic energy gain as a function of the guide magnetic field can be interpreted by the change in the growth rate of the oblique tearing mode.
Comparing the 3D time history in the left panel with the 2D time history in the right panel, we find several small differences: The 2D time history shows clearly several bumps associated with the coalescence of plasmoids, but the 3D time history does not show such a large oscillation in the time history of the energy.The onset of the plasma heating occurs at almost the same time stage t(V A /λ) ∼ 100 for 2D reconnection, but the onset of 3D reconnection with the weak guide magnetic field B G /B 0 < 1/4 occurs earlier at t(V A /λ) ∼ 70, because the drift-kink instability can be excited in the early time stage 18,48,49 .
Although, there are some differences between 3D and 2D, we think that the overall energy history is quite similar and there is no significant difference in the energy gain.

V. MODEL FITTING OF ENERGY SPECTRUM
In Paper I, we discussed the energy partitioning between thermal and nonthermal components during reconnection by using 2D simulation.Here we study it in the 3D system.Shown in Figure 6 is the energy spectra for B G /B 0 = 0, 1 and 2 integrated over the whole simulation domain.In our previous paper of 2D reconnection, we studied the energy spectrum downstream of the separatrix lines sandwiched by the last reconnected magnetic field lines, but in 3D reconnection we cannot define the separatrix surface to distinguish the downstream and the upstream of the reconnected magnetic field line.We then examine the energy spectra including both the preheated plasmas in the upstream and the heated/accelerated plasmas in the downstream.The top three panels are the results of 3D reconnection, while the bottom three are 2D reconnection, and the colored lines illustrate the time evolution of the energy spectra.The left two panels are the case of the anti-parallel reconnection with B G = 0, the middle panels are the case of the guide-field reconnection with B G /B 0 = 1, and the right panels are the case with B G /B 0 = 2.
The horizontal axis represents the particle kinetic energy (γ − 1) normalized by mc 2 and the vertical axis represents the number density N(γ)dγ, where γ = 1/ 1 − (v/c) 2 .The vertical scales of N(γ)dγ are shown in an arbitrary unit, but the ratio of the top to the bottom is set to five orders of magnitude.On the horizontal axis, the number of particles is counted by logarithmic binning in energy , so the one-count level is below the spectra shown here.The colored lines show the time evolution of the energy spectra, whose time stages are indicated by the colored bars on the right.The bluish and reddish colors correspond to the earlier and later stages, respectively.Since we have assumed a uniform plasma temperature The dashed lines in the bottom plot may represent an acceptable residual for reference.
in space for both the Harris plasma sheet and the background components in Equation ( 2), the dark blue spectra in Figure 6 show the initial plasma state with the combination of the non-drifting Maxwellian distribution function and the drifting Maxwellian distribution function.The red line is the final state when one large plasmoid is formed in the simulation box.
Over time, we clearly observe that the Maxwellian plasmas are gradually heated by reconnection heating processes, and these heated plasmas are simultaneously further accelerated to form a nonthermal component for both 3D and 2D reconnection.Roughly speaking, the spectra consist of the lower energy distribution close to the Maxwellian function and the higher energy population with a power-law-like distribution.The lower energy Maxwellianlike part remains to be almost the same temperature, but the higher energy power-law-like part evolves with time.These energy spectra for both 3D and 2D reconnections show similar time evolution.
In the late nonlinear stage at t/(V A /λ) = 1045 indicated by the red curve, we find the well-developed nonthermal energy spectra.We now examine the spectral behavior of the from the simulation data of N data (γ), described as where min() denotes a function that returns the minimum element from N data (γ) and N M +κ (γ).We observed that the residuals in wide energy ranges were less than 2 × 10  represent the plasma situated outside the plasma sheet before heating, but some plasmas may enter into the plasma sheet without any significant heating process 35 .More important is the behavior of the kappa temperature T κ in a strong guide magnetic field regime with B G /B 0 > 3/2.While the 2D kappa temperature T κ,2D increases with increasing the guide magnetic field, the 3D kappa temperature T κ.3D remains almost the same or decreases slightly with increasing the guide magnetic field.
In addition to the thermal plasma property mentioned above, the behavior of the kappa index shown in the middle panel (b) can provide insight into the energy partitioning between the thermal and the nonthermal plasmas.The 3D and 2D results are displayed as the red solid and the blue dashed lines, respectively, as a function of the guide magnetic field.The error bars, which are mostly smaller than the closed squares/circles, represent one-sigma standard deviation.We find that the kappa index for 3D remains around 3-4 regardless of the guide magnetic field.Similarly, for the 2D system, the kappa index also shows a hard spectrum for B G /B 0 < 1, but increases with an increase in the guide magnetic field.Since the thermal temperature of the kappa distribution T κ,2D for 2D in panel (a) increases with an increasing guide magnetic field, we believe that the 2D reconnection can effectively heat the plasma rather than accelerate the nonthermal particle.
To improve the clarity of the spectral behavior in 3D and 2D, Figure 8 has been expanded to show the transition energy from the thermal to the nonthermal regime in Figure 6.
The left panel (a) and the right panel (b) in Figure 8 correspond to panels (c) and (f) in  1) .The thick dashed lines in panels (a) and (b) correspond to those of the high-energy cyan dashed lines below the cutoff energies, which were γ cut = 14.1 in 3D and 22.9 in 2D.In the model-fitting of the 3D result in panel (a), the number density of the kappa distribution is larger than that of the Maxwellian, and the kappa temperature was T κ /mc 2 = 0.79, which is well separated from the the cutoff energy γ cut .Then the slope around 8 (γ − 1) 30 represents the nonthermal particle contribution, and the difference between the thick and thin curves, i.e., the kappa and the power-law distributions is minimal.
However, the slope in panel (b) of the 2D result, the slope below (γ − 1) ∼ 15 is still a part of the thermal contribution of the kappa distribution, because of the kappa temperature T κ /mc 2 = 1.65.Furthermore, since the number density of the kappa distribution is smaller than that of the Maxwellian distribution, we find that the slope indicated by the red line around 5 (γ − 1) 30 comes from both the contributions of the Maxwellian distribution of the lower-energy cyan curve and the kappa distribution of the higher-energy cyan curve.
In addition, the cutoff energy γ cut = 22.9 is not necessarily well separated from the kappa FIG. 8. Expanded energy spectra for panels (c) and (f) in Figure 6, and the format is basically the same as shown in Figure 6.Panels (a) and (b) are 3D and 2D results for B G /B 0 = 2, respectively.
For reference, the thick and thin dashed line show the kappa distribution function and the powerlaw function, respectively.
temperature of T κ /mc 2 = 1.65.These spectra illustrate the challenge of determining the nonthermal power-law index, when the separation between the thermal and the nonthermal populations is poor.
Panel (c) in Figure 7 shows the fraction of the nonthermal number density (ε den ) and the nonthermal energy density (ε ene ) as the function of the guide magnetic field B G /B 0 , where ε den and ε ene are defined as, and where N M κ (γ) represents the fraction of the Maxwellian distribution function in the κ distribution function (see Appendix for the definition of N M κ ).Therefore, the numerator of N κ (γ) − N M κ (γ) corresponds only to the nonthermal power-law component.In Figure 7 the red solid lines show the results of the 3D reconnection, while the blue dashed lines are the results of the 2D reconnection.
In the weak guide magnetic field regime of B B /B 0 < 1/4, the nonthermal energy efficiency is significantly large with ε ene ∼ 0.9, and the nonthermal number density ε den is about 0.5.However, the nonthermal efficiency for both 3D and 2D reconnections decreases with increasing the guide magnetic field B G .For example, in 2D reconnection the nonthermal energy density ε ene,2d decreases from 0.9 at B G /B 0 = 0 to 0.27 at B G /B 0 = 5/2.Although 3D reconnection shows the similar behavior of decreasing nonthermal particle production with increasing the guide magnetic field, the efficiency of the nonthermal production is much higher than that of 2D.The nonthermal energy density of ε ene,3d still holds about 0.6.
Therefore, we find that the strong guide-field reconnection suppresses the conversion of the magnetic energy into the kinetic energy, but the nonthermal production rate against the total plasma heating is not necessarily low in 3D reconnection.
Looking at the time history of the total kinetic energy in Figure 5, the energies increase gradually with time.At the final time stage at t(V A /λ) = 1045, however, we have confirmed that one large magnetic island is formed in the system even for the strong guide magnetic field of B G /B 0 = 2 and 5/2.In addition, from the time history of the energy spectra in Figure 6, we think that the final spectra shown by the red curves seem to reach the almost final stable state.Therefore, the nonthermal efficiencies of ε den and ε ene estimated from the energy spectra are almost the final states in our simulation box.We conclude that 3D reconnection can produce more nonthermal particles compared to 2D reconnection and can maintain a hard energy spectrum even for a strong guide-field reconnection.
We discussed that the patch and turbulent waves generated in 3D reconnections can enhance the efficiency of the nonthermal particle acceleration and that a hard nonthermal energy spectrum can be maintained even in the presence of a strong magnetic field.In order to study the effect of the plasma turbulence for a strong guide-field reconnection, we show the comparison of the ratio of the wave spectra of 3D to 2D reconnection for two cases of the guide-filed B G /B 0 = 1 and 2 in  spectrum is smaller than those of 2D spectrum, because the ratio is smaller than unity.At the edge of the plasma sheet y = 20λ in Panel (c), the 3D amplitudes becomes larger than 2D amplitudes in a wide range of the wavenumber.The behavior is consistent with Panels (a) and (b) in Figure 4, even though the time stages of the spectra are not necessarily same.
Let us look at the case of B G /B 0 = 2 (blue line).We find that the 3D turbulence of the B G /B 0 = 2 is stronger than that of B G /B 0 = 1 from the neutral sheet at y = 0 to the edge of the reconnecting plasma sheet at y = 20λ.This result supports the importance of the turbulence in the 3D guide-field reconnection.

VII. DISCUSSIONS AND CONCLUSION
In the present work, we have investigated the efficiency of the nonthermal particle acceleration during a relativistic reconnection with a guide magnetic field by using 3D particle-incell simulations, and have compared the efficiency with that obtained from 2D simulations.
For simplicity, we assumed a pair plasma with electrons and positrons, and the initial plasma temperature is fixed to be T 0 /mc 2 = 1.Under this assumption, we studied how the nonthermal particle acceleration is modified as the function of the guide magnetic field B G .So far the particle acceleration process during reconnection has been extensively investigated by many authors, but these studies are limited to the cases with the anti-parallel magnetic field or with a weak guide magnetic field of B G /B 0 1 [37][38][39]50,51 . In ths paper, we have studied the strong guide-field reconnection up to B G /B 0 = 5/2.Our simulation study confirmed that the efficiency of the nonthermal particle acceleration decreases with the increase of the guide magnetic field as expected before, but more importantly the energy spectrum with a strong guide magnetic field in the 3D simulation shows a hard energy spectrum with the kappa index of κ ∼ 4.This is contracted to the 2D system where the spectral index increases with increasing the guide magnetic field, and the kappa index becomes about 8 for In this paper, we have studied the case of T 0 /mc 2 = 1 and σ = 2.It is important to study the σ dependence from the nonrelativistic to the relativistic regime, but in our previous studies of two-dimensional reconnection 34,35 , the efficiency of the plasma heating and the particle acceleration does not change in the range of 1 < σ < 10.
The 3D guide-field reconnection inherently has the nature of the patch and turbulent reconnection, as multiple reconnections can occur from the neutral sheet to the outer plasma sheet.Over time, the magnetic islands generated in the different locations can overlap and the coherent island structure seen in the 2D system can be destroyed 37 .The turbulent wave spectrum can form over a wide region of the plasma sheet.These turbulent wave can enhance the efficiency of the nonthermal acceleration by the 2nd-order Fermi acceleration process.If the coherent structure of the magnetic island still remains after the overlapping of the multiple islands, the multiple scattering among the magnetic islands can generate the nonthermal particle through the 1st-order Fermi acceleration process 53 .In any case, the generation of the patchy and turbulent reconnection can improve the efficiency of the particle acceleration.
We have discussed the efficiency of particle acceleration enhanced by the patchy and turbulent nature of 3D reconnection as one of the additional acceleration processes, and we have not necessarily discussed the details of the particle acceleration mechanisms that have been studied by many people before.The main acceleration in relativistic reconnection discussed earlier is the direct acceleration in and around the magnetic diffusion region, because the relative size of the magnetic diffusion region with E > B supported by an inductive electric field is not necessarily small for a relativistic reconnection 15 .In addition to this direct acceleration, the particles ejected from the diffusion region collide with the magnetic field pile-up region formed in front of the magnetic island, and the particles may be accelerated by the gradient B drift and/or the magnetic curvature drift effect 9 .During the nonlinear time evolution, the particles trapped inside the magnetic islands can be accelerated by the 1st-order Fermi process in conversing reconnection flows 55 .During the magnetic coalescence, a secondary diffusion region appears at the interface of two magnetic islands with a stronger reconnection electric field, and then the direct acceleration can occur 52 .
Other various types of plasma instabilities developed in the plasma sheet can contribute to the energization of particles.Recently, the disruption of the magnetic flux tube by the kink instability in the 3D reconnection with a weak guide magnetic field is pointed out as another acceleration process 39 .
We have discussed that 3D reconnection maintains a hard spectrum with κ ∼ 4, even though the kappa index of 2D guide-field reconnection increases with increasing the guide magnetic field.As we can see from the energy spectra in Figure 8, however, the separation of the nonthermal cutoff energy of f cut from the thermal energy is not large enough for the strong guide magnetic field B G /B 0 = 2.We need to carefully analyse the spectral index by paying attention to the separation between the thermal and the nonthermal parts.In our model fitting, we have distinguished the nonthermal kappa distribution function from the thermal Maxwellian function, and we believe that our analysis has been done carefully under the assumption of the model fitting function of N M +κ .We think that it would be better to verify this result by using a more larger simulation system where we can capture the evolution of the nonthermal component.To improve the dynamic energy range, we will probably need an extremely large supercomputer memory with a system size 10 3 times larger in a 3D simulation.
To distinguish thermal plasma heating from nonthermal particle acceleration in our model fitting, we introduce the idea of the Maxwellian part of the kappa distribution function N M κ (γ) = lim κ→∞ N κ (γ).The kappa distribution consists of a thermal Maxwellian distribution at low energies and a nonthermal population approximated by a power-law function at high energies with N κ (γ) ∝ γ −κ+1 , where a power-law index s = κ − 1 for γ ≫ T κ /mc 2 .We defined the thermal part of the kappa distribution by taking the limit of κ → ∞, keeping the other fitted parameters the same, because a kappa distribution function approaches a Maxwellian distribution as κ → ∞, as described below:

FIG. 1 .
FIG. 1. Structure of magnetic reconnection with and without guide magnetic field.(a) 3D reconnection structure with the guide magnetic field B G /B 0 = 1 at t(V A /λ) = 366.The red and blue isosurfaces show high and low plasma densities, respectively.The green lines are magnetic field lines, and two-dimensional contours in xy-, yz-and xz-plane are the slice at z/λ=80.7,x=0, and the neutral sheet y = 0, respectively.(b) 3D reconnection structure without the guide magnetic field.(c) The integrated plasma density along the z-axis of case (a).(d) 2D reconnection with the guide magnetic field B G /B 0 = 1 with the same parameter as case (a), and the white lines are the magnetic field lines.

Figures 1 (
Figures 1 (a) and (b) in the top two panels show the density structures obtained for the 3D system with and without the guide magnetic field B G .The left panel (a) is the case with B G /B 0 = 1, while the right panel (b) is the case without the guide magnetic field.
) and (b) are the plasma density integrated along the x-axis at t(V A /λ) = 365 and at t(V A /λ) = 523, respectively.The red solid and blue dashed curves are the 3D and 2D profiles, respectively.With time, the high-density plasma sheet becomes thicker for both 2D and 3D reconnection, but we can clearly observe that 3D reconnection produces the wider plasma sheet compared to 2D reconnection.The small bump/shoulder appearing around y/λ ∼ ±(25 − 30) for 2D reconnection in panel (a) corresponds to the separatrix of the magnetic island.In the later time stage at t(V A /λ) = 523, the shoulder is

FIG. 2 .
FIG. 2. Comparison of the density profiles of the guide-field reconnection with B G /B 0 = 1 in 3D reconnection (red solid line) and in 2D (blue dashed line).The left panel (a) is taken at time t/(V A /λ) = 365 and the right panel (b) at t(V A /λ) = 523.For reference, the density profile of the 3D reconnection with the anti-parallel reconnection B G = 0 is shown as the red dashed line in the panel (a).
The top three panels (a), (b) and (c) in Figure3show the Fourier spectra of B y (k x , y, k z ) obtained for the guide-field reconnection with B G /B 0 = 1 at three different positions of y/λ = 0, −10 and −20 at the time of t(V a /λ) = 366 , respectively.The amplitudes are normalized to the Harris magnetic field B 0 .The bottom three panels of (d), (e) and (f) are those for the anti-parallel magnetic field reconnection with B G /B 0 = 0. From the top three panels of (a), (b) and (c), we see that the excited waves with long wavelengths kλ < 1 propagate along the x-axis at the neutral sheet at y = 0, and the waves at y/λ = −10 and −20 propagate obliquely away from the x-axis.The propagation angle at y/λ = −20 is more oblique than that at y/λ = −10, and we think that these turbulent and broadband waves are generated to satisfy the resonance condition of k•B = 0 for the oblique tearing mode.In fact, the white lines depicted in Figure3(b) and (c) are the slopes calculated from the average magnetic field magnitudes B x (y) = B x (x, y, z)dxdz and B z (y) = B z (x, y, z)dxdz, indicating that the propagation angle satisfies k • B = 0 37 .In contrast, the anti-parallel magnetic field reconnection in the bottom three panels of (d), (e), and (f) does not show any behavior of the oblique tearing mode.Therefore, we can understand that the patchy density structure of the guide-field reconnection in Figure 2 (a) is due to the oblique tearing mode.

FIG. 4 .
FIG. 4. Comparison of the wave spectra of the guide-field reconnection with B G /B 0 = 1 in 3D reconnection (red solid line) and in 2D (blue dashed line).The left panel (a) is taken at the neutral sheet y = 0 and the right panel (b) at y = −20λ.

FIG. 5 .
FIG. 5. Time evolution of the total kinetic energy in the system as a function of the guide magnetic field B G /B 0 = 0, 1/4, 1/2, 1 and 2. Panel (a) on the left and panel (b) on the right are the 3D and 2D results, respectively.

FIG. 6 .
FIG. 6.Comparison of 2D and 3D energy spectra for the entire reconnection region.The top three panels are the 3D reconnection for (a) B G /B 0 = 0, (b) B G /B 0 = 1 and (c) B G /B 0 = 2, and the bottom three panels are the 2D reconnection for (d) B G /B 0 = 0, (e) B G /B 0 = 1 and (f) B G /B 0 = 2.The colored lines indicate the time evolution of the spectra, whose time stage is indicated in the right side bar.The dark-blue line indicates the initial stage with the Maxwell distribution, and the red line indicates the final stage where a large plasmoid has formed in the system.Two cyan dashed lines represent the model fitting result for the composite function of the Maxwell and kappa distributions.The low-energy region corresponds to the Maxwellian part and the high energy region is the kappa distribution part.The lower panels show the residuals of the model fitting, which indicate the difference between the model fitting and the simulation spectrum.
red curve by using a least-squares fit to a model function.As already discussed in Paper I, the final-stage spectra can be well fitted by a composite model spectrum N M (γ) + N κ (γ) of a Maxwell distribution function N M (γ) and a kappa distribution function N κ (γ).The details of the model fitting function are given in Appendix.In the six panels in Figure 6, two cyan lines are the model fitting results for N M +κ (γ) = N M (γ) + N κ (γ) at the final stage.The low energy part corresponds to the Maxwellian function N M (γ), and the high-energy part is the kappa distribution function N κ (γ).The bottom six panels below the energy spectra show the residual of the model fitting of N M +κ (γ)

− 1
as indicated by the dashed lines, and we think that the model fitting of N M +κ (γ) described a very good model.As the fitting curve of N M +κ (γ) overlaps with the red curve of the final stage, the curve of N M +κ (γ) is not shown here.We confirmed that not only the above six energy spectra shown here but also other simulation results can be well approximated by the model function consisting of the Maxwellian and kappa distributions N M +κ (γ).

FIG. 7 .
FIG. 7. Comparison of 3D and 2D energy spectra.(a) Maxwellian temperature T M (dashed lines) and kappa temperature T κ (solid lines), (b) kappa index for 2D (blue dashed line) and 3D (red solid line), and (c) the nonthermal energy ε ene and the number density ε den as a function of the guide magnetic field B G /B 0 .

Figure 6 .
Figure 6.The red curves of the final state in 3D and 2D look similar, but the spectral behavior of the model fitting above was quite different.To see the difference, the thick black dashed line shows the kappa distribution function neglecting the cutoff contribution, i.e., γ γ 2 − 1 (1 + (γ − 1)/(κT κ /mc 2 )) −(κ+1) , and the thin dashed line represents a power-law function γ −(κ−1) .The thick dashed lines in panels (a) and (b) correspond to those of the

Figure 9 .
The vertical axis in Panels (a), (b), and (c) shows B 3D y (k x )/B 2D y (k x ) for the cases of B G /B 0 = 1 (red line) and B G /B 0 = 2 (blue line).Panels (a), (b), and (c) show the ratio at the neutral sheet y = 0, at y = 10λ, and at y = 20λ, respectively.The ratios of B b /B 0 = 1 and B G /B 0 = 2 are taken at t(V A /λ) = 522 and at t(V A /λ) = 914, when the total electric field energy integrated in the whole system becomes almost maximum.The peak electric field energies normalized by the initial kinetic energy were 0.1 − 0.25.For the case of B G /B 0 = 1 at the neutral sheet y = 0, the amplitudes of 3D

FIG. 9 .
FIG. 9. Comparison of the ratio of the 3D to 2D wave spectra as a function of the wavenumber k x λ.The red solid line is the case of the guide-field B g /B 0 = 1 and the blue solid line is the case of B g /B 0 = 2.The left panel (a) is taken at the neutral sheet y = 0, the middle panel (b) at y = 10λ, and the right panel (c) at y = 20λ.