[PDF] - Detailed dynamics of slow extraction and its influence on transfer PDF Document

SLIDE 1

Le traitement du cancer au moyen d'accélérateurs a une longue histoire avec l'utilisation de linacs, de cyclotrons et, plus récemment, de synchrotrons. Le traitement peut se faire essentiellement de deux manières: par faisceau élargi ou balayage avec faisceau fin. Les hadrons sont particulièrement indiqués pour cette dernière méthode grâce à l'existence du pic de Bragg. La combinaison de petits faisceaux bien focalisées avec la pénétration contrôlable (pic de Bragg) et la haute efficacité biologique (ions légers) ouvre la possibilité de traiter les tumeurs les plus difficiles, radio-résistantes, de forme complexe et situées près d'organes délicats. Le synchrotron est le plus apte à accélérer des ions légers (carbone) avec variation d'énergie d'impulsion à impulsion. Le système de balayage du faisceau est contrôlé par une mesure en temps réel du flux de particules délivré au patient. C'est pourquoi le déversement doit être allongé par l'utilisation d'une éjection lente. La qualité du profil de dose dépend finalement de l’uniformité du déversement du faisceau. C'est le plus grand défi pour le synchrotron du fait de la sensibilité de l'éjection lente. Cette thèse étudie en détails l’éjection lente dans son aspect temporel aussi bien que transversal. Les conclusions de cette étude permettent de choisir, parmi les différentes possibilités, le type d’éjection lente qui procure le déversement le plus régulier. La répartition du faisceau dans l’espace de phase transversal est également importante pour le système de distribution de la dose et la planification du

traitement. A partir de la connaissance détaillée des caractéristiques du faisceau extrait par

éjection lente, un nouveau regard sur la conception des lignes de transfert est proposé, qui permet de mieux résoudre les problèmes dus à l’asymétrie et à la distribution spéciale du faisceau dans l'espace de phase. Des exemples d'applications pratiques de ces idées théoriques sont présentés.

Detailed dynamics of slow extraction and its influence on transfer line design

The treatment of cancer with accelerator beams has a long history with linacs, cyclotrons and now synchrotrons being exploited for this purpose. Treatment techniques can be broadly divided into the use of spread-out beams and scanned ‘pencil’ beams. The Bragg-peak behaviour of charged hadrons makes them ideal candidates for the latter. The combination of precisely focused ‘pencil’ beams with controllable penetration (Bragg peak) and high, radio- biological efficiency (light ions) opens the way to treating the more awkward tumours that are radio-resistant, complex in shape and lodged against critical organs. To accelerate light ions (carbon) with pulse-to-pulse energy variation, a synchrotron is the natural choice. The beam scanning system is controlled via an on-line measurement of the particle flux entering the patient and, for this reason, the beam spill must be extended in time (seconds) by a slow- extraction scheme. The quality of the dose intensity profile ultimately depends on the uniformity of the beam spill. This is the greatest challenge for the synchrotron, since slow- extraction schemes are notoriously sensitive. In this thesis, the resonant slow extraction is studied in detail both in its temporal and transverse aspects. The results from this study indicate which extraction scheme, among the various possibilities, is to be preferred for an application requiring smooth spills. The extracted beam distribution in the transverse phase spaces is also of interest for the dose distribution system and for treatment planning. Armed with a detailed knowledge of the beam characteristics from the slow extraction a novel approach to transfer line design, that copes better with the asymmetry and special phase space distributions, is described and possible implementations of the theoretical ideas are given as examples. KEYWORDS : slow extraction, resonant extraction, hadrontherapy, synchrotron, transfer line Université Claude Bernard - Lyon 1, Institut de Physique Nucléaire

SLIDE 2

N° d’ordre : 254-99 Année 1999

THESE présentée devant l’Université Claude Bernard - Lyon 1 pour l’obtention du DIPLOME DE DOCTORAT (arrêté du 30 mars 1992) présentée et soutenue publiquement le 25 Octobre 1999 par Marco Pullia

Dynamique de l’éjection lente et son influence sur les lignes de transfert

Directeur de thèse: Prof. J. E. Augustin Jury:

U. Amaldi
J. E. Augustin
P. J. Bryant

(rapporteur)

B. Chambon
P. Lefevre

(rapporteur) J Meyer

J. Remillieux
C. Steinbach

(rapporteur)

SLIDE 3

INTRODUCTION 3 BASIC ACCELERATOR PHYSICS 7 1.1 Coordinate system 7 1.2 Accelerator magnets 8 1.3 Transverse optics 10 1.4 Transfer matrix formalism 13 1.5 Off-momentum particles and dispersion function 14 THEORY OF THE THIRD-INTEGER RESONANCE 17 2.1 Sextupole magnets 17 2.1.1 Sign conventions for sextupole magnets 18 2.1.2 Normalised sextupole strength 19 2.2 Basic theory for third-integer resonance 20 2.3 Kobayashi Hamiltonian 22 2.3.1 Properties of the Hamiltonian 23 2.3.2 Stable triangle geometry 24 2.4 Transverse particle dynamics 25 2.4.1 Movement on a separatrix 25 2.4.2 Phase space maps at the sextupole 27 2.5 General first-order Hamiltonian 28 2.5.1 Closed-orbit distortions 31 2.6 Phase-space maps along the machine 32 2.6.1 General description of separatrices 32 2.7 Resonance excitation by sextupoles 34 2.7.1 Combination of sextupoles 35 2.8 Chromatic effects of sextupoles 35 2.8.1 Schemes for chromaticity correction 36 2.9 Sextupole families for small machines 37 RESONANT SLOW EXTRACTION 39 3.1 Set-up of the extraction process 39 3.2 Resonance stopband 40 3.3 Overview of extraction methods 41 3.4 Spiral step 43 3.5 Effect of the electrostatic septum 45 3.6 Separatrix geometry at electrostatic septum 46 3.7 Phase-space representation of beam and resonance 48 3.8 Hardt condition 50 3.9 Transfer between extraction septa 56 3.9.1 Minimisation of chromatic effects 57

SLIDE 4

3.10 Achromatic transfer between extraction septa 58 3.11 Comparison of extraction techniques 62 3.11.1 Varying the sextupole strength 62 3.11.2 Varying the tune distance 62 3.11.3 Varying the amplitude (zero chromaticity extraction) 65 TRANSIT TIME 67 4.1 Translated hamiltonian 67 4.2 Extraction time with constant parameters 69 4.2.1 Time from nearest fixed point 69 4.2.2 Time spent along the side of the stable triangle 72 4.2.3 Time for particles starting along the side of the stable triangle 73 4.3 Extraction time with varying parameters 73 4.3.1 Time from nearest fixed point 74 4.3.2 Time spent along the side of the stable triangle 76 4.3.3 Time for particles starting along the side of the stable triangle 76 TIME PROFILE OF THE EXTRACTED BEAM 79 5.1 “Strip” profile 79 5.2 Linear probability density 80 5.3 Inverting the Transit Time 81 5.4 Spill length 83 5.5 Width of the initial spike 85 5.6 Population of spike and tail 86 5.7 Time profile without ripple 86 5.8 Ripple 87 5.9 Comparison of extraction methods 89 5.10 Elementary "band" for a wide momentum spread 90 5.11 Simplified model 91 5.12 Uniform distribution in phase space 93 5.13 Gaussian distribution in phase space 95 5.14 Simulations of extraction with ripple 97 5.14.1 Amplitude selection 100 5.14.2 Amplitude-momentum selection (uniform distribution) 101 5.14.3 Amplitude-momentum selection (gaussian distribution) 102 5.14.4 Amplitude-momentum selection (gaussian distribution) II 103 5.15 Concluding remarks 104 TRANSVERSE ASPECTS OF THE EXTRACTED BEAM 105 6.1 Simplest model 105 6.2 Equations of motion including the vertical plane 106 6.3 Separatrix width 110 6.4 Numeric check 112 6.5 Separatrix displacement 113 6.6 Effect of vertical, closed-orbit distortion 115 6.7 Conclusions 117

SLIDE 5

A NEW APPROACH TO TRANSFER LINE DESIGN 119 7.1 The “unfilled” ellipse approach 119 7.2 Horizontal beam size control 120 7.3 Vertical beam size control 121 7.4 Initial dispersion 122 7.5 Basic design concepts 122 7.6 Telescopes 123 7.7 Rotator 124 7.8 Example modules 125 7.8.1 Matching section 125 7.8.2 Phase shifter 128 7.8.3 Stepper 128 7.8.4 Phase shifter-stepper 129 7.8.5 Closed-dispersion switching module and extension section. 129 7.8.6 Open dispersion switching module and extension section 130 7.8.7 Rotator 131 7.8.8 “Riesenrad” gantry 132 7.8.9 Resulting lines 133 7.9 Conclusions 136 CONCLUSIONS 137 ACKNOWLEDGEMENTS 139 REFERENCES 141

SLIDE 6

3

Introduction

Hadrontherapy was first proposed by Robert R. Wilson in 1946 [1]. He understood that by exploiting the typical, energy-deposition pattern for charged hadrons (known as the Bragg peak) that is shown in Figure 1, it would be possible to deliver a high radiation dose to deep tumours while limiting the dose received by the surrounding healthy tissue and the entry channel. Figure 1. The energy deposition for charged hadrons, is more favourable than for other particles. By virtue of the Bragg peak behaviour and the low lateral scattering of charged hadrons, it is possible to confine the radiation dose to the target volume of the tumour (conformal treatment) to a much higher accuracy than with conventional radiotherapy (X-rays and electrons). Hadrontherapy is particularly suited to tumours of the eye, brain, salvary glands, spinal chord, for some prostate and uterus tumours and generally for all those tumours that are situated close to a critical organ, which must be spared any radiation dose. In addition to the better conformation of the dose shared by protons and light ions, the latter have the advantage of much higher linear energy transfer (LET). This is due to ionisation and scales with the square of the ion’s charge (Z). Whereas protons mainly cause cell damage by creating free radicals, light ions cause additional damage in the form

f double strand breaks in the DNA via the ionisation electrons. A still more interesting

feature is that the LET function is higher on the Bragg peak than in the entry channel. Too large Z ions cause severe damages in the entry channel as well, resulting in higher morbidity (complications), and show a larger tail due to fragmentation. Thus the charge Z

SLIDE 7

Introduction

________________________________________________________________________________________________________________________________________________________________________________________________________

4

f the ion has to be chosen carefully. Radiobiologists and therapists have agreed that the
ptimum ion is carbon (Z = 6), which delivers a radio-biological effect (RBE) that is

several times that of protons in the Bragg peak region, while still being close to protons in the input channel. This makes carbon ions better suited to oxygen-poor and radio-resistant tumours. Light ions also exhibit fragmentation, which increases with Z. The negative aspect is the creation of an exit dose, while the positive aspect is that for certain ions it creates positron emitting isotopes that can be used with a PET camera to visualise the dose distribution directly after irradiation. Carbon is again the optimum ion that has a weak exit dose, but still creates the useful positron emitters carbon 10 and 11. Since the Bragg peak is rather narrow, a spread-out Bragg peak (SOBP) is necessary to treat thick tumours (see Figure 2). The SOBP is obtained by the superposition of many single Bragg peaks. The range (depth) of the peak can be changed in steps by inserting material in the beam path (degrader) or by the active use of the accelerator producing the

beam. To limit the number of steps, a small energy spread is created in the beam by a

ridge filter, since accelerators are, relatively speaking, mono-energetic. A variable degrader may also be used to create the full energy spread directly. In all cases, the ratio

f the entry channel dose to the peak is less favourable than for the single Bragg peak, but

it is still better than for conventional therapy. Figure 2 also shows the relative importance

f the exit dose due to fragmentation for different ions and protons

Figure 2 The Spread Out Bragg Peak needed for treating thick tumours is still preferable to conventional radiation. The transverse spreading of the beam can be obtained either with a passive system using scatterers and collimators or by actively scanning a pencil beam into the tumour volume.

SLIDE 8

Introduction

________________________________________________________________________________________________________________________________________________________________________________________________________

5 The latter allows a better conformation of the dose and does not require patient-specific

hardware. The two systems are schematically shown in Figure 3.

Figure 3. Active and passive spreading systems. Therapeutic beams can be produced with linacs, cyclotrons and synchrotrons, but when carbon ions are also to be accelerated, the synchrotron is the only practical choice. The first design of a dedicated carbon ion synchrotron was the EULIMA study [26]. To reach a depth of, say, 27 cm, carbon ions have to be accelerated to a total energy of 4800 MeV (400 MeV/u). Such a high energy makes cyclotrons and linacs unsuitable. One of the main clinical specifications is the uniformity of dose that has to be distributed

ver the treatment volume with a maximum tolerable error of r2.5%. This implies that the

dose has to be carefully measured during irradiation and that the flux of particles at the patient has to be low enough to allow such an accurate measurement. This calls for a slow extraction and for a very constant spill to avoid overdose due to fluctuations. With active scanning, the beam is distributed in small volume units named voxels (short for volume pixels). The large number of voxels needed to cover the whole tumour means a short time per voxel and imposes very tough conditions on the spill uniformity. Passive beam delivery systems are less prone to fluctuations but, on the other hand, require patient- specific hardware and the dose conformation is less accurate. The clinical specifications and the consequent physical specifications for the beam are, of course, similar in all the existing and planned facilities. The Proton-Ion Medical Machine Study (PIMMS) at the PS-Division in CERN was set up following an agreement between TERA, MED-AUSTRON, CERN and GSI. The agreed aim of the study was to investigate and design a synchrotron-based medical facility that would allow the direct clinical comparison of protons and ions for cancer treatment. The machine was to be designed primarily for high-precision active beam scanning for both protons and carbon ions, but was also to be capable of delivering proton beams by passive scattering. The main beam specifications agreed for PIMMS are resumed in Table 1. scatterer collimator Patient Patient Scanning magnets

SLIDE 9

Introduction

________________________________________________________________________________________________________________________________________________________________________________________________________

6 PIMMS performance parameters

Active scanning Passive scanning (Pencil beam) (large area beam) Extraction energies for carbon ions 120-400 MeV/u

Extraction energies for protons*

60-220 MeV 60-250 MeV Nominal treatments 60 spills of 1 s +1s to ramp up and down = 2 minute 120 spills of 0.25s +1s to ramp up and down = 2.5 minute Nominal dose delivered 2 Gray in 2 litre 2 Gray in 7.5 litre Number of protons in one spill at patient 1010 2 u 1010 Number of carbon ions in one spill at patient 4 u 108

Spot sizes variation at all energies

(full-width half-height). 4-10 mm

Field size

20 cm u 20 cm 20 cm u 20 cm

Table 1. PIMMS main performance parameters

Within the framework of PIMMS, the aim of this thesis is to study in detail the resonant slow extraction and its consequences on the extracted beam in order to choose the most suitable extraction scheme. The study of the transfer lines to the treatment rooms is also

included. A new approach to beam transport has been used to exploit the particular

transverse characteristics of the extracted beam. Chapters 1 to 3 act as an introduction, summarize and discuss the state of the art and define the main quantities, Chapters 4 to 7 describe the original work of this thesis.

SLIDE 10

7

Chapter 1 BASIC ACCELERATOR PHYSICS

The guiding and focusing of a charged particle beam in a circular accelerator rely on a series of magnetic elements, separated by field-free drift spaces, that form the accelerator

lattice. The design of the lattice is one of the first tasks for the accelerator designer. One

normally starts with a simplified structure, containing only ideal magnetic dipoles and

quadrupoles. This basic structure is referred to as the linear lattice and once it has been

determined the character of the machine is more or less fixed. In a separated function synchrotron, dipole magnets bend the particles onto circular trajectories and quadrupole magnets are needed for the focusing of particles with small deviations compared to the ideal trajectory. The equilibrium orbit can be defined as the orbit of a particle with the design momentum p0 that closes upon itself after one turn and is stable. The motion of a particle in an accelerator is then conveniently described by the deviations of its trajectory with respect to this equilibrium orbit. 1.1 Coordinate system Generally in a synchrotron, there is bending in only one plane referred to as the horizontal plane which contains the equilibrium orbit. For a description of particle motion a right- handed curvilinear coordinate system (x, s, z) following the equilibrium orbit is used. The azimuthal coordinate s is directed along the tangent of the orbit. The radial transverse direction x is defined as normal to the orbit in the horizontal plane and the vertical transverse direction is given by the cross product of the unit vectors x

ˆ and s ˆ . The

vectors s

ˆ and z ˆ then define the vertical plane. The local radius of curvature U0 and the

bending angle s/U are defined as positive for anticlockwise rotation when viewed from positive z. The choice of the coordinate system is shown in Figure 1.1.

Figure 1.1. Curvilinear coordinate system.

s

x ˆ z ˆ s ˆ

local centre

f curvature

x z equilibrium

rbit

particle trajectory U0 horizontal plane

SLIDE 11

Chapter 1 Basic accelerator physics

________________________________________________________________________________________________________________________________________________________________________________________________________

8 1.2 Accelerator magnets For this study, the magnetic elements are considered as purely transverse, two- dimensional fields in the curvilinear coordinate system that follows the equilibrium orbit. It is assumed that each element can be decomposed into a harmonic series that extends only as far as the sextupole component and that each element can be represented by a ‘block’ of field that is discontinuous in the axial direction. This is the so-called hard-edge

approximation. For a final design, end-field corrections will have to be applied to improve

the accuracy of the hard-edged model. The equilibrium orbit corresponds to a particle with the nominal momentum and starting conditions. Particles of the same momentum, but with small spatial deviations will

scillate about this orbit. The equilibrium orbit is variously known as the central orbit, the

reference orbit and the closed orbit. In the hard-edged model, this orbit is a series of straight lines connecting circular arcs of cyclotron motion (see Figure 1.2).

Figure 1.2. Cyclotron motion for a particle of mass Amav with charge ne.

The force experienced by a charged particle, moving in a magnetic field is given by the Lorenz equation as

B v p F u

q

t d d , (1.1)

where v is the particle velocity, q its charge and B the magnetic field. Since the Lorenz force always acts perpendicular to the particle velocity, only the direction of movement changes while the magnitude of the velocity is constant. The nominal trajectory in a circular accelerator is defined by the main bending magnets which provide a dipolar field,

z B ˆ B . (1.2)

Equating the magnetic deflection to the centrifugal force (scalar values) gives,

ne v Am B v Am nev B

z z av 2 av

U

U , (1.3)

where A is the atomic mass number, mav is the average relativistic mass per nucleon, n is the charge state of the particle, e is the elementary charge, and U is the radius of the cyclotron motion. The reluctance of the particle to be deviated is characterised by BU which is known as the magnetic rigidity.

Bz Amavv2/U v U Bznev

SLIDE 12

Basic accelerator physics Chapter 1

________________________________________________________________________________________________________________________________________________________________________________________________________

9 The application of a sign convention for the magnetic rigidity is usually avoided and an ‘engineering’ formula is quoted in which the momentum is always taken as positive.

Magnetic rigidity: ] [GeV/c 299792458 . 1 ] [ ] [

av

p n A m T Bz

U

(1.4)

where pav is the average momentum per nucleon. From (1.4) it can be seen that the higher the particle momentum, the higher the magnetic field needed to keep the particle on the nominal trajectory. Therefore, the magnetic fields in an accelerator have to be continuously adjusted, according to the actual particle momentum during acceleration. In

rder to avoid the momentum dependency when characterising magnetic elements in an

accelerator, the fields are normalised w.r.t. the magnetic rigidity (particle momentum). For a dipole magnet the normalised bending strength h, is given by

U 1 h . (1.5)

The absolute field level required to keep a particle with a certain momentum on the nominal trajectory is then simply found by multiplying (1.5) with the magnetic rigidity (1.4). The equilibrium orbit in a synchrotron is defined by the main dipole magnets. In general particle trajectories will deviate slightly from this ideal design orbit. To maintain stability about the design orbit, focusing forces are needed. These are provided by quadrupole magnets. A quadrupole magnet has four poles with a hyperbolic contour and the magnetic field in the current-free region of the magnet gap can be derived from a scalar potential [2],

z x G z x

)

) , ( , (1.6)

where G is known as the quadrupole gradient, being defined as:

d d

x

B G

z

. (1.7)

The equipotential lines are the hyperbolae x·z = const and the field lines are perpendicular to them. The horizontal and vertical magnetic fields in a quadrupole magnet are linear in the deviation from the magnet centre:

x G B z G B

z x

and

. (1.8)

The transverse forces acting on a particle with a deviation (x,z) from the equilibrium orbit are obtained with (1.1) as

qvGx z x qvB F

z x

) , ( (1.9)

and

qvGz z x qvB F

x z

)

, ( . (1.10)

SLIDE 13

Chapter 1 Basic accelerator physics

________________________________________________________________________________________________________________________________________________________________________________________________________

10 Depending on the sign of G, the force will be focusing in the horizontal direction and defocusing in the vertical direction or vice versa. A quadrupole that focuses horizontally (G < 0), is called a focusing quadrupole and a quadrupole that defocuses horizontally (G > 0), is called a defocusing quadrupole. By alternating focusing and defocusing quadrupoles an overall focusing effect can be obtained in both the horizontal and the vertical plane (alternating gradient focusing). An important property of a quadrupole magnet is that the horizontal force depends

nly on the horizontal and not on the vertical position of the particle trajectory. Similarly,

the vertical component of the Lorentz force depends only on the vertical position. The consequence is that in a linear machine, containing only dipolar and quadrupolar fields, the horizontal and vertical motions are completely uncoupled. In analogy to the normalised bending strength of a dipole magnet (1.5), the quadrupole gradient can be normalised w.r.t. the magnetic rigidity, thus defining the normalised quadrupole strength,

U B G k . (1.11)

1.3 Transverse optics When analysing the transverse motion in an accelerator, it is practical to use the distance s measured along the trajectory, instead of the time as independent variable,

s v t s s t d d d d d d d d . (1.12)

It is usual to consider only dipole and quadrupole fields (linear lattice), which leads to uncoupled motions in the two transverse planes. The equations of motion are obtained from (1.1) while retaining only first order terms in the transverse coordinates [3]. The horizontal and vertical motions can be represented by one simple expression valid for both planes,

) ( ) ( ) (

c

c s y s K s y (1.13)

where, for the horizontal motion,

) ( ) ( ) ( ) ( and

2

s k s h s K s K x y

x

{

(1.14)

and, for the vertical motion,

) ( ) ( ) ( and s k s K s K z y

z

{ . (1.15)

For the hard-edged model, the focusing function Ky(s) is piecewise constant along the trajectory and for a circular machine, it has the periodicity of the lattice,

) ( ) ( s K L s K

y y

,

(1.16)

where the accelerator is composed of N identical sections or cells, with C = N·L and C being the circumference of the machine.

SLIDE 14

Basic accelerator physics Chapter 1

________________________________________________________________________________________________________________________________________________________________________________________________________

11 The motion of mono-energetic particles about their equilibrium orbit is known as the betatron motion and it is usual to parameterise the motion so that its pseudo-harmonic behaviour is brought into evidence. The first stage is to express the motion in a form developed by Courant and Snyder [4],

Betatron motion:

s

B s A s y d cos ) ( V E V E , (1.17)

where y(s) represents either transverse coordinate as a function of the distance s along the equilibrium orbit, A and B are constants depending on the starting conditions and E(s) is the betatron amplitude function (dimension length). The phase P(s) of the pseudo

scillation is given as
s

s d V E V P . (1.18)

To complete this description, the derivative of E(s) is added in the relations,

s s d d 2 1 E D

(1.19)

and

E D J

2

1 s . (1.20)

The expressions (1.17)-(1.20) all depend on a knowledge of E(s). Although the analytic solution for E(s) is more complicated than the original motion equation, it is reasonably easy to evaluate this function numerically and to tabulate it for any lattice. The parameters D(s), E(s) and J(s) are collectively known as the Courant and Snyder parameters, or more usually the Twiss parameters. The above parameterisation is now so commonplace that it is the starting point for nearly all lattice design. The pseudo-harmonic motion can be further transformed into a simple harmonic motion by returning to equation (1.17), introducing P and differentiating, to obtain,

B

A s y

P

E cos ) ( (1.21)

B

A B A s y

c

P E P E D sin 1 cos . (1.22)

The phase terms can be extracted and used to define new coordinates Y(P) and Yc(P) that are known as normalised coordinates,

E

P P 1 cos s y B A Y

(1.23)
E

D E P P s y s y B A Y

c
c

sin . (1.24)

SLIDE 15

Chapter 1 Basic accelerator physics

________________________________________________________________________________________________________________________________________________________________________________________________________

12 It is useful to represent real-space coordinates by lower case symbols and normalised coordinates by upper case symbols (x o X, etc.). Normalised coordinates use the phase advance as independent variable and real-space coordinates use the distance. The transformations between the two systems are conveniently expressed in matrix form as,

Real to normalised:

s

y y Y Y d / d / / 1 d / d E E D E P (1.25) Normalised to real:

P

E E D E d / d / 1 / d / d Y Y s y y . (1.26)

The elimination of the phase advance from equations (1.23) and (1.24) yields an invariant of the motion,

) constant ( 2

2 2 2 2 2

c

c
c
y

y y y Y Y A E D J . (1.27)

Equation (1.27) is in fact the equation of an ellipse in (y, yc) phase space. The constant A2 equals the (area/S) of the ellipse described by the betatron motion in either the normalised (Y , Yc) or the real (y , yc) phase spaces (see Figure 1.3). When referring to a single particle, this area is sometimes called the single-particle emittance. When referring to a beam it is known as the emittance. The situation with a beam is complicated by the definition used for deciding the limiting ellipse that defines the area. This may be related to a number of standard deviations of the beam distribution, or the overall maximum. It is useful to note that the emittances in real and normalised phase spaces are equal and to re- express (1.27) as,

Emittance: S S E S E S

y y y y y

E Y E y Y y E

ˆ

and ˆ ˆ ˆ

2 2

, (1.28)

where y

ˆ and Y ˆ are the maximum excursions in real and normalised phase spaces that

define the beam-size at the observer’s position.

0.0015

0.0015

0.0015

0.0015

0.0015

0.0015

0.0015

0.0015 Figure 1.3. Phase-space ellipses in real and normalised phase spaces.

(EyEy/S) = y

ˆ

(EyJy/S)

Dy(Ey/S Ey)

(Y2+Yc 2) (Ey/S) =Y

ˆ

y [m] yc [rad] Yc [m1/2] Y [m1/2]

SLIDE 16

Basic accelerator physics Chapter 1

________________________________________________________________________________________________________________________________________________________________________________________________________

13 It should be noted that during acceleration the emittance decreases and therefore the phase space density changes. This is not against Liouville’s theorem, but due to the definition of the transverse phase space. The yc-coordinate, describing the divergence of a particle with respect to the nominal orbit,

v v t y v s y y

y

c d d 1 d d , (1.29)

is not a generalised momentum as defined in classical mechanics. As long as the momentum is constant, the difference is just a scaling factor but during acceleration this factor changes and so does the emittance. The physical reason is the increase of the longitudinal momentum during acceleration, whereas the transverse momentum is

constant. Therefore the divergence vy/v decreases and the beam shrinks, which is known

as adiabatic damping. Normalising the emittance w.r.t. the particle momentum gives a constant of the motion, the so called normalised emittance,

y y n

E E

rel rel ,

J E , (1.30)

where Erel and Jrel are the relativistic parameters. 1.4 Transfer matrix formalism Another basic expression that is needed for the present study is the general transfer

matrix. This can be derived by expanding (1.21) into two terms,

P E P E sin cos B A s y

(1.31)

where A and B are new constants. Differentiation of (1.31) with respect to s gives,

P

D P E P P D E sin cos sin cos

c

B A s y . (1.32)

The constants A and B can be replaced using the initial conditions at s = s1, choosing P1 = 0

1 1 1 1 1 1 1

and E D E E y y B y A

c

, (1.33)

to give the general transfer matrix from position s1 to position s2. The phase advance from s1 to s2 is written as 'P,

>

@

'
'

'

'
'

'

'
P

D P E E P D D P D D E E P E E P D P E E sin cos cos sin 1 1 sin sin cos : matrix transfer General

2 2 1 1 2 2 1 2 1 2 1 1 1 2 2 1

s s M . (1.34)

SLIDE 17

Chapter 1 Basic accelerator physics

________________________________________________________________________________________________________________________________________________________________________________________________________

14 The phase-space coordinates at position s2 are then given by

c
c

1 1 2 1 2 2

) ( y y s s y y M . (1.35)

When equation (1.35) is applied to a full turn in a ring, the input conditions equal the

utput conditions (D = D1 = D2, E = E1 = E2, 'P = SQ), so that
Q

Q Q Q Q Q S D S S J S E S D S 2 sin 2 cos 2 sin 2 sin 2 sin 2 cos

turn 1

M , (1.36)

where Q is known as the betatron tune and is equal to the number of transverse

scillations the particle makes during one turn in the machine. Expression (1.36) describes

the evolution of the phase-space coordinates of a particle at a certain position in the

machine. A plot of the coordinates for a large number of turns gives a phase-space
trajectory. In a linear machine phase-space trajectories are always of elliptical shape; the
rientation of the ellipses at any position s in the machine is determined by the local Twiss

parameters and the beam size is found with the emittance according to Figure 1.3. The general transfer matrix for normalised coordinates is simply a 2x2 rotation matrix describing a clockwise rotation by the phase advance 'P between position s1 and s2,

'

'

'

'

P

P P P cos sin sin cos ) (

2 1 N

s s M . (1.37)

The single-turn matrix for normalised coordinates is

Q

Q Q Q S S S S 2 cos 2 sin 2 sin 2 cos

turn N,1

M . (1.38)

The 2x2 transfer matrix formalism is commonly known as Twiss-matrix formalism. 1.5 Off-momentum particles and dispersion function A simple extension of (1.17) allows the motion of particles with different momenta to be described,

motion

Dispersion motion Betatron

) ( ) ( d cos ) ( ) ( p p s D B s s A s y

s

G E V E

(1.39)

where D(s) is known as the dispersion function and Gp/p = (ppart. - p0)/p0 is the relative momentum deviation of the particle. The dispersion is created by the momentum dependency of the bending radius in dipole magnets and appears therefore only in the plane

f bending (generally the horizontal plane). The equilibrium orbit for an off-momentum

particle is to first order displaced from the central orbit by the product of the dispersion function and the momentum deviation,

p p D y p p D y

y y

G G c c

EQ.O EQ.O

and . (1.40)

SLIDE 18

Basic accelerator physics Chapter 1

________________________________________________________________________________________________________________________________________________________________________________________________________

15 An analytic derivation of the dispersion function is given in [4], but it is again common practice to rely on lattice programs to supply numerical listings of D(s) and its derivative with distance, as for the betatron amplitude function. The momentum deviation Gp/p is treated as a quasi-variable and particles are transferred through the lattice with 3u3 matrices of the form,

1 23 22 21 13 12 11 2

/ 1 /

c
c

p p y y m m m m m m p p y y G G (1.41)

where m11, m12, m21 and m22 are the coefficients of the general transfer matrix (1.34) and m13, m23 are additional dispersion coefficients (see Section 3.10). The dispersion vector in the form (D, Dc, Gp/p=1) also propagates through the lattice according to (1.41). The normalised form of the dispersion function (Dn, Dcn), similar to (1.25), will also be frequently used,

s

D D D D d / d / / 1 d / d

n n

E E D E P . (1.42)

It should be mentioned that (1.41) is strictly applicable to only small values of Gp/p and to linear lattices. For trajectories with more than a few per mil momentum deviation,

r for trajectories that pass through non-linear magnetic lenses, it is advisable to perform a

numerical tracking if the orbit position is required to a high precision. Exactly this situation arises when calculating the position and angle of the separatrices for a resonant extraction [5] and it is useful to be able to incorporate the more exact tracking method into the general transfer matrix as,

2

1

s s M

> @

>

@

c
c

c

1

/ / / /

1 EQ, 22 1 EQ, 21 2 EQ, 22 21 1 EQ, 12 1 EQ, 11 2 EQ, 12 11

p p y m y m y m m p p y m y m y m m G G (1.43)

where yEQ and ycEQ are the position and angle of the tracked off-momentum equilibrium

rbit.

A more detailed and complete introduction to accelerator physics can be found in References [3,6,7].

SLIDE 19

17

Chapter 2 THEORY OF THE THIRD-INTEGER RESONANCE

The third-integer resonance can be used to extract particles from a synchrotron during a large number of turns [8]. In this case, the slowly extracted beam is known as the spill. In a medical machine, it is essential to extend the extraction time (spill time) to about one second (~106 turns), in order to facilitate the measurement of the radiation dose delivered to the patient. This is done by making particles of the circulating beam unstable in a controlled way. For this, the betatron tune in the extraction plane (in general the horizontal plane is used) has to be close to the resonance condition, 3Qh = integer. The resonance is excited by sextupole magnets. 2.1 Sextupole magnets In the current-free region of a magnet gap, the field can always be derived from a scalar potential ). Assuming that the magnetic field has only transverse components, the scalar potential of a magnet with 2m poles is given by [2],

magnets

Normal m m magnets Skew m m

) Im( ) Re( iz x B iz x A

)

. (2.1)

The field components are obtained by differentiation according to

x Bx w w )

and

z Bz w w )

.

(2.2)

The transverse fields in a normal sextupole magnet (m = 3) are then found to be

B x z

B xz

x

, 6

3

B x z

B x z

z 3

, ( )

3

2 2 .

(2.3)

The relationship between the coefficient B3 and the magnetic field is found by comparing (2.3) to the Taylor expansion of the magnetic field in the horizontal plane,

B x z B B x x B x x

z z z

d d ( , ) ! ! d d ...

1

1 1 2

2 2 2

(2.4)

so that

B d B dx

z 3 2 2

1 6

.

(2.5)

The vertical and horizontal fields in a sextupole can then be written as

B

B x x z B B x xz

z x

1

2

2 2 2 2 2 2

d d and d d

z z

. (2.6)

SLIDE 20

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

18 The effect of a sextupole on a particle trajectory can be described in a simple way by considering the magnet as a thin lens. The integrated magnetic field seen by the particle is given by the product of the magnet length "S and the field at the average transverse particle

position. As the lens is considered to be infinitely thin, the particle position remains

unchanged during the passage, i.e. 'x = 'z = 0. The sextupole gives a kick to the particle that changes the direction of the trajectory. With the magnetic rigidity |BU| of the particle the thin-lens kicks given by the sextupole are

2

2 S 2 2 2 z 2 S

2 1 d d 2 1 z x k z x x B B B B x

S z

c
c

' " " " U U (2.7)

and

xz k xz x B B z c

c

'

S 2 z S

d d " " U , (2.8)

where kc is the normalised sextupole gradient. 2.1.1 Sign conventions for sextupole magnets To define a sign convention for sextupole magnets, the trajectory of a particle in the horizontal plane (i.e. z = 0) is considered. Figure 2.1 illustrates the sign convention and shows the corresponding magnet geometry. x If the sextupole focuses for positive x, then it can be said to be an F-sextupole. x If the sextupole defocuses for positive x, then it can be said to be a D-sextupole.

Figure 2.1. Sign convention for sextupole magnets.

S N S N z N S B x s z Bz Bz

F-Sextupole (Bz negative, kc negative)

x s z Bz Bz N N S z S N B S

D-Sextupole (Bz positive, kc positive)

SLIDE 21

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

19 2.1.2 Normalised sextupole strength In order to describe the effect of a thin-lens sextupole in normalised coordinates, the kick introduced by the sextupole and the particle position have to be expressed in terms of normalised coordinates. From the above, the effect of a thin lens sextupole in real space is

'x 0 and

2

2 2 z 2 S S

d d 2 1 z x x B B B B x

z

c

' U U " " . (2.9)

Many problems are simplified by a transformation into normalised coordinates that converts the basic betatron motion into a simple harmonic motion. The normalisation matrix and its inverse are quoted below,

y

y y y

E E D E 1

y N,

M

y

y y y

E E D E 1 ) (

1 y N,

M and

c
c

y y Y Y

y N,

M . (2.10)

When applied to (2.9) the relations between normalised and unnormalised coordinates are found to be,

' ' c c x X

x

1 E x X

x

E z Z

z

E . (2.11)

Thus the effect of the sextupole in normalised coordinates appears as,

'X 0

c

'

2 2 2 2 2 z 2 S 2 / 3 x

d d 2 1 Z X S Z X x B B X

x z x z

E E E E U E " (2.12)

where S is the normalised sextupole strength,

k x B B S

x x

c

S

2 / 3 2 z 2 S 2 / 3

2 1 d d 2 1 " " E U E . (2.13)

The expressions describing the effect of the thin-lens sextupole on the vertical trajectory are obtained in the same way as

'Z 0 XZ S XZ x B B Z

x z x z

E E E E U E 2 d d 2 1 2

2 z 2 S 2 / 3 x

c

' " . (2.14)

From the above expressions, it can be seen that a sextupole couples the horizontal and vertical motion of a particle. The strength of the coupling is proportional to the ratio of the vertical and horizontal beta functions (Ez/Ex) at the position of the sextupole. For a horizontal extraction, choosing a small Ez/Ex ratio and provided the vertical tune does not satisfy a resonance condition, the influence of the vertical motion can be neglected to first

rder. For this reason only the horizontal motion is considered in the following analysis.

SLIDE 22

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

20 2.2 Basic theory for third-integer resonance A simple method that gives an immediate insight into the physics of resonant extraction was developed by Kobayashi [9,10]. The method describes the effect of the sextupole as a perturbation to the linear machine. The general transfer matrix Mk for normalised coordinates, describing k turns in the machine is given by:

)

( 2 cos ) ( 2 sin ) ( 2 sin ) ( 2 cos

h h h h

kQ kQ kQ kQ

k

S S S S M . (2.15)

Consider a particle with a horizontal betatron tune close to a third-integer, Qh = n r 1/3 + GQ, (n is integer, |GQ|<<1/3), where GQ is defined as the tune distance of the particle from the resonance,

GQ Q Q

particle resonance

,

(2.16)

the explicit transfer matrices for one, two and three turns in the unperturbed machine can be written as:

>

@

>

@

>

@

>

@

M1 2 1 3 2 1 3 2 1 3 2 1 3 1 2 3 2 3 2 1 2 r

r
r
r
#
r
cos

/ sin / sin / cos / / / / / S G S G S G S G n Q n Q n Q n Q # (2.17)

>

@

>

@

>

@

>

@

M2 4 1 3 4 1 3 4 1 3 4 1 3 1 2 3 2 3 2 1 2 r

r
r
r
#
r
cos

/ sin / sin / cos / / / / / S G S G S G S G n Q n Q n Q n Q # (2.18)

>

@

>

@

>

@

>

@

M3 6 1 3 6 1 3 6 1 3 6 1 3 1 6 6 1 r

r
r
r
#
cos

/ sin / sin / cos / S G S G S G S G SG SG n Q n Q n Q n Q Q Q . (2.19)

For the following calculations, the tune distance GQ is kept only in the transfer matrix for the full three turns. The coordinates of the particle after one turn are then found as:

X X X X c

#
r
c
1

1 2 3 2 3 2 1 2 / / / / # (2.20)

and similarly for 2 turns

X X X X c

#
r
c
2

1 2 3 2 3 2 1 2 / / / / # . (2.21)

By defining the small quantity H = 6SGQ, which will be referred to as the modified tune distance, it can be seen that a particle with exactly resonant tune (i.e. H = 0) will return to its initial position every three turns.

X X X X ' '

#
3

1 1 H H (2.22)

SLIDE 23

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

21 The effect of the sextupole during three turns in the machine is now calculated as a perturbation by linearly adding: (A) The effect of 3 turns with a sextupole placed after the 3rd turn, M Sextupole

3

(B) the effect of 3 turns with a sextupole placed after the 2nd turn, M Sextupole M

2 1

(C) the effect of 3 turns with a sextupole placed after the 1st turn.

M Sextupole M

1 2

as shown schematically in Figure 2.2 below.

(A) (C) (B) Add

Figure 2.2. Perturbative treatment of the sextupolar effect during three turns.

With the effect of the sextupole as derived earlier,

'X 0 ' c X SX 2, (2.23)

the calculation of the three terms gives: (A) 3 turns + sextupole

2

3 3

X X S X X X X X X c

c
c

c

H

H H (2.24)

(B) 2 turns + sextupole+1 turn

X X X X X S X X X X X X X S X X

3 2 3 2

1 2 1 2 3 2 3 2 3 2 1 2 1 2 3 2 3 2 1 2 3 2 1 2 3 2 1 2 1 2 3 2

c
r

r

c
c
c
c
r
c
c
#

# # # # (2.25)

(C) 1 turn + sextupole+2 turns

X X X X X S X X X X X X X S X X

3 2 3 2

1 2 1 2 3 2 3 2 3 2 1 2 1 2 3 2 3 2 1 2 3 2 1 2 3 2 1 2 1 2 3 2

r

c

c
r

c

c r
r

c

c
r

c

#

# # (2.26)

SLIDE 24

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

22 Subtracting the initial coordinates (X0, Xc0) and adding the three terms (A), (B) and (C) while keeping only first-order correction terms in H gives

' ' X X S X X S X X X X SX S X X S X X

3 2 2 3 2 2 2

3 2 1 2 3 2 3 2 1 2 3 2 1 2 1 2 3 2 1 2 1 2 3 2 c r

c
r

c

c
c
r

c

H

H # # # . (2.27)

The cancellation of signs shows that there is no fundamental difference between the 1/3rd and 2/3rd resonances.

' ' X X S X X S X X X X SX S X X S X X

3 2 2 3 2 2 2

3 2 1 2 3 2 3 2 1 2 3 2 1 2 1 2 3 2 1 2 1 2 3 2 c

c
c
c
c
c
H

H (2.28)

Finally, the following expressions for the change of position and divergence of the particle within three revolutions in the machine including the sextupolar effect as a perturbation are

btained:
'

' X X SX X X X S X X

3 3 2 2

3 2 3 4 c c c

c

H H (2.29)

2.3 Kobayashi Hamiltonian The time needed for three revolutions in the machine is short compared to the spill time and can be safely used as the basic time unit. The elementary changes occurring in this time are also the smallest that need to be resolved to understand the physics of the

extraction. Thus the subscripts are no longer needed and (2.29) can be treated as a

continuous function that is derived from a Hamiltonian H, such that

2

2 turn) (3 1 3 turn) (3 1 3

4 3 ' 2 3 X X S X X t X X X SX X X t X X

t t

c

'

c '

'

c

c

c

'

'

'

' '

H w w H w w H H (2.30)

The Hamiltonian is then found by integrating the above partial differentials,

H
c
c

H 2 4 3

2 2 2 3

X X S XX X . (2.31)

It should be noted that in this formulation time is dimensionless.

SLIDE 25

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

23 2.3.1 Properties of the Hamiltonian The Hamiltonian is time independent and therefore a constant of the motion. The plot of this function, for constant values of H, creates maps of particle trajectories in normalised phase space at the sextupole. The above Hamiltonian consists of two parts. The term

H

2

2 2

X X

c

(2.32)

describes particle motion in the linear, unperturbed machine (i.e. S = 0). These trajectories are circles with the radius (2H/H) in normalised phase space. The perturbation term

S

XX X 4 3

2 3

c (2.33)

is proportional to the sextupole strength and distorts the circles. The resulting trajectories are symmetric with respect to the X-axis since Xc appears only in a quadratic form. A change in sign of the modified tune distance H corresponds to a reflection about the Xc axis and a change in sign of the normalised sextupole strength S is equivalent to a rotation

f the phase-space trajectories by 180º as illustrated in Figure 2.3.

All the properties of the system can be derived from the Hamiltonian. In particular, when H has the value [(2H/3)3/S2], it factorises into three straight lines,

S X X X S X X S 4 6 3 4 3 3 4 3

c
c
H

H H . (2.34)

These three lines are referred to as the separatrices that define the boundaries between stable and unstable regions in phase space. The size of the stable region is determined by the ratio |H/S|. For a particle that has exactly the resonance tune, there are no stable

trajectories. Figure 2.4 shows phase-space maps at the sextupole illustrating the

dependency of the stable region on the particle tune, represented by H.

1 0.5 0.5 1 1 0.5 0.5 1

Norm. phase space H/S0 > 0

X’ X

1 0.5 0.5 1 1 0.5 0.5 1

Norm. phase space H/S0 < 0

X’ X

Figure 2.3. Particle trajectories at the sextupole.

SLIDE 26

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

24

1 0.5 0.5 1 1 0.5 0.5 1

Norm. phase space H = 0; S0

X’ X

1 0.5 0.5 1 1 0.5 0.5 1

Norm. phase space 0.5(H/S0) > 0

X’ X

Figure 2.4. Particle trajectories at the sextupole.

2.3.2 Stable triangle geometry The Hamiltonian yields the equations of the three separatrices from (2.34) that divide up the phase space and hence the geometry of the phase space as shown in Figure 2.5. The equations of the separatrices denoted by A, B and C are explicitly given by:

A: X S 2 3 H B: c 3 2 1 2 2 3 X X S H C: 3 2 1 2 2 3 c X X S H . (2.35)

A B C h X’ X P0 Unstable region surrounding stable triangle stable region P1 P3 P2

Figure 2.5. Geometry of separatrices and stable triangle for a third-integer resonance.

SLIDE 27

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

25 Fixed points in phase space are derived from the Hamiltonian by solving the system

w w w w H H c X X 0. (2.36)

The phase-space map for a third-integer resonance contains four fixed points. The three unstable fixed points are the crossing points of the separatrices, the stable fixed point P0 is given by the equilibrium orbit. The coordinates of the unstable fixed points are:

P1 = 0 4 3 , H S

P2 =
2

3 2 3 H H S S , P3 =

2

3 2 3 H H S S , . (2.37)

The geometry of the stable triangle is conveniently described by introducing the apothem h

f the stable triangle,

h S S Q 2 3 4 H S G . (2.38)

x For a positive value of h, the upright separatrix corresponds to a negative value of X. x From the equations of the separatrices and the coordinates of the fixed points it can be seen that a change in sign of h (i.e. either S or H changes sign) is equivalent to a 180q rotation of the stable triangle around the origin or a reflection about the Xc axis. x The area of the stable triangle is given by:

Area of triangle

Stable emittance = 3 3 48 3

2 2 2

h S Q S G S (2.39)

The trajectory of a particle in normalised phase space in the unperturbed machine (i.e. S=0) is a circle. The area of this circle is the so called single particle emittance. When the sextupole is raised to the nominal strength in an adiabatic way (i.e. during a number of turns much larger than the revolution period in phase space), a particle with a certain GQ will remain stable if its emittance in the unperturbed machine was not bigger than the area

f the stable triangle calculated for the nominal sextupole strength. Therefore the area of

the stable triangle is sometimes referred to as the stable emittance. 2.4 Transverse particle dynamics Up to now, only the time-independent geometry of the phase space for a third-integer resonance has been considered, but it is also necessary to analyse the dynamic aspects and to describe particle movement in phase space. 2.4.1 Movement on a separatrix The change of the particle coordinates within three turns in the machine is derived from the Hamiltonian:

'

' X X X SXX X X X S X X

3 3 2 2

3 2 3 4 c c c

c

w w H w w H H H ' (2.40)

SLIDE 28

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

26 All separatrices are essentially equal (see Section 2.6) and there is no loss in generality if

nly the upright separatrix (A) from Figure 2.5 is considered. By inserting the equation of

the separatrix

X S 2 3 H (2.41)

into (2.40), the phase-space velocity for particles moving on the separatrix is obtained as,

' c

c

X S SX H2

2

3 4 and 'X 0 . (2.42)

It can be seen that there are two contributions to the phase-space velocity in (2.42); the constant term H2/S and the quadratic term dependent on the square of Xc. The contributions

f the two terms are of opposite sign summing to zero at the fixed points P2 and P3. The

velocity therefore changes sign when crossing a fixed point. This means that one fixed point acts as an attracting pole with particles on either side moving towards it, while it is the contrary for the other. Thus the separatrix is divided into three sections: the incoming section, the side of the stable triangle and the outgoing section. This is shown schematically in Figure 2.6, the arrows indicate the direction of particle movement.

Figure 2.6. Particle movement on a separatrix

The separatrices divide the phase space into stable and unstable regions. Trajectories within the stable triangle are closed and correspond to stable particles whereas trajectories

utside are not closed and belong to unstable particles. However, the movement of a

particle on any phase-space trajectory follows the same direction as the movement on the closest separatrix. Particle movement on stable trajectories in the sense of clockwise or anticlockwise rotation is determined only by the sign of the modified tune distance H and is independent

f the sign of the sextupole strength S. This follows immediately from (2.41) and (2.42) as

a change of the sign of S reverses the velocity on the separatrix but at the same time rotates the phase-space trajectories by 180q, resulting in an unchanged direction of rotation

f the particles.

fixed point separatrix

utgoing

section incoming section side of stable triangle velocity 'X’ fixed point X’

SLIDE 29

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

27 2.4.2 Phase space maps at the sextupole The geometry of the phase-space map for a third-integer resonance is determined by the normalised sextupole strength S and the modified tune distance H of the particles. Figure 2.7 summarises the possible geometry of the separatrices at the sextupole, the arrows indicate the direction of particle movement. x The size of the stable triangle is proportional toHS . x The orientation of the stable triangle at the sextupole is determined by the sign of H/S. A positive value corresponds to Figure 2.7 (a),(d), with the upright separatrix at a negative value of X. Changing the sign rotates the phase-space map by 180q. x The sign of H determines absolutely the direction of rotation in phase space within the stable triangle. A positive value means that the operational tune is above resonance and particles within the stable triangle will perform a clockwise rotation. x Finally, the direction of the outgoing section of the upright separatrix is always downwards for positive values of S and upwards for negative values.

Figure 2.7. Overview of stable triangle geometry at the resonance sextupole.

X’ X X’ X X’ X X’ X H H = 6SG SGQ > 0 tune above resonance H H = 6SG SGQ < 0 tune below resonance S > 0 S < 0 (a) (b) (c) (d)

SLIDE 30

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

28 2.5 General first-order Hamiltonian In the Kobayashi theory as derived above, it is assumed that the equilibrium orbit is coincident with the centre of the sextupole. However, the equilibrium orbit of particles with a momentum deviation compared to the design momentum of the accelerator is displaced to first order from the origin by the product of the relative momentum deviation (Gp/p) and the dispersion function (D,Dc) in the machine. In normalised coordinates, the equilibrium orbit for an off momentum particle is given by:

p p D X

n

G

EQ.O

p p D X

n

G c c

EQ.O

. (2.43)

where (Dn, Dcn) is the normalised dispersion function. When going round the machine, particles will perform betatron oscillations around the orbit defined by the dispersion

function. For the derivation of a general Hamiltonian including the effect of an off-centred

equilibrium orbit at the sextupole, it is convenient to introduce a second coordinate system (XE, XcE) with its origin on the off-centred orbit. The particle coordinates are then split into a constant term, given by the dispersion function, and a betatron part as shown in Figure 2.8.

(2.44) (2.45) Figure 2.8. Coordinates of the betatron motion.

The general Hamiltonian is derived following exactly the mathematics of the Kobayashi

theory. Particle motion in the unperturbed machine is, as before, described by transfer
matrices. For k revolutions in the machine,

X X X X

k k E E E E

c

c
M

. (2.46)

The main difference appears when calculating the kick of the sextupole on the off-centred equilibrium orbit:

Sextupole in dispersion-free region: ' ' X X SX

E E E

c

2

(2.47) Sextupole in dispersion region: ' ' X X SX S X D p p

n E E E

G c

2

2

(2.48) X X D p p X X D p p

n n

E

E

G G c c c

c

c c X X D p p X X D p p

n n E E

G G

X Xc XE XcE X XE DnGp/p

SLIDE 31

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

29 This yields the more general Hamiltonian:

H
c
c
c

H G

E E E E E E E

2 4 3 3 2

2 2 2 3 2 2

X X S X X X S D p p X X

n

(2.49)

It can be seen that the dispersion-dependent term affects (apart from the shift of origin)

nly the circular part of the trajectories. Reordering (2.49) shows that this can be

considered as a change in tune of the particle.

H
c
c

1 2 3 4 3

2 2 2 3

H G

E E E E E

S D p p X X S X X X

n

(2.50)

The initial tune distance of the particle was given by

G S H Q 1 6 . (2.51)

Assuming that the dispersion-dependent term is of the same form leads to the following expression for an additional change in tune,

G S G ~ Q S D p p

n

1

6 3 . (2.52)

The above expression contains the normalised sextupole strength and dispersion. Rewriting the tune change in terms of unnormalised coordinates with (2.13) gives

p p D k p p D x B B Q

x x

G E S G U E S G c

S

2 z 2 S

4 1 d d 4 1 ~ " " (2.53)

with kc being the normalised sextupole gradient

2 2

1

c

dx B d B k

z

U (2.54)

By introducing the chromaticity Qc as the linear change of the betatron tune with momentum, the chromatic effect of a sextupole can be expressed as

D k p dp Q d Q c

c

S

4 1 / ~ " E S , (2.55)

which corresponds exactly to the well-known linear form for the tune shift introduced by a

sextupole. Although the above derivation is restricted to the region in tune close to the

third-order resonance the result is in fact generally valid.

SLIDE 32

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

30 From the above theory it can be deduced immediately that the resonance sextupole is best positioned in a dispersion-free region, otherwise any change of the sextupole strength will also alter the tune distance of the particles and therefore the extraction phase space. The general Hamiltonian describes the physics of a third-integer resonance to first

rder correctly without any further restrictions on sextupole locations or particle momenta.

Comparing (2.50) with the Kobayashi Hamiltonian (2.31) shows that both expressions are

f the same qualitative shape although they are defined in different coordinate systems.

This means that the considerations about phase-space geometry and dynamics derived in the previous sections can be adopted by simply substituting

H SG S G G

6

6 3 Q Q SD p p

n

(2.56) X X D p p

n

G

(2.57) c o c c X X D p p

n

G (2.58)

in all the relevant expressions. Equation (2.56) represents the chromatic effect of the sextupole, (2.57) and (2.58) are the coordinate transformations from the system (XE, XcE) back to the reference system with its origin at the sextupole centre. The general expressions for the separatrices are given below. Figure 2.9 shows the phase-space geometry at the sextupole for positive values of h, Dn, Dcn and Gp/p.

(2.59) (2.60) (2.61) Figure 2.9. Geometry of the separatrices at the sextupole for an off-momentum beam.

Finally the change of the phase-space map when the sextupole is located in a non-zero- dispersion region can be summarised: x All phase-space trajectories are shifted by the offset of the equilibrium orbit. x The size of the stable triangle is changed due to the tune change caused by the chromatic action of the sextupole.

X’ X A B C h equilibrium orbit shifted by dispersion

'X

D p p

n G /

' c

c X D p p

n G /

A: X

D p p h

n

G

B:

c c

3

2 1 2 X D p p X D p p h

n n

G G

C:

3 2 1 2 c c

X

D p p X D p p h

n n

G G

SLIDE 33

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

31 Figure 2.10 shows phase-space maps for an off-momentum particle (Gp/p > 0) with a modified tune distance H0 in the unperturbed machine. To illustrate the effect of the dispersion, the sextupole position is changed from a zero dispersion region to a region with finite dispersion while keeping the normalised sextupole strength S0 constant.

Figure 2.10. Change of phase-space map for off-momentum beam due to dispersion at the sextupole.

2.5.1 Closed-orbit distortions In a real machine, it has to be expected that there is a distortion of the equilibrium orbit due to magnet imperfections and misalignments. This will give an effect similar to the one caused by a finite dispersion function. The fundamental difference is that the closed orbit distortion will be to first order independent of the particle momentum. It is straight forward to also include orbit distortions at the sextupole in the Hamiltonian. The coordinates of the particle at the sextupole can be split into three parts: x Contribution of the dispersion function x Closed orbit distortion at the sextupole x Betatron oscillation of the particle The Hamiltonian is derived following exactly the same scheme as above:

3

2 2 2 CO.

3 4 3 2 1

E E E E E

G H X X X S X X X p p D S

n

c
c
'
H

. (2.62)

The geometry and dynamics of the phase-space can be analysed as in the previous section. The shift of phase-space trajectories consists of the momentum-dependent dispersion part and the constant offset given by the closed orbit distortion. The chromatic term shows the same dependency. The total change in tune will be:

CO.

S S x CO. Disp.

4 1 4 1 ~ ~ ~ x k p p D k Q Q Q

x

' c

c
"

" E S G E S G G G (2.63)

1 0.5 0.5 1 1 0.5 0.5 1

norm. phase space H/S > 0 Dn = Dn’ = 0

X’ X

1 0.5 0.5 1 1 0.5 0.5 1

norm. phase space H0/S0 > 0 D n = D’n < 0

X’ X

SLIDE 34

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

32 2.6 Phase-space maps along the machine To draw the phase-space map at any position s in the machine, it is sufficient to describe the evolution of the separatrices around the machine. All the relevant physics of the slow extraction process can be obtained from this. The following considerations are based on the general Hamiltonian (2.50), derived earlier. It is assumed that there is only one sextupole , defining the reference position (betatron-phase P = 0) in the machine. 2.6.1 General description of separatrices The description of particle motion in the linear machine (outside the sextupole) is done by using the coordinate systems introduced in Figure 2.8. There are two contributions to the evolution of the separatrices and the stable triangle when tracking around the machine. x The momentum-dependent equilibrium orbit describes the centre of the stable triangle at any position s in the machine

X s D s p p

n EQ.O

G c c X s D s p p

n EQ.O

G , (2.64)

x The phase advance 'P from the sextupole to a position s determines the orientation of the stable triangle, in fact the separatrices are simply rotated clockwise around the equilibrium orbit by 'P. This behaviour is implied by the transfer matrix, describing the betatron motion in normalised phase space, which has the form of a rotation matrix,

P

P

P P P P

' '

'

'

'

' R M cos sin sin cos (2.65)

The size of the stable triangle is determined absolutely by the normalised sextupole strength S and the modified tune distance H (including chromatic effects) and remains unchanged in normalised phase space. A schematic example is given in Figure 2.11.

Figure 2.11. Normalised phase-space maps separated by 'P = 90°.

This is not a restriction on the generality of the theory at the level of approximation being used. It will be shown later that from the

standpoint of resonance excitation (described on a time scale of 3 turns) many sextupoles can be replaced by a single virtual sextupole.

phase advance 'P = 90° clockwise rotation by 90° Normalised phase space at sextupole (a) Xc X A B h C Dn,SextGp/p Xc X Normalised phase space at element with 'P = 90° (b) C A B h Dn(s)Gp/p

SLIDE 35

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

33 From the considerations above, it follows that a general description of the separatrices has to contain the dispersion function Dn(s) and the phase advance from the sextupole 'P(s) as

parameters. However, the earlier derived expressions (2.59)-(2.61), describing the

separatrices at the sextupole, contain only the dispersion function. The constant dispersion function at the sextupole can be simply replaced by Dn(s). The dependency on the betatron phase advance is implemented via geometrical considerations. A standard form for the description of a straight line is shown in Figure 2.12.

(2.66) Figure 2.12. Perpendicular form of a straight line.

Rewriting the expressions of the separatrices in the form of (2.66) gives exactly the required angular dependency. With h being the distance from the side of the stable triangle to its centre, the equations of the separatrix will have the form:

h p p s D X p p s D X

n n

c
c
D

G D G sin cos . (2.67)

By inspection of Figure 2.11.(a), the values for D at the sextupole are: (A) D = 180° anticlockwise (B) D = 300° anticlockwise (C) D = 420° anticlockwise The separatrices rotate clockwise with the betatron phase advance measured from the

sextupole. Thus, the effective angle at an element with a phase advance 'P from the

sextupole will be (D'P). For convenience, it is easier to think in terms of 'P appearing as a positive term, which gives:

(A)

c

c

X

D s p p X D s p p h

n n

' ' ' ' cos sin P P (2.68) (B)

q

c c

q

X D s p p X D s p p h

n n

' ' ' ' cos sin P P 120 120 (2.69) (C)

q

c c

q

X D s p p X D s p p h

n n

' ' ' ' cos sin P P 120 120 (2.70)

Equations (2.68)-(2.70) describe the separatrices and therefore the stable triangle at any position s with a phase advance 'P from the sextupole. It was mentioned earlier that all the separatrices are essentially equal, which is

bvious now by inserting the phase advance for one turn, 'Pturn # 2Sn/3 = ±120° mod

360°, in the above equations. This means that a particle will appear to jump from one separatrix to the next on consecutive turns.

D h x y

x y h cos sin D D

The angle D is measured anticlockwise

from the x-axis to the perpendicular h.

SLIDE 36

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

34 2.7 Resonance excitation by sextupoles The theory derived so far shows how a single sextupole magnet excites the third-integer

resonance. In general, in a real machine there will be several sextupole magnets (e.g. for

the correction of chromaticity). Apart from these magnets there will also be some unwanted sextupole field components introduced mainly by the end fields of the main dipole magnets. The common effect of all these sextupole fields on the resonance is described by the so-called driving term N [11@ For the third-order resonance 3Qh = n, the driving term is,

^ ` s

x B B C

x z C x

d i 3 exp d d 1 24 1

2 2 2 / 3

P U E S N

.

(2.71)

The integral is made around the full machine circumference and includes all sextupole

fields. For short sextupoles, the above can be rewritten as a sum, using the normalised

sextupole strength S, as defined earlier,

N

S P

1

12 3 C Sn

x,n n

i exp (2.72)

An equivalent sextupole can be found by evaluating the above driving term and equating to a single ‘virtual’ sextupole.

S

S

n virt x,virt n x,n

i exp i exp 3 3 P P

.

(2.73)

By separating real and imaginary parts,

S

S

virt x,virt n n x,n

cos cos 3 3 P P

S

S

virt x,virt n n x,n

sin sin 3 3 P P

,

(2.74)

the betatron phase of the equivalent sextupole is found from

tan

sin cos 3 3 3 P P P

x,virt n n x,n n x,n

S

S

n

, (2.75)

and the equivalent sextupole strength is given by

S

S S

virt n x,n n n x,n n 2 2 2

3 3

cos

sin P P . (2.76)

The virtual sextupole describes the resonance excitation of all sextupole fields in the

machine. The stable triangle will have its characteristic shape with one upright separatrix

at the position with the betatron phase Pvirt of the virtual sextupole. The orientation of the separatrices at an element with a betatron phase P can now be evaluated by inserting 'P = (P - Pvirt) into equations (2.68)-(2.70).

SLIDE 37

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

35 However, the size of the stable triangle depends on two parameters: x The strength of the sextupoles driving the resonance which is represented by Svirt. x The modified tune distance H of the particles given by the initial value H0 in the linear machine (without sextupoles) and the chromatic effects of all the sextupoles. For the determination of the size of the stable triangle it is therefore necessary to evaluate, separately from the driving term, the chromatic effect of all the sextupoles in dispersion regions which is done in Section 2.8 below. 2.7.1 Combination of sextupoles In a series of an even number of sextupoles with identical normalised strength, S0, it can be seen from (2.76) that a separation in phase of 'P = S/3 leads to cancellation of the resonance driving term,

^

`

^

`

2

x,1 x,1 2 x,1 x,1 2 Virtual

... 60 3 sin 3 sin ... 60 3 cos 3 cos

q
q
P

P P P S S S S S . (2.77)

Similarly, a phase separation of 'P = 2S/3 leads to a reinforcement of the driving term,

^

`

^

`

2

x,1 x,1 2 x,1 x,1 2 Virtual

... 120 3 sin 3 sin ... 120 3 cos 3 cos

q
q
P

P P P S S S S S . (2.78)

2.8 Chromatic effects of sextupoles The chromatic effect of a thin-lens sextupole was derived earlier (2.53), (2.55). In general the effect of sextupoles (in a region of finite dispersion) on the horizontal and vertical chromaticities of the linear machine,

,

/ d d

c

p p Q Q

x x ,

/ d d

c

p p Q Q

z z

(2.79)

is given by,

c

c

'

C x x x

s s D s s k Q d 4 1 E S

c

c '

C x z z

s s D s s k Q d 4 1 E S . (2.80)

For short sextupoles, using the thin-lens approximation, the above can be rewritten as

c
c

'

i i x x i i x

4 1 D k Q E S "

c

c '

i i x z i i x

4 1 D k Q E S " . (2.81)

Expressions (2.80) and (2.81) describe the chromatic effect of all sextupoles in the

machine. The modified tune distance H of a particle with a momentum deviation Gp/p is

now obtained as

p p Qx G S H H c '

6

. (2.82)

SLIDE 38

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

36 Finally, the concepts for the description of the phase space of a third-integer resonance, derived earlier for a single sextupole magnet, can now be generalised for any sextupole distribution in the machine by using: x Expression (2.82) for the modified tune distance of the particles. x The virtual sextupole strength (2.76) as normalised sextupole strength. x The phase of the virtual sextupole (2.75) as the reference position. 2.8.1 Schemes for chromaticity correction Chromatic effects in accelerators are caused by the momentum dependency of the focal properties of lattice elements. Sextupole magnets, placed in dispersion regions, are used to correct such effects and to adjust the chromaticity of the machine, e.g. to ensure the transverse stability of the beam [12], or to meet the Hardt-condition [13] for resonant extraction (see Section 3.8). For an independent adjustment of the horizontal and vertical chromaticities at least two magnets are needed. It can be seen from (2.80) that 'Qcx and 'Qcz are of opposite signs and therefore focusing (kcSF) and defocusing (kcSD) sextupoles, as defined in Section 2.1, are used. In general, chromaticity correction is done not with one F- and one D- sextupole, but with groups consisting of several magnets, each group powered by a single power converter. Such sets of magnets with identical fields form so called magnet

families. A simple scheme for chromaticity control is given below.

Consider two sextupole families, one focusing with strength kcSF, the other one defocusing with kcSD. Assuming that all the magnets are of the same length and using the thin-lens approximation, gives two expressions for the change of the chromaticities,

' c

c c

Q

k D k D

n n x S SF x x n SD x x n ND NF

" 4

1 1

S E E (2.83)

' c

c c

Q

k D k D

n n z S SF z x n SD z x n ND NF

" 4

1 1

S E E , (2.84)

with NF and ND being the numbers of magnets within the focusing and defocusing

families. To get some degree of orthogonality for the correction, the positions of the

F-sextupoles should be in ‘large Ex -small Ez’ regions and the D-sextupoles should be in ‘small Ex -large Ez’ regions. The unwanted effect of a chromaticity correction is that in general the resonance driving term will also be changed. It should be noted that spacing the chromaticity sextupoles at intervals of 600 in betatron phase as shown in (2.77), is not sufficient for a cancellation of the resonance driving term. It can be seen that (2.83) and (2.84) contain the normalised sextupole gradient kc, but for a cancellation the normalised sextupole strengths have to be identical. This means it is also essential that the horizontal beta functions are identical at all sextupole locations.

SLIDE 39

Theory of the third-integer resonance Chapter 2

________________________________________________________________________________________________________________________________________________________________________________________________________

37 2.9 Sextupole families for small machines Equation (2.71) for the driving term of the third-integer resonance shows that in principle all sextupole magnets in a machine contribute to the resonance excitation. Thus, any correction of the chromaticity affects not only the modified tune distance H for an off- momentum particle, but at the same time changes the strength of the resonance excitation. Ideally one would like to have completely independent handles for the control of the resonance strength and the horizontal and vertical chromaticities, especially in a machine using a third-integer resonance as these parameters determine the extraction set-up. In small medical machines, it is not always possible to have several sextupole families, mainly for the lack of space and for the cost of the additional power supplies. As mentioned above there are three parameters to control, x the resonance excitation represented by Svirt, x the horizontal chromaticity Qcx, x the vertical chromaticity Qcz. Therefore three independent power supplies will be needed and the minimum number of sextupole magnets will be three. From (2.80) it can be seen that sextupoles in dispersion-free regions do not change the chromaticity, for this reason sextupoles that excite or control the resonance are best placed in a zero dispersion region. Unfortunately, for the chromaticity sextupoles no similar situation exists, but by taking advantage of the rather special third-integer tune, an approximate solution that leaves the resonance excitation unchanged, can be found when the tune is near an even number (i.e. Qh | nr1/3, with n even). The scheme is based on a lattice with a periodicity of two, which means the ring can be split in two halves with identical sequences of magnetic elements. Due to the symmetry, the dispersion and beta functions on opposite sides of the ring will be equal. Consider now two sextupoles placed on opposite sides of the ring. If the betatron phase at the first sextupole is taken as the origin, the second sextupole will then be at Px = QS | (nr1/3)S The effect on the resonance and the chromaticities is now evaluated with (2.76), (2.83) and (2.84). Resonance excitation:

S

S S Q S S Q

virt 1 2 1 2 2 2 2

3 3

cos

cos sin sin S S (2.85) S S S

virt 1 2

r (2.86)

Since the tune is very close to (nr1/3)S the sine terms are zero and the cosine terms are

unity. The upper sign is for n odd and the lower sign for n even.

Chromaticity:

2

1 x z S z 2 1 x x S x

4 ; 4 k k D Q k k D Q c

c

c ' c

c
c

' E S E S " " (2.87)

SLIDE 40

Chapter 2 Theory of the third-integer resonance

________________________________________________________________________________________________________________________________________________________________________________________________________

38 Because of the equal lattice functions at the sextupoles the relation between normalised and unnormalised sextupole strength (2.13) can be written as c

k S

1 1

c

and c

k S

2 2

c

, with c = constant. Now (2.87) can be rewritten as

2

1 x z S z 2 1 x x S x

c 4 ; c 4 S S D Q S S D Q

c

'

c

' E S E S " " (2.88)

From (2.86) and (2.88) it can be seen that for n even, the excitation of the resonance is proportional to (S1-S2) whereas (S1+S2) changes the chromaticity. Consider now a family

f two sextupoles in series on one power converter. When connected with the same

polarity (i.e. S1 = S2) the sextupoles change only the chromaticities without affecting the resonance, connecting them with opposite polarity (i.e. S1 = -S2) has the converse effect. For n odd, there is no such fortuitous situation.

SLIDE 41

39

Chapter 3 RESONANT SLOW EXTRACTION

The use of a resonant slow extraction scheme allows the continuous extraction of particles from a synchrotron during a large number of turns. Slow extraction is a dynamic process that can be described by combining the different phase-space properties that characterise particle behaviour in a synchrotron with a horizontal betatron tune close to a third-integer, as derived in the previous chapter. For clarity and simplicity it is assumed that there is

nly one sextupole for the resonance excitation, located in a zero-dispersion region. The

chromaticities of the machine and the tunes are therefore independent of the resonance

strength. It is also assumed that the resonance is positioned in the machine centre, i.e. the

resonance tune (GQ = 0) corresponds to the central orbit momentum (Gp = 0). 3.1 Set-up of the extraction process First, the linear machine is considered. During injection and acceleration of particles, the presence of resonances has to be carefully avoided, as this leads to immediate particle loss. This is done by adjusting the tune values such that they are sufficiently far from dangerous resonances and by compensating for non-linear field errors that could drive resonances. At the end of the acceleration process the so-called flat top is reached (i.e. the fields in the main magnets are kept constant). The resonance sextupole is still switched off and the particle trajectories are circles in normalised phase space. The radius of each circle is determined by the particle position and divergence, (X0, Xc0), and is referred to as normalised amplitude A, with

A X X

c

2 2 .

(3.1)

The circle area is the so-called single particle emittance,

S

S

2 2 2

X X A E c

.

(3.2)

Before starting the extraction process, the circulating beam has to be positioned with respect to the resonance. The horizontal tune is moved closer to the resonance by adjusting the currents in the quadrupole magnets and the resonance sextupole is slowly raised to the nominal strength S in a quasi-adiabatic way. As shown earlier, the circular trajectories are deformed and the phase space is split into stable and unstable regions by the separatrices that define the stable triangle. A particle with a certain tune distance GQ remains stable if its single particle emittance in the unperturbed machine was not larger than the area of the corresponding stable triangle. Figure 3.1 shows the trajectory of a stable particle at the position of the sextupole in normalised phase space. The particle starts on a circle and ends on an almost exactly triangular shaped trajectory but remains stable.

SLIDE 42

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

40

Figure 3.1. Particle trajectory during adiabatic raise of sextupole, normalised phase space.

3.2 Resonance stopband The resonance sextupole strength is one of the main parameters for the slow extraction. For a nominal value of S, the size of the stable triangle for a certain tune distance GQ, is

fixed. The stability of a particle with GQ depends on its initial amplitude in the

unperturbed machine. A criterion is found by comparing (3.2) and expression (2.39) for the size of the stable triangle, as

S

G S S

2 2 2

3 48 Q S A E

stable stable

d . (3.3)

A particle will be stable if its emittance is smaller than or equal to the area of the stable

triangle. Therefore the area of the stable triangle is often called stable emittance.

Inverting the inequality in (3.3) and solving with respect to GQ gives the tune interval for which particles with a certain amplitude will be unstable. This interval is referred to as the stopband of the third-integer resonance,

S A Q S G 3 48 1

.

(3.4)

Equation (3.4) can be rewritten in terms of absolute tune with GQ

Q Q

particle resonance

as

A S Q Q A S Q

res part res

S S 3 48 1 3 48 1

.

(3.5)

The above expressions are best illustrated graphically in the Steinbach diagram as shown in Figure 3.2.

SLIDE 43

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

41

Figure 3.2. Steinbach diagram for third order resonance stopband, amplitude-tune space.

From (3.4) it can be seen that for every tune distance GQ, there is a maximum stable amplitude A(GQ). Particles will be unstable if their amplitude is bigger than A(GQ) and stable if their amplitude is smaller. There are no stable amplitudes for particles on resonance, i.e. all particles with GQ = 0 are unstable. The range of stable amplitudes increases when the tune is moved away from the resonance. Replacing the inequality in (3.4) by an equality yields the equation of the boundary between stable and unstable regions in Figure 3.2, sometimes referred to as resonance line,

S Q A G S 3 48 . (3.6)

Particles that exactly fulfil (3.6) are on the limit of stability, which means they are moving

n a stable triangle in phase space.

3.3 Overview of extraction methods The width of the third-integer stopband is symmetric with respect to the resonance, therefore there is no loss in generality when considering only particles with GQ < 0. It is assumed that the particles are stable after the machine has been set to extraction

conditions. In order to extract particles from the stable region in Figure 3.2, they either

have to be moved into the unstable region or the size of the unstable region has to be

increased. Inspection of (3.6) shows that this can be done by:

x Widening the stopband by increasing the resonance sextupole strength S. x Moving the particles into the stopband by changing their tune values. x Increasing the particle amplitudes until the critical value A(GQ ) is reached. The different methods are explained in more detail below. Extraction by changing the resonance strength For the description of the extraction process two particles with a tune distance GQ0 and amplitudes A1 and A2 are considered. Figure 3.3 (a) shows the situation before the start of

extraction. Both particles are stable, their amplitudes are smaller than the critical

tune, Q Amplitude A , emittance, (

Qres

unstable region stable region stable region A0 stopband for particles with amplitude A0

Qres-|GQ(A0)| Qres+|GQ(A0)|

resonance line

SLIDE 44

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

42 amplitude A(S0) given by (3.6) and their phase space trajectories are deformed circles inside the stable triangle defined by S0 and GQ0.

Figure 3.3. Extraction via resonance sextupole change.

When increasing the sextupole strength, the unstable region widens, the critical amplitude A(S) is decreasing which also means, the size of the corresponding stable triangle

decreases. Figure 3.3 (b) shows the extraction process at the moment when the critical

amplitude has reached the value of the amplitude A1. This means the particle is moving on a stable triangle and will start going outwards on one of the separatrices, finally being

extracted. Basically, the same applies to the particle with the smaller amplitude A2, but

later on, when the sextupole strength has been further increased. The main difference is that this particle will be extracted from a smaller stable triangle. Inspection of equation (3.6) shows that zero-amplitude particles, with a finite tune distance, can never be extracted as this would require an infinite sextupole strength. Extraction by changing the tune distance The substantial difference from the above method is that, instead of increasing the resonance strength, the tune of the particles is slowly moved towards resonance tune. In Figure 3.4 (a) both particles are stable, their tune distance is larger than the width of the resonance stopband.

Figure 3.4. Extraction via tune change.

It can be seen that every amplitude is linked to a critical tune value at which the particle becomes unstable. The smaller the particle’s amplitude, the further away it’s tune from the corresponding critical tune. Therefore, when the tune is moved towards the resonance, the particle with the larger amplitude will be extracted earlier, as indicated in Figure 3.4 (b). As in the extraction method described above, particles that reach the limit of stability are moving on stable triangles. The size of these triangles is again determined by the initial particle amplitudes.

tune Q A(S0) A1 A , E

Qres Qres-GQ0 GQ0

A2 tune Q A1 = A(S1) A , E

Qres Qres-GQ0 GQ0

A2

(a) (b)

tune Q A(GQ1) A1 A , E

Qres Qres-GQ1 GQ1

A2

(b)

tune Q A(GQ0) A1 A , E

Qres Qres-GQ0 GQ0

A2

(a)

SLIDE 45

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

43 Extraction by increasing the particle amplitude The third possibility for extraction is to increase the transverse amplitudes of the particles. The main difference compared to the extraction techniques described earlier is that particles with the same tune distance will be extracted from stable triangles with equal size, independent of the initial particle amplitudes. 3.4 Spiral step After a particle has become unstable, it starts moving on the outgoing part of one of the

separatrices. In Section 2.6, it was shown that all three separatrices are essentially equal,

they simply exchange their position from turn to turn. However, when analysing the extraction process, it is more convenient to consider the separatrices as geometrically fixed in phase space. Then, an unstable particle moves out along a separatrix, one turn later it is

n the next and after every third turn it returns to its initial separatrix. Earlier (2.40), the

change of the particle coordinates within three revolutions in the machine was derived as

2

2 3 3

4 3 2 3 X X S X X X SX X X c

c

' c

c

' H H (3.7)

By inserting (3.7) into (3.1) the amplitude increase during three turns, for a particle that is just on resonance (i.e. H = 6SGQ = 0) is obtained as

2 3

4 3 A S A A A

'

. (3.8)

From (3.8) it can be seen that the growth increases rapidly as the particle progresses along the unstable section of the separatrix. After a certain number of turns, the particle amplitude has increased so much, that the particle ‘jumps’ into the electrostatic septum and is extracted. Depending on the phase advance from the sextupole, 'PES, the electrostatic septum is positioned in the outer or inner side of the vacuum chamber. The phase advance has to be chosen such that only one separatrix is intercepted by the electrostatic septum. Figure 3.5 gives a schematic view of the normalised phase space at the resonance sextupole and the electrostatic septum for a phase advance of 'PES = 225°.

Figure 3.5. Amplitude increase during last three turns before extraction.

X X’ extraction separatrix X X’ extraction separatrix phase advance 'PES = 225° clockwise rotation by 225° electrostatic septum 1 2 3 3 2 1 T XES 'R 'Rc

SLIDE 46

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

44 The maximum possible amplitude, Alast, of a particle, which is not extracted, can be calculated from the normalised position of the electrostatic septum, XES, and the angle T between the extraction separatrix and the horizontal axis as

ES last

X A ) cos( 1

.

(3.9)

Three turns later the particle will be extracted and by inserting (3.9) into (3.8), the increase

f amplitude during these last three turns in the machine is found to be

2 ES 2 last

cos 4 3 4 3

'
X

S A S A . (3.10)

From (3.10) it can be calculated how far the particle jumps into the septum, the so-called spiral step 'R, and the spiral pitch 'Rc (see Figure 3.5),

2 ES last

cos 1 4 3 cos X S A R

'

' , (3.11)

and

2 ES last

cos tan 4 3 sin X S A R

'

c ' . (3.12)

The spiral step and pitch define the size of the extracted beam. It is important to note that the slow-extracted beam differs markedly from the circulating beam. In the plane

f extraction, the phase-space footprint is a narrow rectangle that corresponds to the

segment of the outgoing separatrix that is cut and deflected by the electrostatic septum and the emittance is visibly smaller then the circulating beam’s one. In the orthogonal plane, the phase space distribution and the emittance are to first order the same as those of the circulating beam [14]. This asymmetry of the extracted beam requires special techniques for the design of the beam delivery system towards the patient, discussed in Chapter 7. It should be noted that (3.11) and (3.12) are only valid for particles which are exactly on resonance (GQ = H = 0). The spiral step and pitch for particles with a finite tune distance are in general different and can be evaluated by transforming the coordinates of a particle, that just misses the electrostatic septum, to the position of the resonance sextupole and then applying (3.7). Equation (3.11) gives the maximum distance that a particle can jump into the

septum. The space on the extraction separatrix from the position of the electrostatic

septum to this point will be filled by particles with different starting conditions, some of the particles will hit the septum directly and will be lost. Therefore the extraction efficiency depends critically on the thickness of the septum. This is the main reason for the use of an electrostatic septum, which is in generally built as a wire septum with a wire thickness of the order of 0.1 mm. Electrical fields up to 100 kV/cm can be obtained in such septa. The drawback is that these fields are often not sufficient to kick the beam directly out of the

machine. The small deflection provided by the electrostatic septum translates into a

physical gap between the circulating beam and the extracted particles further downstream, where one or more magnetic septa are positioned to finally extract the beam. Ideally the extraction channel has to be designed such that the only particle loss occurs at the electrostatic septum, which is unavoidable due to its finite thickness. A detailed discussion

f the layout of the extraction channel is given in Sections 3.9 and 3.10.

SLIDE 47

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

45 3.5 Effect of the electrostatic septum The electrostatic septum cuts off particles from the separatrices. Particles that have jumped into the septum are deflected by the electrostatic field, corresponding to an angle M. This kick transforms into a physical gap between the extracted and the circulating part of the separatrix further downstream, where the magnetic septum is positioned.

Figure 3.6. Kick of electrostatic septum in normalised phase space.

The gap size is calculated by comparing the movement of two on-momentum (Gp = 0) particles from the electrostatic to the magnetic septum; particle A starts just inside the electrostatic septum, B just outside, as shown in Figure 3.6. The thickness of the electrostatic septum (typically 0.1 mm) is neglected, both particles are assumed to start from the radial position of the septum xES, and with the same angle xcES but only particle A receives the kick M. With the 2x2 Twiss-matrix (real phase space), the positions and angles of the particles at the magnetic septum are obtained as

A: M

c
12

12 11

m x m x m x

ES ES MS

M

c
c

22 22 21

m x m x m x

ES ES MS

and (3.13) B:

ES ES MS

x m x m x c

12

11 ES ES MS

x m x m x c

c

22 21

, where the suffices ES and MS denote the positions of the electrostatic and magnetic septa. Thus, the effect of the kick appears at the magnetic septum as a difference in position and angle of the particles,

'x m

MS

12 M

(3.14)

and

M

c

'

22

m xMS , (3.15)

where 'xMS is the gap, available for the thicker magnetic septum, which is explicitly given by

'xMS

M

E E P

ES MS sin

. (3.16)

From (3.16) it follows that to relax on the kick (voltage), that has to be provided by the electrostatic septum to create the space for the magnetic septum, the lattice functions of the machine and the positions for the septa have to be chosen such that: x The phase advance between the septa is close to 90° + n180q. x The beta functions at the septa have reasonably large values.

phase advance P Electrostatic Septum X X’ B kick M A B A gap Magnetic Septum X X’

SLIDE 48

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

46 3.6 Separatrix geometry at electrostatic septum So far, in all the figures of the normalised phase space at the electrostatic septum, the extraction separatrix has been shown in the first quadrant, where X and Xc > 0, with a positive kick, M > 0. In fact, only the first and third quadrants are useable. This follows from purely geometrical considerations as shown in Figure 3.7. For the first and fourth quadrant X is always positive, which means that the kick of the electrostatic septum must also be positive, otherwise the extracted particles are driven back into the septum wires. For the first quadrant, a positive kick opens a gap for the magnetic septum further downstream, whereas it is the contrary in the fourth quadrant, as a positive kick drives the extracted beam segment closer to the machine centre. The same considerations apply when comparing second and third quadrant operation.

Figure 3.7. Possible working quadrants for the electrostatic septum.

A further limitation for the separatrix geometry at the electrostatic septum is introduced by the angular spacing of 120° between adjacent separatrices. From equation (3.16) it was found that, for an optimum usage of the electrostatic septum’s kick, the phase advance to the magnetic septum should be 90° + n·180°. Figure 3.8 shows a logical layout for first quadrant operation with the extraction separatrix at 45° at the electrostatic septum and a phase advance of 90° to the magnetic septum.

Figure 3.8. Ideal separatrix geometry at electrostatic and magnetic septa for first quadrant operation.

Kick is positive at a negative position and drives beam back to centre and wires

X’ X

NO O.K. NO

Kick is negative at a positive position and drives beam back to centre and into wires Kick is positive at a positive position and creates a gap for the magnetic septum Kick is negative at a negative position and creates a gap for the magnetic septum

O.K.

Xc X MS ES X Xc kick phase advance q 45°

45°

safety margin safety margin

I IV II III

SLIDE 49

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

47 The extraction separatrix could be rotated by a maximum of 15° anticlockwise before the safety margin in Figure 3.8 is lost and the preceding separatrix hits the electrostatic septum as shown in Figure 3.9 (a), or by 15° clockwise before the following separatrix hits the magnetic septum, Figure 3.9 (b). This limits the possible angles for the extraction separatrix at the electrostatic septum to 45° ± 15° for first quadrant operation and to 225° ± 15° for third quadrant operation.

Figure 3.9. Limitations of separatrix geometry by electrostatic and magnetic septa.

The alternatives to the extraction layout with a phase advance of 90° between the septa, as shown above, are limited. x The first is to accept a much smaller phase advance, which might be imposed by lack of space, but this requires a stronger electrostatic septum which may have implications for the reliability. x The second is to use 270° phase advance, but this has two unfavourable aspects. Firstly, the extracted beam has to be transported for a longer distance in the machine, which for a small synchrotron means that crossing of non-linear magnetic elements (e.g. chromaticity or resonance sextupoles) is almost unavoidable. Any change in these elements would result in a change of the extraction geometry. Secondly, a phase advance of 270° (more generally 270° + n·360°) means that the electrostatic and magnetic septa are on opposite sides of the vacuum chamber. This may have the drawback that aperture is lost for the circulating beam by ‘encasing’ it between the septa.

Xc X X Xc phase advance q

30°

60°

no safety margin

MS ES kick (a) Limited at electrostatic septum. Xc X MS ES X Xc kick phase advance q 30°

60°

no safety margin

(b) Limited at magnetic septum.

SLIDE 50

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

48 3.7 Phase-space representation of beam and resonance Up to now only single particles have been considered, but the beam in an accelerator always has a certain momentum spread Gp/p and, due to the chromaticity in the machine, this momentum spread translates into a tune spread according to

p p Q Q G G c . (3.17)

The size of the beam in an accelerator is characterised by the emittance, which is defined as an area in phase space that contains a certain percentage of the beam particles. In general, a beam will contain particles with amplitudes between zero and a maximum amplitude, corresponding to a maximum single particle emittance, which is numerically equal to the total emittance of the beam. At a given position s in the machine, the beam can be represented in phase space by a series of ellipses (circles in normalised phase space), centred around the dispersion vector D(s)·Gp/p, as shown in Figure 3.10.

Figure 3.10. Representation of a beam in phase space and normalised phase space.

The circles that represent the beam emittance in normalised phase space become triangles under the influence of the resonance. The resonance for particles of each momentum is represented by a triangle, corresponding to the last stable orbit, and the extensions along the outward separatrices. According to (3.17) each momentum corresponds to a certain tune distance from the resonance and the corresponding triangles are of different sizes as shown in Figure 3.11.

Figure 3.11. Normalised phase space representation of a beam under the influence of the resonance.

(Gp/p)2 (Gp/p)1 (Gp/p)0 Real phase space Areas are emittance of betatron motion for each momentum. x xc Xc X Normalised phase space Areas in normalised phase space are still equal to the emittance. Normalised phase space Size of the stable triangles depends on the tune distance from the resonance. Orientation of the triangles depends on phase advance from resonance sextupole. X’ X (Gp/p)2 (Gp/p)1 (Gp/p)res

SLIDE 51

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

49 When considering the interactions of the beam and the resonance, it is very convenient to use the Steinbach diagram that represents the stopband of the third-integer resonance. With (3.17) the expression for the stopband (3.6) can be rewritten with the momentum spread and the chromaticity as

p p S Q A G S

c

3 48 . (3.18)

It is sometimes more convenient to use the momentum spread rather than the tune as abscissa, since the momentum spread is an independent beam parameter (changing the chromaticity will change the tune but not the momentum). Figure 3.12 (a) shows the circulating beam before extraction, corresponding to the phase space representation in Figure 3.10. Figure 3.12 (b) shows the beam during the extraction process, corresponding to the phase-space diagram 3.11.

(a)

(b)

Figure 3.12. Steinbach diagram of the beam before and during extraction.

It can be seen that during the extraction a quasi-static situation is reached and particles with all amplitudes in the beam enter the resonance. These particles have different momenta and, due to the chromaticity, different tunes, which means the corresponding stable triangles are of different sizes. Figure 3.13 shows the separatrix geometry at the resonance sextupole and the electrostatic septum in normalised phase space.

Figure 3.13. Separatrices at the resonance sextupole and the electrostatic septum.

The dotted lines are the separatrices corresponding to zero amplitude particles P0 in Figure 3.12 (b), the full lines are the separatrices of particles with the maximum amplitude, P2.

'p/p A , E

Gp/p) res Gp/p) stack

stack A max 'p/p A , E

Gp/p) res=Gp/p)0

stack A max P2 P0

band of unstable particles

X X’ Extraction separatrices X X’ Extraction separatrices phase advance 'PES = 225° clockwise rotation by 225° Electrostatic septum M XES

Gp/p)2 Gp/p)inst

SLIDE 52

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

50 The separatrices of all other particles are found in between these two extreme cases. The instantaneous momentum spread of the extracted beam can be calculated from (3.18) with the maximum amplitude Amax as

Q S A p p c S G 3 48 1

max inst.

. (3.19)

3.8 Hardt condition During the extraction process, there will in general be particles with all amplitudes becoming unstable at any given instant and being extracted along separatrices from stable triangles with different sizes. As shown in Figure 3.13, the separatrices at the electrostatic septum are in general not superimposed, which means that particles, moving outwards along different separatrices, reach the septum with different angles. Therefore, not only particles that hit the septum cross-section directly will be lost, but due to the finite septum length, also particles that are moving close to the septum wall with a certain divergence.

(a) (b) Figure 3.14. Particle losses at the electrostatic septum due to angular spread of separatrices.

Figure 3.14 (a) shows the trajectories of two particles in real space, one starting just short

f the electrostatic septum at x < xES, and the other just beyond the septum wires at

x > xES. From geometrical considerations it follows that particles starting from x < xES will be lost along the septum wall if their relative divergence, compared to the septum angle, is larger than (xES - x)/LES. A similar situation exists for particles inside the septum with x > xES, except that the electric field curves their orbits and modifies the forbidden

angles. Figure 3.14 (b) shows the so called ‘shadow’ regions (grey shaded) in normalised

phase space. All the particles on the segments of the separatrices that cross the shadow regions are lost along the septum wall. It was mentioned earlier that stable triangles of different sizes correspond to different

momenta. At a position in the machine where the dispersion is non-zero, the triangles will

be shifted according to their momentum (see Section 2.6). With a finite dispersion at the electrostatic septum, this effect can be used to superimpose the extraction separatrices in

rder to avoid particle losses in the shadow regions. It should be noted that in

Figures 3.13 and 3.14 (b) all the triangles are centred which means the dispersion was assumed to be zero at the electrostatic septum. The condition, to arrange the separatrix geometry and to set the optics at the electrostatic septum such that all separatrices are superimposed, is known as the ‘Hardt condition’ [13]. The requirements on the lattice functions can be derived with a

xES x<xES Septum length LES s E0 x particle trajectories XES X X’ XcES x>xES

SLIDE 53

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

51 purely mathematical approach from the general expressions for the separatrices (2.67)- (2.70), given in Section 2.6,

h

p p s D X p p s D X

n n

'

c
c
'
P

D G P D G sin cos . (3.20)

The angle D describes the orientation of the separatrices at the sextupole, 'P is the phase advance from the sextupole to the considered lattice element and h is the distance from the side of the stable triangle to its centre,

h S S Q 2 3 4 H S G . (3.21)

With (3.17) the above expression can be rewritten as

p p Q S S h G S H c 4 3 2 . (3.22)

Inserting (3.22) into (3.20) gives the general expression for the separatrices as a function

f particle momentum and chromaticity of the machine,
p

p Q S p p s D X p p s D X

n n

G S P D G P D G c '

c
c
'
4

sin cos . (3.23)

It can be seen that the size of the stable triangles and the displacement of the separatrices due to the dispersion is determined by the momentum of the corresponding particles. In

rder to superimpose the extraction separatrices the momentum dependency has to be

removed from the equation of the separatrices, giving the Hardt condition,

Q

S D D

n n

c

'
c
'
S

P D P D 4 sin cos . (3.24)

The above equation may appear very flexible, but there are some boundary conditions that restrict the parameters. x The choice of the angle (D-'P) is restrained by the geometry of the extraction. For

ptimised operation in the first or third quadrant, as discussed earlier, it can be shown

that the phase (D-'P) is either 135° for particles with tune values below resonance (GQ < 0), or 315° for particles above resonance (GQ > 0). The different possibilities for the extraction layout are summarised in Figures 3.17 and 3.18. x The sextupole strength cannot be used as a variable since it determines the spiral step and spiral pitch and therefore the horizontal size of the extracted beam. x The horizontal chromaticity Qc is the main variable, but for small, low-energy machines which are working below transition, the chromaticity should be negative to ensure the stability of the coasting beam [12]. x The normalised dispersion function depends on the lattice. The task of the designer is to foresee a position for the electrostatic septum with suitable values of Dn and Dcn.

SLIDE 54

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

52 The left-hand side (LHS) of (3.24) can be seen as a scalar product of two vectors,

'
'
c

) sin( ) cos( = LHS

n n

P D P D D D . (3.25)

where the (cos,sin) vector is the unit vector, perpendicular to the extraction separatrices as was shown in Section 2.6. The value of the LHS is therefore fixed by the extraction geometry and the dispersion vector at the position of the electrostatic septum. In order to superimpose the separatrices and to fulfil the Hardt condition, the right hand side of (3.24) is then adjusted by changing the horizontal chromaticity of the machine. For a more physical insight, the expressions for the stopband of the resonance (3.18), the momentum spread of the extracted beam (3.19) and the Hardt condition (3.24) have to be considered

together. A change of the chromaticity changes the width of the stopband (the slope of the

resonance line) and therefore the extracted momentum spread. This means that the momenta corresponding to different-size triangles change and the triangles are shifted along the dispersion vector, according to their new momenta. In general, the extraction should be arranged such, that the momentum spread of the extracted beam is small for two main reasons: x The optics of the beam delivery system towards the patient (Chapter 7) is rather

complicated. In a gantry it is likely that the dispersion takes large values which,

together with a large momentum spread, leads to enlarged magnet apertures and therefore bigger and more expensive magnets. x The transfer between electrostatic and magnetic septum is in general chromatic (see Section 3.9). These chromatic effects can be partially corrected with a higher voltage

f the electrostatic septum, but are less important for a small momentum spread.

Inspection of (3.19) shows that the extracted momentum spread will be small for large absolute values of chromaticity and therefore a large absolute right-hand side in (3.24). The LHS (3.25) is zero for perpendicular vectors, which is the case when the normalised dispersion is parallel to the extraction separatrix. For this configuration it is obvious that the separatrices cannot be superimposed for any finite value of chromaticity. (The zero chromaticity extraction is considered in Section 3.11.) The absolute value of the LHS will have a maximum for parallel vectors in (3.25), this means, the normalised dispersion is directed at right angle to the extraction separatrices and the shift of triangles with different momenta, relative to each other, is most effective. The following Figures 3.15 and 3.16 show the separatrix geometry at the electrostatic septum for first quadrant operation. The full lines represent the extraction separatrices corresponding to zero amplitude (zero-size triangle) and maximum amplitude (maximum-size triangle) particles. The dotted lines indicate the position of the maximum- size triangles if the dispersion function was zero. In 3.15 the normalised dispersion vector is almost parallel to the separatrices, thus the absolute value of the LHS is small, whereas in 3.16 the normalised dispersion is perpendicular to the separatrices and the LHS is large. In 3.15 (a) and 3.16 (a) the chromaticity is not adjusted to fulfil the Hardt condition, the RHS of (3.24) is equal for both diagrams. In 3.15 (b) and 3.16 (b) the chromaticity has been adjusted in order to superimpose the separatrices. In 3.15 (b) the dispersion vector is not well suited, the final absolute chromaticity is small and the extracted momentum spread is large, as can be seen from the Steinbach diagram. In 3.16 (b) the normalised dispersion vector is perpendicular to the separatrices that are superimposed for a small momentum spread of the extracted beam and a large absolute chromaticity.

SLIDE 55

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

53

Figure 3.15. Hardt condition for small chromaticity and large momentum spread of extracted beam. Figure 3.16. Hardt condition for large chromaticity and small momentum spread of extracted beam.

'p/p A, E

Gp/p)

stack

band of unstable particles

extracted

momentum spread

Dn Dn(Gp/p) X, Dn shift Dn(Gp/p) Xc Dcn 'p/p A, E

Gp/p)

stack

band of unstable particles

large extracted

momentum spread

Dn Dn(Gp/p) X, Dn shift Dn(Gp/p) Xc Dcn

(a) (b)

Dn Dn(Gp/p) X, Dn shift Dn(Gp/p) Xc Dcn 'p/p A, E

Gp/p)

stack

band of unstable particles

extracted

momentum spread

'p/p A, E

Gp/p)

stack

band of unstable particles

small extracted

momentum spread

Dn Dn(Gp/p) X, Dn shift Dn(Gp/p) Xc Dcn

(a) (b)

SLIDE 56

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

54 Figures 3.17 and 3.18 summarise all the possible extraction layouts for optimised first and third quadrant operation (derived according to Section 2.6). The arrows indicate the required direction of the normalised dispersion vector to fulfil the Hardt condition with the smallest possible momentum spread of the extracted beam. The dotted lines are the extraction separatrices for the zero-amplitude particles which are exactly on resonance (GQ = 0, Gp = 0), the full lines correspond to maximum amplitude particles with a momentum deviation Gp and a tune deviation GQ = QcGp/p.

Figure 3.17. Separatrix geometry for first quadrant operation, ideal normalised dispersion vector. Figure 3.18. Separatrix geometry for third quadrant operation, ideal normalised dispersion vector.

Qc c > 0 G Gp/p > 0 (a) (b) (c) (d) Qc c < 0 G Gp/p < 0 X’ X

Dn Dc

cn

ES

GQ > 0 X’ X

Dn Dc

cn

ES

GQ < 0 X’ X

Dn Dc

cn

ES

GQ < 0 X’ X

Dn Dc

cn

ES

GQ > 0 Qc c > 0 G Gp/p > 0 (a) (b) (c) (d) Qc c < 0 G Gp/p < 0 X’ X

Dn Dc

cn

ES

GQ > 0 X’ X

Dn Dc

cn

ES

GQ < 0 X’ X

Dn Dc

cn

ES

GQ < 0 X’ X

Dn Dc

cn

ES

GQ > 0

SLIDE 57

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

55 It was mentioned earlier that for a small, low-energy machine that works below transition the chromaticity should be negative in order to ensure the transverse stability of the beam [12]. This constraint leaves only four possibilities for the extraction geometry at the electrostatic septum. In Figures 3.17 (b) and (d) the electrostatic septum is on the outside of the vacuum chamber, the dispersion required for the Hardt condition is Dn > 0 and Dcn < 0. In Figures 3.18 (b) and (d) the electrostatic septum is on the inside of the vacuum chamber, the dispersion required for the Hardt condition is Dn < 0 and Dcn > 0. In general, in small machines the dispersion Dn is positive. This fact leaves only the two layouts in Figures 3.17 (b) and (d), for the extraction. Figure 3.19 shows the qualitative shape of the normalised dispersion for a dispersion bump, an arc with regular distributed dipoles and a regular cell structure. The favoured position for the electrostatic septum is on the downward slope of the dispersion (shaded areas), where Dn > 0 and Dcn < 0.

Figure 3.19. Qualitative shape of the normalised dispersion for typical lattice structures.

In Figures 3.17 (b) and (d) the electrostatic septum is in the outer half of the chamber. Earlier, in (3.16) it was shown that the ideal phase separation between electrostatic and magnetic septa is P = 90° + n·180°. An odd value of n (i.e. P = 270°, 630° etc.) means that the septa are on opposite sides of the chamber and, with the electrostatic septum being

n the outside, the magnetic septum has to go to the inside. One feels intuitively that the

extraction should not be done to the inside of the machine because of the high magnetic rigidity of the beam. Therefore, the phase separation should be around 90° to have both septa on the outside of the chamber. A larger phase advance of 90°+ n·360°is less convenient as the extracted beam has to be transported for a longer distance in the machine. The final decision for the extraction layout of Figure 3.17 (d) is based on aperture

considerations. The resonance should be positioned in the centre of the chamber, (i.e. the

resonance tune should correspond to the central orbit momentum) to ensure a balanced growth of the separatrices. During the extraction set-up the beam has to be kept sufficiently far from the resonance. With both septa being on the outside of the vacuum chamber it is natural to position the beam in the inner half of the chamber to avoid aperture

limitations. In this case the beam is below the resonance in momentum and, due to the

negative chromaticity, the resonance is reached from above in tune.

Dn P P P Dn Dn Correct orientation for (Dn,Dn’) in a dispersion bump Correct orientation for (Dn,Dn’) in an arc with spread out dipoles Correct orientation for (Dn,Dn’) in a regular cell lattice

SLIDE 58

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

56 3.9 Transfer between extraction septa Generally, during the extraction process, particles with different momenta are extracted at the same time. When the Hardt Condition is fulfilled, all the extraction separatrices are superimposed at the electrostatic septum and all particles reach the septum on the same separatrix, independent of their momentum. In Section 3.5, the transfer of on-momentum particles between electrostatic and magnetic septa was considered. However, the extracted beam has a certain momentum spread, and in general, the trajectories of particles with a momentum deviation are different to those of on-momentum particles. This may cause a reduction of the space available for the magnetic septum, as calculated in (3.16). It is therefore also necessary to analyse the transfer of particles with the maximum momentum deviation. This is done in a similar way to above, by considering the trajectories of two particles, C and D, starting from the radial position of the electrostatic septum, xES with the septum angle xcES, but only particle C receives the kick M of the

septum. In a linear lattice, the movement of particles with a momentum deviation can be

described with the 3x3 transfer matrix formalism. The horizontal coordinates of the particles at the magnetic septum are derived as

C:

p

p m m x m x m x /

13 12 12 11

G M

c
ES

ES MS

,

p

p m m x m x m x /

23 22 22 21

G M

c
c

ES ES MS

,

and (3.26) D:

p

p m x m x m x /

13 12 11

G

c
ES

ES MS

,

p

p m x m x m x /

23 22 21

G

c
c

ES ES MS

. Comparison of (3.26) and (3.13) shows that the separation between circulating and extracted particles 'xMS is the same as calculated for on-momentum particles, but the gap appears at a different position and angle. The shift in position reduces the effective gap width for the magnetic septum. Figure 3.20 shows the transfer between the septa for particles representing the full momentum spread of the extracted beam. The grey shaded areas represent particles with momenta in between these limits.

Figure 3.20. Transfer from electrostatic to magnetic septum.

phase advance P Electrostatic Septum X X’ B, D kick M A, C B A X X’ C D gapeff Magnetic Septum B A C D gapeff gap for Gp = 0 '"extr 'divextr

SLIDE 59

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

57 The effects of the chromatic terms, m13 and m23 in the transfer matrix, on the phase space geometry at the magnetic septum and the extracted beam are described below. x A non-zero m13 causes a loss of space for the magnetic septum and has to be corrected with a stronger kick of the electrostatic septum. The effective gap width is

gap

m m p p

eff

12

13

M G / . (3.27)

x The horizontal width of the extracted beam is increased, and, to avoid losses inside the magnetic septum, the horizontal aperture has to be enlarged.

'l

m p p

extr

13 G /

(3.28)

x A non-zero m23 is leading to a larger overall divergence of the extracted beam at the magnetic septum and also requires an enlarged horizontal aperture to avoid losses.

'div

m p p

extr

23 G /

(3.29)

It should be mentioned that at the electrostatic septum any angle error of the extraction separatrices will increase particle losses. At the magnetic septum it is usual to foresee a small clearance of say a few millimetres between beam and septum and therefore angular spreads up to 1mrad (approx.) will not lead to losses. For this reason, only the m13 term will be considered in the further discussion. 3.9.1 Minimisation of chromatic effects One way of reducing the above mentioned disadvantages is to arrange the extraction such, that the momentum spread of the extracted beam is small. However, the momentum spread is not a free parameter as it is fixed when applying the Hardt condition to avoid losses at the electrostatic septum. It is therefore of prime importance to optimise the lattice in such a way that the superposition of the extraction separatrices results in a small momentum spread of the extracted beam. The explicit form of the element m13, expressed in terms of normalised dispersion,

P

P E sin cos

ES , ES , MS , MS 13 n n n

D D D m c

,

(3.30)

shows that the loss of space for the magnetic septum is proportional to EMS. Decreasing EMS reduces the influence of m13, but the gap created by the electrostatic septum (3.16) is also proportional to EMS and becomes smaller. The effective gap (3.27) for the magnetic septum is reduced by decreasing EMS.

gap

m m p p D D D p p

eff n n n

c
12

13

M G E E P P P G

MS ES MS ES ES

sin cos sin

, , ,

. (3.31)

The only effective approach to reduce chromatic effects is to minimise m13 directly.

SLIDE 60

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

58 3.10 Achromatic transfer between extraction septa A momentum-independent transfer between electrostatic and magnetic septa requires the chromatic term m13 in the general 3x3 transfer matrix to vanish. Inspection of (3.30) shows that the size of m13 is determined by the normalised dispersion functions at the

septa. In order to reduce m13, the lattice needs to provide suitable values of the dispersion

at the positions of the septa. The shape of the dispersion function in a lattice is determined by the distribution of the dipole magnets. In the analysis below, the effect of the dipoles is approximated as a point kick. In this model, every passage of a dipole adds a kick Dcn,0 to the actual value of the normalised dispersion function, according to

, D D D D D D

n n entrance to dipole n po kick

f dipole

n n n exit of dipole

c

c
c

c

,0

int ,0

.

(3.32)

In a bending-free region, the dispersion function acts like the betatron oscillation of a particle and can be described with the 2x2 Twiss-matrix formalism. The transformation of the dispersion function between two lattice elements, denoted by the suffices 1 and 2, without crossing dipoles is given by

1 2

c
c

D D m m m m D D

22 21 12 11

. (3.33)

In normalised phase space, the normalised dispersion vector simply rotates by the phase advance 'P between the elements, according to

1 2

c
'

'

'

'

c

n n n n

D D D D P P P P cos sin sin cos . (3.34)

Combining (3.34) and (3.32) allows an approximate analysis of the transfer between the extraction septa. Some particular cases are demonstrated below. Both septa in a bending-free dispersion region

Figure 3.21. 180°dispersion bump. (a) Normalised coordinates (P, Dn), (b) normalised phase space (Dn, Dcn).

Figure 3.21 shows a dispersion bump as created by two dipole point kicks of identical strength, spaced by 180° in phase. The electrostatic septum is positioned in the first half of the bump, the magnetic septum 90° later, in the second half.

phase 180° 90° MS ES

Dn Dcn Dn

kick dipole1 kick dipole2

(Dn, Dcn)ES (Dn, Dcn)MS

(a) (b)

SLIDE 61

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

59 To prove that in such a case the transfer is always achromatic, the closed orbit for particles with a momentum deviation Gp is considered. The closed orbit at any position s in the machine is to first order defined by the product of dispersion function and momentum deviation,

p

p s D p s x / ,

EQ.O.

G G

and
p

p s D p s x / ,

EQ.O.

G G

c

c . (3.35)

The 3x3 transfer matrix links the orbit coordinates between different positions (1 and 2) in the machine according to

c
c

p p p x p x m m m m m m p p p x p x / ) ( ) ( 1 / ) ( ) (

1 1 23 22 21 13 12 11 2 2

G G G G G G . (3.36)

By inserting (3.35) into (3.36), general expressions for m13 and m23 in terms of the dispersion functions and the 2x2 Twiss-matrix elements are found,

1 12 1 11 2 13

D m D m D m c

and

1 22 1 21 2 23

D m D m D m c

c

. (3.37)

In regions without bending the dispersion function transforms according to (3.33) and by inserting this into (3.37) it follows directly that m13 and m23 are zero and therefore: x The transfer via a dispersion region without crossing bending magnets is always achromatic with respect to position and angle. It should be noted that this result is exact and valid not only inside a dispersion bump, but in any lattice section without bending. A typical structure where this result can be applied is the so called ‘square’ lattice. The dipole magnets are grouped in the corners, to create dispersion bumps on two

pposite sides and dispersion-free straight sections on the remaining two sides. The fact of

the fully achromatic transfer between the septa may give the impression that such a structure is best suited for a slow extraction scheme, but there are two major problems with a square lattice: x The natural position for the electrostatic septum (i.e. the first half of the dispersion bump) is the worst place to put it when adjusting for the Hardt condition. It was shown in Section 3.8, that for ideal operation Dn > 0 and Dcn < 0, as in the second half of a bump is required, whereas in the first half Dn > 0 and Dcn > 0, as can be seen from Figure 3.21. x The second problem is that in general there is too little useable phase advance in the dispersion straight-section. An important fraction of the 180° is lost because of the length of the corner dipoles. To get a reasonable phase separation between the extraction septa, they have to be positioned close to the corners on either side of the straight section. Such a layout requires a very strong kick from the magnetic septum in

rder to clear the dipole magnets that close the dispersion bump. Moving the magnetic

septum closer to the centre of the drift space reduces the phase separation from the electrostatic septum. In this case, the electrostatic septum has to kick very hard, to create the space needed for the magnetic septum.

SLIDE 62

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

60 Both septa in regions with dispersion and bending

Figure 3.22. Extended dispersion bump, (a) normalised coordinates (P, Dn), (b) normalised phase space (Dn,Dcn).

Figure 3.22 shows a so-called ‘extended-dispersion-bump’. The free space within the dispersion bump is increased by adding a central dipole in order to overcome the problems

f the ‘square’ lattice. The disadvantage of such a structure is, that in general, the transfer

between the extraction septa is chromatic, as a bending magnet has to be crossed. Some design guidelines can be found for an extraction layout with optimised phase advances of P = 90q + n180q between the septa. Using expression (3.30) for the chromatic term,

P

P E sin cos

ES n, ES n, MS n, MS 13

D D D m c

,

(3.38)

and assuming a phase advance P = 90q + n360° between electrostatic and magnetic septa gives

ES

n, MS n, MS 13

D D m c

E

. (3.39)

To make the transfer achromatic (m13 = 0),

ES n, MS n,

D D c (3.40)

is required, but as shown earlier, one needs to work with a negative Dcn,ES for adjusting the Hardt condition, therefore m13 can only be made zero by having negative dispersion at the magnetic septum. For a phase advance of P = 270q + n360q (septa on opposite sides of the vacuum chamber) it follows

ES

n, MS n, MS 13

D D m c

E

(3.41)

and therefore to make the transfer achromatic,

ES n, MS n,

D D c

(3.42)

is required. In this case, m13 can be made zero by having a positive Dn,MS and a negative Dcn,ES just as required for the Hardt Condition. A disadvantage of this solution might be that the particles which are extracted, have to be transported for a longer distance in the machine (e.g. crossing of sextupoles between the septa would be more difficult to avoid).

If a sextupole is crossed (either resonance or chromaticity) between the ES and the MS, then there is a variable optical element in the

extraction channel. Any change in the Qc or resonance-strength alters the extraction geometry.

phase >180° P MS ES

Dn Dcn Dn

kick dipole1 kick dipole2

(Dn, Dcn)ES (Dn, Dcn)MS

kick dipole3

(a) (b)

SLIDE 63

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

61 Electrostatic septum in dispersion region and magnetic septum in zero-dispersion region

Figure 3.23. 180°dispersion bump. (a) Normalised coordinates (P, Dn), (b) normalised phase space (Dn, Dcn).

Figure 3.23 shows a 180° dispersion bump; contrary to the extraction layout in Figure 3.21, the electrostatic septum is positioned in the second half of the bump where Dn > 0 and Dcn < 0, as required for the Hardt condition. The magnetic septum is positioned in the dispersion-free straight section after the dipoles that close the bump. As for the extended dispersion bump, the transfer between the septa is in general chromatic due to the passage of dipoles. Assuming the dispersion bump is created by dipole point kicks, the normalised dispersion function can be described with (3.32) and (3.34) as

sin

, n n

c

D D and

cos

, n n

c

D D , (3.43)

where - is the phase advance counted from the first dipole kick and Dcn,0 is the strength of the kick. With Dn,MS = 0 in (3.30), the chromatic term m13 can be written as

P

P E sin cos

ES n, ES n, MS 13

c
D

D m , (3.44)

and by inserting (3.43) into (3.44) a simple expression for m13 is obtained,

P
E
sin

n MS 13 ,0

D m . (3.45)

It follows immediately that the transfer between the septa will be achromatic( m

13

) for

P
n 180

(3.46)

Obviously it is impossible to use exactly n = 1, since this gives the position of the dipole which is closing the dispersion bump. To keep m13 small, the magnetic septum has to be positioned as close to the dipole as possible. This solution has been adopted for the lattice of the medical synchrotron. The electrostatic septum is positioned in the second half of the dispersion bump, which allows the Hardt condition to be fulfilled with negative chromaticity and a small extracted momentum spread while keeping the chromatic effects small. The only resulting compromise is a reduced phase advance of less than 90° between the septa.

For larger n, there is again the problem of transporting the extracted part of the beam through a large distance in the machine.

Dcn Dn

kick dipole1 kick dipole2

(Dn, Dcn)ES

phase 180° P MS ES

Dn

T (a) (b)

SLIDE 64

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

62 Transfer for un-fulfilled Hardt Condition Adjusting for the Hardt Condition fixes the chromaticity and therefore the slope of the resonance line (3.18) and the momentum spread of the extracted beam (3.19). Instead of superimposing the extraction separatrices at the electrostatic septum, the chromaticity can be used to arrange the extraction geometry such that particles with different momenta arrive at the septum with different angles. With a proper choice of these starting conditions a reduction of the space available for the magnetic septum can be avoided even though the transfer between the septa is chromatic. This method is used in the present CERN-PS slow extraction scheme [15]. Transfers between zero dispersion regions By inspection of the general expression (3.30) for m13, it follows immediately that transfers between zero-dispersion regions are always fully achromatic. However, with the electrostatic septum located in a zero-dispersion region, the separatrices can be superimposed only for zero chromaticity, which will be discussed later. 3.11 Comparison of extraction techniques In order to slowly extract a beam, either the particles have to be moved progressively into the resonance stopband or the stopband has to be moved onto the beam. This means that some physical quantities have to be changed. The varying parameter, the way of changing it and the initial geometry of the beam in the amplitude-momentum (amplitude-tune) space distinguish the extraction method. In Section 3.3, a general overview of the different extraction techniques was given. 3.11.1 Varying the sextupole strength Extraction by changing the sextupole strength is generally not used. The main argument against a change of the resonance strength is found by considering expressions (3.11) and (3.12) for the spiral step and pitch that determine the size of the extracted beam. As both terms are proportional to the resonance strength S, any change of S during the extraction leads to a change of the beam size and therefore, (in the case of a medical machine), the beam spot at the patient, which is not acceptable. 3.11.2 Varying the tune distance The most common extraction method is to slowly change the tune distance GQ, of the circulating particles. There are two fundamentally different approaches. x To move the resonance through the beam by changing the betatron tune of the machine with quadrupole magnets. x To move the beam into the stationary resonance by acceleration (deceleration) of the circulating particles. The change in momentum translates via the chromaticity into a change of the tune (3.17). The main problem, when moving the resonance through the beam, is that the Hardt condition is ineffective. The extraction separatrices can be superimposed at any time but during the extraction they change their position. This means that the divergence of particles that arrive at the electrostatic septum changes from the beginning to the end of the extraction process. This variation in angles also affects the transfer between the extraction septa and reduces the space available for the magnetic septum. For the

SLIDE 65

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

63 compensation of these effects, e.g. a dipole correction bump, that changes the orbit around the electrostatic septum, can be used. However, such a dynamic correction complicates the operation of the machine and should be avoided, especially for a medical machine where simplicity and reliability of operation are of major concern. Figure 3.24 (a) shows the principle of this method for a beam with a large horizontal emittance Ex and a small momentum spread. The first particles that are extracted have large amplitudes and a large tune distance (GQ = QcGp/p), the small amplitude particles are extracted at the end of the spill with a small GQ. The variation of the mean tune distance

f the extracted particles results in the problem of a changing beam size during extraction

as can be seen from equations (3.7). The extracted momentum spread is always equal to the total momentum spread of the beam. Figure 3.24 (b) shows the same extraction technique for a beam with a small transverse emittance and a large momentum spread. During the whole spill particles of all amplitudes are extracted, the beam size remains therefore unchanged. The instantaneous momentum spread of the extracted beam is small but the mean momentum varies. This has two undesirable features. x Firstly, it could require dynamic changes of the septa strengths and the magnets in the extraction line within a single spill. x Secondly, for a medical machine, a change of the mean momentum of the beam during the spill causes a variation in the depth of the Bragg-peak in the target volume.

Figure 3.24. Slow extraction via changing the tune by moving the resonance. Extraction via tune change by moving the resonance

Amplitude Amplitude Amplitude 'p/p 'p/p 'p/p

unstable particles unstable particles

Amplitude Amplitude Amplitude 'p/p 'p/p 'p/p

unstable particles unstable particles

Large horizontal emittance Small momentum spread Small horizontal emittance Large momentum spread

x Size of extracted beam varies

x Momentum of extracted beam is constant

x Momentum of extracted beam varies

x Size of extracted beam is constant

(a) (b)

SLIDE 66

Chapter 3 Resonant slow extraction

________________________________________________________________________________________________________________________________________________________________________________________________________

64 For the second method, the particles are accelerated and their tune changes due to the chromaticity according to (3.17). In this case the resonance is stationary, the Hardt condition can be applied, the transfer between the septa is unchanged during the

extraction. A further advantage is that all lattice parameters are kept constant during the

extraction and all lattice power supplies have constant settings. Figure 3.25 (a) shows the extraction process for a beam with a large horizontal emittance and a small momentum spread. As in Figure 3.24 (a) the first particles that are extracted have large amplitudes, the main problem is that the stationary resonance transforms the horizontal emittance into a momentum spread. It can be seen that the mean momentum of the extracted beam varies during a single spill, resulting in the same difficulties as discussed for Figure 3.24 (b). The most preferable situation is shown in Figure 3.25 (b). The beam has a large momentum spread and a small horizontal emittance. During the whole spill time particles

f all amplitudes are extracted. The momentum spread of the extracted beam is small and

the mean momentum is constant. With such a geometry, all machine parameters can be kept constant and no dynamic correction systems are needed.

Figure 3.25. Slow extraction via changing the tune by acceleration (deceleration) of the beam. Extraction via tune change by acceleration (deceleration ) of the beam

Amplitude Amplitude Amplitude 'p/p 'p/p 'p/p

unstable particles unstable particles

Amplitude Amplitude Amplitude 'p/p 'p/p 'p/p

unstable particles

Large horizontal emittance Small momentum spread Small horizontal emittance Large momentum spread

x Size of extracted beam varies x Momentum of extracted beam varies x Size of extracted beam is constant x Momentum of extracted beam is constant

unstable particles

(a) (b)

SLIDE 67

Resonant slow extraction Chapter 3

________________________________________________________________________________________________________________________________________________________________________________________________________

65 Acceleration (deceleration) of the beam into the resonance There are a number of different methods for the acceleration of the beam into the stationary resonance. The best known are listed below, this is a partial list and a complete discussion of the different techniques would go beyond the scope of this thesis. x Acceleration of the beam with the rf cavity. x Use of longitudinal stochastic noise [16]. x Use of a betatron core, inductive acceleration [17,18]. 3.11.3 Varying the amplitude (zero chromaticity extraction) The third method for slow extraction is to increase the particle amplitudes until they enter the resonance stopband. This method is particularly well suited to an extraction set-up based on a machine with zero chromaticity. In Section 3.1 it was shown that the stability of a particle, for a given sextupole strength, is determined by its tune distance GQ and its amplitude A. From (3.17) it follows with Qc = 0 that all particles have the same tune distance and their stability is therefore

nly a function of the particle amplitude. This in turn means that all particles that become

unstable will be extracted from stable triangles of equal size, independent of their

momentum. From the expression for the Hardt condition (3.24), it can be seen that the

right hand side will be zero for Qc = 0. Inspection of (3.25) shows that in this case, the separatrices can be superimposed by simply positioning the electrostatic septum in a zero- dispersion region. In order to further optimise the extraction geometry, the transfer between electrostatic and magnetic septa should be achromatic. In Section 3.10, it was shown that transfers between zero-dispersion regions are always fully achromatic. Therefore the magnetic septum should also be positioned in a dispersion-free region, 90° downstream from the electrostatic septum. Finally, the extraction process is performed by blowing up the beam with transverse rf excitation. Rf kicks applied at the betatron frequency (product of fractional horizontal tune and revolution frequency) cause the betatron amplitudes to grow and particles to enter the stopband. Alternatively, transverse stochastic noise can be used to blow-up the beam.

Figure 3.26. Zero chromaticity extraction by transverse beam blow up.

Figure 3.26 shows the amplitude momentum space and the resonance stopband for the zero-chromaticity extraction. During the extraction all machine parameters are kept constant; as for the extraction technique with the stationary resonance, no dynamic correction systems are needed. The extracted beam has the momentum spread of the circulating beam, the mean momentum and the beam size are constant during the

extraction. The only disadvantage of this extraction method is the lack of intrinsic

transverse stability of the coasting beam, due to the zero-chromaticity.

Amplitude 'p/p Resonance Stopband Amplitude 'p/p

unstable particles

SLIDE 68

67

Chapter 4 TRANSIT TIME

The precise control of the slow-extracted beam is a key point in designing a synchrotron for hadrontherapy. The flux of extracted particles, in particular, should be as smooth as possible in order to avoid overdosing. A more detailed study of the dynamics of the extraction is thus needed to get a deeper understanding and to choose the extraction scheme that suits the need for a constant spill best. The first step is to evaluate the time that a particle spends in the accelerator after becoming unstable and before being extracted, which will be referred hereafter as transit time [19]. 4.1 Translated hamiltonian The motion of the particles is described by the Kobayashi hamiltonian (2.31). This is an approximate hamiltonian, chosen because it can be manipulated relatively easily and because numerical simulations show that the final expressions describe the physics to an accuracy of a few percent. In Figure 4.1 the geometry of phase space for positive H and S is shown with the direction of motion of a particle just out of the stable triangle illustrated by dashed arrows. The equivalent pictures when H and/or S change sign can be determined by mirroring/rotating Figure 4.1 as shown in Chapter 2. Furthermore the position in which

ne of the separatrices intercepts the electrostatic septum downstream in the lattice, is

shown.

A B C h X’ X P0 Unstable region surrounding stable triangle stable region P1 P3 P2 Electrostatic Septum Intercept

Figure 4.1. Separatrices geometry at the sextupole.

Since, as already mentioned in Chapter 2, after each revolution the separatrices exchange their position and since the infinitesimal time considered is three revolution periods, all the separatrices are equivalent and there is no loss of generality in choosing a particular one.

SLIDE 69

Chapter 4 Transit Time

________________________________________________________________________________________________________________________________________________________________________________________________________

68 When making an extraction, the stable area is shrunk over a period of the order of 1 s, whereas the revolution time in the machine is typically of the order of 1 Ps and the transit time of the order of some hundreds of revolutions. Thus, particles that cross the limiting separatrix of the stable area stay very close to the separatrices during their

extraction. In the following, the particles that become unstable along separatrix A in

Figure 4.1 will be considered. The proximity of the extraction trajectories to the separatrices A and C means that the influence of the stable point P3 will be strong and the particles will spend most of their time in the vicinity of this point. It is convenient therefore to change the coordinate system and translate the origin to P3, as illustrated in figure 4.2.

A B C h

' X X

P0 Unstable region surrounding stable triangle Stable region P1 P3 P2 Electrostatic Septum Position

Figure 4.2. Translation of the reference system to stable point P3.

The Kobayashi hamiltonian H [m], in the new coordinates, reads

) ' 3 ' 3 6 6 4 ( 4

3 2 2 3

X X X X X h X h h S H

(4.1)

where

S h H 3 2

[m1/2] is the apothem of the triangle and S [m-1/2] is the normalised sextupole strength. To simplify the notation, from now on we will omit the bar over H, X and X’. The Kobayashi hamiltonian has been derived by approximating the phase-space displacement of a particle every third revolution of the accelerator, and considering the time of 3 turns as the infinitesimal increment for time. Thus, in the following, the unit of measurement of the time is 3 revolution periods of the machine, dt { 3TREV but is dimensionless, as already mentioned in Chapter 2. The times derived from the following formulae, multiplied by 3, will thus give the number of turns necessary to reach the electrostatic septum.

SLIDE 70

Transit Time Chapter 4

________________________________________________________________________________________________________________________________________________________________________________________________________

69 The equations of motion in the new reference frame are given by:

).

( 3 ' 3 ' 3 6 12 4 ' ) ( ' 6 3 6 4

2 2

b X X hX hX S X H dt dX a XX hX S Y H dt dX

w

w w w (4.2)

The transit time for a particle will be worked out integrating the equations of motion between the starting position and the electrostatic septum position. 4.2 Extraction time with constant parameters If H and S are constant during the extraction time, the particle will follow a trajectory with constant H. Thus it is possible to work out the trajectory equation by equating the general hamiltonian to the value that it assumes with the initial conditions (X0, Xc0). From (4.1)

). ' 3 ' 3 6 6 4 ( 4 ) ' 3 ' 3 6 6 4 ( 4

3 2 2 3 3 2 2 3

X X X X hX hX h S X X X X hX hX h S

(4.3)

The total transit time is calculated in the following by evaluating and summing the time needed by the particle to travel from P3 to the electrostatic septum and the time from the position where the particle jumps across the separatrix to P3. As an example the expressions derived, will be evaluated for the following static parameters: x TREV = 0.5 Ps x Spill length = 0.5 s = 106 TREV x GQ = 3.782 10-3 x H = 6SGQ = 7.13 10-2 x O† = 15 10-6 x S = 27.7897 m-1/2 x h = (2/3)(H/S) = 1.71 10-3 m1/2 x n‡ = 4.79 4.2.1 Time from nearest fixed point In the vicinity of the fixed point P3, that is |X|, |Xc| << h, the third order terms in X and Xc can be neglected. Simplifying, the trajectory equation (4.3) results in

X X X X X X 3 ' 3 '

2 2

(4.4)

† The value of O is chosen to correspond to the distance moved by the separatrix during 3 turns in the example in next sections, when

acceleration is taken into account. Its meaning will be explained later.

‡ The value of n is a measurement of the ES position. Its definition is given later in this paragraph.

SLIDE 71

Chapter 4 Transit Time

________________________________________________________________________________________________________________________________________________________________________________________________________

70 where X0 and Xc0 are the coordinates of the starting point. The substitution of this trajectory into the equation of motion for X, equation (4.2) (a), yields:

. 3 6 ' 6 3 6 3 6 4

2 2

X

X X X hX S dt dX (4.5)

This equation can be integrated by variable separation. By definition, the trajectory equation (4.4) is valid only close to P3, but as the particle approaches asymptotically the separatrix

X X 3 1 '

the third order terms in the hamiltonian cancel out. With this

consideration they can also be neglected far from P3 along the outgoing separatrix and integration can be extended to the electrostatic septum. Using the standard integral formula:

ac b b ax ac b b ax ac b dx c bx ax 4 2 4 2 ln 4 1 1

2 2 2 2

(4.6)

valid for b2 > 4ac, the transit time Tc,s can be written as

ES

c

X X T t t s c

dX X X X hX X S dt T

2 2 ,

' 3 3 1 3 2

.

' 3 4 9 3 2 ' 3 4 9 3 2 ln ' 3 4 9 1 3 2

2 2 2 2 2 2

ES

X X

X X X h h X X X X h h X X X X h S

Since |X0|, | Xc0| << h, we can expand the square roots to first order:
|

ES

X X s c

h X X X h h X h X X X h h X Sh T

2 2 2 2 ,

' 3 9 2 1 3 3 2 ' 3 9 2 1 3 3 2 ln 3 3 2 (4.7) . 3 ' 3 3 3 ' 3 ln 3 3 2

2 2

ES

X X

h X X X h X h X X X X Sh

|

It is useful to express all the distances in units of h, thus we define

SLIDE 72

Transit Time Chapter 4

________________________________________________________________________________________________________________________________________________________________________________________________________

71

XES -nh, X0 -Oh, h X /

3

2 '0

. These definitions are illustrated in Figure 4.3

A B C h

' X X

P0 P1 P3 P2 ES X0 = -Oh Xc0 = -23/h XES = -nh

Figure 4.3. Definition of O, / and n.

For O = 0, / = 0 corresponds to the stable point P3 and / = 1 to stable point P2. Using the new parameters O, / and n in (4.7) gives:

/
/
/
/
|

h h h h h h h nh h nh Sh T s

c

O O O O O O O O O O 3 3 1 3 3 1 3 3 3 1 3 3 3 1 ln 3 3 2

2 2 2 2 ,

O 3 3 ln 3 3 2

|

n n Sh

. Inside the logarithm, the term

h

/

O

O 3 3 2

2

has been neglected with respect to nh in the first fraction and with respect to Oh in the denominator of the second fraction (since / is small near to P3). In the numerator of the second fraction only the term 3h has been kept. Remembering that

S h H 3 2

, it results that the transit time from the nearest fixed point in the static case is

SLIDE 73

Chapter 4 Transit Time

________________________________________________________________________________________________________________________________________________________________________________________________________

72

. 3 3 ln 3 1

,

O H

|

n n T s

c

(4.8)

Note that, provided / << 1, the time needed for a particle to be extracted is independent

f the initial Xc coordinate. For the numerical data given,

Tc,s = 95 which means that the particle needs 285 revolutions to reach the ES. A numerical simulation of the extraction process, yields 312 revolutions for / = 0.1 and 306 revolutions for / = 0. The simulation just advances the particle, applying the displacement that it undergoes each third revolution, that is the Kobayashi equations (2.29), until it reaches the ES. The number of iterations multiplied by 3 gives the number of revolutions. This simplified tracking is necessary in order for the used quantities h, O etc to be

meaningful. An exact tracking, as will be seen better in Chapter 6, yields to a distortion of

the separatrices that destroys the simple parameterisation used herein. The validity of this simple description will be confirmed by the results described in Chapter 5, which are checked against an exact tracking. 4.2.2 Time spent along the side of the stable triangle For a particle starting somewhere along the side of the stable triangle, far from the unstable fixed point P3, the extraction time is given by the sum of the time needed to travel to P3 and the time needed to go from P3 to the ES. The motion along the side of the triangle is determined by the second equation of motion, (4.2) (b). The trajectory of a particle starting near the side of the stable triangle is not a straight line but it remains close to the side for most of its length. Until |X| << |Xc| it is possible to approximate equation (4.2) (b), by neglecting the terms in X and X2, so that

2

' 3 ' 3 6 4 ' X hX S dt dX

|

(4.9)

which, similar to the previous case, can be integrated by variable separation, yielding:

A A T

s s

/
/

| 3 2 1 ln 3 1

, 1

H (4.10)

where the integration has been done between the starting point

h X /

3

2 '0

and

Ah X F

'

a point near the corner where the curve will be joined to the trajectory calculated in the previous Section 4.2.1. As the junction point has to be near P3 to allow the use of (4.8), A has to be smaller or of the order of 0.1 3 2 , which corresponds to / = 0.1.

SLIDE 74

Transit Time Chapter 4

________________________________________________________________________________________________________________________________________________________________________________________________________

73 4.2.3 Time for particles starting along the side of the stable triangle To correctly join the two trajectories, it is necessary to evaluate the position of the start of the second. To do this we can evaluate the hamiltonian (4.1) at the starting point

h

h /

3

2 , O

and work out the XF coordinate at XcF = -Ah. From

)

3 2 3 3 2 3 6 6 4 ( 4 ) 3 3 6 6 4 ( 4

3 2 2 3 3 2 2 3

h h h h h h h h h S X Ah X Ah hX hX h S

F F F F

O O O O

/
/
neglecting the higher order terms in O and XF, XF is obtained
.

3 2 1 12 3 2 1 12 h h A h A A X

F F

O O O

{

/

/
|
/
/
(4.11)

Starting from this point it is possible to evaluate the total time needed to reach the electrostatic septum:

. 3 3 3 2 1 ln 3 1 3 3 3 2 1 ln 3 1

, , 1 , F F s c s s s s

n n A n n A A T T T O H O H

/
/

|

/
/

|

And finally, by substituting the definition of OF (4.11), the transit time for a particle starting

along the side of the stable triangle in the static case can be evaluated as

O

H 3 3 1 1 ln 3 1

2 ,

/
|

n n T s

s

(4.12)

which does not depend on A, that is on the exact position of the junction point between the two trajectories. This expression is valid until (1-/) >> O, such that the velocity in Xc does not depend on X. A comparison of (4.12) with numerical simulations is summarised hereafter / Formula (4.12) (revolutions) Numerical Simulation (revolutions) 0.1 285 312 0.5 312 339 0.9 393 423 0.999 615 675 4.3 Extraction time with varying parameters Resonant slow extraction can be “activated” in many ways, but the most common is to vary the tune Q of the machine, by varying the focusing quadrupoles, to bring the beam into resonance. Alternatively it is possible to vary the momentum of the particles, as GQ = Qc 'p/p. Varying GQ causes the triangle to shrink, in fact H SG 6 Q is directly

SLIDE 75

Chapter 4 Transit Time

________________________________________________________________________________________________________________________________________________________________________________________________________

74 proportional to the apothem h of the stable triangle. Reducing the size of the triangle causes some particles to pass from the stable region inside the triangle to the unstable region outside and thus to be extracted. In this Section a linear variation of the tune with time, that is Q = constant, will be considered. The main approximations used in the following are: x the relative variation in the size of the triangle during the transit time will be small, that is 'h << h; x instead of considering the movement of the separatrices as the triangle shrinks, the relative movement of the separatrix to the particle, will be considered as an additional contribution to the particle velocity and the triangle will be considered fixed during the extraction time. This approximation is illustrated in Figure 4.4, where the situation is sketched at two different times t1 and t2, with t2 > t1.

Approximation Particle position at t2 Real situation Particle position at t1 Stable region at t2 Stable region at t1 Particle position at t2 Particle position at t1 Stable region at t1 & t2 'h 'h

Figure 4.4. Relative motion of particle and stable triangle.

For numerical checks few more data are needed beyond the ones used in the static case. It will be assumed: x 'GQ per turn = -1.891 10-8 x

Q = 'GQ in 3 turns = -5.673 10-8

x O = -'h in 3 turns/h = - Q/(GQ)0 = 15 10-6 4.3.1 Time from nearest fixed point It was shown in the static case, that in this region the extraction time does not depend on the initial value /. Following this hint, / is set to zero and the particle is assumed to move

n the outgoing separatrix

X X 3 1 '

. With this assumption, equation (4.2) (a) becomes

SLIDE 76

Transit Time Chapter 4

________________________________________________________________________________________________________________________________________________________________________________________________________

75

2

3 2 3 X hX S dt dX

.

When considering acceleration, it is necessary to add the velocity of the separatrix as it recedes from the particle. It is easy to evaluate it in the non-translated frame and to note that the velocity has to be the same in both frames. In this frame, the separatrix equation is given by X = -h. Thus:

dt dQ S dt dh S 4

.

(4.13)

As the separatrix moves to the right, the particle moves to the left with respect to it. Thus, the relative velocity of the particle with respect to the separatrix, is given by:

Q

S X hX S dt dX

S

4 3 2 3

2

.

(4.14)

This equation can again be integrated by separation of variables and the transit time results:

|
2

2 2 2 2 2 2 2 2 ,

4 3 2 4 3 2 3 2 4 3 2 3 2 4 3 2 ln 3 1 3 3 8 2 1 1 3 3 3 3 3 8 2 1 1 3 3 3 ln 3 1 3 8 3 3 3 3 8 3 3 3 ln 3 8 3 1 H S O H S O H S H S H H S H H H S H H H S H H S H H S H

O O

Q Q Q n Q n Q SX Q SX Q SX Q SX Q T

nh h nh h d c

whence the transit time from the nearest fixed point in the dynamic case is
|

2 2 ,

3 1 3 3 ln 3 1 4 3 2 3 2 3 ln 3 1 H H O H H S O H

n

n Q n n T d

c

(4.15)

where the assumption 'h << h which translates into:

SLIDE 77

Chapter 4 Transit Time

________________________________________________________________________________________________________________________________________________________________________________________________________

76

2

3 2 1 3 2 3 1 4 H S H H S

Q

S S Q

has been used and O and

2

H Q

have been neglected with respect to unity. In principle,

O H H is also negligible with respect to H H 2 . This meams that particles start within one step of the separatrix and the separatrix moves many times during extraction. The term in O has not been suppressed in (4.15) in order to recover (4.8) when

Q

.

In the numerical example the transit time results Tc,d = 231 machine revolutions, to be compared with the numerical simulation result of 249 turns for / = 0, and 252 turns for / = 0.1 4.3.2 Time spent along the side of the stable triangle The next step is to evaluate the time needed to travel along the side of the stable triangle. The same approximations will be used as before, with the only difference that the velocity

f the separatrix has now to be multiplied by

3 2 , to obtain the relevant component of the

motion of the extraction separatrix C in Figure 4.2. The side of the stable triangle is also moving, but the motion in Xc is quasi-independent of O. Separatrix C moves with the same speed as separatrix A in the previous section, but only the Xc component of its motion is of

interest. Equation (4.2) (b), provided that X remains negligible during the particle

movement, becomes

Q

S X hX S dt dX

3

8 ' 3 ' 3 6 4 '

2

S

(4.16)

which, integrated between

h X /

3

2 '0

and

Ah X

d F

,

'

gives:

A A T

d s

/
/

| 3 2 1 ln 3 1

, 1

H (4.17)

which is exactly the same as in the case of the static situation, equation (4.10). This is reasonable considering the approximation used, in which the velocity in Xc does not depend on X and in which the variation in the length of the triangle side is negligible in the time considered. 4.3.3 Time for particles starting along the side of the stable triangle It only remains to evaluate the point where the second integral begins. The assumption is made that this point can be calculated by adding the position calculated under stationary conditions using (4.11) and the movement of the separatrix, evaluated as the product of the time needed to move along the side of the triangle times the velocity at which the separatrix moves. The starting point for the second integral is then:

SLIDE 78

Transit Time Chapter 4

________________________________________________________________________________________________________________________________________________________________________________________________________

77

Ah

X h Q S A A h A A X

d F d F d F

/
/
/
/
,

, ,

' 4 3 2 1 ln 3 1 3 2 1 12 O S H O

(4.18)

The total time to reach the ES for a particle starting from

H

H /

3

2 , O

with / greater

r of the order of 0.1, taking into account the shrinking of the stable triangle, is given by

the sum of the time to move along the side and the time to reach the ES from the junction point near P3, that is:

Q Q n n A A T T T

d F d F d c d s d s

2

, 2 , , , 1 ,

4 3 2 4 3 2 3 2 3 ln 3 1 3 2 1 ln 3 1 H S O H S O H H

/
/

|

And finally the transit time for a particle starting along the side of the stable triangle in the

dynamic case can be evaluated to

/
/
/
/
|

. 3 1 3 1 3 2 3 ln 3 1 3 2 3 1 3 2 3 ln 3 1

2 , 2 , ,

H H O H H S O H

d

F d F d s

A A n n Q A A n n T (4.19)

This time, A does not fortuitously disappear, as in the static case, thus a value must be given for it. As the two regions (near and far from P3) have been previously divided at / = 0.1, it seems reasonable to choose A such that the junction point corresponds to / = 0.1. This means A

01 2 3

. . Evaluating (4.18) and (4.19) with the numerical data used earlier, and comparing with the simulation, yields: / Formula (4.19) (revolutions) Numerical Simulation (revolutions) 0.1 231 252 0.5 258 267 0.9 300 306 0.999 399 432 For an initial tune distance GQ = 3.782 10-4, which corresponds to an amplitude 10 times smaller and thus a single particle emittance 100 times smaller, one would have the following transit times:

SLIDE 79

Chapter 4 Transit Time

________________________________________________________________________________________________________________________________________________________________________________________________________

78 / Formula (4.19) (revolutions) Numerical Simulation (revolutions) 0.1 1323 1431 0.5 1575 1584 0.9 1986 2061 0.95 2136 2349 Particles with / = 0.999, as discussed in the next chapter, are extracted along the following separatrix and are thus not reported.

SLIDE 80

79

Chapter 5 TIME PROFILE OF THE EXTRACTED BEAM

In Chapter 4 the time needed for a particle to be extracted after becoming unstable has been evaluated. The transit time formulae can be used to estimate the temporal structure

f the extracted beam, that is the number of particles reaching the electrostatic septum as a

function of time, when the beam is slowly brought into the resonance. This will be expressed as the density in time P(t), defined such that P(t) dt is the number of particles that reach the ES between t and t+dt. This function depends on the extraction method and can thus be used to compare different extraction techniques. 5.1 “Strip” profile During the extraction process, the stable region is shrunk slowly, such that a narrow strip*

f particles becomes unstable and is thus extracted. This situation is represented in Figure

5.1.

A B C

' X X

P0 P1 P3 P2 Electrostatic Septum Position

Figure 5.1. Shrinking the stable region leaves a narrow strip of particles unstable. In the figure

nly one of the three sides of the triangle is considered and the motion of the particles

is shown by the arrows.

The first particle will reach the ES at t = Tc,d and the last one at t = Ts,d(/=/F). Note that the maximum time does not correspond to / = 1. This is due to the fact that the particles that start very near to the stable point P2 will be overtaken by the inward movement of separatrix B. These particles will be extracted along the following separatrix. This is

* The name “strip” will be reserved for the narrow region of monoenergetic particles on the side of a stable triangle. The name “band” will

be introduced later for the series of strips in a beam with a momentum spread.

SLIDE 81

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

80 equivalent to considering them as starting near P3 (they are in fact replaced by particles from the separatrix C that are overtaken by the movement of separatrix A). The delay with which they appear near P3 will not be considered, for the moment. The value of /F will be given later in this Chapter. As noticed in Chapter 4, all the particles starting with 0</<0.1, reach the electrostatic septum at virtually the same moment, so the spill will start with a spike at t = Tc,d. If U(/) is the linear probability density of particles in the strip, and Nstrip is the total number of particles contained in the strip, the spike will contain Nspike particles, where:

separatrix

by the

vertaken

Particles 1 3 P

f

vicinity in the Particles 1 .

) ( ) (

/

/ /

/

/

F

d N d N N

strip strip spike

U U (5.1)

where

1 ) (

1

/ /

d

U

. Thus the initial spike can be approximately described by:

Nstrip Pspike(t)dt = Nspike G(t-Tc,d)dt. (5.2)

After the initial spike, the spill shape can be evaluated by noticing that the particles coming

ut between T(/) and T(/)+dt are the ones which started between / and /+d/. This

means:

dt dt d N d N dt T P N dt t P N

strip strip tail strip tail strip

/ / / / / ) ( ) ( )) ( ( ) ( U U (5.3)

for t0 < t < tF, which corresponds to 0.1 < / < /F. The time profile for the elementary strip is then given by the sum of Pspike and Ptail.

Pstrip (t)dt = Pspike (t)dt + Ptail (t)dt. (5.4)

Thus to evaluate the time profile of the elementary strip Pstrip, /F, U(/) and /(t) have to be

determined. The shape of this elementary spill is the key to calculating the spill profile for

the whole beam and for including the influence of ripple (an explanation of ripple is given later in Section 5.8). 5.2 Linear probability density Assuming that the initial beam is smoothly distributed in phase space and that the resonance is applied adiabatically, then the density probability U(/) in the strip will be proportional to 1/v(/), where v(/) = dXc/dt is the velocity in normalised phase space of the particles. With the Kobayashi Hamiltonian translated to P3, v(/) is given by (4.2) (b)

2

2 3

' 3 ' 3 6 12 4 ' X X hX hX S X H dt dX

w

w . (4.2) (b)

SLIDE 82

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

81 Restricting (4.2) (b) to trajectories close to separatrix A, substituting the definitions of O and / neglecting the term in O2, yields

2

2

3 3 3 ) ( ' /

/
/

O Sh v dt dX

which yields

2

2

3 3 3 1 ) ( 1 ) ( /

/
/

/ O U Sh k v k

where k is a normalisation constant defined by

1 ) (

1

/ /

d

U

. Thus

3 ln 2 9

2

O Sh k

.

Finally, the line density is given by

3

ln 3 3 2 3 ) (

2

O O U

/
/

/ . (5.5)

5.3 Inverting the Transit Time The second element that is necessary for the evaluation of Ptail, is d//dt. The time needed for a particle to reach the electrostatic septum starting from X = -Oh and

h X /

3

2 '

, near separatrix A, is:

2 , ,

3 1 3 1 3 2 3 ln 3 1 H H O H

/
/
d

F d s

A A n n T . (4.19)

Re-arrangement of (4.19) to solve for /, and the use of t instead of Ts,d, to stress that it is now the independent variable, yields

/
/
H

H H O

H

3

1 3 1 3 2 3

, 3 d F t

A A n n e

SLIDE 83

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

82 and, using (4.18),

/
/
/
/

/

/
1

3 2 ln 1 3 1 3 2 1 12 1 1 3 2 3 3

3

A A A A A A e n n

t

H H H O

H

.

(5.6)

As O is of the order of

H H

, and H << 1, the first term in the denominator of the last

fraction in (5.6) can be neglected with respect to the second. Equation (5.6) can then be rewritten as

R K K R ln

(5.7)

where

t

e n n K A A R

H

H H

3 2

3 3 3 ; 1 3 2

/
/
.

(5.8)

In order to invert T(/), it is necessary to solve (5.7) with respect to R. This equation has real solutions only when K t 1, which corresponds to t t Tc,d. This is to be expected since no particle is extracted prior to Tc,d. Let R = R(K) be the solution of equation (5.7), then:

) ( 1 ) ( K g K g

/

(5.9)

where

) ( 3 2 ) ( K R A A K g

. Referring back to equation (5.3), the aim is to evaluate
dt

d dt d dt d / /

/

| /

/
/

/ / 3 ln 3 3 2 3 3 ln 3 3 2 3

2 2

O O O U . (5.10)

Substitution of expression (5.9), yields:

K

R K R K g K g K g K g K g K g dt d

3

ln 2 1 3 ln 2 1 1 3 ln 2 1

2 2

O O O U

/

/ (5.11)

which no longer depends on A. Note that neglecting O with respect to (/2 - / in equation (5.10), means that / >> O, which is always true for / > 0.1, and (1 - /) >> O. This will be shown later to be also true for / < /F. Derivation of (5.7) with respect to time, gives

SLIDE 84

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

83

R R K R K K R

)

ln( K R K R R

1

) ln( 1

which, substituted into (5.11), gives

t e n n Q K K R K R K K K R K R dt d

H

H S O O U

3

3 ) ( )) ( ln( 1 3 ln 1 ) ( )) ( ln( 1 3 ln 2 1

/

/

.

(5.12)

Note that neglecting the derivatives of n and H was already included in the model when it was assumed that the stable region does not change during the extraction time. A plot of expression (5.12) is shown in Figure 5.2 compared with a numerical simulation.

100 200 300 400 500 600 10 20 30 40 50

Figure 5.2. Strip Profile (blue line) shown with simulation results (points).

When K >> 1, then R >> K and K can be neglected in the denominator. From (5.7),

K R R 1 ) ln( 1

. Thus, for K >> 1, the asymptotic value of the spill density is

3 ln 2 3 3 ln 2 1 ) ( O H O U

|

/ /

f

K

K dt d

t

(5.13)

5.4 Spill length The end of the spill (and thus the spill length) can be evaluated by noticing that separatrix B in its movement overtakes some particles with / sufficiently close to unity that their velocity is slower than the velocity of the separatrix itself. Those particles will then be extracted along the following separatrix. The effect is that some particle near P2 will

Ptail (t) Turns

SLIDE 85

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

84 disappear to reappear near P3. The spill therefore ends with the particle whose velocity is equal to the velocity of the separatrix. The Xc coordinate of this particle is given by

Q

S X X hX hX S dt dX

3

8 3 ' 3 ' 3 6 12 4 '

2 2

S

.

(5.14)

Neglecting the term in X and X2 with respect to the right-hand side, and substituting

h X /

3

2 '

, yields

3 8 9 9

2 2 2

/
/

Q S Sh Sh

S

whose solutions are

2 3 9 16 1 1

2 2 2 , 1

/

h S Q

#

S

. The solution of interest is the one close to the stable point P2,

2 2 2 2

3 3 1 3 2 1 3 9 8 1 H H H S S

/

Q h S Q

F

(5.15)

Note that 1-/F >> O. This justifies the omission of O in equation (5.10) in the previous

section. The spill will thus end at
|
/
/
/
/
|

/

F F F F F d s

A A A A n n T 1 3 2 ln 1 3 3 1 3 2 3 ln 3 1

2 ,

H H H

1

2 2 2

3 3 3 2 ln 1 3 3 3 3 3 2 3 ln 3 1

|

H H H H H H H

A

A A A n n

where the term in O in the last fraction has been neglected as in (5.6). Expressing this time in units of Tc,d , yields

|

/

H

H H H H H H H H H

2

1 2 2 2 , ,

3 3 3 ln 3 1 3 3 3 2 ln 1 3 3 3 3 3 2 3 ln 3 1 n n A A A A n n T T

d c F d s

SLIDE 86

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

85

2 3 3 ln 3 ln 3 3 3 2 ln 1 ln 3 3 ln 2 3 2 ln 3 ln

2 2 2

|

H

H H H H H

n

n A A A A n n

. Where all the addenda have been considered negligible with respect to

H H

2

ln

. Hence this complicated derivation converges to the simple result,

d

c F d s

T T

, ,

2 / . (5.16)

Thus the first particle reaches the electrostatic septum at Tc,d, the last one at 2 Tc,d, and the spill length is Tc,d. 5.5 Width of the initial spike It is now possible to evaluate the time needed for the last particle to be overtaken by the

separatrix. Consider the velocity of the particle along the side of the stable triangle and,

similar to Chapter 4, add the velocity of the separatrix:

Q

S Sh dt dX

3

8 3 3 3 '

2 2

S O

/
/
.

This yields

Q
H

S O H 2 3 3 3

2

/
/

/

. Finally, the time for a particle starting at /i to be overtaken by separatrix B is

H

H H H H H H S O H H H

3

2 1 3 2 3 2 ln 3 1 2 3 3 3 1

1 2

/

/ | /

/
/
/

i i

i

d Q T

For /i -> /F, T goes to infinity, as expected since this is the time needed for the separatrix to reach a particle which moves with the same velocity. So it is necessary to consider /i = /F +

3 2O

, that is one step of the separatrix away from /F. Then

H H 6 1 ln 3 1 T (5.17)

which may be a large fraction of Tc,d for small H. However, for the sake of simplicity, and considering that most of the particles are near the stable points and thus will be overcome in a few turns, this delay will be neglected and all the particles starting near P2 will be considered as if they were starting near P3.

SLIDE 87

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

86 5.6 Population of spike and tail It is interesting to evaluate the fraction of the beam contained in Pspike and Ptail. This can be done in a straightforward way by integrating the density U(/). It is easier to integrate Ptail between / = 0.1 and / = /F, as in this range O can be neglected in the density. It results,

2 1 3 ln 2 3 27 ln

1 .

| / /

/

O H O U

F

d (5.18)

neglecting ln(273 H) and ln(3) with respect to ln(O) yields the result that 1/2 of the beam is in Ptail and 1/2 in the initial spike. Thus the two components of the strip profile are equally important. This result is of consequence for the response to ripple and the efficiency of feedback systems. For perturbations up to frequencies corresponding to the width of the initial spike (of the order of 100 kHz for revolution periods of the order of 1 Ps), half of the beam behaves coherently with a definite delay while the other half is extended over a period that will cause overlap with ripple frequencies as low as 1 kHz. 5.7 Time profile without ripple Up to now, only the microscopic behaviour has been considered. When performing an extraction, the whole beam is present and the total spill is the important feature. In the following the aim is to derive the global behaviour from the elementary strip profile. In the Steinbach diagram, which shows the stable and unstable regions in the amplitude vs tune (or momentum) space, the strip corresponds to a fixed amplitude and tune (momentum), as shown in Figure 5.3.

Figure 5.3. A strip in the Steinbach diagram corresponds to a certain amplitude and momentum.

In general the beam may have any shape in the amplitude-momentum space, and the beam distribution shown in the figure is meant to represent a generic distribution. The particular cases of wide and narrow beam, with respect to the resonance, will be analysed later. 'p/p Amplitude BEAM Strip Stable region Resonance region

SLIDE 88

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

87 The total time profile of the beam is obtained by summing the contributions of the different

strips. In a schematic way this is shown in Figure 5.4:

100 200 300 400 500 600 0.0025 0.005 0.0075 0.01 0.0125 0.015 0.0175 0.02

Figure 5.4. Strip profiles for different strips. The total spill is the sum of all the strips.

where Pstrip,1, Pstrip,2, Pstrip,3 represent the profiles in time for strips that become unstable at different times or with different amplitudes and so on. In principle they are all different and the extracted current is given by: I(t) = 6i Pstrip,i (t). The evaluation of I(t), can be greatly simplified when the total extraction time is much longer than the duration of a strip. In this case, indeed, the relative variation of H, and of the other parameters, during the strip extraction is small and thus all the terms Pstrip,i(t), relative to the same H, effectively contributing to the total extracted current are equal. Then, the sum of all these terms just gives the integral over one of them. Since the same argument applies to all the amplitudes, or modified tune distances H, involved then the extracted current equals the number of particles that become unstable as a function of time NT(t)! This means that: x if the extraction is performed in amplitude the time profile is the amplitude distribution, x if the extraction is performed in momentum the time profile is the momentum distribution, x combinations of amplitude and momentum yield the combined distribution. 5.8 Ripple As it has just been shown, under ideal extraction conditions, apart from an initial and a final transient, the intuitive result is found, and the effort of this complicated analysis seems

worthless. Unfortunately, ideal conditions are never met in reality and the fields in the

magnets, and thus all the parameters that depend on them, are not perfectly stable, but

scillate around their value. This effect is called ripple and in general many ripple

Pstrip,3 Pstrip,1 Pstrip,2

t dN/dt

SLIDE 89

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

88 frequencies with different amplitudes are present at once. In the following only one ripple frequency will be considered. The effect of a ripple in tune is to move the “V” of the resonance back and forth while a ripple in the resonant sextupole strength causes the slope of the “V” to increase and decrease, as shown in Figure 5.5. In the very simple model in which a particle is extracted as soon as it enters the resonance, the profile in time dN/dt of the extracted beam is given by the distribution of particles that become unstable per unit tune variation F = dN/dQ and by the relative motion between the resonance and the stack:

t

Q Q Q Q dt dQ dQ dN dt dN

R r

Z ZG F F cos

(5.19)

where a sinusoidal ripple has been assumed. If the product of the ripple frequency Z and ripple amplitude GQR equals or exceeds the stack velocity

Q , the extracted beam is

completely or over modulated, respectively, and there are instants during which no beam is extracted because the resonance is “escaping” faster then the beam.

Figure 5.5. Beam and resonance with ripples.

This approximation holds until the transit time of the particles is short compared to the period of the ripple. As soon as the transit time is no longer negligible, a more accurate analysis is needed. As it has already been shown, the transit time is not the same for all the particles, thus a mixed concept considering both the transit time and the number of particles having a certain fixed delay has to be developed. In practice, this means evaluating the time profile for the set of particles that become unstable in one step of the ejection process, taking into account all the possible amplitudes. The “width” of this profile is the quantity that has to be compared to the ripple period to have an indication of the overlapping of these elementary profiles. If these overlap, the continuity of the extracted beam will not be broken even if the product ZGQR exceeds Q0 . 'p/p Amplitude Resonance excitation ripple Tune ripple

SLIDE 90

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

89 5.9 Comparison of extraction methods To extract the beam, either the particles have to be progressively moved into the resonance stopband or the stopband must be moved onto the beam. This means that some physical quantities have to be changed. The varying parameter, the way in which it is changed and the initial geometry in the amplitude momentum space distinguish the extraction method. Of the many possibilities illustrated in Section 3.11, the three cases in which the average extracted momentum does not vary during extraction are considered: x Amplitude selection by moving the resonance: the beam is narrow in momentum and the tune of the particles is changed, by varying the focusing of the machine, see Figure 5.6 (a). x Amplitude-momentum selection by moving the beam [13]: the beam is wide in momentum and the tune changes, due to chromaticity, when accelerating the particles, see Figure 5.6 (b). x RF knock-out [20]: the chromaticity is zero (or quasi-zero), so that all particles have the same tune. The particles are made unstable by blowing up the beam with transverse RF excitation. Thus the resonance is reached in amplitude, see Figure 5.6 (c).

Figure 5.6. Comparison of the main extraction methods: (a) moving the resonance, (b) moving the beam, (c) increasing the particle amplitude.

In the amplitude selection case, the movement of the resonance maroons the large amplitude particles first in the unstable region. Since the momentum spread is small, the basic element of spill marooned in the elementary extraction process (one elementary movement of the resonance) is nothing more than the strip profile just described. So the “building block” of this particular spill starts with a very narrow peak containing half of the particles involved. Moreover in the initial phase the transit time is very short and sensitivity to the ripples is high. At the end of the spill the transit times are longer and the particles in the tail are distributed over a longer period. In the RF-knockout method the blow-up velocity is fixed by the spill length. When looked at in tune, the Steinbach diagram is similar to the amplitude selection case one. The beam 'p/p Amplitude 'p/p Amplitude 'p/p Amplitude (a) (b) (c)

SLIDE 91

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

90 in this case has null width and the resonance is reached in amplitude. A ripple in tune would move the beam right and left, which causes the resonance line to move up and down in the “'p/p diagram”. As the particles enter the resonance at high amplitudes, and thus large tune distances, the transit time is short and the sensitivity to ripples is high. In the amplitude-momentum selection case, particles of all amplitudes become unstable at the same moment. This leads to an enlarged leading peak, which can fill the time interval during which no beam enters the resonance in a much more efficient way. A more detailed study is presented to estimate the smoothing of the ripples with this method in the following sections. 5.10 Elementary "band" for a wide momentum spread Consider the amplitude-momentum selection case. The beam is wide in tune and a stationary situation develops in which all betatron amplitudes are extracted simultaneously. In Figure 5.7 an amplitude-momentum selection extraction is shown and the elementary band* of beam that becomes unstable in one step is highlighted.

Figure 5.7. Amplitude-momentum selection extraction. A complete band of particles becomes unstable in one step.

A hollow beam is considered to set a minimum H which can be considered large with respect to its variation in the transit time. If the number of particles that would normally be in the hole is small, then a good representation for the beam with no hole is found. This approximation is not so restricting since firstly the radial particle density across a beam rises from zero at the centre and secondly particles with sufficiently small amplitude will cross the resonance without being extracted. Let the band that is extracted contain particles with tune shifts between GQmin and GQmax. The time profile of this band, is given by considering the time profile (5.4) for each H value between HF = 6SGQmin, corresponding to TF = Tc,d(HF), and H0 = 6SGQmax, corresponding to T0 = Tc,d(H0), and summing at each instant all these contributions, as shown in Figure 5.8.

* Here the word “band” is used to indicate the series of strips corresponding to different momentum and emittance that become unstable

simultaneously.

Qres GQmax No particles with small amplitude (hollow beam) H0, T0 HF, TF GQmin Amplitude

SLIDE 92

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

91 The resultant spill can be written as,

F

F

d t P t P N d t P N t P N

tail spike strip strip strip band B H H H H

H H H H H H H ) , ( ) , ( ) ( ) , ( ) ( ) ( (5.20)

where NB is the total number of particles in the band and Nstrip(H)dH is the number of particles marooned with modified tune distance between H and H + dH.

100 200 300 400 500 600 0.0025 0.005 0.0075 0.01 0.0125 0.015 0.0175 0.02

Figure 5.8. The total flux of particles is obtained summing the “strip profiles” for all the amplitudes present.

5.11 Simplified model The exact integral is far too complicated to be evaluated analytically. So, to facilitate the task, each Pstrip(t) is approximated by one delta function plus a rectangle representing Pspike and Ptail respectively. The width of the rectangles, as mentioned in section 5.4, is equal to the time needed for the first particle to reach the ES, that is Tc,d, and hence the height is 1/Tc,d. Not all the Pstrip(t, H) give their contribution at any given instant. Most of them are zero, either because they have not yet “started” or because they have already “finished”. At time t, only the H for which t/2 < Tc,d(H) < t give a non zero contribution. Thus the part of the integral arising from Ptail, becomes

)

2 / ( ) ( , ,

) ( 1 2 ) ( ) (

t t d c T tail band B

d T N t P N

H H

H H H

. To be correct the integration limits should be max(H(t), HF), which corresponds to min(t, TF), and min(H(t/2), H0) which corresponds to max(t/2, T0).

Pstrip,3

H3

Pstrip,1

H1

Pstrip,2

H2 t dN/dt Only half of the particles are in the rectangle

SLIDE 93

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

92 After a change of the integration variable,

)

), 2 / ( min( ) ), ( max( ) , 2 / max( ) , min( ,

1 2 )) ( ( ) ( 1 2 ) (

H H H H

H H H H H

t t T t T t T d c T

F F

dT dT d T T N d T N . (5.21)

It is necessary to find H(T), that is to invert

H

H H O H

3

1 3 3 ln 3 1

,

n n T d

c

. To simplify the calculation, it is assumed that the ES is far away so that n >> 1, thus n/(n+3) | 1. Neglecting O as usual, gives

H H H H H H H

2

,

3 3 ln 3 1 3 1 3 ln 3 1

d

c

T

which can be written

H H H H H

3

3 ln 3 3 1 3 2 1

,

d

c

T

It is necessary to invert the equation Y = ln(X)/X. A simple fit to the inversion, valid to within 10% in the range 4 < X < 30 000, is X = -1.42 ln(Y)/Y.

0. 05 0. 1 0. 15 0. 2 0. 25 0. 3 0. 35

Figure 5.9. Fit to the inversion (Logarithmic scale).

SLIDE 94

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

93 Note that X = 10 corresponds to

2 2

10 5

|

H H

and values of X smaller than this should be

avoided to stay within the limits of the transit time formulae (

2

H H

measures the relative

variation of the stable region during the transit time). Thus

T T T T

3

2 1 ln 3 2 42 . 1 3 2 1 3 2 1 ln 3 3 42 . 1 H H H H H

(5.22)

from which

T

T dT d 3 2 1 ln 1 1 3 2 42 . 1

2

H H

.

(5.23)

It is now necessary to estimate NT(H)dH. This is the number of particles that are in the border of thickness dH of a triangle of surface

2 2

3 2 3 3 3 3

S

H H

. This border corresponds to the annulus of width dR around a circle of the same surface SR2 in the initial beam. Thus,

H

H H S U H U H H d S S N H d R N d N

B B T 2 2

3 3 8 ) 3 2 3 3 ( 3 3 )) ( ( ) ( (5.24)

where NB is the total number of particles in the band. 5.12 Uniform distribution in phase space Assuming a uniform phase space distribution in the initial beam, of value U0 = 1/Triangle Area =

2 2

3 4 3 H S

, yields for the integration of the tails from all the differential strips in the band,

)

, 2 / max( ) , min( ,

1 2 )) ( ( 1 ) (

T t T t T B tail band

F

dT dT d T T N N t P H H

)

, 2 / max( ) , min( 4 2 2 2

3 2 1 ln 3 2 1 ln 2 1 3 4 42 . 1 3 3 8

T t T t

F

dT T T T S H H U

)

, 2 / max( ) , min( 4 2 2 2

3 2 1 ln 3 2 1 ln 3 42 . 1 4

T t T t

F

dT T T T H H H

.

SLIDE 95

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

94 Using the standard formulae

1

2 1

1 1 1 ) ln( ) ln(

p p p

x p x p ax dx x ax

and

1

3 1 2 1 2 2

1 2 1 ) ln( 2 1 ) ln( ) ln(

p p p p

x p x p x x p x dx x x

, the integral can be evaluated. It results that

)

, 2 / max( ) , min( 3 2 2 2 ) , 2 / max( ) , min( 3 3 3 2 3 3 2 2 ,

27 1 3 2 1 ln 3 3 2 1 ln 9 3 42 . 1 4 27 2 9 3 2 1 ln 2 3 3 2 1 ln 9 1 3 3 2 1 ln 3 42 . 1 4 ) (

T t T t T t T t tail band

F F

T T T T T T T T T T T t P

u
H

H H H H H H

and finally, for times 2T0 < t < TF, that is apart from the initial and final parts
u
3

2 2 2 2 ,

27 7 ) 2 ln( 24 ) 2 ln( 72 3 2 1 ln 21 ) 2 ln( 144 3 2 1 ln 63 3 42 . 1 4 ) ( t t t t P

tail band

H H H

(5.25)

To correctly evaluate the initial and the final part of the beam, the correct integration limits have to be considered. In fact for times t < T0, no particle has reached the septum and Pband,tail = 0. For times T0 < T < 2T0, the integral has to be performed between T and T0 (there are no particles with H = H(T/2)). For the same reason for times TF < T < 2TF, the integral has to be performed between TF and T/2. For times greater than 2TF, Pband,tail = 0 again.

SLIDE 96

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

95 The contribution from the initial spikes of the elementary strips in the band, gives

.

3 2 1 ln 3 2 1 ln 3 42 . 1 4 2 )) ( ( ) (

3 2 2 2 ,

t t t dT dT d T t T N t P

T T T spike band

F

H

H H H G H

The graph of the spill shape for a particular case is drawn in Figure 5.10.

20 40 60 80 0. 00 2 0. 00 4 0. 00 6 0. 00 8 0. 01

Figure 5.10. Band profile for GQmax = 3.782 10-3, GQmin = 3.782 10-4,

Q =-5.673 10-8,

S = 27.7897 m-1/2. The two contributions and the sum are shown.

5.13 Gaussian distribution in phase space Beams are often assumed to have a gaussian distribution in phase space

x

E R x

e E R

2

2 1

S

U (5.26)

where Ex is the horizontal one-V emittance*. Assuming that the area of the biggest stable triangle corresponds to an emittance of n-V, the one-V emittance is

2 2 2 2

3 3 4 n S n Area Ex S H S (5.27)

Thus

2 2 2

2 2 2 2

3 8 3

H H

H H U

n

e n S

(5.28)

* The definition of emittance used is phase space area divided by S

Sum of the leading peaks Sum of the tails Total flux Pband(t) t

SLIDE 97

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

96

H H H H H

H H

d e n N d N

n T T

2 2 2

(5.29)

and

u
)

, 2 / max( ) , min( 4 2 3 3 2 1 ln 2 42 . 1 ) , 2 / max( ) , min( 2 2 2 ,

2 2 2 2 2

. 3 2 1 ln 3 2 1 ln 3 42 . 1 2 1 2 )) ( ( ) (

T t T t T T n T t T t T tail band

F F

dT T T T e n dT dT d T T N t P H H H H H

H H

(5.30)

The integral can be evaluated numerically and plotted, as shown in Figure 5.11.

20 40 60 80 0. 00 1 0. 00 2 0. 00 3 0. 00 4 0. 00 5

Figure 5.11. Band profile for a gaussian beam for GQmax = 3.782 10-3, GQmin = 3.782 10-4,

Q =-

5.673 10-8, S = 27.7897 m-1/2, Etot = 9 Erms. The two contributions and the sum are shown.

If the gaussian and the uniform cases are plotted on the same graph, the comparison in Figure 5.12 is obtained

20 40 60 80 0. 00 25 0. 00 5 0. 00 75 0. 01 0. 01 25 0. 01 5

Figure 5.12. Band profile for gaussian and uniform beam. A 3V emittance of the gaussian beam has been considered as total emittance.

Sum of the leading peaks Sum of the tails Total current Pband t Pband t

SLIDE 98

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

97 From the plots above, it appears that the width of the “band profile” is of the order of T0. Thus it can be expected that when the ripple period is of the order of T0, the width of the profile partially fills the time interval during which no beam enters the resonance and smoothes the modulation of the extracted beam. The band profiles have been worked out in this Chapter considering a constant Q . Thus considering them as “filling the gap” in the presence of ripple is an approximation which neglects the instantaneous variation and uses an average Q ( Q = Q 0 + Z GQR cos(Zt) | Q 0). This can be accepted by considering that

Q appears in the expression for the transit time inside a logarithm and visually by

comparing the profiles in time for a variation in Q of 2 orders of magnitude, as shown in figure 5.13:

20 40 60 80 0. 00 2 0. 00 4 0. 00 6 0. 00 8 0. 01 0. 01 2 20 40 60 80 0. 00 5 0. 01 0. 01 5 0. 02

Figure 5.13. Time profiles for gaussian and uniform beams for different values of

Q . Even a large variation in Q does not vary the orders of magnitude and, to a first approximation, the band profile can be considered when analysing the spill with ripple. The enlarged leading peak quoted in Section 5.9 can be thus estimated in hundreds of turns to be compared with the tens of turns of the strip profile, demonstrating the superior ripple tolerance of the amplitude-momentum selection scheme. 5.14 Simulations of extraction with ripple All the theory and approximations made so far, call for a check to verify the correctness of the conclusions drawn. With this in mind, three numerical simulations have been made: 1) amplitude selection extraction with a mono-energetic beam; 2) amplitude-momentum selection with a uniform distribution in phase space; 3) amplitude-momentum selection with a gaussian distribution in phase space. In order to be consistent with the numerical examples shown so far, the following parameters have been used for the simulations: Trev = 0.734 Ps

Q = -5.673 10-8 =>

Q = -0.0257 s-1. It should be noted that these simulation parameters correspond to the extraction of protons at 60 MeV from the PIMMS ring and are therefore representative of a practical situation. t t

Q = -5.673 10-9

.

Q = -5.673 10-7

. Pband Pband

SLIDE 99

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

98 From inspection of Figure 5.12, it results that the width of the “band profile” should be able to fill gaps of the order of 'T = 100. To demonstrate this effect a round number for the tune ripple frequency has been chosen as fripple = 4000 Hz => 'T = 113. The corresponding amplitude for 100% modulation of the extracted beam in the “instantaneous transfer” model is Aripple = 10-6. The tune ripple is excited by a quadrupole gradient ripple on a single quadrupole family. The corresponding normalised quadrupole gradient ripple amplitude is easily obtained with a lattice program. It should be noticed that the relative ripple amplitude, in this case, is just 1.3 10-6, which is a very tough constraint on the quality of the power supplies! Since the binning in the analysis of the results cannot be infinitely large, a slightly larger (20% larger) amplitude has been used in the simulations. The program used to simulate the extraction is MAD [21]. Note that the particle tracking does not follow the approximated hamiltonian used in the theoretical derivation, but is an exact tracking through the PIMMS synchrotron. In order to limit the computing time and to have a sufficiently large number of extracted particles, only a small fraction of a beam is generated which occupies the region of interest. This is schematically shown in Figure 5.14.

Figure 5.14. The initial distributions are generated only in the region of interest.

The resonant sextupole strength is raised linearly from zero to its final value in 4000 turns. The beams are generated such that a small part of them is unstable once the sextupole is fully raised. In order to extract those particles that are already unstable, after the sextupole has reached its full strength, the program lets the particles circulate for 6000 turns before extraction starts. As an example the main steps in the preparation of the beam are shown in Figures 5.15 to 5.17 where the initial distribution, the distribution after 4000 turns (sextupole reaches its final value) and 10000 turns (just before extraction starts) are shown for the uniform, band-like distribution. The time profile is analysed using bins 30 turns wide, which correspond to an integration time of 22 Ps for an hypothetical extracted current measurement.

“band-like” beam monochromatic beam

SLIDE 100

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

99

Figure 5.15. Initial beam distribution for the uniform band case. Figure 5.16. Beam distribution after raising the sextupole. Figure 5.17. Beam distribution just before extraction starts.

SLIDE 101

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

100 5.14.1 Amplitude selection As is evident from inspection of Figure 5.18, in this case the extracted beam is 100%

modulated. This is to be expected since the basic element for building the time profile is

the strip profile, whose width is small with respect to the 4 kHz period.

10 20 30 40 50 60 70 80 4020 4260 4500 4740 4980 5220 5460 5700 5940 6180 6420 6660 6900 7140 7380 7620 7860 20 40 60 80 100 120 4020 4260 4500 4740 4980 5220 5460 5700 5940 6180 6420 6660 6900 7140 7380 7620 7860

Figure 5.18. Time profile without (upper) and with (lower) ripple.

SLIDE 102

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

101 5.14.2 Amplitude-momentum selection (uniform distribution) For the “uniform phase-space distribution” – “band profile”, the basic element for the time profile construction is the band profile shown in Figure 5.10. It is clear that 113 time units after the band profile start, there is still some beam in the resonance that fills the “gap”. The simulation, shown in Figure 5.19, confirms this intuition, showing that the continuity in the extracted flux is not broken, although the modulation is visible.

10 20 30 40 50 60 2010 2250 2490 2730 2970 3210 3450 3690 3930 4170 4410 4650 4890 5130 5370 5610 5850 10 20 30 40 50 60 70 80 2010 2250 2490 2730 2970 3210 3450 3690 3930 4170 4410 4650 4890 5130 5370 5610 5850

Figure 5.19. Time profile without (upper) and with (lower) ripple.

SLIDE 103

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

102 5.14.3 Amplitude-momentum selection (gaussian distribution) For the “gaussian phase-space distribution” – “band profile”, the basic element for the time profile construction is the band profile shown in Figure 5.11. By direct inspection, it is clear that 113 time units after the band profile start, the particle flux is still close to maximum and, consequently, a very good smoothing of the ripple is expected. The simulation of Figure 5.20 shows that in this case the extracted spill is almost indistinguishable with or without ripple.

10 20 30 40 50 60 4020 4260 4500 4740 4980 5220 5460 5700 5940 6180 6420 6660 6900 7140 7380 7620 7860 10 20 30 40 50 60 70 4020 4260 4500 4740 4980 5220 5460 5700 5940 6180 6420 6660 6900 7140 7380 7620 7860

Figure 5.20. Time profile without (upper) and with (lower) ripple.

SLIDE 104

Time profile of the extracted beam Chapter 5

________________________________________________________________________________________________________________________________________________________________________________________________________

103 5.14.4 Amplitude-momentum selection (gaussian distribution) II In order to have a second check for the “gaussian phase-space distribution” – “band profile”, the same simulation has been run with a 10 times larger ripple amplitude. The simulation of Figure 5.21 shows that also in this case the extracted spill is very well behaved.

10 20 30 40 50 60 4020 4260 4500 4740 4980 5220 5460 5700 5940 6180 6420 6660 6900 7140 7380 7620 7860 10 20 30 40 50 60 70 4020 4260 4500 4740 4980 5220 5460 5700 5940 6180 6420 6660 6900 7140 7380 7620 7860

Figure 5.21. Time profile without (upper) and with (lower) ripple.

SLIDE 105

Chapter 5 Time profile of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

104 5.15 Concluding remarks At kHz frequencies and above the amplitude-momentum selection extraction has been shown to offer a very useful smoothing effect on the spill. This allows a relaxation of the tolerances on the power supplies in a region where the constraints are very tough. It should be noticed that the time profile is expressed in number of turns and the revolution period varies with machine size and extraction energy. Full suppression of ripple effects starts for ripple periods of the order of T0 = Tc,d(dQmax), that is the transit time for the largest amplitude particles. For PIMMS, a typical proton-carbon ion medical synchrotron, this is of the order of 4 kHz for the protons at the minimum extraction energy. The smoothing diminishes to near zero at around 1 kHz. At higher extraction energies, the revolution period diminishes and the full smoothing threshold moves up in frequency, which reduces the efficiency at a given frequency. However, at higher extraction energies, the magnet power supplies work at a higher current and the relative ripple is lower. Fortuitously, this gain outweighs the loss in “band smoothing” due to the higher revolution frequency. The suppression of ripple in the spill is of paramount importance for raster scanning. Other smoothing techniques are available: x front-end acceleration by RF noise or empty bucket channelling (efficient up to 2 kHz) x compensation via scanning speed (efficient up to 1 kHz) x feedback (efficient up to 1 kHz) x geometric smoothing from beam size (efficient at 10 kHz and above) x eddy currents smoothing (efficiency increases from 1 kHz upwards) However the intrinsic smoothing from the band profile is embedded in the extraction design and is, as such, totally reliable not requiring any active components. Moreover it covers the difficult region of the few kHz.

SLIDE 106

105

Chapter 6 TRANSVERSE ASPECTS OF THE EXTRACTED BEAM

For a passive beam delivery system, which scatters the beam into a much larger emittance than the incoming beam and then cuts away everything except the central region, the beam distribution at the entrance to the scatterer is likely to be of little concern. On the contrary, when active scanning is considered, the extracted beam distribution is directly delivered to the patient. A good understanding of the transverse distribution from a slow extraction is thus of prime importance for treatment planning and for treatment quality and is the subject of the present Chapter. 6.1 Simplest model In Chapter 3, the resonant extraction has been depicted using the stable triangle and the

separatrices. The situation at any position in the machine, can be obtained by rotating the

pattern at the resonant sextupole clockwise by the phase advance to the considered

position. At the electrostatic septum, when the Hardt condition is fulfilled, the horizontal

phase space is as shown in Figure 6.1.

Figure 6.1 Stable region and extraction separatrices at the electrostatic septum (ES)

The extracted beam is the part of the separatrices which is cut by the electrostatic septum. In this very simplified model, its emittance is zero, since the separatrix has no “thickness”. An already better approximation in this model would be considering a “thickness” for the separatrix equal to the unit separatrix movement Oh. In the case of amplitude selection this would be approximately correct, since the extraction separatrix is effectively moving, while for momentum-amplitude selection this would be an upper limit. With the numbers used in the previous Chapters and assuming a spiral step of 0.01 m, EES = 16 m and an angle between the separatrix and the X axis of 45º, this simple model would yield:

X Xc ES

SLIDE 107

Chapter 6 Transverse aspects of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

106 (thickness) Oh = 15 10-6 · 1.7 10-3 m1/2 = 2.55 10-8 m1/2 and (length) 2 ·0.01 m/ 16 m1/2 = 3.54 10-3 m1/2 that is an emittance of Ex = 2.9 10-5 S mm mrad. Alternatively, the Liouville theorem can be invoked and the emittance of the extracted beam can be estimated by the circulating beam emittance, the spill length and the extracted momentum spread as:

spill rev ext circ circ x

T T p p p p E E ) / ( ) / ( ' '

which is of the same order of magnitude as before. Numbers of this order of magnitude are unrealistic and a deeper study is thus needed. 6.2 Equations of motion including the vertical plane The horizontal emittance of the extracted beam can arise from many sources such as ripples “shaking” the separatrix, non-adiabaticity of the process and other effects neglected in the approximations. The first striking quantity which has been neglected and which might give rise to an emittance is the coupling of the vertical motion into the horizontal plane via the strong resonance sextupole. As was shown in Chapter 2, the field seen by a particle when it traverses a sextupole depends both on x and on z. The kick in normalised phase space was shown in (2.12) and (2.14), which are reported hereafter for simplicity:

c

'

2 2

Z X S X

x z

E E XZ S Z

x z

E E 2

c

'

. Using these expressions together with the single-turn transfer matrix and with the help of a powerful calculation tool such as Mathematica [22] it is possible to evaluate the exact transfer function for three turns. This yields a very large expression, whose size renders it

useless. As an example the term 'X3 is shown in Figure 6.2. Similar expressions can be

found for 'Xc3, 'Z3 and 'Zc3.

SLIDE 108

Transverse aspects of the extracted beam Chapter 6

________________________________________________________________________________________________________________________________________________________________________________________________________

107

Out[5]= X0 X0Cos

#2dQxS 2S Qx '3 3XP0Cos #2dQxS 2S Qx'2Sin #2dQxS 2S Qx ' 2SX02Cos#2dQxS 2S Qx '3Sin#2dQxS 2S Qx'

SX02Cos#2dQxS 2S Qx

'4Sin #2dQxS 2S Qx' 2BRSZ0 2 Cos #2dQxS 2S Qx' Cos #2S Qz'2Sin #2dQxS 2 S Qx '

BRSZ02Cos#2 S Qz

'4Sin#2dQxS 2S Qx ' 3X0Cos #2dQxS 2S Qx ' Sin #2dQxS 2S Qx'2

4SX0XP0Cos

#2dQxS 2S Qx '2 Sin #2dQxS 2S Qx'2 4SX0XP0Cos#2dQxS 2S Qx '3 Sin #2dQxS 2S Qx'2

2S2X0

3 Cos #2dQxS 2S Qx'4 Sin #2dQxS 2 S Qx '2 2BRS2X0Z0 2Cos #2dQxS 2 S Qx '2Cos #2S Qz'2Sin#2dQxS 2S Qx '2

XP0Sin

#2dQxS 2S Qx'3 2SXP0 2 Cos #2dQxS 2S Qx' Sin #2dQxS 2 S Qx '3 2SX02Cos #2dQxS 2 S Qx '2Sin#2dQxS 2S Qx'3

4SXP02Cos#2dQxS 2S Qx

'2Sin #2dQxS 2S Qx'3 8S2X0 2 XP0Cos#2dQxS 2S Qx '3Sin #2dQxS 2S Qx'3

S3X04Cos#2dQxS 2S Qx

'4 Sin #2dQxS 2S Qx'3 4BRS 2 XP0Z02Cos #2dQxS 2S Qx' Cos #2S Qz'2 Sin #2dQxS 2S Qx'3

2BRS3X0

2Z02Cos #2dQxS 2S Qx'2Cos #2S Qz'2 Sin #2dQxS 2 S Qx '3 BR 2 S 3Z0 4 Cos #2S Qz'4Sin #2dQxS 2S Qx'3

4SX0XP0Cos

#2dQxS 2S Qx ' Sin #2dQxS 2S Qx'4 2S2X0 3 Cos #2dQxS 2S Qx'2Sin #2dQxS 2S Qx '4

10S2X0XP02Cos

#2dQxS 2 S Qx '2Sin#2dQxS 2S Qx'4 4S3X03XP0Cos#2dQxS 2S Qx '3Sin#2dQxS 2S Qx'4

2BRS2X0Z0

2Cos#2 S Qz '2Sin #2dQxS 2S Qx'4 4BRS3X0XP0Z02Cos#2dQxS 2S Qx' Cos #2S Qz'2 Sin #2dQxS 2S Qx'4

SX02Sin#2dQxS 2S Qx

'5 4S2X02XP0Cos #2dQxS 2S Qx ' Sin #2dQxS 2S Qx'5 4S2 XP0 3Cos#2dQxS 2S Qx ' Sin #2dQxS 2S Qx'5

6S3X0

2 XP0 2Cos#2dQxS 2S Qx '2Sin #2dQxS 2S Qx'5 2BRS3XP0 2Z02Cos #2S Qz '2 Sin #2dQxS 2S Qx'5

2S2X0XP0

2Sin#2dQxS 2S Qx '6 4S3X0XP03Cos #2dQxS 2 S Qx ' Sin#2dQxS 2S Qx '6 S3XP04Sin #2dQxS 2 S Qx '7

4BRSZ0ZP0Cos

#2dQxS 2S Qx ' Cos #2S Qz ' Sin #2dQxS 2S Qx' Sin #2S Qz' 4BRSZ0ZP0Cos#2 S Qz '3Sin#2dQxS 2S Qx ' Sin #2S Qz '

4BR2S2 X0Z02Cos

#2dQxS 2S Qx' Cos #2S Qz'3Sin #2dQxS 2S Qx' Sin #2S Qz' 4BRS2X0Z0ZP0

Cos

#2dQxS 2S Qx'2Cos #2S Qz' Sin #2dQxS 2S Qx '2 Sin #2S Qz ' 4BR2S2XP0Z02Cos #2S Qz'3Sin #2dQxS 2 S Qx '2Sin#2 S Qz '

8BRS2XP0Z0ZP0Cos

#2dQxS 2S Qx ' Cos#2 S Qz ' Sin#2dQxS 2S Qx '3 Sin #2S Qz '

4BRS3X0

2Z0ZP0Cos #2dQxS 2S Qx '2 Cos #2S Qz ' Sin #2dQxS 2S Qx'3Sin #2S Qz'

4BR2S3 Z0

3ZP0Cos#2 S Qz '3Sin#2dQxS 2S Qx '3 Sin #2S Qz ' 4BRS2X0Z0ZP0Cos #2S Qz' Sin #2dQxS 2S Qx '4Sin#2 S Qz ' 8BRS3X0XP0Z0ZP0

Cos

#2dQxS 2S Qx' Cos #2S Qz ' Sin#2dQxS 2 S Qx '4Sin #2S Qz ' 4BRS3XP02Z0ZP0Cos #2S Qz' Sin #2dQxS 2S Qx'5Sin #2S Qz'

2BRSZP02Cos#2dQxS 2S Qx

' Sin #2dQxS 2S Qx' Sin #2S Qz'2 2BRSZ0 2 Cos #2S Qz '2Sin #2dQxS 2S Qx' Sin #2S Qz'2

4BRSZP02Cos#2 S Qz

'2Sin#2dQxS 2S Qx ' Sin #2S Qz '2 12BR2S2X0Z0ZP0Cos #2dQxS 2S Qx' Cos #2S Qz'2Sin#2dQxS 2S Qx ' Sin #2S Qz '2

4BR3S3 X0

2Z0 2 Cos #2dQxS 2S Qx'2 Cos #2S Qz'2Sin #2dQxS 2S Qx ' Sin#2S Qz '2

2BRS2X0ZP0

2Cos #2dQxS 2S Qx '2 Sin #2dQxS 2S Qx'2Sin #2S Qz'2 12BR2S2 XP0Z0ZP0Cos #2S Qz'2Sin #2dQxS 2 S Qx '2Sin#2S Qz '2

8BR3S3 X0XP0Z02Cos#2dQxS 2S Qx

' Cos #2S Qz '2Sin #2dQxS 2S Qx'2Sin #2S Qz'2

4BRS2XP0ZP02Cos

#2dQxS 2 S Qx ' Sin#2dQxS 2S Qx '3Sin #2S Qz'2

2BRS3X0

2ZP02Cos#2dQxS 2S Qx '2Sin #2dQxS 2S Qx'3Sin #2S Qz'2 4BR 3 S 3XP0 2 Z0 2Cos #2S Qz '2 Sin #2dQxS 2S Qx'3Sin #2S Qz'2

6BR2S3 Z0

2ZP0 2Cos#2 S Qz '2Sin#2dQxS 2S Qx '3Sin #2S Qz'2 2BRS2X0ZP02Sin#2dQxS 2S Qx '4Sin #2S Qz '2

4BRS3X0XP0ZP02Cos

#2dQxS 2 S Qx ' Sin#2dQxS 2S Qx '4Sin #2S Qz'2 2BRS3XP02ZP02Sin #2dQxS 2S Qx'5Sin#2 S Qz '2

4BRSZ0ZP0Cos

#2S Qz' Sin #2dQxS 2S Qx ' Sin #2S Qz '3 4BR2S2X0Z02Cos#2dQxS 2S Qx' Cos #2S Qz' Sin #2dQxS 2S Qx' Sin #2S Qz'3

8BR2S2 X0ZP02Cos#2dQxS 2S Qx

' Cos #2S Qz ' Sin #2dQxS 2S Qx' Sin #2S Qz'3 8BR 3 S 3X0 2 Z0ZP0

Cos

#2dQxS 2S Qx'2Cos #2S Qz' Sin #2dQxS 2S Qx ' Sin #2S Qz '3 4BR2S2XP0Z02Cos #2S Qz' Sin #2dQxS 2 S Qx '2Sin#2 S Qz '3

8BR2S2 XP0ZP02Cos

#2S Qz' Sin #2dQxS 2 S Qx '2Sin#2 S Qz '3 16BR 3 S 3X0XP0Z0ZP0Cos#2dQxS 2S Qx '

Cos

#2S Qz' Sin #2dQxS 2S Qx '2Sin #2S Qz'3 8BR 3S3XP0 2Z0ZP0Cos #2S Qz ' Sin #2dQxS 2S Qx'3Sin #2S Qz'3

4BR2S3 Z0ZP03Cos#2 S Qz

' Sin#2dQxS 2S Qx '3Sin #2S Qz'3 BRSZ02Sin #2dQxS 2S Qx ' Sin #2S Qz '4

4BR2S2 X0Z0ZP0Cos

#2dQxS 2S Qx ' Sin #2dQxS 2S Qx' Sin #2S Qz'4 4BR3S 3 X0 2ZP0 2Cos#2dQxS 2S Qx '2Sin #2dQxS 2S Qx' Sin #2S Qz'4

4BR2S2 XP0Z0ZP0Sin

#2dQxS 2 S Qx '2Sin#2 S Qz '4 8BR 3S3X0XP0ZP0 2Cos #2dQxS 2S Qx' Sin #2dQxS 2S Qx'2Sin#2 S Qz '4

4BR3S3 XP0

2ZP02Sin #2dQxS 2 S Qx '3Sin#2 S Qz '4 BR 2S3ZP0 4Sin #2dQxS 2S Qx'3Sin #2S Qz'4

Figure 6.2 Exact calculation of the displacement 'X3 in three turns.

In order to compare with the Kobayashi equations, only terms up to second order in H, X0, Xc0, Z0 and Zc0 have to be kept. Thus it is necessary to expand the expression for 'X3 in power series in dQx up to second order and then neglect the higher order terms. This results in

> @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @

2

2 2 2 2 3 3 2 4 2 2 2 2 4 2 2 3

' 2 sin 2 cos 3 2 ' 2 sin 2 3 ' 2 sin 2 cos 3 2 ' 2 sin 2 cos 3 2 ' 2 sin 2 cos 3 2 sin 2 3 2 sin 2 cos 3 2 cos 2 3 2 cos 2 3 ' 2 3 ' Z Q Q S Z Q S Z Z Q Q S Z Z Q Q S Z Z Q Q S Z Q S Z Q Q S Z Q S Z Q S X SX X X

z z z z z z z z z z z z z z x z

S S S S S S S S S S S S S S E E H

'

(6.1)

for tunes Qx = n+1/3 and

SLIDE 109

Chapter 6 Transverse aspects of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

108

> @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @

2

2 2 2 2 3 3 2 4 2 2 2 2 4 2 2 3

' 2 sin 2 cos 3 2 ' 2 sin 2 3 ' 2 sin 2 cos 3 2 ' 2 sin 2 cos 3 2 ' 2 sin 2 cos 3 2 sin 2 3 2 sin 2 cos 3 2 cos 2 3 2 cos 2 3 ' 2 3 ' Z Q Q S Z Q S Z Z Q Q S Z Z Q Q S Z Z Q Q S Z Q S Z Q Q S Z Q S Z Q S X SX X X

z z z z z z z z z z z z z z x z

S S S S S S S S S S S S S S E E H

'

(6.2)

for tunes Qx = n-1/3. In both cases the Kobayashi equation is clearly recognisable in the first two terms. All the other terms depend only on Z and Zc and only an overall sign changes for the two families nr1/3. In the following Qx = n-1/3 will be considered, the results for Qx = n+1/3 being easily derivable. Proceeding in the same way as for 'X3 an expression can be found for 'Xc3. This results in a single expression for both Qx = n-1/3 and Qx = n+1/3,

>

@ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @

2

6 2 4 2 2 2 4 2 2 2 2 2 5 3 3 3 5 3 2 4 2 2 4 2 2 4 2 2 2 2 6 2 4 2 2 2 2 3

' 2 sin ' 2 sin 2 cos 6 ' 2 sin 2 cos 9 ' 2 sin 2 cos 2 ' 2 sin 2 1 ' 2 sin 2 cos 6 ' 2 sin 2 cos 20 ' 2 sin 2 cos 2 ' 2 sin 2 cos 6 ' 2 sin 2 cos 2 ' 2 sin 2 cos 2 sin 2 cos 9 2 sin 2 1 2 sin 2 cos 6 2 sin 2 cos 2 cos 2 cos 2 1 2 cos 2 1 ' 4 3 ' Z Q S Z Q Q S Z Q Q S Z Q Q S Z Q S Z Z Q Q S Z Z Q Q S Z Z Q Q S Z Z Q Q S Z Z Q Q S Z Z Q Q S Z Q Q S Z Q S Z Q Q S Z Q Q S Z Q S Z Q S Z Q S X X S X X

z z z z z z z z z z z z z z z z z z z z z z z z z z z z z z x z

S S S S S S S S S S S S S S S S S S S S S S S S S S S S S S E E H

'

(6.3)

Again the first three terms correspond exactly to the Kobayashi equation while the rest is independent of X0, Xc0 and H. Collecting terms and using standard trigonometric formulae to simplify the resulting coefficients, expression 6.2 can be rewritten as

SLIDE 110

Transverse aspects of the extracted beam Chapter 6

________________________________________________________________________________________________________________________________________________________________________________________________________

109

> @ > @ > @

>

@

'

2 2 2 2 3

' 8 cos 4 3 ' 4 cos 4 3 ' 8 sin 2 3 ' 4 sin 2 3 ' 2 3 ' Z Z Q S Z Z Q S Z Z Q S Z Z Q S X SX X X

z z z z x z

S S S S E E H

A better understanding can be achieved by expressing Z and Zc in terms of single particle emittance and angle in normalized phase space, Z0 = Hz cos(Iz) and Z0c = Hz sin(Iz). This results in

> @ > @ > @ > @ > @ > @ > @ > @ > @

>

@ > @ > @

'

2 2 2 2 3

sin cos 8 cos 4 3 sin cos 4 cos 4 3 sin cos 8 sin 2 3 sin cos 4 sin 2 3 ' 2 3 '

z z z z z z z z z z z z z z z z x z

Q SE Q SE Q SE Q SE X SX X X I I S I I S I I S I I S E E H

A further simplification, using standard trigonometric formulae, is then possible for 'X3 which gives

> @ >

@

z z z z x z

Q Q S E X SX X X S I S E E H 6 2 sin 2 sin 2 3 ' 2 3 '

, 3

'

. (6.4)

Applying the same procedure to the expression for 'Xc3, results in

>

@

> @ > @

.

6 2 sin 10 2 sin 2 2 sin 2 1 ' 4 3 6 '

, , 2 2 3 z z z z z z x z x

Q Q Q S E X X S X Q X S I S I S E E SG

'

Using the standard trigonometric formulae

2 sin[a-b] + sin[a-c] = (2 cos[b] + cos[c]) sin[a] - (2 sin[b] + sin[c]) cos[a] = = - {(2 cos[b] + cos[c])2 + (2 sin[b] + sin[c])2}1/2 cos[a+k] = = - {5 + 4 cos[b-c]}1/2 cos[a+k]

where

cos[k]= (2 sin[b] + sin[c]) / {(2 cos[b] + cos[c])2 + (2 sin[b] + sin[c])2}1/2 sin[k]= (2 cos[b] + cos[c]) / {(2 cos[b] + cos[c])2 + (2 sin[b] + sin[c])2}1/2

finally yields

SLIDE 111

Chapter 6 Transverse aspects of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

110

k

Q Q S E X X S X X

z z z z x z

'

, 2 / 1 2 2 3

2 Cos 4 Cos 4 5 2 Sin 2 1 ' 4 3 ' I S S E E H

(6.5)

Equations (6.4) and (6.5) model the second-order description of the horizontal motion of a particle at the sextupole. For null initial amplitude in z, Ez = 0, or for Ez/Ex << 1 they coincide with the Kobayashi equations. When Ez cannot be neglected, the effect of the vertical plane is shown. 6.3 Separatrix width Consider a particle starting with X0 = -2H/3S. If the Z plane is neglected (or the amplitude in Z is null), as in the Kobayashi model, the equations of motion give

2

2 2 2 3 3

' 4 3 ' 4 3 ' ' 2 3 ' SX S X X S X X X SX X X

'
'

H H H

which shows that a particle starting on the separatrix X = -2H/3S, see Figure 6.3, will stay

n the separatrix forever, that is until it reaches the electrostatic septum. This is equivalent

to saying that the “separatrix width” is null.

Figure 6.3 In the Kobayashi model, the separatrix width is null.

Consider now equations (6.4) and (6.5). The influence of the Z plane is considered at the same order as for the X plane. This time a particle starting with X0 = -2H/3S no longer stays on the separatrix as an additional displacement in X is introduced by the coupling. After three turns, a particle starting with X0 = -2H/3S, Z0 = Ez cos(Iz,0), Zc0 = Ez sin(Iz,0) will be “out” of the separatrix by

> @ >

@

z z z z x z

Q Q S E S I S E E 6 2 sin 2 sin 2 3

,

Initial position After 3 turns

SLIDE 112

Transverse aspects of the extracted beam Chapter 6

________________________________________________________________________________________________________________________________________________________________________________________________________

111 After 3n turns the particle position in the Z plane is given by the angle Iz,n = Iz,0 - n 6SQz. And the additional displacement in the three following turns is

> @

>

@

z z z z x z

Q n Q S E S I S E E 6 1 2 2 sin 2 sin 2 3

,

.

The total displacement after 3N turns is thus

> @

>

@

1

,

6 1 2 2 sin 2 sin 2 3

N n z z z z x z

Q n Q S E S I S E E

. The sum is a standard form

> @

2

sin 1 2 sin 2 1 sin sin

1

y y N y N x ny x

N n

and yields

> @

> @ > @

> @

> @ > @ .

6 sin 6 sin 6 2 sin 6 sin 6 sin 12 2 1 6 2 sin 12 6 2 sin

, , 1 , z z z z z z z z z N n z z z

Q Q N Q N Q Q N Q N Q Q n Q S S S I S S S S I S S I

Since sin[2Iz,0 - N 6SQz] sin[N 6SQz] varies between -1 and +1, the amplitude of
scillation of the particles around the separatrix is proportional to cosec[6SQz]. When

Qz -> n/6, n integer, the oscillation amplitude grows linearly with the number of turns needed to reach the ES. Including the coefficient, the amplitude of oscillation of the particles around the separatrix is

> @ > @

z z z x z

Q Q S E S S E E 6 sin 2 sin 2 3 . (6.6)

In the derivation of this amplitude, the “recall force” towards the separatrix (the trajectories in the Kobayashi approximation approach asymptotically the separatrix) has been neglected. On the other hand, if a “null width” separatrix is considered at the last position before the electrostatic septum, the particles will spread in a band of amplitude given by the last addendum in (6.4) that is, dropping the term sin(2Iz,0 - 6SQz) which ranges between -1 and +1,

> @

z z x z

Q S E S E E 2 sin 2 3 . (6.7)

Expressions (6.6) and (6.7) give an upper and a lower limit to the half-width of the

separatrix. The distance of the trajectory to the separatrix as foreseen by the Kobayashi

(Upper Limit) (Lower Limit)

SLIDE 113

Chapter 6 Transverse aspects of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

112 approximation should be added, but in a slow extraction this is generally negligible. When the vertical tune is far from n/6, the difference between the upper and lower limit is small: for a tune Qz = 1.72, for example (PIMMS), 1/sin(6SQz) = 1.18, that is the upper limit is

nly 20% larger than the lower limit

6.4 Numeric check The formulae derived for the separatrix width can be checked by simulating the extraction with MAD. Particles with a given unique amplitude in x and z as illustrated in Figure 6.4, are generated. Once the sextupole has been raised, 3000 turns are added for the unstable particles to be extracted and finally the remaining particles are accelerated into the resonance.

Figure 6.4 Particle distributions at ES during extraction simulation: (a) initial distribution in (x, xc), (b) initial distribution in (z, zc), (c) after raising the sextupole in (x, xc), (d) after raising the sextupole in (z, zc), (e) extracted beam (x, xc), (f) extracted beam (z, zc). (a) (b) (c) (d) (e) (f)

SLIDE 114

Transverse aspects of the extracted beam Chapter 6

________________________________________________________________________________________________________________________________________________________________________________________________________

113 In Table 6.1 the calculated, expression 6.7, and simulated separatrix “thicknesses” are reported for Ex = 5.7 S mm mrad (2 sigmas for PIMMS) and for Ez ranging 0 to 6 in steps

f 1 S mm mrad.

Separatrix “thickness” Ez [S mm mrad] Simulation [P radians] Theory [P radians] 4 1 7 5 2 10 9 3 14 14 4 17 19 5 20 23 6 24 28

Table 6.1 Comparison of numeric and analytic results for the “thickness” of the separatrix

The analytic expression and the results from simulations agree reasonably well. Other unforeseen effects are present which are visible at small vertical amplitudes giving a finite “thickness” where the theory would foresee zero effect. The theoretical “thickness” increase is larger than the simulations, but the order of magnitude is correct. These differences are caused by the higher-order terms which have been neglected in the equations of motion (6.4) and (6.5). Since expanding the equations to third order in H, X0, Xc0, Z0 and Zc0 leads to very complicated, and thus useless, expressions and since the correct orders of magnitude are given by the simple expressions derived above, no further effort is made in this thesis to refine the formulae. 6.5 Separatrix displacement When analysing the results of a real tracking, the first characteristic that is observed is that the separatrices are not straight lines, as they have been depicted, but they are curved. This, of course, depends on the higher order terms in x which have been neglected in the analytic approach. This in turn partially invalidates many of the results reported in the general part, such as the expression for the Hardt condition for example. This is, however, a minor concern because, when designing the extraction, the equations and the results of Chapters 1 to 3 are used only as a guideline and the Hardt condition, for example, is set exactly by tracking the separatrices and then adjusting the parameters until the separatrices are overlapped at the electrostatic septum. In Figure 6.5 the extraction plot at the electrostatic septum is shown in un-normalised coordinates. The curvature of the separatrices is evident and the part of the separatrix cut by the electrostatic septum, that is the extracted beam, is clearly visible in the right hand part of the figure. The extracted beam looks denser than the separatrix because, to make it clearly visible, all the particles extracted up to the moment of the snapshot have been plotted (equivalent to a time integration). The second effect which is noticed is that the separatrices coresponding to different amplitudes in z are not centered on top of each other as would be expected from expression 6.7.

SLIDE 115

Chapter 6 Transverse aspects of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

114 Figure 6.6, shows the four extracted particle distributions when the initial amplitude in (x, z) is (0,0), (2,0), (0,2) and (2,2) sigmas. The two long lines correspond to zero amplitude in x while the two short lines are associated with 2 sigma amplitude in x. The two thick lines are, as explained in the previous paragraph, associated with 2 sigma amplitude in z and finally the two thin lines correspond to zero amplitude in z. Since the separatrices are overlapped at the electrostatic septum the spirit of the Hardt condition is fulfilled, if not the equation.

Figure 6.5. Extraction plot (exact tracking)

0.0003 0.00032 0.00034 0.00036 0.00038 0.0004 0.00042 0.00044 0.00046 0.00048 0.0005 0.035 0.037 0.039 0.041 0.043 0.045 0.047

Figure 6.6. The four extraction separatrices for amplitude in x and z of 0 and 2 sigma.

As anticipated, the separatrix besides getting thicker in proportion to Ez, also gets

displaced. Such a displacement can be explained if the tune dependance on emittance is
considered. The presence of sextupoles, indeed, creates such a dependance since the

focusing force is not proportional to the particle position. An expression for those tune shifts can be found in [23] and is reported hereafter:

x

SLIDE 116

Transverse aspects of the extracted beam Chapter 6

________________________________________________________________________________________________________________________________________________________________________________________________________

115

>

@

>

@

>

@ > @

>

@

>

@

>

@

'
z

x z x k z j z k x j x z kz jz z x z x k z j z k x j x z kz jz x x k x j x x kx jx x x k x j x z kz jx jz kx x kx jx K j K k kx jx k j x

Q Q Q Q N s s s s J Q Q Q Q N s s s s J Q Q N s s J Q Q N s s J J l k l k Q 2 sin 2 1 2 2 cos 2 2 sin 2 1 2 2 cos 2 3 sin 1 2 3 cos sin 1 2 cos 3 3 ' ' 64 1

1 1 2

S S \ \ \ \ E E S S \ \ \ \ E E S S \ \ E E S S \ \ E E E E E E E E S

>

@

>

@

>

@

>

@

>

@

>

@

u
u
'
z

x z x k z j z k x j x z x kz jz z x z x k z j z k x j x z x kz jz x x k x j x z kz jz x kx jz K j K k kx jx k j z

Q Q Q Q N s s s s J J Q Q Q Q N s s s s J J Q Q N s s J J l k l k Q 2 sin 2 1 2 2 cos 2 2 sin 2 1 2 2 cos 2 sin 1 2 cos 4 ' ' 64 1

1 1 2

S S \ \ \ \ E E S S \ \ \ \ E E S S \ \ E E E E E E S

which can be summarised in

'

'

z x z x

J J a a a a Q Q

22 12 12 11 2 2

where J0x,z =(1/2) Ex,z and the aij, are complicated expressions depending on the sextupole strengths and on the Twiss parameters at the sextupoles. The interesting term for the separatrix displacement is a12. A tune shift, via chromaticity, corresponds to a momentum error. In turn a momentum error becomes a displacement in phase space because of the dispersion function at the electrostatic septum. If the separatrix displacement is calculated according to the expression above, a displacement of 4.3 P rad per 1 S mm mrad is found to be compared with the 4.9 found by the simulations. 6.6 Effect of vertical, closed-orbit distortion In Chapter 3 it has been shown that the effect of a horizontal closed orbit distortion at the sextupole can be considered as a tune shift. Since the vertical plane has an effect on the

SLIDE 117

Chapter 6 Transverse aspects of the extracted beam

________________________________________________________________________________________________________________________________________________________________________________________________________

116 “thickness” of the separatrix, it can be expected that a vertical, closed-orbit distortion at the sextupole has also a widening effect. Including a vertical offset ZCO, calculating the exact 3-turn transfer function, expanding to second order in H and keeping only terms up to second order in H, X0, Xc0, Z0, Zc0 and ZCO yields

> @ > @ > @ > @

z z z CO z x z z z z z x z

Q Q Z E S Q Q SE X SX X X I S S E E I S S E E H

'

3 sin sin 3 2 2 6 sin 2 sin 2 3 ' 2 3 '

3

(6.8)

in which the Kobayashi terms and the “old” separatrix width are clearly recognisable and which shows an additional term in the separatrix width due to the orbit distortion. Expression (6.8) can be re-written as

> @ > @ > @

'

x x Q E Z E Q S X SX X X

z z CO z z x z

2 sin sin ) cos( 2 2 sin 2 3 ' 2 3 '

3

S S E E H

where the coefficient of the last term is the separatrix width without orbit distortion, equation 6.7. The maximum for f(x) = a sin(x) + sin(2x) is

2 2 2

32 64 1 1 32 3 4 1

a

a a a . (6.9)

Equation (6.9) can be plotted resulting in Figure 6.6

0. 5 1 1. 5 2 2. 5 3 0. 5 1 1. 5 2 2. 5 3 3. 5 4

Figure 6.7. Separatrix widening due to closed orbit distortion, f(a).

Equation (6.9) represents the ratio between the separatrix width with orbit distortion to the separatrix width for zero closed orbit distortion. The separatrix widening can be checked against numerical simulations. The results for a few cases are summarised in the Table 6.2:

SLIDE 118

Transverse aspects of the extracted beam Chapter 6

________________________________________________________________________________________________________________________________________________________________________________________________________

117 zco [mm]

) cos( 2

z z CO

Q E Z S

Theory Simulation 1 0.7 1.5 1.4 2 1.4 2.1 1.9 4 2.8 3.3 3.4

Table 6.2 Comparison of numeric and analytic results for the separatrix widening as a function of the closed-orbit distortion.

6.7 Conclusions A more detailed study of the slow extraction process, which includes the effect of the vertical plane, has been developed. Approximate but simple expressions have been found which foresee a “thickness” for the extracted beam distribution. An explanation for the displacement of the separatrix due to the vertical oscillation amplitude has been found based on existing literature. The results are confirmed by numerical simulations and result in a horizontal emittance for the extracted beam of the order of 0.2 S mm mrad which is much more realistic than the estimate based on the Kobayashi model or on the Liouville theorem approach. The beam is distributed in the horizontal phase space along the extraction separatrix, while in the vertical phase space it shows the usual “elliptical” phase space trajectories. The strong asymmetry between the vertical and horizontal emittance and beam distribution leads to difficulties in describing the beam in the transfer lines, which will be treated and solved in the next chapter.

SLIDE 119

119

Chapter 7 A NEW APPROACH TO TRANSFER LINE DESIGN

As was shown in the previous Chapter, the slowly extracted beam has a markedly different distribution in the two transverse phase spaces. In the vertical one, the beam looks “normal”, meaning that the phase-space distribution is approximately the same as in the synchrotron with the particles filling the self-ellipse. In the horizontal phase-space, the particles are distributed, for all practical purposes, in a long and thin rectangle, hereafter referred to as the “bar” of charge, which is the part of the extraction separatrix cut by the electrostatic septum. Such a strong asymmetry calls for a special design of the transfer lines to the treatment rooms. The extraction lines have to be designed such that the beam size at the patient can be chosen in the range 4 to 10 mm at all the extraction energies and that the dispersion function at the end is null with its derivative. Since in a modern facility rotating gantries are to be foreseen the additional requirement of a rotational optics has to be considered. 7.1 The “unfilled” ellipse approach In a transfer line the beta function is not uniquely determined by the periodicity as in a circular accelerator. Once the quadrupole gradients, edge angles, drift lengths etc. in a lattice have been defined, the transfer function is fixed. If a set of initial values for the Twiss parameters is chosen, the values along the line will adjust themselves in order to maintain the transfer function unchanged. Thus the trajectory of a particle is independent

f whatever set of beta values is used. Since the initial values for the Twiss parameters are

free, it is reasonable to choose them such that they give a useful description of the beam. Choosing the initial values for the vertical plane equal to those of the synchrotron, or even better matching them to the extracted beam ellipse, gives a good description of the beam in the transfer lines and thus matched values been chosen in this study. In a real facility, a measurement of the Twiss parameters would be, of course, the best beam description. In the horizontal plane, the betatron function from the synchrotron does not represent the beam distribution and a direct fitting of an ellipse to the extracted beam would not be a good description for a series of reasons: 1) it would lead to very large Ex values because the emittance is small while the beam size is large; 2) it would lead to strong chromatic effects because emittance and Ex vary strongly within a small momentum deviation; 3) the thickness of the “bar” of charge, due to adiabatic damping in the vertical plane, depends on the extraction energy and a different Ex value has to be foreseen for every extraction energy; 4) due to unavoidable variations in the closed orbit or other relevant parameters, the thickness of the “bar” of charge may change requiring a new matching of the lines; briefly a direct matching of the “bar” of charge would be source of never-ending troubles.

SLIDE 120

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

120 The approach proposed in this thesis is that of using a “large” emittance and a “small” beta to describe an “unfilled” ellipse containing the beam, as shown in Figure 7.1.

Figure 7.1. An unfilled ellipse is used in the horizontal phase space (a), while a “normal” filled ellipse is used in the vertical plane (b).

The unfilled ellipse is chosen such that the “bar” of charge is a diameter of the ellipse itself. The use of a much larger ellipse solves all the problems listed above, because for off momentum particles the ellipse is just less filled and because the same ellipse can be used to describe the beam also when the thickness varies. The price that is paid for such a great simplification is that the usual expression EE describing the beam size has now to be considered as an upper value for the beam dimensions, since the orientation of the “bar” of charge along the transfer line varies with the phase advance. The beam size corresponds to the projection of the beam distribution onto the x-axis and thus a variation of the

rientation of the “bar” of charge corresponds to a variation of the beam size.

7.2 Horizontal beam size control Despite requiring additional care, the proposed “unfilled ellipse” scheme offers an additional possibility for controlling the horizontal beam dimension. By a proper variation

f the phase advance the “bar” of charge rotation can be used to set the beam size at the

end of the line, as shown in Figure 7.2.

Figure 7.2. Rotating the “bar” of charge allows an independent horizontal beam size control.

z zc xline xcline xsync xcsync (a) (b)

10 mm 4 mm

x xc x xc

SLIDE 121

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

121 For this method to work an insertion must be foreseen that is capable of varying the horizontal phase advance without changing the beta function at its exit. Such a device, hereafter called a “phase-shifter” will be described later in this Chapter. The maximum required beam size at the patient determines the (minimum) product of the betatron amplitude function at the end of the line and the emittance of the unfilled ellipse. A criterion for deciding how much larger the unfilled ellipse should be than the “bar” of charge, arises naturally from the minimum beam size required at the patient. The width of the “bar” must be no bigger than this dimension when the “bar” is upright in the ellipse. This sets a minimum height for the ellipse, since the phase-space area of the “bar” is conserved through the line. As has been shown in the previous Chapter, the width of the separatrix can vary due to closed orbit distortions or because of other effects, thus it is preferable to choose a larger unfilled ellipse emittance than the minimum one to keep some safety margin. Choosing to have Dx equal to zero and choosing that the required maximum beam size is

btained with an horizontal “bar” of charge fixes the betatron amplitude function at the

end of the line. Similarly, choosing Dx equal to zero at the entrance to the line fixes the initial value of the betatron amplitude function*. 7.3 Vertical beam size control In the vertical plane, the phase-space ellipse of the beam is filled, the distribution is near-gaussian and a conventional Twiss description of the line is valid. Thus, the beam size can be controlled by adjusting the vertical betatron amplitude function at the patient. In the present design, a module, called a “stepper” and described later in this Chapter, is introduced in the common section of the line. The stepper changes the vertical betatron amplitude function while leaving the horizontal betatron amplitude function and the horizontal phase advance unchanged. In order to transmit the action of the stepper to each

f the treatment rooms, the downstream lines consist of modules with integer-S phase

advances and fixed magnification in the vertical plane (telescope modules, see later). The stepping range depends on the specific design of the line (total magnification from the stepper to the end of the line) and on the vertical emittance which, since the extraction energy is not unique, has to consider the adiabatic damping as well. However the ratio Emax/Emin is fixed and can be estimated by comparing the betatron amplitude necessary to have the maximum beam size with the minimum emittance (maximum extraction energy) to the betatron amplitude necessary to have the minimum beam size with the maximum

emittance. Let En be the normalised emittance, which does not change with the extraction

energy, then

n r r

E R

min , min , 2 min min

E J E

;

n r r

E R

max , max , 2 max max

E J E

* The choice of D = 0 is a natural one. At extraction, in PIMMS, the bar of charge is near horizontal and it is customary to choose D = 0 at

the patient. Moreover in the following Sections the matching condition between modules in the extraction lines will have D = 0.

SLIDE 122

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

122 whence, with the requirements Rmin = 2 mm, Rmax = 5 mm and considering the energy extraction range (60-250 MeV for protons and 120-400 MeV/u for carbon ions), it results

carbon

for 2 . 12 protons for 4 . 13

min , min , 2 min max , max , 2 max min max r r r r

R R E J E J E E

The normalised emittances for protons and carbon ions in the synchrotron are better chosen such that the beta range for carbon ions falls inside the beta range for protons, to avoid a further enlargement of the total range. In PIMMS, the choice of having the same geometric emittance at the minimum extraction energy has been made. 7.4 Initial dispersion As for the betatron amplitude function, the synchrotron values for Dx and Dxc are also not representative of dispersion for the extracted beam. Again, in principle, any value can be chosen as an initial value, but a useful one should be chosen. Since the dispersion function represents the displacement of the off-momentum particles, it is natural to choose as initial value the one which corresponds to the distance between the center of gravity of the

n-momentum particles and the center of gravity of the off-momentum ones divided by the

relative momentum difference, as shown in Figure 7.3.

0.0003 0.00032 0.00034 0.00036 0.00038 0.0004 0.00042 0.00044 0.00046 0.00048 0.0005 0.035 0.037 0.039 0.041 0.043 0.045 0.047

Figure 7.3. Definition of the initial value for D and Dc.

Although to a close inspection neither is the spiral step exactly linear with the initial particle amplitude nor is the separatrix a straight line in the amplitude vs momentum Steinbach diagram is, such a definition is good for all practical purposes. 7.5 Basic design concepts The principal steps in the design of the transfer lines can be summarised as: x Exploiting the “bar” of charge to create an independent control of the beam size by rotating the “bar” in an unfilled phase-space ellipse using a “phase shifter” in the common part of the lines.

Center of on momentum particles Center of off momentum particles

SLIDE 123

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

123 x Controlling the vertical beam size by a “stepper” that provides a variable magnification

f the vertical betatron amplitude function while leaving the horizontal betatron

amplitude function and the horizontal phase advance unchanged. The variable vertical beta function is ‘handed’ through telescope modules all the way to the different treatment rooms. x Placing the phase shifter and the stepper in the common initial part of the line so that they can act for all treatment rooms in the complex. x Using telescope modules (extension modules, switch modules and rotators) with integer-S phase advances in both planes that have one-to-one, one-to-minus one or fixed magnification properties. x Matching the strongly asymmetric slowly extracted beam and non-zero dispersion functions to a rotating gantry by the use of a “rotator”. The modular character of the extraction lines is illustrated schematically in Figure 7.4. The lay-out starts with a matching module that makes the liaison from the synchrotron to the main extraction line. This module, which must be custom designed for a specific centre, delivers zero dispersion and standard lattice functions to the modular line (Ex = Ez = 3 m and Dx = Dz = 0 for the present example). In the horizontal plane, the standard values are ‘handed’ from module to module throughout the line to the patient. In the vertical plane, the standard values are modified by the ‘stepper’ and then ‘handed’ through the modules to the patient. In this nomenclature, the planes between the modules are known as ‘hand-over’ planes.

Figure 7.4. Schematic layout of extraction lines.

This general strategy has been adapted to a specific example for PIMMS. 7.6 Telescopes The basic matching scheme uses one-to-one and one-to-minus one modules with integer-S phase advances, but in certain cases there is a practical advantage in using the more general telescope module. Consider the transfer matrix,

>

@

'
'

'

'
'

'

'
P

D P E E P DD P D D EE P EE P D P E E sin cos / sin 1 cos sin sin cos / ' '

2 / 1 2 / 1 2 / 1 2 / 1

S C S C

(7.1)

Phase shifter ‘Stepper’ Extension module Switching/extension modules ‘Hand-over’ planes between modules Rotator Gantry or Fixed beam Matching section

SLIDE 124

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

124 where C and S are known as the principal trajectories and the other symbols have their usual meanings. If 'P = nS, then sin('P) = 0 and S = 0. The transfer matrix is always independent of the choice of the initial Twiss parameters, therefore S must be zero for any initial D0 and E0. This shows the first result: A lattice with integer-S phase advance in one plane, has the same phase advance for any incoming lattice functions in that plane. Now consider the transfer matrix for the Twiss parameters with S = 0:

2

2 2 2 2 2 2

' ' ' 2 ' ' ' ' ' ' 2 ' ' ' ' ' 2 J D E J D E J D E S S C C CS CC C S S C C SS SC CS CC S CS C (7.2)

If the lattice is matched for D = D0 in (7.2), then Cc = 0 must be true since CSc = 1, being the determinant of the transfer matrix (7.1). Thus it remains

2

2

' 1 ' J D E J D E S CS C (7.3)

and for any incoming set of lattice functions it follows that

ion) magnificat (constant ;

2

E D D E E C

Thus a telescope module provides a fixed magnification between the incoming and

utgoing betatron amplitude functions while maintaining the alpha-functions equal.

Normally, the extension and switching modules in Figure 4 would use C 2 = 1 and the action of the ‘stepper’ would be transmitted directly to the patient. For practical reasons, it may be necessary to reduce beam sizes in certain regions in which case a demagnification followed by a magnification could be used. An overall magnification can be used to move the range of the stepper. 7.7 Rotator Another key module in the extraction line is the rotator that transcribes the normal modes and dispersion functions from the fixed line into those of the rotating gantry. The ability of this module to also map the dispersion function to the gantry system opens the possibility

f closing the gantry’s dispersion bump in the fixed part of the line; a possibility that is not

available for conventional gantries that use either the ‘symmetric’ or the ‘round’ beam method for matching [24]. The rotator is a quadrupole lattice with 2S and S phase advances in the two planes and unit magnification. The elements are mounted so that they can be rotated in proportion to the gantry angle. The rotator’s transfer matrix has the simple form,

SLIDE 125

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

125

z'

z x' x z' z x' x 1 1 1 1

If the rotator is rotated by half of the gantry angle (2T), then it maps the incoming beam into the co-ordinates of the gantry as shown by the matrix multiplication

1

1 1 1 cos sin cos sin sin cos sin cos 1 1 1 1 cos sin cos sin sin cos sin cos T T T T T T T T T T T T T T T T

(7.4)

where on the left-hand side the three matrices represent (right to left) the rotation by an angle T before entering the rotator (to refer the beam to the proper axes of the rotator), the passage through the rotator and then a further rotation by T before entering the gantry (to refer the beam to the proper axes of the gantry). The right-hand side shows the result, that the beam is mapped one-to-one or one-to-minus one from the fixed line to the gantry, making the rotation completely transparent to the optics. Note that since the rotator is dipole free, any dispersion function entering the lattice behaves as a betatron oscillation. For this reason, the dispersion function is rotated in the same way and matched one-to-one from the fixed beam line to the rotating gantry. 7.8 Example modules In the following sub-sections, design examples for the various modules are presented. As explained above, the slow-extracted beam is described with an unfilled phase-space ellipse in the horizontal plane. This makes it possible to use the standard Twiss-formalism but gives only an upper limit for the beam size. The effective beam size depends on the

rientation of the “bar” inside the ellipse. In order to see the true excursions of the beam,

for a given setting of the phase shifter, it is necessary to track the four ‘corners’ of the “bar” of charge. For designing the aperture, however, the unfilled ellipse is sufficient, since the beam will in principle move over the full ellipse area at different times, due to the action of the phase shifter (0 < 'P < S/2). 7.8.1 Matching section The first part of the line is the matching section, whose task is to match the initial dispersion, as defined in Section 7.4, to zero. Moreover the matching section has to deliver at its end, the standard lattice functions Ex = Ez = 3 m and Dx = Dz = 0 to the modular line. The set of initial values chosen at the entry to the electrostatic septum is reported in Table 7.1 and the resulting ellipses together with the extracted beam distribution are shown in figure 7.5. The particle distribution is the same as in Figure 6.6, except for the coordinate translation to the extraction line reference system. On a “normal” scale it shows how thin the “bar” of charge is. The argument about not having to change the ellipse parameters when the thickness changes becomes evident.

SLIDE 126

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

126 Twiss parameter Chosen value Ex [m] 5.000 Dx 0.000 Hx [S mm mrad] 5.000 Ez [m] 6.864 Dz 0.087 D [m] 1.942 Dc

0.026

Table 7.1 Initial set of Twiss parameters for the matching section.

0.0015
0.001
0.0005

0.0005 0.001 0.0015

0.008
0.006
0.004
0.002

0.002 0.004 0.006 0.008

0.0015
0.001
0.0005

0.0005 0.001 0.0015

0.008
0.006
0.004
0.002

0.002 0.004 0.006 0.008

Figure 7.5. Twiss ellipses and beam distribution at the entrance to the extraction line.

The initial part of the matching section is the part of the ring which is traversed by the beam before reaching the magnetic septa and being effectively extracted of the vacuum pipe of the synchrotron. Since the extracted beam is off axis, it crosses the quadrupoles

SLIDE 127

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

127

ff-center and sees consequently a dipole component. This implies that the quadrupoles

have to be considered as combined-function magnets. Moreover the orbit crosses the faces of the magnets at an angle, which means that the magnets have to be modelled with appropriate edge angle. The calculation of the distorted-orbit equivalent lattice can be done easily with a powerful calculation tool like WinAgile [25]. The geometry of the extraction line is illustrated in Figure 7.6 and the optical functions are shown in Figure 7.7.

Grid size 3.5829 [m] Horizontal plan view [X-Y plane]

Figure 7.6. Initial part of the extraction line (filled magnets) and synchrotron.

Betatron amplitude functions [m] versus distance [m] Dispersion functions [m] versus distance [m] 0.000 0.000 34.250 34.250 35.0000 6.0000

6.0000

Horizontal Vertical

Figure 7.7. Optical functions in the matching section

SLIDE 128

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

128 7.8.2 Phase shifter The horizontal beam size is controlled by varying the horizontal phase advance in the extraction line, so that the “bar” of charge rotates in the unfilled ellipse. Figure 7.8 shows a dedicated insertion that varies the horizontal phase advance while keeping the horizontal and vertical Twiss functions at the exit constant, as well as the vertical phase advance (although this is not strictly necessary).

Betatron amplitude functions [m] versus distance [m] 0.000 7.400 15.0000 Horizontal Vertical

Figure 7.8. Optical functions in the phase shifter. ['Px = S/2 to S in steps of S/10, 'Pz = 2.5 radian, Ex = Ez = 3 m and Dx = Dz = 0 at entry and exit.]

7.8.3 Stepper For the control of the vertical beam size, a dedicated module, called the ‘stepper’, is used. In the design example shown in Figure 7.9, Ez covers the range 2 to 27 m at the exit for an incoming value of 3 m with Dz = 0 at entry and exit. In the horizontal plane, Ex = 3 m at entry and exit, Dx = 0 and 'Px = S/2. The horizontal phase advance has to be kept constant during the stepping in order not to vary the orientation of the “bar” of charge.

Betatron amplitude functions [m] versus distance [m] 0.000 7.400 30.0000 Horizontal Vertical

Figure 7.9. Optical functions in the stepper. [Ez = 3 m at the entry and steps from 2, 3 to 27 m in steps of 6 m at the exit. 'Px = S/2, Ex = 3 m and Dx = Dz = 0 at entry and exit.]

It should be recalled that the important feature is the ratio in Ez provided by the stepper and not the absolute values. These can, according to the vertical emittance in a specific

SLIDE 129

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

129 machine design, be magnified in one of the downstream telescope modules to create whatever spot sizes are needed. 7.8.4 Phase shifter-stepper The modules shown in Figures 7.8 and 7.9 are in fact identical and it is possible to combine their functions into a single unit. This is inconvenient inasmuch as a single module has to span over a two dimensional parameter space, which makes the operation more complicated and may reduce the global ranges, but it represents a considerable saving in

space. Figure 7.10, shows the betatron amplitude functions in the combined phase shifter-

stepper for four extreme cases in the parameter range.

Betatron amplitude functions [m] versus distance [m] 0.000 7.400 40.0000 Horizontal Vertical

Figure 7.10. Extreme optical functions in the combined phase shifter-stepper.

7.8.5 Closed-dispersion switching module and extension section. When designing a centre, it is necessary to have a number of utility modules as indicated in Figure 7.4, a switch with a closed dispersion bend and its corresponding extension section are needed for fixed beam lines and for closed dispersion gantries. Figure 7.11 shows the optical functions for the extension module, Figure 7.12(a) shows the Twiss functions of a switching module with a closed dispersion bend and 7.12(b) shows the corresponding geometry. The lattice functions are shown with the range of 2 to 27 m in the vertical plane corresponding to the stepper in Figure 7.9.

Betatron amplitude functions [m] versus distance [m] 0.000 15.465 40.0000 Horizontal Vertical

Figure 7.11. Optical functions for the extension module.

SLIDE 130

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

130

Betatron amplitude functions [m] versus distance [m] Dispersion functions [m] versus distance [m] 0.000 0.000 15.420 15.420 40.0000 0.8000

0.8000

Horizontal Vertical Grid size 1.7000 [m] Horizontal plan view [X-Y plane]

(a) (b) Figure 7.12. Optical functions (a) and geometry (b) for the closed-dispersion switching module.

7.8.6 Open dispersion switching module and extension section The possibility introduced by a rotator to also map the dispersion function to the gantry system opens the possibility of closing the gantry’s dispersion bump in the fixed part of the

line. In Subsection 7.8.8 a gantry of this type is described. The present module has to

create a dispersion suitable for the following gantry. In a practical case, the open dispersion bend and the gantry must be designed as an integral unit. As an example and as anticipated in Section 7.5, a magnification factor of 1/3 in the vertical plane has been chosen for the deflection module. Thus the module delivers a range

f 0.67 to 9.0 m in vertical betatron amplitudes for the gantry. The reason for this choice

is to limit the Twiss functions in the rotator. In the horizontal plane, the lattice is one-to

ne for Ex = 3 m Dx = 0 at entrance and exit. Figure 7.13 shows the optical functions for

the extension module, Figure 7.14(a) shows the Twiss functions of the switching module and 7.14(b) shows the corresponding geometry

Betatron amplitude functions [m] versus distance [m] 0.000 16.503 40.0000 Horizontal Vertical

Figure 7.13. Optical functions for the extension module.

SLIDE 131

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

131

Betatron amplitude functions [m] versus distance [m] Dispersion functions [m] versus distance [m] 0.000 0.000 20.576 20.576 45.0000 6.5000

6.5000

Horizontal Vertical Grid size 2.0000 [m] Horizontal plan view [X-Y plane]

(a) (b) Figure 7.14. Optical functions (a) and geometry (b) for the closed-dispersion switching module.

7.8.7 Rotator As explained in the general part, the rotator is a telescope with unit magnification and phase advances S and 2S. Its function is to rotate the beam by twice as much as the module rotation. Since the gantry which follows that insertion can be freely rotated by ±180º, the rotator must be able to rotate by ±90º on its longitudinal axis. Since the beam is coupled inside the module, the beam sigma matrix description should be used rather than the Twiss-functions, however, at the exit, the distributions are uncoupled and the Twiss- formalism is again valid. In Figure 7.15 the Twiss functions are shown in the rotator, for a rotator at 0º and at 90º, where the x and z plane exchange their roles.

Betatron amplitude functions [m] versus distance [m] 0.000 10.071 30.0000 Horizontal Vertical

Figure 7.15. Twiss functions for a rotator at 0º and 90º.

Figure 7.16 shows an example of a rotator module, when rotating the beam by q (physical angles: gantry q, rotator q). Note that the range of incoming beam sizes is transferred to the orthogonal plane when going through the rotator. For illustration the incoming dispersion is non-zero (as would be the case for an open-dispersion module) and

SLIDE 132

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

132 also exchanges planes, like the beam sizes. Since the beam is coupled inside the module, the beam size is quoted rather than the Twiss-functions.

Horizontal beam envelope [m] versus distance [m] Vertical beam envelope [m] versus distance [m] Note: Linear optics 0.000 0.000 10.071 10.071 0.0100 0.0100

0.0100
0.0100

Betatron envelope Momentum envelope

Figure 7.16. Beam sizes inside the rotator for beam rotation and gantry angle of q.

7.8.8 “Riesenrad” gantry In conventional gantries, the beam is deflected away from the rotation axis and then is re-deflected towards the patient. This implies that the last bending magnet, which is a very big and heavy element, rotates on a mechanical structure at the maximum radius. Without a rotator there is, practically, no other solution since the dispersion at the entrance to the gantry must be zero and at least two magnets are then needed to have zero dispersion at the patient. The use of a rotator allows a, technically, more clever solution: a unique 90º dipole starting on the rotational axis together with a patient positioning system that positions the patient at the maximum radius. Such a device is referred to as “riesenrad”

gantry. In this example, the initial dispersion created by the open dispersion switching

module and mapped into the rotated gantry system by the rotator, closes inside the “riesenrad” gantry, as shown in Figure 7.17. In the present set of modules, the deflection section is a telescope with magnification 1/3 in the vertical plane. Assuming that the Ez range needed at the patient is the one spanned by the stepper, the gantry has to be a telescope with magnification 3. In the horizontal plane a final beta value has to be chosen that satisfies the spot size

requirements. The value Ex = 7.2 m has been chosen in this example that, recalling that the

unfilled ellipse emittance is 5 S mm mrad, corresponds to a 6 mm beam half size. With this assumption the large spiral step particles are distributed in 12 mm and the short spiral step

nes are distributed in 8 mm. The average, roughly assumed as FWHM, corresponds to

the 10 mm required in the specifications.

SLIDE 133

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

133

Betatron amplitude functions [m] versus distance [m] Dispersion functions [m] versus distance [m] 0.000 0.000 15.722 15.722 45.0000 15.0000

15.0000

Horizontal Vertical

Figure 7.17. Optical functions in the “riesenrad” gantry.

7.8.9 Resulting lines A possible combination of the modules is shown in Figure 7.18. To demonstrate the behaviour of the beam along the line the initial distribution has been tracked for two different settings of the phase shifter ('Px = S/2 and S) and the phase space footprints at different locations along the line are shown in Figures 7.19 and 7.20. Finally the Twiss functions for the first gantry line are shown in Figure 7.21.

Grid size 10.0000 [m] Horizontal plan view [X-Y plane]

Figure 7.18. Geometry of a possible set of lines

SLIDE 134

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

134

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

0.006
0.004
0.002

0.002 0.004 0.006

0.0015
0.001
0.0005

0.0005 0.001 0.0015

0.008
0.006
0.004
0.002

0.002 0.004 0.006 0.008

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

2.50E-03
2.00E-03
1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03 2.00E-03 2.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

Figure 7.19. Beam distributions at entry to: (a) line, (b) phase shifter, (c) switch module, (d) rotator, (e) patient. (a) (b) (e) (c) (d)

SLIDE 135

A new approach to transfer lines design Chapter 7

________________________________________________________________________________________________________________________________________________________________________________________________________

135

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

0.006
0.004
0.002

0.002 0.004 0.006

0.0015
0.001
0.0005

0.0005 0.001 0.0015

0.008
0.006
0.004
0.002

0.002 0.004 0.006 0.008

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

2.50E-03
2.00E-03
1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03 2.00E-03 2.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

1.50E-03
1.00E-03
5.00E-04

0.00E+00 5.00E-04 1.00E-03 1.50E-03

8.00E-03
6.00E-03
4.00E-03
2.00E-03

0.00E+00 2.00E-03 4.00E-03 6.00E-03 8.00E-03

Figure 7.20. Beam distributions at entry to: (a) line, (b) phase shifter, (c) switch module, (d) rotator, (e) patient. (a) (b) (e) (c) (d)

SLIDE 136

Chapter 7 A new approach to transfer lines design

________________________________________________________________________________________________________________________________________________________________________________________________________

136

Betatron amplitude functions [m] versus distance [m] Dispersion functions [m] versus distance [m] 0.000 0.000 95.424 95.424 35.0000 15.0000

15.0000

Horizontal Vertical

Figure 7.21. Optical functions in the line to gantry 1.

7.9 Conclusions The strong asymmetry in the distribution of the extracted beam creates a problem in the design of the transfer lines to the treatment rooms. A standard description of the beam by matching the optical functions to the beam is not suitable and would yield exotic and complicated lattices. A novel approach to the transfer line design has been proposed that solves those problems. The proposed scheme is based on the use of an ‘unfilled’ ellipse and on a modular design that allows to extend the center by just joining the different modules according to the

needs. The functionality of the scheme has been demonstrated by a series of example

modules and by transporting the extracted beam through the resulting line.

SLIDE 137

137

Conclusions

Although there is a wealth of expertise on slow extraction in the high-energy accelerator laboratories and an equally extensive experience of cancer treatment using conventional techniques, the field of conformal tumour treatment by active scanning in a synchrotron- based facility is a relatively new field. The slow extraction and accelerator design need to be optimised for very different boundary conditions to those normally applied in high-energy laboratories and, on the medical side, the usage of ion beams presents an entirely new field of clinical research. The precision of active scanning opens new possibilities for treatment but, on the technical side, treatment programs, scanning systems

etc. all require complete redevelopment.

This thesis has given a contribution to this development by improving the understanding of the slow extraction from a synchrotron, the extracted-beam characteristics and the techniques needed to master the machine operation. The changes in philosophy for the matching to the gantry and for the beam size control are direct results of this optimisation and understanding of the synchrotron extraction. The introduction of the ‘rotator’ is an example of an academic idea that has languished in a ‘drawer’ until the need was fully understood. This thesis is also an example of technology transfer from high-energy physics to a practical application in medical physics. It shows that technology transfers are not

automatic. The problems need to be identified and studied. It does show that such

developments furnish interesting topics and provide the possibility for novel development. The main subject of the thesis is a theoretical analysis of the slow extraction from a synchrotron, but that analysis opens practical problems which are faced and solved both in the extraction technique and in the beam handling field. Many other problems arise from the results of this thesis, as for example the need for a new approach to treatment planning due to the extracted beam distribution, which will have to be studied if a “centre-of-excellence” in hadrontherapy is to be built.

SLIDE 138

139

Acknowledgements

This thesis has been made within the framework of the Proton-Ion Medical Machine Study (PIMMS) at CERN. PIMMS was initially a collaboration between the TERA Foundation and Med-AUSTRON that was hosted and supported by CERN. It was later joined by Onkologie 2000 and has benefited from frequent contacts with GSI. I should like to thank the founders of PIMMS, Professor U. Amaldi, Professor M. Regler and the CERN Directorate, for making this study a reality and making it possible for me to participate. I thank the CERN PS Division for their warm hospitality during the study and, in particular, I thank Dr J. Bosser, Dr M. Chanel, Dr P. Lefevre and Dr C. Steinbach who were always ready to help me whatever the problem. I gratefully acknowledge the help, constructive criticism and corrections that I have received from my university supervisor Professor J.E. Augustin and the kindness and support that Professor J. Meyer gave me concerning the administrative procedures. I am indebted to the members of the jury for accepting to devote their valuable time to me. To my PIMMS colleagues, I should like to convey my sincere appreciation for their friendship, support and stimulating discussions. Finally, my special thanks go to the person who pushed me and directed me into the subject without sparing time or effort every day in these last few years, the PIMMS leader Dr P. Bryant.

SLIDE 139

141

REFERENCES

[1] R.R. Wilson, Radiological Use of Fast Protons, Radiobiology, 47, p487-491 (1946) [2] W.H. Lamb, R.J. Lari, Proceedings of international symposium on magnet technology, editors H. Brechna and H.S. Gordon, SLAC, (1965), p.487-96. [3]

J. Rossbach, P. Schmüser, Basic course on accelerator optics, Proceedings of the

1992 CERN Accelerator School, CERN 94-01, (1994), p.17-88. [4] E.D. Courant, H.S. Snyder, Theory of alternating gradient synchrotrons, Annals of Physics, 3 (1958), p.1-48. [5]

C. Steinbach, K.Endo, Simultaneous sharing of slow extraction and internal target,

KEK, KEK-77-18, (1978). [6] P.J. Bryant, K. Johnsen, The principles of circular accelerators and storage rings, (Cambridge university press), Cambridge, (1993). [7] D.C. Carey, The optics of charged particle beams, (harwood academic publishers), New York, (1987). [8] H.G. Hereward, Proceedings Vth Int. Conference on High Energy Accelerators, Dubna, (1963). [9]

Y. Kobayashi, H. Takahashi, Improvement of the emittance in the resonant ejection,

Proceedings VIth Int. Conference on High Energy Accelerators, Massachusetts, (1967), p.347-51. [10] M.Q. Barton, Beam extraction from synchrotrons, Proceedings VIIIth Int. Conference on High Energy Accelerators, CERN, Geneva, (1971), p.85-8. [11] G. Guignard, A general treatment of resonances in accelerators, CERN, CERN-78-11, (1978). [12] W. Schnell, B. Zotter, A simplified criterion for transverse stability of a coasting beam, and its application to the ISR, CERN, CERN-ISR-GS-RF/76-26, (1976). [13] W. Hardt, Ultraslow extraction out of LEAR (transverse aspects), CERN, PS/DL/LEAR Note 81-6, (1981). [14] M. Pullia, Transverse aspects of the slow extracted beam, to be published. [15] C. Steinbach, Beam optics at resonant extraction septa, European Particle Accelerator Conference 1992, Berlin, (Editions Frontières), (1992), p.857-9. [16] S. van der Meer, Stochastic extraction, a low-ripple version of resonant extraction, CERN, CERN/PS/AA 78-6, (1978). [17] J.C Ciret, Extraction du faisceau de Saturne II par acceleration betatronique. Le Gephyrotron, GERMA 76.02/IE-117, (1976). [18] L. Badano, S. Rossi, Characteristics of a betatron core for extraction in a proton- ion medical synchrotron, CERN, CERN/PS 97-19 (DI), (1997). [19] M. Pullia, Transit time for third order resonance extraction, CERN, CERN/PS 96- 36 (DI), (1996).

SLIDE 140

References

________________________________________________________________________________________________________________________________________________________________________________________________________

142 [20] K. Hiramoto, M. Nishi, Resonant beam extraction scheme with constant separatrix, NIM A 322(1992), 154-160. [21] H. Grote, F.C. Iselin, The MAD Program, CERN/SL 90-13 (AP) (Rev.4) [22] S. Wolfram, Mathematica 3.0, Wolfram Research Inc., http://www.wolfram.com/ [23] J. Bengtsson, Non linear dynamics for storage rings with applications to the Low-Energy Antiproton Ring (LEAR) at CERN, CERN 88-05, (1988), p. 47. [24] M. Benedikt, C. Carli, Matching to gantries for medical synchrotrons, Part. Accel.

Conf. 1997, Vancouver, (1997).

[25] P.J. Bryant, WinAgile: Windows Alternating Gradient Interactive Lattice dEsign, http://nicewww.cern.ch/~bryant. [26] P Mandrillon et al, The EULIMA Project, 11th International Conference on Cyclotrons and their Applications Tokyo, Japan; 13-17 Oct 1986.

Detailed dynamics of slow extraction and its influence on transfer line design

N° d’ordre : 254-99 Année 1999

THESE présentée devant l’Université Claude Bernard - Lyon 1 pour l’obtention du DIPLOME DE DOCTORAT (arrêté du 30 mars 1992) présentée et soutenue publiquement le 25 Octobre 1999 par Marco Pullia

Dynamique de l’éjection lente et son influence sur les lignes de transfert

Directeur de thèse: Prof. J. E. Augustin Jury:

(rapporteur)

(rapporteur) J Meyer

(rapporteur)

Table of contents

Introduction

Chapter 1 BASIC ACCELERATOR PHYSICS

@

> @

@

Chapter 2 THEORY OF THE THIRD-INTEGER RESONANCE

@

@

@

@

@

@

@

@

@

@

@

@

^ ` s

`

`

`

`

c

c

Chapter 3 RESONANT SLOW EXTRACTION

Chapter 4 TRANSIT TIME

h

Chapter 5 TIME PROFILE OF THE EXTRACTED BEAM

Chapter 6 TRANSVERSE ASPECTS OF THE EXTRACTED BEAM

> @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @

> @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @

@ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @ > @

> @ > @ > @

@

> @ > @ > @ > @ > @ > @ > @ > @ > @

@ > @ > @

> @ >

@

@

> @ > @

> @ >

@

> @

@

> @

@

> @

> @

> @ > @

> @

> @ > @ .

> @ > @

> @

@

@

@ > @

@

@

@

@

@

@

@

@

@

> @ > @ > @ > @

> @ > @ > @

Chapter 7 A NEW APPROACH TO TRANSFER LINE DESIGN

@

Conclusions