We give a systematic theoretical analysis of trapped nonadiabatic charged particle motion in two-dimensional taillike magnetic field reversals. Particle dynamics is shown to be controlled by the curvature parameter κ, i.e., the ratio κ² = Rmin/ρmax between the minimum radius of curvature of the magnetic field and the maximum Larmor radius in it for a particle of given energy. κ≫1 corresponds to the usual adiabatic case with the magnetic moment μ as a first-order invariant of motion. As κ decreases toward unity, the particle motion becomes stochastic due to deterministic chaos, caused by the overlapping of nonlinear resonances between the bounce- and the gyro-motion. We determine the threshold of deterministic chaos and derive the related pitch angle diffusion coefficient which describes statistically the particle behavior in the limit κ → 1. Such behavior, which for κ ≅ 1 becomes strongly chaotic, applies, e.g., to thermal electrons in Earth's magnetotail and makes its collisionless tearing mode instability possible. We also show that in sharply curved field reversals, i.e., for κ 1. Both types of trapped particle motion in sharply curved magnetic field reversals κ<1 are closely connected with fast oscillations perpendicular to the reversal plane. However, the trajectories are adiabatic only in the case that they permanently remain crossing the reversal plane. The adiabatic are of a ring type, i.e., they resemble rings in phase space and also in real physical space. For ring-type orbits the action integral over the fast oscillations is an adiabatic invariant in the usual sense. On the other hand, the most common particle trajectories in a sharply curved field reversal with κ<1 are essentially of a cucumberlike quasi-adiabatic type. For quasi-adiabatic cucumberlike orbits the action integral over the fast oscillations is an adiabatic invariant only in a piecemeal way between successive traversals in the phase space of the fast motion of a separatrix between orbits which do and those, which do not cross the reversal plane. Due to the effect of separatrix traversals the slow motion shifts between different cucumber orbits with a conservation of the action integral on average but with its chaotic phase space diffusion even for very small perturbation parameters κ. The case κ<1 is applicable, e.g., to thermal ions and high-energy electrons in Earth's magnetotail. Our findings lead to a systematic interpretation of particle observations in Earth's magnetotail and of numerous numerical calculations, carried out in the past. They also explain rather well, e.g., the pitch angle diffusion of plasma sheet particles, the isotropization of the plasma sheet electron distribution immediately before a substorm and provide with the transition to chaos a mechanism for the onset of a large-scale tail instability and the explosion of isolated substorms. Further implications for magnetotail physics, such as acceleration processes and the influence of the particle escape from the field reversal will be discussed in a second related paper.

Regular and chaotic charged particle motion in magnetotaillike field reversals: 1. Basic theory of trapped motion

Given the solution (Xt ) of a Stochastic Differential System, two situat,ions are considered: computat,ion of Ef(Xt ) by a Monte–Carlo method and, in the ergodic case, integration of a function f w.r.t. the invariant probability law of (Xt ) by simulating a simple t,rajectory. For each case it is proved the expansion of the global approximat,ion error—for a class of discret,isat,ion schemes and of funct,ions f—in powers of the discretisation step size, extending in the fist case a result of Gragg for deterministic O.D.E. Some nn~nerical examples are shown to illust,rate the applicat,ion of extrapolation methods, justified by the foregoing expansion, in order to improve the approximation accuracy

/pdf/expansion-of-the-global-error-for-numerical-schemes-solving-35p3lm5ogz.pdf

Expansion of the global error for numerical schemes solving stochastic differential equations

The paper introduces a new semi-implicit extrapolation method especially designed for the numerical solution of stiff systems of ordinary differential equations. The existence of a quadratic asymptotic expansion in terms of the stepsize is shown. Moreover, the new discretization is analyzed in the light of well-known stability models. The efficiency of the new integrator is clearly demonstrated by solving a series of challenging test problems including real life examples.

A semi-implicit mid-point rule for stiff systems of ordinary differential equations

Chemical reaction systems are certainly one of the most challenging scientific fields in which numerical and analytical methods for ordinary differential equations are used.

Numerical Treatment of Inverse Problems in Chemical Reaction Kinetics

A theory of unimolecular reactions induced by coherent, monochromatic infrared radiation (URIMIR) in the absence of collisions is presented. It is shown that the set of first order linear differential equations for the amplitudes of molecular states (Schrodinger equation) can be reduced, under specified conditions, to a much smaller set of first order linear differential equations for the coarse grained populations of levels for polyatomic molecules (master equation). Four limiting situations are identified in which such linear rate equations provide a reasonable approximation. Rate coefficients are obtained as a function of spectroscopic parameters (energy levels and transition moments). Solutions for the master equations are given as a function of time and at steady‐state. Simple limiting cases (Strong Field Limit, Weak Field Limit, Reaction Threshold Bottleneck, etc.) are identified and very simple rate expressions are obtained for these cases. A complete statistical mechanical theory of URIMIR is form...

Theory of unimolecular reactions induced by monochromatic infrared radiation

Stoer, J. / Bulirsch, R., Einführung in die Numerische Mathematik II, IX, 286 S., 1973. DM 14,80, US $ 5.50. Berlin-Heidelberg-New York. Springer-Verlag

H. Schwetlick

Papers

Stoer, J. / Bulirsch, R., Einführung in die Numerische Mathematik II, IX, 286 S., 1973. DM 14,80, US $ 5.50. Berlin-Heidelberg-New York. Springer-Verlag