Staility of off-axis motion for intense particle eams in periodically focusing channels J.S. Moraes a,, R. Pakter a, and F.B. Rizzato a a Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 15051, 91501-970 Porto Alegre, Rio Grande do Sul, Brasil. and Centro Universitário La Salle, Av. Victor Barreto, 2288, 92010-000, Canoas, RS, Brasil arxiv:physics/0410156v1 [physics.plasm-ph] 19 Oct 2004 A general equation for the centroid motion of free, continuous, intense eams propagating off-axis in solenoidal periodic focusing fields is derived. The centroid equation is found to e independent of the specific eam distriution and may exhiit unstale solutions. A new Vlasov equilirium for off-axis eam propagation is also otained. Properties of the equilirium and the relevance of centroid motion to eam confinement are discussed. PACS numers: 41.85.Ja,05.45.-a A fundamental understanding of the kinetic equilirium and staility properties of high-current, lowemittance eams in periodically focusing systems is crucial for the development of a wide range of advanced particle accelerator and coherent radiation source applications. For a long time, the Kapchinskij-Vladimirskij (KV) distriution [1] was the only Vlasov equilirium distriution known for the propagation of periodically focused intense particle eams. Equilirium and staility analysis ased on the KV eam have een critical to the development and understanding of the physics of intense eams [2, 3, 4, 5, 6, 7, 8, 9, 10]. More recently, it has een shown that the KV distriution can e generalized to allow for rigid eam rotation with respect to the Larmor frame in periodic solenoidal focusing fields [11]. Studies indicate that rotation may have an important role in particle eam staility [12]. In the derivation of these Vlasov equiliria it is always assumed that the eam is perfectly aligned with the symmetry axis of the focusing field [1, 5, 11]. Actually, this simplifying assumption is generally used in the analysis of intense eams [7] ecause the axis is an equilirium for the eam centroid, and the equilirium is stale if smooth-eam approximations are employed where the periodic fluctuations of the focusing field are averaged out [13]. In some cases, however, we may expect the onset of parametric resonances involving the centroid motion and the focusing field oscillations, which would destailize the centroid motion and heavily affect the overall eam dynamics. In such conditions the averaging procedure is no longer valid and a detailed description of the centroid dynamics ecomes mandatory. In this paper, we derive from a kinetic Vlasov-Maxwell description a general equation for the centroid motion of free, continuous, intense eams propagating off-axis in solenoidal periodic focusing fields. It is shown that the centroid oeys a Mathieu type equation. The equation is independent of the specific eam distriution and ecomes unstale whenever the oscillatory frequency of the centroid, which is related to the rms focusing field strength per lattice, is commensurale with the focusing field periodicity itself. In the particular case of a uniform eam density around the eam centroid, we show that there exists a self-consistent Vlasov equilirium distriution for the eam dynamics. The eam envelope that determines de outer radius of the equilirium eam around the centroid is shown to oey the familiar envelope equation [5, 7, 14, 15], eing independent of the centroid motion. An example of the Vlasov equilirium is discussed in detail to show the possiility of finding eam solutions for which the extensively studied envelope equation [3, 4, 8, 9, 10, 16] is stale, whereas the centroid motion is unstale, revealing the importance of the centroid motion to overall eam confinement properties. We consider a free, continuous charged-particle eam propagating with average axial velocity β cê z through a periodic solenoidal focusing magnetic field descried y B(r,s) = B z (s)ê z r 2 B z(s)ê r, (1) where r = xê x + yê y, r = (x 2 + y 2 ) 1/2 is the radial distance from the field symmetry axis, s = z = β ct is the axial coordinate, B z (s+s) = B z (s) is the magnetic field on the axis, the prime denotes derivative with respect to s, c is the speed of light in vacuo, and S is the periodicity length of the magnetic focusing field. Since we are dealing with solenoidal focusing, it is convenient to work in the Larmor frame of reference [7], which rotates with respect to the laoratory frame with angular velocity Ω L (s) = qb z (s)/2γ mc, where q, m and γ = (1 β 2 ) 1/2 are, respectively, the charge, mass and relativistic factor of the eam particles. The Larmor frame is used throughout the paper, such that ê x and ê y are assumed to e versors along the rotating axes. In the paraxial approximation, the eam distriution function f (r,v,s) evolves according to the Vlasov-Maxwell system [7] e-mail: jsmoraes@if.ufrgs.r e-mail: pakter@if.ufrgs.r e-mail: rizzato@if.ufrgs.r f s +v f +( κ z r ψ) v f = 0, (2) 2 ψ = (2πK/N )n (r,s), (3)
n = f dv, (4) where n (r,s) is the eam density profile, κ z (s) = q 2 Bz 2(s)/4γ2 β2 m2 c 4 is the focusing field parameter, K = 2q 2 N /γ 3β2 mc2 is the eam perveance, N = f drdv =const. is the numer of particles per unit axial length, and v r [17]. In Eqs. (2)-(4), ψ is a normalized potential that incorporates oth selfelectric and self-magnetic fields, E s and B s. It is related to the self-scalar and self-vector potentials y φ s = β 1 A s z = γ3 mβ2 c2 ψ(r,s)/q, where A s (r,s) = A s z(r,s)ê z, E s (r,s) = φ s (r,s), and B s (r,s) = A s (r,s). Our first task here is to determine the evolution of the eam centroid located at r(s) N 1 rf (r,v,s)drdv. (5) In order to do that one multiplies Eq. (2) y r and integrates over phase-space to get r = v, (6) where v N 1 vfdrdv. If one now multiplies Eq. (2) y v and integrates over phase-space, one otains v = κ z r ψ, (7) where ψ N 1 ψfdrdv is otained y integration y parts of the v - term in velocity space. Using Eqs. (3) and (4) we can rewrite ψ as ψ = (2πK) 1 ψ 2 ψdr. (8) Then we note that the integrand of Eq. (8) can e cast into the more suitale form ψ 2 ψ = [ ψ ψ I( ψ) 2 /2] (9) where the unit dyadic I reads I ê x ê x + ê y ê y. Now, employing Gauss theorem we otain ψ = (2πK) 1 ê n [ ψ ψ I( ψ) 2 /2 ] da = 0, (10) ecause ψ 0 as r for eams in free space. In Eq. (10), da and ê n are, respectively, the oundary differential element and the unit vector normal to the oundary of integration located at r. In fact, the result ψ = 0 is expected ased on the action-reaction law; since ψ corresponds to the self-force exerted on the eam particles y themselves, its average throughout the eam distriution has to vanish due to the pairwise structure of the interparticle electromagnetic interaction. Using Eqs. (7) and (10) in Eq. (6), we finally otain the centroid equation of motion r +κ z (s) r = 0. (11) Let us stress that we have not made any assumption on the particular form of the eam distriution function so far. Thus, the centroid equation aove is always valid as long as the eam evolves according to the Vlasov-Maxwell system, Eqs. (2)-(4). In the laoratory frame, comined to the oscillatory motion descried y Eq. (11) the centroid also rotates with angular velocity Ω L (s) around the center r = 0. Taking into account that κ z (s) is periodic, Eq. (11) is of the Mathieu type which is known to present unstale solutions related to parametric resonances in the r motion. If we conveniently write the average of κ z (s) over one lattice periodas(1/s) S 0 κ z(s)ds σ0 2/S2,whereσ 0 isadimensionless parameter proportional to the rms focusing field measuring the vacuum phase advance in the small field, smooth-eam approximation, the instailities in the centroid motion are expected when one approachesσ 0 nπ; this condition corresponds to parametric resonances etween the oscillation periodicity of r in the average(rms) focusing field and the periodicity of the focusing field itself. Depending on the exact profile of κ z (s) the size of the unstale regions surrounding σ 0 nπ vary significantly. If the aim is eam confinement, these regions are to e avoided. It is worth mentioning that although Eq. (11) is strictly valid for free eams only, it is expected to provide a good description of the centroid motion in ounded systems if the eam is nearly symmetric and is not excessively displaced from a pipe center located at r = 0. The reason is ecause in this case ψ = ±ê n ψ at the pipe walls, where ê n is now the unit vector normal to the wall, and the surface integral in Eq. (10), performed along the oundary, still vanishes since ψ is approximately constant there. Note also that the presence of a pipe would generally not suppress the centroid instailities discussed in connection with Eq. (11); in fact, it would even enhance it ecause the image charges induced are of opposite sign, attracting the eam to the wall. OurnexttaskistoshowthatwecanconstructaVlasov equilirium for off-axis eam transport. In particular, we assume a eam with a uniform radial density distriuted aroundacenter located at r o (s) = x o (s)ê x +y o (s)ê y, i.e., { N /πr n (r,s) = 2(s), r δ < r (s), (12) 0, r δ > r (s), where r (s) is the equilirium eam envelope and r δ r r o. A schematic of the eam distriution of Eq. (12) and corresponding vectors is shown in Fig. 1. For such eam we can easily recognize r o (s) as eing the centroid coordinate. According to what was shown previously, its evolution must then oey r o +κ z (s)r o = 0. (13) Using the prescried n (r,s) in Eq. (3) we find for the normalized self-potential 2 ψ(r,s) = Krδ 2 /2r2 (s) (14)
3 in the eam interior (r δ < r ). Therefore, a single particle of the eam located at r(s) sujected to the external focusing field force κ z (s)r and the self-field force ψ(r, s) will evolve according to r +κ z r (K/r 2 )r δ = 0. (15) If we now sutract Eq. (13) from Eq. (15) we otain r δ +κ z r δ (K/r 2 )r δ = 0, (16) which descries the motion of the eam particle with respect to the center r o. Equation (16) can e solved with known techniques of physics of eams [5, 7]. Considering the motion along the x-axis, we write x δ = A xδ w(s)sin[ s ζ(s)ds+ζxδ 0]withA xδ andζ xδ 0 constants. Sustituting this expression into Eq. (16) we otain w +κ(s)w = w 3, (17) ζ(s) = w 2 (s), where κ(s) κ z (s) K/r 2 (s), and the constant of motion A xδ can e expressed in the form A 2 x δ = (x δ /w) 2 +(wx δ w x δ ) 2. (18) Performing an equivalent calculation for the motion along the y-axis, one shows that A yδ given y A 2 y δ = (y δ /w) 2 +(wy δ w y δ ) 2. (19) is also a constant of motion. From Eq. (16) one sees that all the forces are central with respect to the centroid r o. Thus, one readily demonstrates that the canonical angular momentum P Θδ given y P Θδ = x δ y δ y δx δ (20) is a constant of motion as well. Because A 2 x δ, A 2 y δ, and P Θδ are exact single-particle constants of motion, a possile choice of Vlasov equilirium distriution function is f EQ (r,v,s) = N [ π 2 δ A 2 x ǫ δ +A 2 y δ T 2ω P Θδ ( 1 ω 2 ) ] ǫt, (21) where df EQ /ds = 0, ǫ T =const is an effective emittance, and the rotation parameter ω =const is in the range 1 < ω < 1 for radially confined eams. Using f EQ in Eq. (4), it is readily shown that the uniform density profile centered at r o of Eq. (12) is consistently otained, provided r (s) = ǫ 1/2 T w(s). Hence, r (s) oeys the familiar envelope equation r +κ z(s)r K r ǫ2 T r 3 = 0. (22) Performing the appropriate averages over the equilirium distriution, Eq. (21), we can show that the eam rigidly rotates around its centroid r o (s) with angular velocity FIG. 1: Schematic of the eam distriution of Eq. (12) and corresponding vectors. C corresponds to the centroid position. Ω δ (s) = ω ǫ T β c/r 2 (s). Also, the rms emittance calculated with respect to the centroid position is given y ǫ χ 4[ χ 2 χ 2 χ χ 2] 1/2 = ǫt, (23) where χ = x δ,y δ, and the rackets indicate averagesover the eam distriution. One thus sees that a Vlasov equilirium distriution can e formed for which the eam envelope oeys Eq. (22) with constant emittance even when the centroid moves off-axis, r o 0, following the dynamics dictated y Eq. (13). We refer to this equilirium as a periodically focused off-axis Vlasov equilirium. Let us call attention to the interesting fact that the centroid motion and the envelope dynamics are uncoupled in this case. In other words, centroid dynamics does not affect the known staility results for the envelope dynamics [3, 4, 8, 10, 16] and is not affected y the latter as well. Oneshouldkeepin mind that forgoodeamconfinement oth centroid and envelope have to e stale. We now illustrate our results with an example of periodically focused off-axis Vlasov equilirium. We consider a particular set of parameters for which the envelope equation (22) is known to e stale, whereas the centroid motion of Eq. (13) was found to e unstale. We investigate eam transport with the aid of self-consistent numerical simulations, where a large numer N = 8000 of macroparticles interact via pairwise electromagnetic interactions [9]. In the simulation we used SK/ǫ T = 5.0 and S 2 κ z (s) = σ 2 0 [1+cos(2π)], with σ 0 = 155 o, over 20 lattice periods. The macroparticles were launched at s = 0 according to the equilirium distriution, Eq. (21), with ω = 0, r o = 0 = r o, and r corresponding to the matched solution with r (s + S) = r (s) of the envelope equation (22). The finite numer of macroparticles in the initial condition acts as a seed for any possile instaility to develop. Simulation results are presented in Fig. 2. The evolution of the centroid displacement r o r o calculated from the macroparticles positions r as r o =< r >, where the rackets indicate average over macroparticles, is shown in Fig. 2(a) (circles). It reveals that the centroid motion develops the typical exponential growth of unstale dynamics that agrees with the fact that the set of parameters considered leads to
4 EMITTANCE r 6.0 3.0 0.0 0.0 5.0 10.0 15.0 20.0 CENTROID DISPLACEMENT9.0 2.0 1.5 1.0 0.5 centroid equation simulation ε x /ε Τ ε y /ε Τ (a) () 0.0 0.0 5.0 10.0 15.0 20.0 1.4 1.2 1.0 0.8 0.6 0.4 envelope equation simulation (c) 17.0 18.0 19.0 20.0 FIG. 2: Multiparticle self-consistent simulation results. (a) The centroid motion; () rms emittance; and (c) the envelope dynamics. Centroid displacement and envelope are normalized to (Sǫ T) 1/2. an unstale solution of Eq. (13). The solid line corresponds to the solution otained from Eq. (13). Despite the centroid instaility, the eam equilirium distriution is preserved as verified in Fig. 2() that shows that rms emittance is well conserved as the eam evolves. RMS emittance is calculated according to Eq. (23), considering averages over macroparticles. Finally, Fig. 2(c) compares the envelope otained from the envelope equation (22) with that otained from the simulation, r = [2 < (r r o ) 2 >] 1/2, for the last 3 periods of the focusing channel. The perfect agreement proves once more the preservation of the equilirium distriution. Moreover, we see that in spite of the unstale centroid the envelope is stale, as predicted. To conclude, ased on kinetic grounds we have derived a general equation for the centroid motion of free, continuous, intense eams propagating off-axis in solenoidal periodic focusing fields. It was shown that the centroid equation is independent of the specific eam distriution and may exhiit unstale solutions. In the particular case of a uniform eam density around the eam centroid, we have shown the existence of a periodically focused off-axis Vlasov equilirium distriution descriing a eam that rigidly rotates with a prescried angular velocity around a moving centroid. The eam envelope around the centroid was shown to oey the familiar envelope equation, eing independent of the centroid motion. An example of periodically focused off-axis Vlasov equilirium was discussed in detail to show the possiility of finding eam solutions for which the envelope equation is stale, whereas the centroid motion is unstale, revealing the importance of centroid motion to the overall eam confinement properties. We acknowledge partial support from CNPq, Brazil. [1] I. M. Kapchinskij and V. V. Vladimirskij, in Proceedings of the International Conference on High Energy Accelerators (CERN, Geneva, 1959), p. 274. [2] I. Hofmann, L. J. Laslett, L. Smith, and I. Haer, Part. Accel., 13, 145 (1983). [3] C. J. Struckmeier and M. Reiser, Part. Accel., 14, 227, (1984). [4] C. Chen and R. C. Davidson, Phys. Rev. Lett., 72, 2195 (1994). [5] C. Chen and R. C. Davidson, Phys. Rev. E., 49, 5679 (1994). [6] R.L. Gluckstern, W.-H. Cheng, and H. Ye, Phys. Rev. Lett. 75 2835 (1995). [7] R.C. Davidson and H. Qin, Physics of Intense Charged Particle Beams in High Energy Accelerators (World Scientific, Singapore, 2001). [8] R. Pakter and F.B. Rizzato, Phys. Rev. Lett. 87, 044801 (2001).
5 [9] R. Pakter and F.B. Rizzato, Phys. Rev. E, 65, 056503 (2002). [10] S.M. Lund and B. Bukh, Phys. Rev. ST Accel. Beams 7, 024801 (2004). [11] C. Chen, R. Pakter, and R.C. Davidson, Phys. Rev. Lett. 79 225 (1997). [12] C. Chen, R. Pakter, and R.C. Davidson, Phys. Plasmas. 6 3647 (1999). [13] M. Hess and C. Chen, Phys. Plasmas 7, 5206 (2000). [14] P. M. Lapostolle, IEEE Trans. Nucl. Sci. NS-18, 1101 (1971) [15] F. J. Sacherer, IEEE Trans. Nucl. Sci., NS-18, 1105 (1971). [16] J.S. Moraes, F.B. Rizzato, and R. Pakter, Phys. Plasmas 10, 4811 (2003). [17] A note on the notation used. The vector derivatives with respect to s in the paper represent derivations directly in the Larmor frame, where ê x and ê y are considered fixed versors.