[en] We study the dynamics of the gyrophase of a charged particle in a magnetic field which is uniform in space but changes slowly with time. As the magnetic field evolves slowly with time, the changing of the gyrophase is composed of two parts. The rst part is the dynamical phase, which is the time integral of the instantaneous gyrofrequency. The second part, called geometric gyrophase, is more interesting, and it is an example of the geometric phase which has found many important applications in different branches of physics. If the magnetic field returns to the initial value after a loop in the parameter space, then the geometric gyrophase equals the solid angle spanned by the loop in the parameter space. This classical geometric gyrophase is compared with the geometric phase (the Berry phase) of the spin wave function of an electron placed in the same adiabatically changing magnetic field. Even though gyromotion is not the classical counterpart of the quantum spin, the similarities between the geometric phases of the two cases nevertheless reveal the similar geometric nature of the different physics laws governing these two physics phenomena.

[en] The Courant-Snyder (CS) theory and the Kapchinskij-Vladimirskij (KV) distribution for high-intensity beams in a uncoupled focusing lattice are generalized to the case of coupled transverse dynamics. The envelope function is generalized to an envelope matrix, and the envelope equation becomes a matrix envelope equation with matrix operations that are non-commutative. In an uncoupled lattice, the KV distribution function, first analyzed in 1959, is the only known exact solution of the nonlinear Vlasov-Maxwell equations for high-intensity beams including self-fields in a self-consistent manner. The KV solution is generalized to high-intensity beams in a coupled transverse lattice using the generalized CS invariant. This solution projects to a rotating, pulsating elliptical beam in transverse configuration space. The fully self-consistent solution reduces the nonlinear Vlasov-Maxwell equations to a nonlinear matrix ordinary differential equation for the envelope matrix, which determines the geometry of the pulsating and rotating beam ellipse. These results provide us with a new theoretical tool to investigate the dynamics of high-intensity beams in a coupled transverse lattice. A strongly coupled lattice, a so-called N-rolling lattice, is studied as an example. It is found that strong coupling does not deteriorate the beam quality. Instead, the coupling induces beam rotation, and reduces beam pulsation.

[en] In a linear trap confining a one-component nonneutral plasma, the external focusing force is a linear function of the configuration coordinates and/or the velocity coordinates. Linear traps include the classical Paul trap and the Penning trap, as well as the newly proposed rotating-radio- frequency traps and the Mobius accelerator. This paper describes a class of self-similar nonlinear solutions of nonneutral plasma in general time-dependent linear focusing devices, with self-consistent electrostatic field. This class of nonlinear solutions includes many known solutions as special cases.

[en] In an uncoupled lattice, the Kapchinskij-Vladimirskij (KV) distribution function first analyzed in 1959 is the only known exact solution of the nonlinear Vlasov-Maxwell equations for high- intensity beams including self-fields in a self-consistent manner. The KV solution is generalized here to high-intensity beams in a coupled transverse lattice using the recently developed generalized Courant-Snyder invariant for coupled transverse dynamics. This solution projects to a rotating, pulsating elliptical beam in transverse configuration space, determined by the generalized matrix envelope equation.

[en] The phase-space dynamics of runaway electrons is studied, including the influence of loop voltage, radiation damping, and collisions. A theoretical model and a numerical algorithm for the runaway dynamics in phase space are developed. Instead of standard integrators, such as the Runge-Kutta method, a variational symplectic integrator is applied to simulate the long-term dynamics of a runaway electron. The variational symplectic integrator is able to globally bound the numerical error for arbitrary number of time-steps, and thus accurately track the runaway trajectory in phase space. Simulation results show that the circulating orbits of runaway electrons drift outward toward the wall, which is consistent with experimental observations. The physics of the outward drift is analyzed. It is found that the outward drift is caused by the imbalance between the increase of mechanical angular momentum and the input of toroidal angular momentum due to the parallel acceleration. An analytical expression of the outward drift velocity is derived. The knowledge of trajectory of runaway electrons in configuration space sheds light on how the electrons hit the first wall, and thus provides clues for possible remedies.

[en] In tokamaks, Ware pinch is a well known neoclassical effect for trapped particles in response to a toroidal electric field. It is generally believed that there exists no similar neoclassical effect for circulating particles without collisions. However, this belief is erroneous, and misses an important effect. We show both analytically and numerically that under the influence of a toroidal electric field parallel to the current, the circulating orbits drift outward toward the outer wall with a characteristic velocity O ((var_epsilon)^-1) larger than the E x B velocity, where (var_epsilon) is the inverse aspect-ratio of a tokamak. During a RF overdrive, the toroidal electric field is anti-parallel to the current. As a consequence, all charged particles, including backward runaway electrons, will drift inward towards the inner wall.

[en] In this paper, a 3-D nonlinear perturbative particle simulation code (BEST) [H. Qin, R.C. Davidson and W.W. Lee, Physical Review Special Topics on Accelerators and Beams 3 (2000) 084401] is used to systematically study the stability properties of intense nonneutral charged particle beams with large temperature anisotropy (T_{perpendicularb} >> T_parallelb). The most unstable modes are identified, and their eigenfrequencies, radial mode structure, and nonlinear dynamics are determined for axisymmetric perturbations with ∂/∂θ = 0

[en] The centroid and envelope dynamics of a high-intensity charged particle beam are investigated as a beam smoothing technique to achieve uniform illumination over a suitably chosen region of the target for applications to ion-beam-driven high energy density physics and heavy ion fusion. The motion of the beam centroid projected onto the target follows a smooth pattern to achieve the desired illumination, for improved stability properties during the beam-target interaction. The centroid dynamics is controlled by an oscillating 'wobbler', a set of electrically-biased plates driven by RF voltage.

[en] A multi-species perturbative nonlinear (δf) electromagnetic particle simulation scheme has been developed for studying the propagation of intense charged particle beams in high-intensity accelerators and transport systems. The scheme is based on the Darwin approximation of Ampere's law, in which the transverse inductive electric field is neglected, resulting in the elimination of high-frequency transverse electromagnetic effects and, consequently, the associated numerical restrictions from the simulation. However, as noted previously, the presence of the time derivative of the vector potential in the equations of motion for the Darwin model can cause numerical instability. To circumvent this difficulty, we have adopted an approach by replacing the mechanical momentum, p_z, in the direction of beam propagation, by the canonical momentum, P_z=p_z+qA_z/c, as the phase-space variable. The resulting Ampere's law is then modified by the presence of an additional shielding term associated with the skin depth of the species. In order to minimize the numerical noise and to easily access both linear and nonlinear regimes for the physics of interest, we have also adopted the δf formalism for the Darwin model. The absence of unwanted high-frequency waves also enables us to use the adiabatic particle pushing scheme to compensate for the mass-ratio disparities for the various species of charge. The scheme is ideal for studying two-stream and filamentation instabilities, which may cause deterioration of the beam quality in the heavy ion fusion driver and fusion chamber

[en] The longitudinal and transverse dynamics of a heavy ion fusion beam during the drift compression and final focus phase is studied. A lattice design with four time-dependent magnets is described that focuses the entire beam pulse onto a single focal point with the same spot size