No Arabic abstract
The spreading of the X-line out of the reconnection plane under a strong guide field is investigated using large-scale three-dimensional (3D) particle-in-cell (PIC) simulations in asymmetric magnetic reconnection. A simulation with a thick, ion-scale equilibrium current sheet (CS) reveals that the X-line spreads at the ambient ion/electron drift speeds, significantly slower than the Alfven speed based on the guide field $V_{Ag}$. Additional simulations with a thinner, sub-ion-scale CS show that the X-line spreads at $V_{Ag}$ (Alfvenic spreading), much higher than the ambient species drifts. An Alfvenic signal consistent with kinetic Alfven waves develops and propagates, leading to CS thinning and extending, which then results in reconnection onset. The continuous onset of reconnection in the signal propagation direction manifests as Alfvenic X-line spreading. The strong dependence on the CS thickness of the spreading speeds, and the X-line orientation are consistent with the collisionless tearing instability. Our simulations indicate that when the collisionless tearing growth is sufficiently strong in a thinner CS such that $gamma/Omega_{ci}gtrsimmathcal{O}(1)$, Alfvenic X-line spreading can take place. Our results compare favorably with a number of numerical simulations and recent magnetopause observations. A key implications is that the magnetopause CS is typically too thick for Alfvenic X-line spreading to effectively take place.
The orientation and stability of the reconnection x-line in asymmetric geometry is studied using three-dimensional (3D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x-line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2D simulations indicate that this x-line orientation maximizes the reconnection rate. The divergence of the non-gyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2D counterpart. We then design 3D simulations with one dimension being short to fix the x-line orientation, but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. The development of oblique structure easily leads to turbulence inside small periodic systems.
A new look at the structure of the electron diffusion region in collisionless magnetic reconnection is presented. The research is based on a particle-in-cell simulation of asymmetric magnetic reconnection, which include a temperature gradient across the current layer in addition to density and magnetic field gradient. We find that none of X-point, flow stagnation point, and local current density peak coincide. Current and energy balance analyses around the flow stagnation point and current density peak show consistently that current dissipation is associated with the divergence of nongyrotropic electron pressure. Furthermore, the same pressure terms, when combined with shear-type gradients of the electron flow velocity, also serve to maintain local thermal energy against convective losses. These effects are similar to those found also in symmetric magnetic reconnection. In addition, we find here significant effects related to the convection of current, which we can relate to a generalized diamagnetic drift by the nongyrotropic pressure divergence. Therefore, only part of the pressure force serves to dissipate the current density. However, the prior conclusion that the role of the reconnection electric field is to maintain the current density, which was obtained for a symmetric system, applies here as well. Finally, we discuss related features of electron distribution function in the EDR.
Electron dynamics surrounding the X-line in magnetopause-type asymmetric reconnection is investigated using a two-dimensional particle-in-cell simulation. We study electron properties of three characteristic regions in the vicinity of the X-line. The fluid properties, velocity distribution functions (VDFs), and orbits are studied and cross-compared. On the magnetospheric side of the X-line, the normal electric field enhances the electron meandering motion from the magnetosheath side. The motion leads to a crescent-shaped component in the electron VDF, in agreement with recent studies. On the magnetosheath side of the X-line, the magnetic field line is so stretched in the third dimension that its curvature radius is comparable with typical electron Larmor radius. The electron motion becomes nonadiabatic, and therefore the electron idealness is no longer expected to hold. Around the middle of the outflow regions, the electron nonidealness is coincident with the region of the nonadiabatic motion. Finally, we introduce a finite-time mixing fraction (FTMF) to evaluate electron mixing. The FTMF marks the magnetospheric side of the X-line, where the nonideal energy dissipation occurs.
We perform a theoretical and numerical study of anti-parallel 2D magnetic reconnection with asymmetries in the density and reconnecting magnetic field strength in addition to a bulk flow shear across the reconnection site in the plane of the reconnecting fields, which commonly occurs at planetary magnetospheres. We predict the speed at which an isolated X-line is convected by the flow, the reconnection rate, and the critical flow speed at which reconnection no longer takes place for arbitrary reconnecting magnetic field strengths, densities, and upstream flow speeds, and confirm the results with two-fluid numerical simulations. The predictions and simulation results counter the prevailing model of reconnection at Earths dayside magnetopause which says reconnection occurs with a stationary X-line for sub-Alfvenic magnetosheath flow, reconnection occurs but the X-line convects for magnetosheath flows between the Alfven speed and double the Alfven speed, and reconnection does not occur for magnetosheath flows greater than double the Alfven speed. We find that X-line motion is governed by momentum conservation from the upstream flows, which are weighted differently in asymmetric systems, so the X-line convects for generic conditions including sub-Alfvenic upstream speeds. For the reconnection rate, while the cutoff condition for symmetric reconnection is that the difference in flows on the two sides of the reconnection site is twice the Alfven speed, we find asymmetries cause the cutoff speed for asymmetric reconnection to be higher than twice the asymmetric form of the Alfven speed. The results compare favorably with an observation of reconnection at Earths polar cusps during a period of northward interplanetary magnetic field, where reconnection occurs despite the magnetosheath flow speed being more than twice the magnetosheath Alfven speed, the previously proposed suppression condition.
Using fully kinetic simulations, we study the x-line orientation of magnetic reconnection in an asymmetric configuration. A spatially localized perturbation is employed to induce a single x-line, that has sufficient freedom to choose its orientation in three-dimensional systems. The effect of ion to electron mass ratio is investigated, and the x-line appears to bisect the magnetic shear angle across the current sheet in the large mass ratio limit. The orientation can generally be deduced by scanning through corresponding 2D simulations to find the reconnection plane that maximizes the peak reconnection electric field. The deviation from the bisection angle in the lower mass ratio limit can be explained by the physics of tearing instability.