Particle drift and loss rates under strong pitch angle diffusion in Dungey's model magnetosphere

Strong pitch angle diffusion so thoroughly randomizes the first two adiabatic invariants of charged-particle motion that it conserves (as a first approximation) the surrounded phase space volume Lambda = p(3) Psi(0), where Psi(0) = closed integral(ds/B) is the flux tube volume per unit magnetic flux and p is the common scalar momentum of the particles involved, (This conservation law is widely used in plasma physics. Its nonrelativistic form is usually derived from the adiabatic gas law but is not really a "fluid" result, Since the conservation law does not couple particles with different energies in the same flux tube, such particles can have different values of Lambda.) Here the flux tube volume Psi is computed for selected field lines in Dungey's model magnetosphere (dipole moment mu(E) approximate to -30.2 mu T-R-E(3) plus uniform southward Delta B-z approximate to -9.0 cos(6) Lambda* mu T, Lambda* being the invariant latitude of the quiet time auroral oval, which maps outward to a circular neutral line at radial distance b = (mu(E)/Delta B-z)(1/3) in the equatorial plane). The integral closed integral(ds/B) is fitted within 0.2% at all L values by the expression Psi(0) approximate to (L(4)a(4)/mu(E)){(32/35) - [2.045 + 1.045(r(0)/b)(3) + 0.095(r(0)/b)(6) + 0.075(r(0)/b)(9)] ln[1 - (r(0)/b)(3)]}, where r(0) (equatorial radius of the field line) is obtained by solving the cubic equation (r(0)/b) = (La/b)[1 + (1/2)(r(0)/b)(3)]. The particle lifetime tau against strong pitch angle diffusion is given in turn by tau approximate to [4 Psi BnBs/(B-n + B-s)(1 - eta)](m/p), where B-n and B-s (which can differ in an offset-dipole field model) are the field intensities at the northern and southern foot points of the field line, m is the relativistic mass, and eta is the backscatter coefficient. The same expression for Psi thus enters both the quasi-adiabatic Hamiltonian (to whose derivative with respect to 1/L the ensemble-averaged gradient-curvature drift rate is directly proportional) and the particle loss rate 1/tau (by which the phase space density f = J(4 pi)/4 pi p(2) is exponentially attenuated along the ensemble averaged drift trajectory). This construction provides a means of numerically simulating kinematical aspects of electron motion in the outer magnetosphere and precipitation in the diffuse aurora, as well as ion motion in part of the plasma sheet (regarded as a source of the ring current). Its use in magnetospheric plasma simulations should save computer time and reveal analytical structure relevant to magnetospheric dynamics.