lattice design and study tools regarding the super-b projectlattice design and study tools regarding...

3
LATTICE DESIGN AND STUDY TOOLS REGARDING THE SUPER-B PROJECT F. M´ eot, CEA/DSM/IRFU & CNRS N. Monseu, IN2P3 LPSC INPG UJF, Grenoble Abstract Lattice design and study tools are being developed, and related beam and spin dynamics simulations are being per- formed, in the framework of the international collaboration regarding the super-B project. The present contribution re- ports on this work. INTRODUCTION The assessment and qualification of the super-B lattice parameters as produced by the on-going design studies [1] require high precision dynamic simulations and ad hoc tools. Preliminary investigations have been undertaken us- ing the stepwise ray-tracing code Zgoubi [2]. That code is based on a numerical resolution of F = qv × B, thus allowing (i) high precision magnetic fields modelling, (ii) high precision particle and spin tracking; in addition, it pro- vides most of the tools relevant to detailed numerical sim- ulation of beam and spin dynamics in super-B high-energy (HER) and low-energy (LER) rings, as spin tracking [3, 4], synchrotron radiation [5], together with a number of optical elements and other dynamics numerical tools. This paper reports on series of preliminary numerical ex- periments performed in the recent months, of which the aim is twofold, (i) evaluating the capabilities of the code as applied to super-B HER and LER lattices, (ii) starting delivering precision analysis of lattice behavior. Table 1: Super-B Parameters MAD8 Ray-tracing HER Energy /GeV 6.7 Orbit length /m 1258.3581 1258.3582 Qx, Qy 40.5750, 17.5950 [40].5750, [17].5950 Q x ,Q y 0.042, -0.0038 0.062, -0.0019 α, 14.361 10 4 , 47.88 4.371 10 4 , 47.83 Max βx, βy /m 388.57, 1126.40 388.55, 1126.35 Max Dx /m 0.6346 0.6346 LER Energy /GeV 4.18 Orbit length /m 1258.3581 1258.3582 Qx, Qy 42.5744, 18.6019 [42].5749, [18].5949 Q x ,Q y -0.620, -0.678 -0.624, -0.676 α, 14.049 10 4 , 49.69 4.053 10 4 , 49.67 Max βx, βy /m 396.54, 1013.17 387.25, 1146.77 Max Dx /m 0.5118 0.5118 0 50 100 150 200 250 300 350 400 0 200 400 600 800 1000 1200 beta x (m) s (m) ray-tracing MAD 0 200 400 600 800 1000 1200 0 200 400 600 800 1000 1200 beta y (m) s (m) ray-tracing MAD Figure 1: Optical functions in LER. Ray-tracing and MAD8 curves are superimposed. PARAXIAL MOTION The input data files to Zgoubi are obtained from transla- tion from MAD8 survey output. A preliminary step prior to further ray-tracing simulations is to get agreement between the two methods as to first and main second order param- eters. This is illustrated in Tab. 1 and Fig. 1. All optical elements in the ray-tracing methods are hard-edge but for the eight spin rotator solenoids in LER. Field in bends and multipoles is derived from the scalar potential V n (s, x, y)=(n!) 2 q=0 () q α (2q) n,0 (s) 4 q q!(n+q)! (x 2 + y 2 ) q n m=0 sin(m π 2 )x n-m y m m!(nm)! with n the multipole order, α(s) the longitudinal form fac- tor, a step function in the present study, however liable to simulate extended fringe fields [2]. Solenoids (in LER) are simulated using for the on-axis field s/ s 2 + r 2 +(L s)/ (L s) 2 + r 2 followed by Taylor expansion off-axis, with L, r respec- tively length and radius. LARGE EXCURSION MOTION, DAS Ray-tracing allows computing large amplitude motion with great precision on magnetic fields and, inherent to the numerical integration, including effects in relation with non-linearities of kinematical origin. Sample results of maximum stable amplitudes tracking are given in Fig. 2. TUPEB001 Proceedings of IPAC’10, Kyoto, Japan 1512 01 Circular Colliders A02 Lepton Colliders

Upload: others

Post on 11-Feb-2021

4 views

Category:

Documents


0 download

TRANSCRIPT

  • LATTICE DESIGN AND STUDY TOOLS REGARDING THE SUPER-BPROJECT

    F. Méot, CEA/DSM/IRFU & CNRSN. Monseu, IN2P3 LPSC INPG UJF, Grenoble

    Abstract

    Lattice design and study tools are being developed, andrelated beam and spin dynamics simulations are being per-formed, in the framework of the international collaborationregarding the super-B project. The present contribution re-ports on this work.

    INTRODUCTION

    The assessment and qualification of the super-B latticeparameters as produced by the on-going design studies [1]require high precision dynamic simulations and ad hoctools. Preliminary investigations have been undertaken us-ing the stepwise ray-tracing code Zgoubi [2]. That codeis based on a numerical resolution of �F = q �v × �B, thusallowing (i) high precision magnetic fields modelling, (ii)high precision particle and spin tracking; in addition, it pro-vides most of the tools relevant to detailed numerical sim-ulation of beam and spin dynamics in super-B high-energy(HER) and low-energy (LER) rings, as spin tracking [3, 4],synchrotron radiation [5], together with a number of opticalelements and other dynamics numerical tools.

    This paper reports on series of preliminary numerical ex-periments performed in the recent months, of which theaim is twofold, (i) evaluating the capabilities of the codeas applied to super-B HER and LER lattices, (ii) startingdelivering precision analysis of lattice behavior.

    Table 1: Super-B Parameters

    MAD8 Ray-tracingHEREnergy /GeV 6.7Orbit length /m 1258.3581 1258.3582Qx, Qy 40.5750, 17.5950 [40].5750, [17].5950Q′x,Q

    ′y 0.042, -0.0038 0.062, -0.0019

    α,√

    1/α 4.361 10−4, 47.88 4.371 10−4, 47.83Max βx, βy /m 388.57, 1126.40 388.55, 1126.35Max Dx /m 0.6346 0.6346

    LEREnergy /GeV 4.18Orbit length /m 1258.3581 1258.3582Qx, Qy 42.5744, 18.6019 [42].5749, [18].5949Q′x,Q

    ′y -0.620, -0.678 -0.624, -0.676

    α,√

    1/α 4.049 10−4, 49.69 4.053 10−4, 49.67Max βx, βy /m 396.54, 1013.17 387.25, 1146.77Max Dx /m 0.5118 0.5118

    0 50

    100 150 200 250 300 350 400

    0 200 400 600 800 1000 1200

    beta

    x (m

    )

    s (m)

    ray-tracingMAD

    0

    200

    400

    600

    800

    1000

    1200

    0 200 400 600 800 1000 1200

    beta

    y (m

    )

    s (m)

    ray-tracingMAD

    Figure 1: Optical functions in LER. Ray-tracing andMAD8 curves are superimposed.

    PARAXIAL MOTION

    The input data files to Zgoubi are obtained from transla-tion from MAD8 survey output. A preliminary step prior tofurther ray-tracing simulations is to get agreement betweenthe two methods as to first and main second order param-eters. This is illustrated in Tab. 1 and Fig. 1. All opticalelements in the ray-tracing methods are hard-edge but forthe eight spin rotator solenoids in LER. Field in bends andmultipoles is derived from the scalar potential

    Vn(s, x, y) = (n!)2{∑∞

    q=0(−)qα

    (2q)n,0 (s)

    4qq!(n+q)! (x2 + y2)q

    }

    {∑nm=0

    sin(mπ2 )xn−mym

    m!(n−m)!}

    with n the multipole order, α(s) the longitudinal form fac-tor, a step function in the present study, however liable tosimulate extended fringe fields [2]. Solenoids (in LER) aresimulated using for the on-axis field

    s/√s2 + r2 + (L− s)/√(L− s)2 + r2

    followed by Taylor expansion off-axis, with L, r respec-tively length and radius.

    LARGE EXCURSION MOTION, DAS

    Ray-tracing allows computing large amplitude motionwith great precision on magnetic fields and, inherent tothe numerical integration, including effects in relation withnon-linearities of kinematical origin. Sample results ofmaximum stable amplitudes tracking are given in Fig. 2.

    TUPEB001 Proceedings of IPAC’10, Kyoto, Japan

    1512

    01 Circular Colliders

    A02 Lepton Colliders

  • Horizontal Vertical

    -.0005

    0.0

    0.0005

    0.001

    Zgoubi|Zpop 21-Mar-10 Y’ (rad) vs. Y (m)

    -.00015

    -.0001

    -0.5E-4

    0.0

    0.5E-4

    0.0001

    0.00015

    0.0002

    Zgoubi|Zpop 17-Mar-10 Z’ (rad) vs. Z (m)

    -.0005

    0.0

    0.0005

    0.00117 Mar 10

    -.00015

    -.0001

    -0.5E-4

    0.0

    0.5E-4

    0.0001

    0.00015

    0.000217 Mar 10

    -.004 -.002 0.0 0.002 0.004

    -.0005

    0.0

    0.0005

    0.00121 Mar 10

    -.006 -.004 -.002 0.0 0.002 0.004 0.006

    -.00015

    -.0001

    -0.5E-4

    0.0

    0.5E-4

    0.0001

    0.00015

    0.000217 Mar 10

    Figure 2: Maximum stable motion in HER, from ray-tracing, dp/p = +1%, 0,−1% from top to bottom row.

    Comparisons with MAD8 have been undertaken [6] as itis often taken as a reference in the literature, it has beenobserved that strong discrepancies arise; given this, com-parisons with other numerical methods are in order.

    Dynamic apertures are obtained by repeating maximumstable H and V amplitude search above, that yields DAs inHER and LER are shown in Fig.3, as observed at centerof QDI, RF section, over 1000 turns - the damping time is≈5000 turns, see next Section.

    0

    0.1

    0.2

    0.3

    0.4

    0.5

    0.6

    0.7

    -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5

    y (c

    m)

    x (cm)

    dp +1%dp +0.5%

    dp 0dp -0.5%

    dp -1%

    0

    0.2

    0.4

    0.6

    0.8

    1

    1.2

    -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8

    y (c

    m)

    x (cm)

    dp +1%dp +0.5%

    dp 0dp -0.5%

    dp -1%

    Figure 3: 1000-turn DA in HER (top) and LER (bottom).

    EQUILIBRIUM EMITTANCES

    Stochastic synchrotron radiation and its effects on par-ticle dynamics are simulated in a very regular manner, asfollows [2, 5]. The number of photons, k, emitted in a stepΔs is given by

    Table 2: Theoretical Parameters Entering SR Simulationsin LER

    Synchronous energy Es GeV 4.18Orbit length m 1258Revolution time μs 4.20Average energy loss /particle/turn, Eloss keV 865Equivalent ρ, 88.5E4s [GeV ]/Eloss[keV ] m 31.2V̂ sin(φs) kV 865.13

    p(k) = λ−k exp(−k) / k!with λ = 20er0 β2 BρΔs / 8h̄

    √3 the mean number of

    photons per radian, and r0 = e2/4π�0m0c2 the classicalradius of the particle, e the elementary charge, h̄ = h/2π, hthe Planck constant, β = v/c, Bρ the particle stiffness. λis evaluated at each integration step from the current valuesβ, Bρ and Δs, then a value of k is drawn by a rejectionmethod. These k photons are assigned energies � = hν atrandom, in the following way. The cumulative distributionof the energy probability law p(�/�c)d�/�c writes

    P(�/�c) = 35π∫ �/�c0

    ∫∞�/�c

    K5/3(x)dx

    with K5/3 the modified Bessel function, �c = h̄ωc, ωc =2π3γ3c/2ρ the critical frequency of the radiation in con-stant field with bending radius ρ; ωc is evaluated at eachintegration step from the current values γ and ρ (assumesconstant magnetic field over Δs). Quantities in the “theo-retical” column in Tab. 3 derive from the following regularexpressions :

    Let ei[keV ] = 88.5E4[GeV ]αi/2π be the energy lossin dipole i (of N in the ring) characterized by its deviationαi with local radius ρi, and �i = 0.683E3[GeV ]/ρi[m]the average photon energy, then, energy loss per particleEloss[keV ] =

    ∑Ni ei, a sum over the N dipoles; num-

    ber of photons per particle, Nphot =∑N

    i ei/�i; beamσE/E = γ

    √Cq / 2ρ with Cq = 55h̄/32

    √3m0c ≈

    3.84 10−13; damping time τE [turns] = Es/Eloss,τE [ms] = 4.2 10

    −3Es/Eloss.The stepwise SR simulation method yields equilibrium

    emittances and other parameters as damping times, H̄, etc.Simulations have been realized with the hypotheses gath-ered in Tab. 2. SR in the sole dipoles is assumed here,whereas the code allows switching SR in other magnets.Note that the quantities in Tab. 2, Eloss = 865.13 keV,Nphot = 541, Ephot = 3.45 keV as obtained from the-

    Table 3: SR in LER, Comparing Theoretical and MonteCarlo Data

    theoretical ZgoubiEloss keV 865.1 878.4Nphot 541.4 541.3�x/π 1.8 10

    −9 1.4 10−9 (a)

    τx ms ≈ 28 (b,a)σE/E 0.67 10

    −3 0.68 10−3 (c)

    τE ms 20.3 ≈ 15 (b,c)(a) After Fig. 4-top.

    (b) Taken at 1/e of initial amplitude.

    (c) After Fig. 4-bottom.

    Proceedings of IPAC’10, Kyoto, Japan TUPEB001

    01 Circular Colliders

    A02 Lepton Colliders 1513

  • 1.E+4

    2.E+4

    3.E+4

    4.E+4

    5.E+4

    6.E+4

    7.E+4

    -.0006

    -.0004

    -.0002

    0.0

    0.0002

    0.0004

    0.0006

    0.0008 x (m) vs. Time (mu_s)

    1.E+4

    2.E+4

    3.E+4

    4.E+4

    5.E+4

    6.E+4

    7.E+4

    -0.4E-4

    -0.2E-4

    0.0

    0.2E-4

    0.4E-4

    0.6E-4

    y (m) vs. Time (mu_s)

    1.E+4

    2.E+4

    3.E+4

    4.E+4

    5.E+4

    6.E+4

    7.E+4

    -.004

    -.002

    0.0

    0.002

    0.004

    0.006 dp/p vs. Time (mu_s)

    Figure 4: 20000-turn tracking, reaching equilibrium inpresence of synchrotron radiation. Top and middle : re-spectively x-motion and y-motion, starting with pure x;bottom: starting with pure dp/p.

    450 500 550 600 650 700 750 800-1.

    -.5

    0.0

    0.5

    1.

    S s

    S x

    S y

    s (m)

    Figure 5: Spin components in IR region (spin is vertical,Sy ≡ 1, elsewhere).

    oretical radiation formulas are identical to those obtainedfrom the stochastic method. Sample behavior is displayedin Fig. 4 : a particle is launched with either x motion ornon-zero dp/p, equilibrium of respectively �x/π ≈ 1.5 nm,σdE/E ≈ 6.5 10−4 are reached.

    SPIN MOTION

    Thomas-BMT equation for spin motion is solved simul-taneously with particle motion, if requested [3]. In LER,the spin vector of the electron beam has to be longitudinalat the interaction point. This is ensured by means a rotatorcomprised of two sets of four solenoids, located on bothsides of the IR. One way to find the �n0(s) spin vector is touse the built-in fitting procedure in Zgoubi, with constraintof identical spin coordinates after one turn (Fig 6). �n0(s)is vertical in the arcs at appropriate momentum; it is thenclear that the method can be used for adjusting the rotatorand IR bends so to reach longitudinal polarization at IP.

    Typical spin behavior for the present tuning of V12 op-tics is displayed in Fig. 5; improvement is being carriedon, in particular it has to be noted that the conditions ofvertical �n0(s) in the arcs as displayed in Fig. 5 are for

    δp/p = +1.252% chromatic closed orbit.An other way to find the �n0(s) is the stroboscopic av-

    eraging [7] method, This second method permit to get ridof closed orbit, and allow full phase space study, especiallywith synchrotron radiation allowing spin diffusion, work isbeing carried on.

    -1

    -0.5

    0

    0.5

    1

    -0.015 -0.01 -0.005 0 0.005 0.01 0.015

    spin n

    0 com

    pone

    nts

    dp/p0 (%)

    Sx (longitudinal)Sy (transverse horizontal)

    Sz (transverse vertical)

    -1

    -0.5

    0

    0.5

    1

    -0.015 -0.01 -0.005 0 0.005 0.01 0.015

    spin n

    0 com

    pone

    ntsdp/p0 (%)

    Sx (longitudinal)Sy (transverse horizontal)

    Sz (transverse vertical)

    Figure 6: Spin components at the IP (up) and RF (down)section in function of momentum.

    COMMENTS

    Given the acute insight in the physics which the ray-tracing tool allows, checks as well as tight comparisonswith theory and other codes are considered worth beingcarried on. It is foreseen to explore further analysis meth-ods as frequency map techniques, based on Zgoubi outputs.Sokholov-Ternov effect will be implemented - however thetime scale in super-B lattice is a priori well beyond storagetime scale. Implementation of beam-beam effects is beingcontemplated.

    REFERENCES

    [1] http://www.pi.infn.it/SuperB

    [2] The ray-tracing code Zgoubi, F. Méot, NIM-A 427 (1999)353-356; http://sourceforge.net/projects/zgoubi/

    [3] A numerical method for combined spin tracking and ray trac-ing of charged particles, F. Méot, NIM A 313 (1992) 492-500.

    [4] Spin dynamics simulations at AGS, H. Huang et al., theseproceedings; Commissioning RHIC spin flipper, M. Baı̈ etal., these proceedings.

    [5] Numerical tools for the simulation of SR loss and induceddynamical effects in high energy transport lines, F. Méot,J. Payet, Report CEA DSM DAPNIA/SEA-00-01 (2001).

    [6] Optical parameters in the HER and LER super-B V12 rings,from ray-tracing, F. Méot, N. Monseu, Note IN2P3/LPSC-10-30 (2010).

    [7] A tracking algorithm for the stable spin polarisation field instorage rings using stroboscopic averaging, K. Heinemannand G. H. Hoffstätter, DESY Report 96-078.

    TUPEB001 Proceedings of IPAC’10, Kyoto, Japan

    1514

    01 Circular Colliders

    A02 Lepton Colliders