Understanding complex astrophysical and laboratory plasmas begins with the dynamics of a solitary charged particle in prescribed electric and magnetic fields. In strong magnetic fields, rapid cyclotron gyration decouples from the slower secular motion of the orbit’s orbit center. This separation of scales forms the bedrock of Guiding Center Theory, which governs particle trapping in planetary magnetospheres, cosmic ray diffusion, and tokamak plasma confinement.
1. The Lorentz Force and Cyclotron Gyration
A particle of mass (m) and charge (q) in arbitrary electromagnetic fields (mathbf{E}(mathbf{r}, t)) and (mathbf{B}(mathbf{r}, t)) obeys the relativistic Lorentz force law:
[
rac{dmathbf{p}}{dt} = q(mathbf{E} + mathbf{v} imes mathbf{B})
]
In a uniform, static magnetic field (mathbf{B} = B hat{mathbf{z}}) with (mathbf{E} = 0), the parallel velocity (v_parallel) remains invariant, while the perpendicular velocity vector (mathbf{v}_perp) rotates uniformly at the Cyclotron (Gyro) Frequency (omega_c):
[
omega_c = rac{|q| B}{m}
]
The radius of gyration, or Larmor Radius (Gyroradius), is determined by the balance of Lorentz and centripetal forces:
[
r_L = rac{v_perp}{omega_c} = rac{m v_perp}{|q| B}
]
2. Drift Velocities in Inhomogeneous Fields
When external forces or field gradients perturb the circular gyromotion, the guiding center (mathbf{R}) drifts across magnetic field lines with distinct drift velocities:
- (mathbf{E} imes mathbf{B}) Drift: Under an electric field perpendicular to (mathbf{B}), balancing (q(mathbf{E} + mathbf{v}_E imes mathbf{B}) = 0) yields:
[
mathbf{v}_E = rac{mathbf{E} imes mathbf{B}}{B^2}
]
Crucially, (mathbf{v}_E) is independent of particle charge (q), mass (m), and energy. Electrons and ions drift collectively with identical velocities, generating zero net electric current. - Grad-(B) Drift: A gradient in magnetic field strength perpendicular to (mathbf{B}) alters the Larmor orbit curvature between regions of stronger and weaker field:
[
mathbf{v}_{
abla B} = rac{m v_perp^2}{2 q B^3} (mathbf{B} imes
abla B) = rac{mu}{q B^2} (mathbf{B} imes
abla B)
]
where (mu quiv rac{m v_perp^2}{2 B}) is the magnetic moment. This drift is charge-dependent, inducing current and charge separation. - Curvature Drift: A charged particle streaming along a curved magnetic field line with radius of curvature (mathbf{R}_c) experiences a centrifugal force (mathbf{F}_{ ext{cent}} = rac{m v_parallel^2}{R_c^2} mathbf{R}_c):
[
mathbf{v}_R = rac{m v_parallel^2}{q B^2 R_c^2} (mathbf{R}_c imes mathbf{B})
]
In vacuum magnetic fields ((
abla imes mathbf{B} = 0)), curvature and grad-(B) drifts sum naturally to (mathbf{v}_{ ext{curv}+
abla B} = rac{m}{q B^4}left(v_parallel^2 + rac{1}{2}v_perp^2
ight)(mathbf{B} imes
abla B^2 / 2)). - Polarization Drift: A time-dependent electric field (partial_t mathbf{E}) induces an inertial lag:
[
mathbf{v}_p = rac{m}{q B^2} rac{partial mathbf{E}_perp}{partial t}
]
producing a polarization current density (mathbf{J}_p = sum n_s q_s mathbf{v}_{p,s} = rac{
ho_m}{B^2} rac{partial mathbf{E}_perp}{partial t}) dominated entirely by heavy ions.
3. The Three Canonical Adiabatic Invariants
When the background magnetic field varies slowly compared to the characteristic orbital period of a periodic motion, the classical action integral (mathcal{J} = oint mathbf{p} cdot dmathbf{q}) constitutes an Adiabatic Invariant:
- First Invariant ((mu)): Associated with cyclotron gyration on timescale ( au_c sim omega_c^{-1}):
[
mu = rac{p_perp^2}{2 m B} = rac{m v_perp^2}{2 B} = ext{const}
]
As a particle moves into a region of converging magnetic field lines ((B uparrow)), (v_perp) increases at the expense of (v_parallel) to conserve total kinetic energy (E = rac{1}{2}m v_parallel^2 + rac{1}{2}m v_perp^2). When (v_parallel = 0), the particle is reflected, creating Magnetic Mirrors in planetary Van Allen radiation belts and mirror confinement machines. - Second Invariant ((J)): Associated with the longitudinal bounce motion between mirror points:
[
J = oint v_parallel ds = ext{const}
] - Third Invariant ((Phi)): Associated with the slow azimuthal drift orbit around the geomagnetic dipole:
[
Phi_M = oint mathbf{A} cdot dmathbf{l} = iint mathbf{B} cdot dmathbf{S} = ext{const}
]
conserving the total magnetic flux enclosed by the guiding center drift shell.
