Adiabatic formation and properties of a quasi-equilibrium beam distribution matched to a periodic focusing lattice *
M. Dorf,R. Davidson,E. Startsev,H. Qin
2010-01-01
Abstract:This paper reports on recent advances in the development of a numerical scheme for describing the quiescent loading of a quasi-equilibrium beam distribution matched to a periodic focusing lattice [1]. The scheme allows for matched-beam distribution formation by means of the adiabatic turn-on of the oscillating focusing field, and it is examined here for the cases of alternating-gradient quadrupole and periodic solenoidal lattices. Furthermore, various distributions are considered for the initial beam equilibrium. The self-similar evolution of the matched-beam density profile is observed for arbitrary choice of initial distribution function and lattice type. The numerical simulations are performed using the WARP particle-in-cell code. INTRODUCTION Periodic-focusing accelerators and transport systems have a wide range of applications ranging from basic scientific research to industrial applications [2]. It is therefore important to develop an improved theoretical understanding of intense beam transport, stability, and equilibrium properties. However, the effects of the intense self-fields produced by the beam space-charge introduce a significant challenge for analytical studies, and the only known distribution function that provides an exact beam equilibrium in an periodic-focusing lattice is the Kapchinskij-Vladimirskij (KV) distribution [2]. This distribution, however, is unstable at sufficiently high beam intensities, and therefore approximate analytical analysis or numerical models have to be developed to describe other (stable) intense beam quasi-equilibrium solutions. Recently, a numerical approach for describing the quiescent formation of a quasi-equilibrium beam distribution matched to an alternating-gradient quadrupole lattice by means of adiabatic turn-on of the oscillating focusing field has been proposed in [1]. In this approach, an equilibrium beam distribution is first loaded into the equivalent smooth-focusing uniform channel [2]; then the oscillating feature of the lattice is adiabatically turned-on, and is accompanied by a corresponding decrease in the uniform focusing component. It was demonstrated that the method allows for the quiescent formation of a quasiequilibrium beam distribution matched to a quadrupole lattice for a broad range of beam intensities and vacuum phase advances describing the strength of the oscillating focusing field. Furthermore, properties of the quasiequilibrium matched-beam distribution have been investigated and compared to the predictions of analytical theory that applies Hamiltonian averaging techniques to the nonlinear Vlasov-Maxwell equations [2], assuming sufficiently small vacuum phase advance, σv. In particular, correct to order ε, where ε≡σv/2π, the theory predicts the self-similar evolution of the beam density profile. The evolution of the quasi-equilibrium beam density profile obtained in the numerical simulations making use of the adiabatic formation scheme has been investigated. The self-similarity feature was preserved to good accuracy for σv<66, and for a wide range of beam intensity ranging from an emittance-dominated beam with σ/σv≈0.26, to a space-charge-dominated beam with σ/σv≈0.91. Here, σ denotes the depressed phase advance. Furthermore, it was shown that for higher values of vacuum phase advance (for instance, σv=87.5) the selfsimilarity feature becomes less accurate, which demonstrates the validity limits of the theory. In this paper, we continue developing this method for adiabatic formation of a matched-beam quasi-equilibrium, and generalize it to the case of a periodic focusing solenoidal lattice. Furthermore, other initial smoothfocusing equilibrium distributions, namely a thermal equilibrium and a waterbag equilibrium, are considered. MODEL DESCRIPTION In this section, following [1], we summarize the method for adiabatic formation of a quasi-equilibrium beam distribution matched to an alternating-gradient quadrupole lattice, and generalize it to the case of a periodic solenoidal focusing lattice. The method works as follows. First, the oscillating focusing field of the lattice is replaced with the smooth-focusing force, and an equilibrium beam distribution is loaded into a uniform focusing channel. Then the oscillating focusing field is adiabatically turned on as the amplitude of the uniform component is adjusted to maintain the average (smoothfocusing) effects of the total focusing field fixed. For the case of a quadrupole lattice, the total focusing force acting on the beam particles is specified by ( ) ( ) ( ) ( )( ) 2 1 q q sf x y q x y F V s x y V s s x y κ κ = − + − − ⎡ ⎤ ⎣ ⎦ e e e e , (1) where V(s) is a function describing the transition of the applied lattice force and varying from zero to unity, κq(s) is the quadrupole lattice coupling coefficient, and q sf κ is the corresponding smooth-focusing value. For the choice of the quadrupole lattice model shown in Fig. 1 it is straightforward to demonstrate to leading order that ___________________________________________ *Research supported by the U.S. Department of Energy. TH6PFP090 Proceedings of PAC09, Vancouver, BC, Canada 3916 Beam Dynamics and Electromagnetic Fields D02 Non-Linear Dynamics Resonances, Tracking, Higher Order ( ) ( ) 3 2 1 ˆ 16 1 2 2 2 q q q q sf S η κ η κ − = . To model the adiabatic transition of the lattice, we adopt a simple model in which V(s) varies according to [1] ( ) ( ) [ ] [ ] ( ) [ ] 1 2 1 1 2 1 exp 1 exp 1 − − + − − + = tr tr L L L s L s V , (2) where 2L1/2 is the length of the transition section, and Ltr is the characteristic length scale for variation of V(s). The present approach for adiabatic formation of a beam quasi-equilibrium matched to a quadrupole lattice, can be generalized to the case of a periodic solenodial lattice. For the case of a solenoidal lattice, to maintain the average (smooth-focusing) effects of the total focusing field fixed, the transition of the applied lattice force is specified by ( ) ( ) ( ) ( ) { }( ) 2 1 1 sol sol sol sol sf sf sol x y F V s V s s x y κ κ κ δκ = − ⎡ − − ⎤ + + ⎣ ⎦ e e , (3) where κsol(s) is the solenoidal lattice coupling coefficient,