Search

Separatrices of magnetic reconnection host intense perpendicular Hall electric fields. The fields are produced by the decoupling of the ion and electron components and are associated with the in-plane electrostatic potential drop between the inflow and outflow regions. The width of these structures is typically less than the ion inertial length, which is small enough to demagnetize ions as they cross the layer. We investigate ion acceleration at separatrices by means of 2D particle-in-cell simulations of magnetic reconnection for two limiting cases: (1) a "GEM-like" setup (here GEM stands for geospace environmental modeling reconnection challenge) with the lobe ion thermal velocity equal to the thermal velocity of the initial current sheet ions, which is comparable to the Alfven velocity and (2) a "cold" ion setup, in which the temperature of the background lobe ions is 1/100 of the initial current sheet temperature. The separatrix Hall electric field is balanced by the ion inertia term in the cold background simulations. The effect is indicative of the quasi-steady local perpendicular acceleration. The electric field introduces a cross field beam of unmagnetized particles, which makes the ion distribution function strongly non-gyrotropic and susceptible to sub-ion scale instabilities. This acceleration mechanism nearly vanishes in the hot ion background simulations. Our particle-in-cell simulations are complemented by one-dimensional test particle calculations. They show that the hot ion particles experience energy-scattering after crossing the accelerating layer, whereas cold ions are uniformly energized up to the energies comparable to the electrostatic potential drop between the inflow and outflow regions.

We study inflow density dependence of substructures within electron diffusion region (EDR) of collisionless symmetric magnetic reconnection. We perform a set of 2.5D particle-in-cell simulations which start from a Harris current layer with a uniform background density n(b). A scan of n(b) ranging from 0:02 n(0) to 2 n(0) of the peak current layer density (n(0)) is studied keeping other plasma parameters the same. Various quantities measuring reconnection rate, EDR spatial scales, and characteristic velocities are introduced. We analyze EDR properties during quasisteady stage when the EDR length measures saturate. Consistent with past kinetic simulations, electrons are heated parallel to the B field in the inflow region. The presence of the strong parallel anisotropy acts twofold: (1) electron pressure anisotropy drift gets important at the EDR upstream edge in addition to the E x B drift speed and (2) the pressure anisotropy term -del.P-(e)/(ne) modifies the force balance there. We find that the width of the EDR demagnetization region and EDR current are proportional to the electron inertial length similar to d(e) and similar to d(e)n(b)(0.22), respectively. Magnetic reconnection is fast with a rate of similar to 0.1 but depends weakly on density as similar to n(b)(-1/8). Such reconnection rate proxies as EDR geometrical aspect or the inflow-to-outflow electron velocity ratio are shown to have different density trends, making electric field the only reliable measure of the reconnection rate. Published under license by AIP Publishing., QC 20200214

We present a new model to describe the electron pressure nongyrotropy inside the electron diffusion region (EDR) in an antiparallel magnetic reconnection scenario. A combination of particle-in-cell simulations and analytical estimates is used to identify such a component of the electron pressure tensor in the rotated coordinates, which is nearly invariant along the outflow direction between the X line and the electron remagnetization points in the outer EDR. It is shown that the EDR two-scale structure (inner and outer parts) is formed due to superposition of the nongyrotropic meandering electron population and gyrotropic electron population with large anisotropy parallel to the magnetic field upstream of the EDR. Inside the inner EDR the influence of the pressure anisotropy can largely be ignored. In the outer EDR, a thin electron layer with electron flow speed exceeding the E x B drift velocity is supported by large-momentum flux produced by the electron pressure anisotropy upstream of the EDR. We find that this fast electron exhaust flow with |V(e)xB|>|E| is in fact a constituent part of the EDR, a finding which will steer the interpretation of the Magnetospheric Multiscale Mission (MMS) data., QC 20161219

The role of the dawn-dusk magnetic field component in stabilizing of the magnetotail flapping oscillations is investigated in the double-gradient model framework (Erkaev et al., Phys. Rev. Lett. 99, 235003 (2007)), extended for the magnetotail-like configurations with non-zero guide field By. Contribution of the guide field is examined both analytically and by means of linearized 2-dimensional (2D) and non-linear 3-dimensional (3D) MHD modeling. All three approaches demonstrate the same properties of the instability: stabilization of current sheet oscillations for short wavelength modes, appearing of the typical (fastest growing) wavelength lambda(peak) of the order of the current sheet width, decrease of the peak growth rate with increasing B-y value, and total decay of the mode for B-y similar to 0: 5 in the lobe magnetic field units. Analytical solution and 2D numerical simulations claim also the shift of lambda(peak) toward the longer wavelengths with increasing guide field. This result is barely visible in 3D simulations. It may be accounted for the specific background magnetic configuration, the pattern of tail-like equilibrium provided by approximated solution of the conventional Grad-Shafranov equation. The configuration demonstrates drastically changing radius of curvature of magnetic field lines, R-c. This, in turn, favors the "double-gradient" mode (lambda > R-c) in one part of the sheet and classical "ballooning" instability (lambda < R-c) in another part, which may result in generation of a "combined" unstable mode.

We present a numerical investigation of the double-gradient mode, which is believed to be responsible for the magnetotail flapping oscillations the fast vertical oscillations of the Earth's magnetotail plasma sheet (quasiperiod similar to 100 - 200 s). It is known that this mode has an unstable solution in the region of the tailward-growing normal magnetic field component. The kink branch of the mode is the focus of our study. The instability is studied using the magnetotail near-equilibrium configuration, fixed by the approximate solution of the Grad-Shafranov equation. The linear three-dimensional numerical analysis is complemented with full 3-D MUD simulations. The results of our linearized MHD code agree with the theory, and the growth rate is found to be close to the peak value provided by an analytical estimate. Also, the eigenfunctions, calculated analytically, are very similar to the perturbations obtained numerically. The full 3D MHD simulations are initialized with the numerically relaxed magnetotail equilibrium, similar to the linear code initial condition. The calculations show that the double-gradient mode is excited in a region of small radii of the magnetic field lines curvature, which is in accordance with the analytical predictions. In contrast to the linearized MHD simulations, non-local interactions are involved; hence, the overall growth rate turns out to be close to the theoretical estimate averaged over the computational domain., QC 20190205

The paper presents the detailed numerical investigation of the "double-gradient mode," which is believed to be responsible for the magnetotail flapping oscillations-the fast vertical (normal to the layer) oscillations of the Earth's magnetotail plasma sheet with a quasiperiod similar to 100-200 s. The instability is studied using the magnetotail near-equilibrium configuration. For the first time, linear three-dimensional numerical analysis is complemented with full 3-D MHD simulations. It is known that the "double-gradient mode" has unstable solutions in the region of the tailward growth of the magnetic field component, normal to the current sheet. The unstable kink branch of the mode is the focus of our study. Linear MHD code results agree with the theory, and the growth rate is found to be close to the peak value, provided by the analytical estimates. Full 3-D simulations are initialized with the numerically relaxed magnetotail equilibrium, similar to the linear code initial condition. The calculations show that current layer with tailward gradient of the normal component of the magnetic field is unstable to wavelengths longer than the curvature radius of the field line. The segment of the current sheet with the earthward gradient of the normal component makes some stabilizing effect (the same effect is registered in the linearized MHD simulations) due to the minimum of the total pressure localized in the center of the sheet. The overall growth rate is close to the theoretical double-gradient estimate averaged over the computational domain.