Non-evolutionarity of a reconnecting current sheet as a cause of its splitting into MHD shocks

Numerical simulations of the magnetic reconnection process in a current sheet show that, in some cases, MHD shocks appear to be attached to edges of the sheet. The appearance of the shocks may be considered to be a result of splitting of the sheet. In the present paper we suppose that this splitting takes place in consequence of non-evolutionarity of the reconnecting current sheet as a discontinuity. The problem of time evolution of small perturbations does not have a unique solution for a non-evolutionary discontinuity, and it splits into other (evolutionary) discontinuities. Such an approach allows us to determine conditions under which the splitting of the-sheet occurs. The main difficulty of this approach is that a current sheet is not reduced to a classified 1D discontinuity, because inhomogeneity of flow velocity inside the sheet is two-dimensional. To formulate the non-evolutionarity problem, we solve the linear MHD equations inside and outside the sheet and deduce linearized 1D boundary conditions at its surface. We show that for large enough conductivity, small perturbations exist which interact with the sheet as with a discontinuity. Then we obtain a non-evolutionarity criterion, with respect to these perturbations, in the form of a restriction on the flow velocity across the surface of the sheet.     

Solar wind access to the plasma sheet along the flanks of the magnetotail

Low-energy particle trajectories in an idealized magnetotail magnetic field are investigated to determine the accessibility of magnetosheath protons and electrons to the plasma sheet along the flanks of the tail magnetopause. The drift motion of the positively (negatively) charged particles incident on the dawn (dusk) magnetotail flank causes such particles to penetrate deeper into the magnetotail. For certain combinations of particle energy, incident velocity vector and initial penetration point on the tail magnetopause, the incident particles can become trapped in the plasma sheet, after which their net drift motion then provides a current capable of supporting the entire observed magnetotail field. The results further indicate that the bulk of the solar wind plasma just outside the distant tail boundary, which streams preferentially in a direction along the magnetopause away from the Earth at velocities around 400 km s^?1, can be caught up in the tail if the initial penetration point is within about 2R_E, of the quasi-neutral sheet. It is suggested that a large fraction of the magnetotail plasma is composed of former solar wind particles which have penetrated the magnetospheric boundary at the tail flanks.     

A note on the motion of charged particles in one-dimensional magnetic current sheets

By considering the integrals of the motion of charged particles moving in one-dimensional current sheets, a simple and exact proof is given that particles which are either magnetically or electrostatically trapped about such a current sheet exhibit zero net drift. The transition to the special case of a strictly neutral sheet, a limit remaining unclear from previous studies, is also elucidated. Finally, the relationship between the results and existing self-consistent current sheet solutions is discussed.     

Origin of Mercury’s double magnetopause: 3D hybrid simulation study with A.I.K.E.F.

During the first and second Mercury flyby the MESSENGER spacecraft detected a dawn side double-current sheet inside the Hermean magnetosphere that was labeled the “double magnetopause” (Slavin, J.A. et al. [2008]. Science 321, 85). This double current sheet confines a region of decreased magnetic field that is referred to as Mercury’s “dayside boundary layer” (Anderson, M., Slavin, J., Horth, H. [2011]. Planet. Space Sci.). Up to the present day the double current sheet, the boundary layer and the key processes leading to their formation are not well understood. In order to advance the understanding of this region we have carried out self-consistent plasma simulations of the Hermean magnetosphere by means of the hybrid simulation code A.I.K.E.F. (Müller, J., Simon, S., Motschmann, U., Schüle, J., Glassmeier, K., Pringle, G.J. [2011]. Comput. Phys. Commun. 182, 946–966). Magnetic field and plasma results are in excellent agreement with the MESSENGER observations. In contrast to former speculations our results prove this double current sheet may exist in a pure solar wind hydrogen plasma, i.e. in the absence of any exospheric ions like sodium. Both currents are similar in orientation but the outer is stronger in intensity. While the outer current sheet can be considered the “classical” magnetopause, the inner current sheet between the magnetopause and Mercury’s surface reveals to be sustained by a diamagnetic current that originates from proton pressure gradients at Mercury’s inner magnetosphere. The pressure gradients in turn exist due to protons that are trapped on closed magnetic field lines and mirrored between north and south pole. Both, the dayside and nightside diamagnetic decreases that have been observed during the MESSENGER mission show to be direct consequences of this diamagnetic current that we label Mercury’s “boundary-layer-current“.     

Generalized adiabatic theory applied to the magnetotail current sheet

We use the generalized first adiabatic invariant, an extension of the magnetic moment for regions of large field gradients, to treat particles in the magnetotail current sheet. The equations of motion can be expressed in terms of drift parameters which vary slowly and smoothly at the drift rate, not at the gyration rate. The analysis leads to boundaries in phase space which form a generalized loss cone and separate particles drifting into and out of the layer from particles trapped within the layer. These boundaries can be used in the moment integrals for densities and currents when the drifting particles differ in temperature, or in other properties, from the trapped population, as has been suggested by observations. We give examples of how different kinds of particle orbits contribute to the spatial profiles of density and current and thus to the field structure of the current sheet. We find that the parallel pressure of the drifting particles must exceed the transverse pressure for self-consistent solutions to exist, and based on this result, we give examples of fully self-consistent solutions using bi-Maxwellian ion and Maxwellian electron distributions. We give a proof, using generalized adiabatic theory, of Cowley's (1978a) theorem that particles trapped in the current layer experience zero net drift.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

The ducting of wave energy by field-aligned current sheets

A theoretical treatment of the modes of oscillation of an idealized current sheet (in which there is no perpendicular temperature) is given. For the simple case of a monoenergetic current sheet computational results are presented. These results indicate the existence of two types of ducted mode which may have relevance to observation. The first of these is a ducted “whistler” mode and the second occurs at somewhat higher frequencies and at phase velocities comparable with the streaming velocity of the current carriers in the sheet region. A simple explanation of how the ducted “whistler” mode can arise is given.     

Modeling of Proton Acceleration in a Magnetic Island Inside the Ripple of the Heliospheric Current Sheet

Crossings of the heliospheric current sheet (HCS) at the Earth’s orbit are often associated with observations of anisotropic beams of energetic protons accelerated to energies from hundreds of keV to several MeV and above. A connection between this phenomenon and the occurrence of small-scale magnetic islands (SMIs) near reconnecting current sheets has recently been found. This study shows how pre-accelerated protons can be energized additionally due to oscillations of multiple SMIs inside the ripple of the reconnecting HCS. A model of the electromagnetic field of an oscillating 3D SMI with a characteristic size of ~0.001 AU is developed. A SMI is supposed to be bombarded by protons accelerated by magnetic reconnection at the HCS to energies from ~1keV to tens of keV. Numerical simulations have demonstrated that the resulting longitudinal inductive electric fields can additionally reaccelerate protons injected into a SMI. It is shown that there is a local “acceleration” region within the island in which particles gain energy most effectively. As a result, their average escape energies range from hundreds of keV to 2 MeV and above. There is almost no particle acceleration outside the region. It is shown that energies gained by protons significantly depend on the initial phase and the place of their entry into a SMI but weakly depend on the initial energy. Therefore, low-energy particles can be accelerated more efficiently than high-energy particles, and all particles can reach the total energy limit upon their escape from a SMI. It is also found that the escape velocity possesses a strong directional anisotropy. The results are consistent with observations in the solar wind plasma.

     

Plasma-sheet dynamics and magnetospheric substorms

We present a conceptual model of the formation of the plasma sheet and of its dynamical behavior in association with magnetospheric substorms. We assume that plasma mantle particles $\text{E}\text{～}\text{×}\text{B}\text{～}$ drift toward the current sheet in the center of the tail where they are accelerated by magnetic-field annihilation to form the plasma sheet. Because of the velocity-dependent access of mantle particles to the current sheet, we argue that the convection electric field and the corresponding rate of field annihilation decrease with increasing radial distance. As a consequence, there exists no steady-state configuration for the plasma sheet, which must instead shrink continuously in thickness until the near-earth portion of the current sheet is disrupted by the formation of a magnetic neutral line. The current-sheet disruption launches a large-amplitude hydromagnetic wave which is largely reflected from the ionosphere. The reflected wave sets the neutral line in motion away from the earth; the neutral line comes to rest at a distance (which we estimate to be a few hundred earth radii) where the incoming mantle particles enter the current sheet at the local Alfvén velocity. At this “Alfvén point” reconnection ceases and the thinning of the plasma sheet begins again. Within this model, the magnetospheric substorm (which is associated with the current-sheet disruption) is a cyclical phenomenon whose frequency is proportional to the rate of convection in the magnetospheric tail.     

Regular and chaotic dynamics in 3D reconnecting current sheets

We consider the possibility of particles being injected at the interior of a reconnecting current sheet (RCS), and study their orbits by dynamical systems methods. As an example we consider orbits in a 3D Harris type RCS. We find that, despite the presence of a strong electric field, a 'mirror' trapping effect persists, to a certain extent, for orbits with appropriate initial conditions within the sheet. The mirror effect is stronger for electrons than for protons. In summary, three types of orbits are distinguished: (i) chaotic orbits leading to escape by stochastic acceleration, (ii) regular orbits leading to escape along the field lines of the reconnecting magnetic component, and (iii) mirror-type regular orbits that are trapped in the sheet, making mirror oscillations. Dynamically, the latter orbits lie on a set of invariant KAM tori that occupy a considerable amount of the phase space of the motion of the particles. We also observe the phenomenon of 'stickiness', namely chaotic orbits that remain trapped in the sheet for a considerable time. A trapping domain, related to the boundary of mirror motions in velocity space, is calculated analytically. Analytical formulae are derived for the kinetic energy gain in regular or chaotic escaping orbits. The analytical results are compared with numerical simulations.     

POLAR 5—an electron accelerator experiment within an aurora. 3. Evidence for significant spacecraft charging by an electron accelerator at ionospheric altitudes

The POLAR 5 rocket experiment carried an electron accelerator on a “daughter” payload which injected a 0,1 A beam of 10 keV electrons in a pulsed mode every 410ms. With spin and precession, injections were made over a wide range of pitch angles. Measurements from a double probe electric field instrument and from particle detectors on the “mother” payload and from a crude R.P.A. on the “daughter” payload are interpreted to indicate that the “daughter” charges to a potential between several hundred volts and 1 kV. The neutralizing return current to the “daughter” is shown to be assymetrically distributed with the majority being collected from the direction of the beam. The additional electrons necessary to neutralize the daughter are thought to be produced and heated through beam-plasma interactions postulated by Maehlum et al. (1980b) and Grandal et al. (1980) to explain the particle and optical measurements. Significant electric fields emanating from the charged “daughter” and the beam are seen at distances exceeding 100 m at the “mother” payload.     

Generalized analytical models of Syrovatskii’s current sheet

Two-dimensional stationary magnetic reconnection models that include a thin Syrovatskii-type current sheet and four discontinuous magnetohydrodynamic flows of finite length attached to its endpoints are considered. The flow pattern is not specified but is determined from a self-consistent solution of the problem in the approximation of a strong magnetic field. Generalized analytical solutions that take into account the possibility of a current sheet discontinuity in the region of anomalous plasma resistivity have been found. The global structure of the magnetic field in the reconnection region and its local properties near the current sheet and attached discontinuities are studied. In the reconnection regime in which reverse currents are present in the current sheet, the attached discontinuities are trans-Alfvénic shock waves near the current sheet endpoints. Two types of transitions from nonevolutionary shocks to evolutionary ones along discontinuous flows are shown to be possible, depending on the geometrical model parameters. The relationship between the results obtained and numerical magnetic reconnection experiments is discussed.     

Electron surfing acceleration by electrostatic waves in current sheet

The electron surfing acceleration in the current sheet with perpendicular propagating electrostatic waves is studied using analytical theories and test particle simulations. The trapped electron moving with the phase velocity v _p of wave may be accelerated effectively in the outflow direction by force until the electron is de-trapped from the wave potential. A criterion K>0 for the electron surfing acceleration is obtained. The electron will escape from the boundary of current sheet quickly, if this criterion does not hold. The maximum velocity of surfing acceleration is about the same as the electric drift velocity. Superposed longitudinal magnetic field along the wave propagation is favorable for the electron surfing acceleration in the current sheet.      

Solar wind energy transfer regions inside the dayside magnetopause: Accelerated heavy ions as tracers for MHD-processes in the dayside boundary layer

Plasma and magnetic field data from PROGNOZ-7 have revealed that solar wind (magnetosheath) plasma elements may penetrate the dayside magnetopause surface and form high density regions with enhanced cross-field flow in the boundary layer.The injected magnetosheath plasma is observed to have an excess drift velocity as compared to the local boundary layer plasma, comprising both “cold” plasma of terrestrial origin and a hot ring current component. A differential drift between two plasma components can be understood in terms of a momentum transfer process driven by an injected magnetosheath plasma population. The braking action of the injected plasma may be described as a dynamo process where particle kinetic energy is transferred into electromagnetic energy (electric field). The generated electric field will force the local plasma to $\text{ε×}\text{B}\text{-drift}$, and the dynamo region therefore also constitutes an accelerator region for the local plasma. Whenever energy is dissipated from the energy transfer process (a net current is flowing through a load), there will also be a difference between the induced electric field and the $\text{v}\text{×}\text{B}$ term of the generator plasma. Thus, the local plasma will drift more slowly than the injected generator plasma.We will present observations showing that a relation between the momentum transferred, the injected plasma and the momentum taken up by the local plasma exists. For instance, if the local plasma density is sufficiently high, the differential drift velocity of the injected and local plasma will be small. A large fraction of the excess momentum is then transferred to the local plasma. Conversely, a low local plasma density results in a high velocity difference and a low fraction of local momentum transfer.In our study cases the “cold” plasma component was frequently found to dominate the local magnetospheric plasma density in the boundary layer. Accordingly, this component may have the largest influence on the local momentum transfer process. We will demonstrate that this also seems to be the case. Moreover we show that the accelerated “cold” plasma component may be used as a tracer element reflecting both the momentum and energy transfer and the penetration process in the dayside boundary layer.The high He⁺ percentage of the accelerated “cold” plasma indicates a plasmaspheric origin. Considering the quite high densities of energetic He⁺ found in the boundary layer, the overall low abundance of He⁺ (as compared to e.g. O⁺) found in the plasma sheet and outer ring current evidently reduces the importance of the dayside boundary layer as a plasma source in the large scale magnetospheric circulation system.     

Large-scale transport of plasma in the Earth's plasma sheet: comparative analysis for adiabatic and non-adiabatic cases

A system of multi-fluid MHD-equations is used to compare adiabatic and non-adiabatic transport of the energetic particles in the magnetospheric plasma sheet. A “slow-flow” approximation is considered to study large-scale transport of the anisotropic plasma consisting of energetic electrons and protons. Non-adiabatic transport of the energetic plasma is caused by scattering of the particles in the presence of both wave turbulence and arbitrary time-varying electric fields penetrating from the solar wind into the magnetosphere. The plasma components are devided into particle populations defined by their given initial effective values of the magnetic moment per particle. The spatial scales are also given to estimate the non-uniformity of the geomagnetic field along the chosen mean path of a particle. The latters are used to integrate approximately the system of MHD-equations along each of these paths. The behaviour of the magnetic moment mentioned above and of the parameter which characterizes the pitch-angle distribution of the particles are studied self-consistently in dependence on the intensity of non-adiabatic scattering of the particles. It is shown that, in the inner magnetosphere, this scattering influences the particles in the same manner as pitch-angle diffusion does. It reduces the pitch-angle anisotropy in the plasa. The state of the plasma may be unstable in the current sheet of the magnetotail. If the initial state of the plasma does not correspond to the equilibrium one, then, in this case, scattering influences the particles so as to remove the plasma further from the equilibrium state. The coefficient of the particle diffusion across the geomagnetic field lines is evaluated. This is done by employing the Langevin approach to take the stochastic electric forces acting on the energetic particles in the turbulent plasma into account. The behaviour of the energy density of electrostatic fluctuations in the magnetosphere is estimated.     

Estimates of the wind speeds required for particle motion on Mars

We have obtained estimates of the threshold wind speed V_gt near the top of the atmospheric boundary layer on Mars and of the rotation angle α between this wind velocity and the direction of the surface stress. this calculation has been accomplished by combining wind tunnel determinations of the friction velocity with semi-empirical theories of the Earth's atmospheric boundary layer. Calculations have been performed for a variety of values of the surface pressure, ground temperature, roughness height, boundary layer height, atmospheric composition atmospheric stability, particle density, particle diameter, and strength of the cohesive force between the particles.The curve of threshold wind speed as a function of particle diameter monotonically decreases with decreasing particle diameter for a cohesionless soil but has the classical $\text{U}$ shape for a soil with cohesion. Observational data indicate that the latter condition holds on Mars. Under “favorable” conditions minimum threshold wind speeds between about 50 and 100m/sec are required to cause particle motion. These minimum values lie close to the highest wind speeds predicted by general circulation models. Hence, particle motion should be an infrequent occurence and should be strongly correlated with nearness to small topographic features. The latter prediction is in accord with the correlation found between albedo markings and topographic obstacles such as craters. For equal wind speeds at the midpoint of the boundary layer, particle movement occurs more readily in general at night than during the day, more readily in the winter polar areas than the equatorial areas noon, and more readily for ice particles than for silicate particles.The boundary between saltating and suspendable particles is located at a particle diameter of about 100 μm. This value is close to the diameter at which the V_gt curve has its minimum. Hence, the wind can set directly into motion both saltating and larger-sized suspendable particles, but dust-storm-sized particles usually require impact by a saltating particle for motion to be initiated. Albedo changes occur most often in regions containing a mixture of dust-stoorm-sized particles and saltating particles. The threshold wind speed for surfaces containing large, nonerodible roughness elements can either be larger or smaller than the value for surfaces with only erodible material. The former condition for V_gt holds when the roughness height z₀ is less than about 1 cm and may be illustrated by craters that have experienced less erosion than their environs. The latter condition for V_gt may be partly responsible for albedo changes detected on the elevated shield volcano, Pavonis Mons. Values of the angle α generally lie between 10 and 30°. These figures place a modest limitation on the utility of surface albedo streaks as wind direction indicators.     

Stationary distribution function for non-adiabatic particle motion in a one dimensional current sheet

The problem of pitch angle scattering in field configurations similar to those found in the geomagnetic tail has been studied previously by Tsyganenko (1982). Tsyganenko used a scattering matrix to map pitch angle distributions through the current sheet. By using numerical solutions of the resulting integral equations he showed for weakly non-adiabatic particles the Stationary Distribution Function (SDF) was isotropic. Using his procedure the SDF was found to develop anisotropies with increasing non-adiabaticity. The work presented here shows analytically that for any degree of scattering the SDF must be isotropic for a general planar field reversal. Computations of particle trajectories have been used to verify some aspects of the analytic work.     

Azimuthal structure of the solar wind velocity andK-corona brightness in the heliospheric current sheet

The solar wind velocity near Earth shows systematic structure in and around the heliospheric current sheet. The solar wind velocity measurements at IMF sector boundary crossings at 1 AU during 1972–1977 have been used to infer the azimuthal structure of the solar wind velocity in the current sheet. We found that the solar wind velocity in the in-ecliptic portion of the current sheet varies from longitude to longitude, where it originates from the corona. Also, the yearly average value of solar wind velocity in the HCS is found to vary with the phase of the solar cycle; with a maximum value around 1974. TheK-corona brightness on the source surface corresponding to the IMF sector boundary crossings during the period of study also shows a similar but opposite pattern of variation when the data are averaged over a long period. However, this relation is not observed when we considered them individually. So, we conclude that there exists a longitudinal variation of solar wind velocity in the heliospheric current sheet.     

Pitch-angle scattering of energetic particles in the current sheet of the magnetospheric tail and stationary distribution functions

An investigation of pitch-angle scattering of energetic particles in magnetic field configurations with a current sheet similar to that observed in the geomagnetotail has been performed. The magnetic field model is specified by two parameters which are the current sheet thickness in units of particle gyroradius and the angle between the magnetic field lines and the sheet plane. Computations of a considerable number of trajectories (about 20,000 for each model case) has provided the possibility of obtaining the matrix of pitch-angle scattering and the corresponding kernel function of the integral equation for the stationary particle distribution function. Solution of this equation shows that isotropic distributions are formed only in the case of a sufficiently thick current sheet. Particle scattering in a thin field reversal region leads to the formation of an anisotropic stationary distribution. The results can be used for interpretation of the data on the spatial distribution of energetic particle fluxes in the near part of the magnetospheric tail and in the vicinity of the outer boundary of the radiation belt.     

Precipitation of charged particles by a parallel electric field

A time-dependent model of the effect of a parallel electric field on particle precipitation from a closed field-line has been constructed and the results are presented. A pattern of field-aligned pitch-angle distributions and energy peaks develops rapidly and then persists unchanged in shape while the intensity decreases for a time of the order of the bounce period of the energetic particles. It is shown that the structures in velocity space are created by the juxtaposition of particles from different source populations. Four sources are found to be sufficient to reproduce the principal features observed frequently by rockets and satellites. They are, a trapped plasma sheet distribution, a loss-cone partially filled by pitch-angle diffusion at the equator, cold ionospheric plasma which has flowed outward along the field line and particles backscattered from the precipitation into the atmosphere.The model develops density gradients and discontinuities far sharper than any observed, so that any parallel electric field actually occurring in an aurora must be accompanied by strong wave-particle interactions either as part of the accelerating mechanism or as a result of the density gradients produced by it.