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.     

The self-consistent geomagnetic tail under static conditions

A self-consistent calculation of the magnetic field and plasma distribution in the magnetotail has been undertaken for static conditions. We find the best agreement with experimental observations by satellite in the tail for an isotropic particle pitch angle distribution, a slow decrease of magnetic field intensity as a function of distance from the Earth |x|^−0.3, and a northward field in the equatorial plane of about 1 gamma at the position of the lunar orbit. Knowing the field and plasma distribution we calculate the sources of the electrical current and find that the magnetic field in the tail can be supported entirely by solar wind drifting into the curved and diminishing magnetic field of the magnetotail. Furthermore the plasma present in the tail at any one instant is swept out in a very short time due to its large curvature drift velocity; the constant plasma sheet is maintained by constant renewal of entering solar wind plasma.     

A note on adiabatic solutions of the one-dimensional current sheet problem

It has recently been shown that adiabatic solutions of the one-dimensional current sheet problem exist provided that magnetically trapped particles are included in the model together with the current-carrying untrapped “beam” particles. We show here that a formulation of the problem in terms of particle velocity and pitch angle is advantageous, and we derive some general properties of the solutions. In particular it is shown that there is, in general, no discontinuity in the value of the particle distribution function ? across the boundary in velocity space between “beam” and trapped particles, but that there will be a discontinuity in the gradients of ?. An example is given in which the beam population is of bi-Maxwellian form at the outer boundary of the current sheet.     

Can gravitational effects damp Alfvén waves?

We show that Alfvén-gravity waves propagating in a gravitationally stratified atmosphere do not suffer damping as a result of the rate of working of the gravity drift current on the electric field of the waves. A self-consistent treatment involving conservation of total energy, Poynting's theorem, and the rate of working of the various drift currents on the electric field demonstrates that dissipation only arises from real dissipative processes such as Ohmic heating or viscous effects, otherwise the system is adiabatic.     

Interaction of energetic particles with HM-waves in the magnetosphere

Energetic particle response in electromagnetic fields of ULF HM-waves in the magnetosphere is reviewed. Pc4–5 geomagnetic pulsations observed at the synchronous altitude are classified into three types, in respect to their major magnetic field polarization in different directions, local time dependence, and different characteristics of accompanied flux modulations of energetic particles, i.e., two nearly transverse waves with the azimuthal and the radial polarization, and the compressional stormtime pulsations. Firstly, we formulate the drift kinetic theory of particle flux modulations under the constraint of the magnetic moment conservation. A generalized energy integral of the particle motion interacting with a ULF-wave with the three-dimensional structure propagating to the azimuthal direction is obtained in the L-shell coordinate of a mirror magnetic field. Its linearized form is reduced to the same form as the previously derived energy change, including the bounce-drift resonant interaction. It is shown that the perturbed guiding center distribution function of energetic particles consists of four contributions, the adiabatic mirror effect corresponding to pitch-angle change, the kinetic effects due to energy change and the accompanying L-shell displacement, and the bounceaveraged drift phase bunching. Secondly, the basic HM-wave modes constitutingcoupling ULF oscillations in non-uniform plasmas are discussed in different models of approach for different plasma states. The diamagnetic drift Alfvén wave and the compressional drift wave with a larger azimuthal mode number in a high-beta plasma are candidates for the stormtimes pulsations. The former is intrinsically a guided localized mode, while the latter is a non-localized mode. By making use of the above preparation, we apply the developed drift kinetic theory to interpret the phase relationships between the ion flux modulation and the geomagnetic pulsation in some selected examples of observations, demonstrating a fair agreement in theoretical results with the observations.     

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.     

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“.     

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.     

The effect of pressure anisotropy on the equilibrium structure of magnetic current sheets

The equilibrium structure of two-dimensional magnetic current sheets is investigated for systems in which the plasma pressure dominates the bulk flow energy, as appears appropriate for the quiet time plasmasheet in the geomagnetic tail. A simple model is studied in which the field is contained between plane parallel boundaries and varies exponentially along the system, while the plasma pressure is anisotropic, the anisotropy being arbitrary but constant along the centre plane. When the field is highly inflated by the plasma current it is found that adiabatic solutions exist only when the plasma pressure is close to isotropic. For the case P_∥ > P_⊥ it is argued that a thin, non-adiabatic current layer will in general form at the sheet centre, usually embedded within a much broader adiabatic current distribution. When P_⊥ > P_∥, a broad region of very depressed fields develops about the centre of the current sheet, terminated at its outer boundary by a spike in the current density. This central region becomes unstable to the mirror mode well before the limiting adiabatic solution is reached.     

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.     

Stability of Electrostatic Electron Cyclotron Waves in a Multi-Ion Plasma

We have studied the stability of the electrostatic electron cyclotron wave in a plasma composed of hydrogen, oxygen and electrons. To conform to satellite observations in the low latitude boundary layer we model both the ionic components as drifting perpendicular to the magnetic field. Expressions for the frequency and the growth rate of the wave have been derived. We find that the plasma can support electron cyclotron waves with a frequency slightly greater than the electron cyclotron frequency ω _ce ; these waves can be driven unstable when the drift velocities of both the ions are greater than the phase velocity of the wave. We thus introduce another source of instability for these waves namely multiple ion beams drifting perpendicular to the magnetic field.     

Acceleration of Electrons and Protons in Reconnecting Current Sheets Including Single or Multiple X-points

Kinematic characteristics of electrons and protons in the magnetic reconnecting current sheet in the presence of a guide field are investigated. Particle trajectories are calculated for different values of the guide field by a test-particle calculation. The relationship between the final energy and the initial position has also been studied. We found that the addition of a guide field not only allows particles to get more energy and not only results in the separation of electrons and protons, but also causes the reconnecting electric field to selectively accelerate electrons and protons for different initial positions. The energy spectrum eventually obtained is the common power-law spectrum, and as the guide field increases, the index for the spectrum of electrons decreases rapidly. However, for a weak background magnetic field, proton spectra are not very sensitive to the guide field; but for a strong background field, the dependence of the spectrum index is similar to the electron spectrum. Meanwhile, kinematic characteristics of the accelerated particles in the current sheet including multiple X-points and O-points were also investigated. The result indicates that the existence of the multiple X- and O-points helps particles trapped in the accelerating region to gain more energy, and yields the double or multiple power-law feature.     

2 1/2-Dimensional Reconnection Model and Energy Release in Solar Flares

We employ a 2 1/2-dimensional reconnection model to analyse different aspects of the energy release in two-ribbon flares. In particular, we investigate in which way the systematic change of inflow region variables, associated with the vertical elongation of current sheet, affects the flare evolution. It is assumed that as the transversal magnetic field decreases, the ambient plasma-to-magnetic pressure ratio increases, and the reconnection rate diminishes. As the transversal field decreases due to the arcade stretching, the energy release enhances and the temperature rises. Furthermore, the magnetosonic Mach number of the reconnection outflow increases, providing the formation of fast mode standing shocks above the flare loops and below the erupting flux rope. Eventually, in the limit of a very small transversal field the reconnection becomes turbulent due to a highly non-linear response of the system to small fluctuations of the transversal field. The turbulence results in the energy release fragmentation which increases the release efficiency, and is likely to be responsible for the impulsive phase of the flare. On the other hand, as the current sheet stretches to larger heights, the ambient plasma-to-magnetic pressure ratio increases which causes a gradual decrease of the reconnection rate, energy release rate, and temperature in the late phase of flare. The described magnetohydrodynamical changes affect also the electron distribution function in space and time. At large reconnection rates (impulsive phase of the flare) the ratio of the inflow-to-outflow magnetic field strength is much smaller than at lower reconnection rates (late phase of the flare), i.e., the corresponding loss-cone angle becomes narrower. Consequently, in the impulsive phase a larger fraction of energized electrons can escape from the current sheet downwards to the chromosphere and upwards into the corona – the dominant flare features are the foot-point hard X-ray sources and type III radio bursts. On the other hand, at low reconnection rates, more particles stay trapped in the outflow region, and the thermal conduction flux becomes strongly reduced. As a result, a superhot loop-top, and above-the-loop plasma appears, as sometimes observed, to be a dominant feature of the gradual phase.     

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.     

Motion of charged particles in the magnetosphere

The adiabatic motion of charged particles in the magnetosphere has been investigated using Mead-Fairfield magnetospheric field model (Mead and Fairfield, 1975). Since the motion of charged particles in a dipolar field geometry is well understood, we bring out in this paper some important features in characteristic motion due to non-dipolar distortions in the field geometry. We look at the tilt averaged picture of the field configuration and estimate theoretically the parameters like bounce period, longitudinal invariant and the bounce averaged drift velocities of the charged particle in the Mead-Fairfield field geometry. These parameters are evaluated as a function of pitch angle and azimuthal position in the region of ring current (5 to 7 Earth radii from the centre of the Earth) for four ranges of magnetic activity. At different longitudes the non-dipolar contribution as a percentage of dipole value in bounce period and longitudinal invariant show maximum variation for particles close to 90° pitch angles. For any low pitch angle, these effects maximize at the midnight meridian. The radial component of the bounce averaged drift velocity is found to be greatest at the dawn-dusk meridians and the contribution vanishes at the day and midnight meridians for all pitch angles. In the absence of tilt-dependent terms in the model, the latitudinal component of the drift velocity vanishes. On the other hand, the relative non-dipolar contribution to bounce averaged azimuthal drift velocity is very high as compared to similar contribution in other characteristic parameters of particle motion. It is also shown that non-dipolar contribution in bounce period, longitudinal invariant and bounce averaged drift velocities increases in magnitude with increase in distance and magnetic activity.     

Two-dimensional equilibrium of a neutral sheet in a gravitational field

We solve the two-dimensional MHD-equations to find the equilibrium structure of a neutral sheet having its axis parallel to the gravitational field in an atmosphere. The evolution of an initially plane sheet which is a self-consistent, non-equilibrium solution of the MHD-equations, is followed until static equilibrium of forces is obtained. The effect of field line anchoring in a cool dense layer at the bottom of the sheet is taken into account.Presently at The Auroral Observatory, Boks 953, 9001 Tromsø, Norway.     

First and second order approximations of the first adiabatic invariant for a charged particle interacting with a linearly polarized hydromagnetic plane wave

In the present paper the effect of a sinusoidal modulation of an electromagnetic field on the invariance of the magnetic moment is studied. Such a generalized invariant plays an important role in problems concerning the motion of charged particles in the non-uniform magnetic field of the magnetosphere or the solar wind. In order to find an adiabatic invariant J, a canonical transformation is introduced, and J is expanded in an asymptotic series in the relative modulation amplitude. We are studying the first and second order terms of this expansion. It is further shown that the curves J = constant closely fit the results obtained by a numerical integration of the system of differential equations governing the motion of the particles.     

磁重联电流片的参数变化对电子加速的影响

通过采用试验粒子的方法,研究了在有引导磁场Bz存在的磁重联电流片中,电子被super-Dreicer电场Ez加速后的运动特征.首先,考虑了引导磁场恒定且与电场有不同方向时对粒子加速的影响.在这种情况下,Bz方向的改变直接改变了电子的运动轨迹,使其沿着不同的路径离开电流片.在Bz和Ez同向时,高能电子的pitch-angle接近于180°.然而,当2者反向时,高能电子的pitch-angle接近0°.引导磁场的取向只是使电场有选择地对不同区域的电子进行加速,不会最终影响电子的能量分布,最终得到的能谱是普遍的幂率谱E-γ.在典型的日冕条件下, γ大约等于2.9.进一步的研究表明γ的大小依赖于引导磁场及磁重联电场的强弱,以及电流片的尺度.随后,也研究了包含多个X-点和O-点电流片中被加速粒子的运动特征.结果表明X-点和O-点的存在使得粒子被束缚在加速区并获得最大的加速,而且最终的能谱具有多幂率谱的特征.