A 2 2 5 - dimensional particle-mesh computer model for the simulation of the current-sheet region of the geomagnetic tail is described. Important features are (a) the use of Fast Fourier techniques for the efficient solution of Ampere's equation, (b) the incorporation of sources and sinks of particles, (c) facilities for simulating finite width effects and (d) the option of including a normal magnetic field component linking through the sheet. Simulations carried out using this model indicate that current sheets with a non-zero normal magnetic field component and an infinite width are stable. The particles trace out Speiser-like orbits in such a case. Sheets with B normal = 0 and a finite width are unstable with respect to the ion tearing-mode instability. However the presence of a normal magnetic field stabilises the system provided ρ 0<2 L y where ρ 0 is the characteristic length associated with the normal field and where L y is the width of the sheet. On the basis of these results it is suggested that a geomagnetic substorm occurs when the normal magnetic field drops below the critical value needed for stability.
Read full abstract