No Arabic abstract
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.
Data from the NASA Magnetospheric Multiscale (MMS) mission are used to investigate asymmetric magnetic reconnection at the dayside boundary between the Earths magnetosphere and the solar wind (the magnetopause). High-resolution measurements of plasmas, electric and magnetic fields, and waves are used to identify highly localized (~15 electron Debye lengths) standing wave structures with large electric-field amplitudes (up to 100 mV/m). These wave structures are associated with spatially oscillatory dissipation, which appears as alternatingly positive and negative values of J dot E (dissipation). For small guide magnetic fields the wave structures occur in the electron stagnation region at the magnetosphere edge of the EDR. For larger guide fields the structures also occur near the reconnection x-line. This difference is explained in terms of channels for the out-of-plane current (agyrotropic electrons at the stagnation point and guide-field-aligned electrons at the x-line).
Based on global hybrid simulation results, we predict that foreshock turbulence can reach the magnetopause and lead to reconnection as well as Earth-sized indents. Both the interplanetary magnetic field (IMF) and solar wind are constant in our simulation, and hence all dynamics are generated by foreshock instabilities. The IMF in the simulation is mostly Sun-Earth aligned with a weak northward and zero dawn-dusk component, such that subsolar magnetopause reconnection is not expected without foreshock turbulence modifying the magnetosheath fields. We show a reconnection example to illustrate that the turbulence can create large magnetic shear angles across the magnetopause to induce local bursty reconnection. Magnetopause reconnection and indents developed from the impact of foreshock turbulence can potentially contribute to dayside loss of planetary plasmas.
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.
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.