Search

Pitch angle scattering is a factor that helps determine the dawn-to-dusk current, controls particle energization, and it has also been used as a remote probe of the current sheet structure. Previous studies have interpreted their results under the exception that randomization will be greatest when the ratio of the two timescales of motion (gyration parallel to and perpendicular to the current sheet) is closet to one. Recently, the average expotential divergence rate (AEDR) has been calculated for particle motion in a hyperbolic current sheet (Chen, 1992). It is claimed that this AEDR measures the degree of chaos and therefore may be thought to measure the randomization. In contrast to previous expectations, the AEDR is not maximized when Kappa is approximately equal to 1 but instead increases with decreasing Kappa. Also contrary to previous expectations, the AEDR is dependent upon the parameter b(sub z). In response to the challenge to previous expectations that has been raised by this calculation of the AEDR, we have investigated the dependence of a measure of particle pitch angle scattering on both the parameters Kappa and b(sub z). We find that, as was previously expected, particle pitch angle scattering is maximized near Kappa = 1 provided that Kappa/b(sub z) greater than 1. In the opposite regime, Kappa/b(sub z) less than 1, we find that particle pitch angle scattering is still largest when the two timescales are equal, but the ratio of the timescales is proportional to b(sub z). In this second regime, particle pitch angle scattering is not due to randomization, but is instead due to a systematic pitch angle change. This result shows that particle pitch angle scattering need not be due to randomization and indicates how a measure of pitch angle scattering can exhibit a different behavior than a measure of chaos.

A one-dimensional, hybrid simulation code is used to study current sheets with a nonzero normal magnetic field B(sub z) and a dawn-to-dusk electric field E(sub y). Such configurations are dependent upon only two parameters: we use the normalized normal magnetic field B-normalized (sub z) = B(sub z)/(4(pi)(n(sub b)) (v(exp 2 sub T))(exp 1/2) and normalized electric field V-normalized (sub D) = (1/V(sub T)(cE(sub y)/B(sub z)), where V(sub T) is the thermal velocity of ions prior to their interaction with the current sheet and n(sub b) is the number density outside the current sheet (at the simulation boundary). A third parameter that is relevant to the motion of particles in current sheets is kappa(sub A), the value of kappa = (R(sub min)/rho(sub max))(exp 1/2) for particles of average energy. We find that if either B-normalized (sub z) is close to or greater than 1, or if kappa(sub A) is close to 1, a rotational mode develops in which the z = 0 current rotates with the ion sense about the normal magnetic field, while for small values of both B-normalized (sub z) or kappa(sub A), the configuration is quasi-steady. To achieve values of kappa(sub A) of the order of or larger than 1, we decrease the value of V-normalized (sub D) uniformly. We find that the magnetic field fluctuations and particle distribution functions are similar in many respects to what was observed in the day 240, 1986, Active Magnetospheric Particle Tracer Explorer (AMPTE)/CCE current disruption event, an event that appears to be located at the site of initiation of current disruption and related particle energization.

The ion tearing instability is investigated in a magnetotail configuration that consists of a diffuse plasma sheet current and an embedded, thin current sheet with a strong current. For historical reasons, the thin embedded current sheet will be called a 'neutral sheet', even though the normal component of the magnetic field, Bn, is nonzero. In particular, we assume that the current within the thin current sheet is due to the acceleration of 'Speiserlike' ion trajectories by a cross-tail electric field Ey. It is found that the strong current within the neutral sheet is essentially unimportant to the growth rate of the tearing instability, and that the growth rate scales as (lambda(0)/Lz) squared, where Lz is the overall half thickness of the plasma sheet and lambda (0) is the ion inertial length. In the absence of the current outside the neutral sheet, current filamentation is stable.

A particle model of the magnetotail neutral sheet self-consistent structure was developed by assuming that the thickness of the neutral sheet is much less than the thickness of the plasma sheet. This made it possible to approximate the neutral sheet as an x-independent structure. The role of nonlinear particle motions in determining the conductivity and the global current sheet structure is investigated.

The principles of the current sheet catastrophe models are briefly reviewed, and observations of some of the signatures predicted by the theory are presented. The data considered here include AMPTE/CCE observations of fifteen current sheet disruption events. According to the model proposed here, the root cause of the current disruption is some process, as yet unknown, that leads to an increase in the k sub A parameter. Possible causes for the increase in k sub A are discussed.

Analytic calculations and simulations of the steady-state structure of the dissipation region during collisionless reconnection are described. Ions are treated as particles and electrons as a dissipative fluid. A limiting reconnection rate results as the ions 'blowout' the magnetic field in the dissipation area. The implications for comprehending substorm dynamics are discussed.

The modified Harris field model is used to investigate the process of pitch angle scattering by a current sheet. The relationship between the incoming asymptotic pitch angle alpha(in) and the outgoing asymptotic pitch angle alpha(out) is studied from first principles by numerically integrating the equation of motion. Evidence that charged particles undergo chaotic scattering by the current sheet is found. For fixed alpha(in), it is shown that alpha(out) exhibits sensitive dependence on the energy parameter in certain energy ranges. For a fixed energy parameter value in the same energy ranges, alpha(out) sensitively depends on alpha(in). For other energy values, alpha(out) does not show sensitive dependence on alpha(in) for most phase angles. A distribution of alpha(in) is mapped from the asymptotic region to the midplane, and it is found that the resulting particle distribution should have beam structures with well-collimated pitch angles near each resonance energy value. Implications for the particle distribution functions in the earth's magnetotail are discussed.

The process of 'differential memory' in the earth's magnetotail is studied in the framework of the modified Harris magnetotail geometry. It is verified that differential memory can generate non-Maxwellian features in the modified Harris field model. The time scales and the potentially observable distribution functions associated with the process of differential memory are investigated, and it is shown that non-Maxwelllian distributions can evolve as a test particle response to distribution function boundary conditions in a Harris field magnetotail model. The non-Maxwellian features which arise from distribution function mapping have definite time scales associated with them, which are generally shorter than the earthward convection time scale but longer than the typical Alfven crossing time.

The results are presented of an analytic and numerical investigation of the structure of the X line during steady state magnetic reconnection in collisionless plasma. The structure of the X line essentially depends on a single dimensionless parameter F, which is a measure of the influx of plasma into the reconnection region. For small F the self-consistent plasma current driven at the X line is small, and the magnetic fields are nearly unchanged from the initial vacuum state. With increasing F the current driven at the X line becomes large and the dissipation region collapses in the direction of the inflow and elongates along the outflow. For sufficiently large F the velocity of the plasma ejected form the X line exceeds the local Alfven velocity. In this regime a fast mode shock forms at the outflow end of the dissipation region which slows the high-velocity outflow plasma to the subsonic flow characteristic of the broader outflow region. Finally, at a critical plasma flux the dissipation region collapses to zero thickness.

The structure of the dissipation region during magnetic reconnection in collisionless plasma is investigated by examining a prescribed two-dimensional magnetic x line configuration with an imposed inductive electric field E(y). The calculations represent an extension of recent MHD simulations of steady state reconnection (Biskamp, 1986; Lee and Fu, 1986) to the collisionless kinetic regime. It is shown that the structure of the x line reconnection configuration depends on only two parameters: a normalized inductive field and a parameter R which represents the opening angle of the magnetic x lines.

The integrodifferential equation describing the linear tearing instability in the bi-Maxwellian neutral sheet is solved without approximating the particle orbits or the eigenfunction psi. Results of this calculation are presented. Comparison between the exact solution and the three-region approximation motivates the piecewise-straight-line approximation, a simplification that allows faster solution of the integrodifferential equation, yet retains the important features of the exact solution.

The problem of linear bi-Maxwellian neutral-sheet tearing instability has been solved, using neither an approximation of the orbits nor the constant-psi approximation. The present method has a larger range of validity than the best previous approximation, the three-region approximation. It is found that as the temperature anisotropy is increased, the growth rate increases and the wavelength associated with the peak growth decreases. It is suggested that if small-island formation is experimentally observed, then temperature anisotropy may be present.