Syllabus and slides

1. Syllabus and slides Lecture 1: Overview and history of Particle accelerators (EW) Lecture 2: Beam optics I (transverse) (EW) Lecture 3: Beam optics II (longitudinal) (EW) Lecture 4: Liouville's theorem and Emittance (RB) Lecture 5: Beam Optics and Imperfections (RB) Lecture 6: Beam Optics in linac (Compression) (RB) Lecture 7: Synchrotron radiation (RB) Lecture 8: Beam instabilities (RB) Lecture 9: Space charge (RB) Lecture 10: RF (ET) Lecture 11: Beam diagnostics (ET) Lecture 12: Accelerator Applications (Particle Physics) (ET) Visit of Diamond Light Source/ ISIS / (some hospital if possible) The slides of the lectures are available at http://www.adams-institute.ac.uk/training/undergraduate Dr. Riccardo Bartolini (DWB room 622) r.bartolini1@physics.ox.ac.uk

2. Beam Instabilities High current operation in storage rings Wakefields Single bunch effects Multi-bunch effects Numerical analysis of instabilities with tracking codes “Good wakefields” Plasma Wakefield accelerators Lecture 8: wakefields and basic principles of beam instability R. Bartolini, John Adams Institute, 9 May 2013

3. To have high luminosity in colliders and high brilliance in synchrotron light sources, storage rings operate with high circulating currents. The bunches are intense and short ≤ 1 cm. They can generate a strong e.m. field as they travel down the vacuum pipe. The electric fields generated by a bunch act back on the bunch itself or on subsequent bunches via the chamber (wakefields), giving rise to current dependent collective phenomena (collective effects). Wakefields and collective effects R. Bartolini, John Adams Institute, 9 May 2013

4. Collective effects • Main issues • wakefields interact with the stored beam • RF heating of the vacuum chamber components • The beam can become unstable • beam loss: suddenpartial or total loss of the beam • saturated injection: difficulty accumulating and storing beam • and/or its properties are compromised • beam oscillations • emittance blow up • Increase effective emittance (via transverse or longitudinal jitter) • Increased bunch length and energy spread • Jitter in arrival times of bunches • These reduce the luminosity of a collider and the brilliance of a light source R. Bartolini, John Adams Institute, 9 May 2013

5. The trailing electric field is termed a wake field. Discontinuitiesalong the vacuum chamber and finite resistivity of the vacuum chamber walls are the main sources for the generation of wake fields in colliders or light sources Small discontinuities in the chamber act like cavities, where the electron bunch can deposit energy in the form of e.m. fields (the wakefields) The trapped fields can have a long decay time. The energy goes into the heating of the chamber and the em fields are sensed by following bunches over many turns Examples of wakefields Trapped modes in a cavity like structure

6. 10 mA ALS Data 29 mA Bolometer signal (V) 40 mA Time (msec) Some instabilities increase gradually with increasing current whilst others have sharp thresholds. Instabilities are of two types, transverse and longitudinal, for both single and multi-bunch filling of the ring Collective effects: examples Bunch shape distortion R. Bartolini, John Adams Institute, 9 May 2013

7. Example: CSR in storage rings Diamond will operate with short electron bunches for generation of Coherent synchrotron radiation (in the THz region) This operating regime is severely limited by the onset of the microbunch instabilities. Sub-THz radiation bursts appeared periodically while the beam was circulating in the ring A ultra-fast Schottky Barrier Diode sensitive to the radiation with 3.33-5mm wavelength range was installed in a dipole beamport; 1.9 mA 3.0 mA 5.2 mA

8. Transverse beam instabilities Vertical multi bunch instability seen at a pinhole camera (Diamond) Stable beam Vertically unstable beam R. Bartolini, John Adams Institute, 9 May 2013

9. Broad band parasitic losses These wake fields decay quickly, with a time scale of less than one turn. These fields are produced by: • vacuum chamber transitions (bellows, tapers, flanges, insertion devices) • vacuum ports • beam position monitors, strip lines • fluorescent screens • injection elements (kickers and their vacuum chamber) The fields affect charges in the same bunch leading to single bunch instabilities, i.e. electrons at the head of the bunch act on the tail (microwave instability (L), head-tail instabilities (H/V)) Narrow band parasitic losses These wake fields ring for a long time affecting the same bunch or another for many turns. The fields are produced by cavity like objects. Multi-bunch instabilities (H/V)arise from such fields. Short range and long range wakefields R. Bartolini, John Adams Institute, 9 May 2013

10. In its full generality, the computation of the wakefields generated by a bunch of charged particles is a very complex electromagnetic problem. • It requires the solution of the Maxwell’s equations with given source terms and boundary conditions imposed by the vacuum chamber • It is a 3D, time dependent problem that should be solved self consistently, i.e. the equation of the e.m. field are coupled to the equation of motion of the charged particles, so that the driving terms (charge and current densities) are also changing with time. • PIC codes runs, for complex EM problems, may take many weeks of CPU time !!! Finite elements codes also require an intensive computational effort despite some simplifying approximationscustomarily used to deal with this problem: • the wakefields are computed assuming that the bunch distribution does not change under the action of the wakefields themselves in the structure • the charge particle beam is assumed to be ultra-relativistic (v ~ c) • the net effect of the wakefield on each particle is assumed to be small and is computed with the momentum transfer of the particle in the structure • where possible, symmetric geometry allows the simplification of the problem. Computation of Wakefields (I)

11. The electromagnetic field generated by the leading particle with charge q1 is computed at a distance Δz where the trailing particle of charge q is located. The distance Δz is usually considered fixed in the computation of the wakefield. This helps simplifying considerably the task Computation of wakefields (II) q1 is the leading particle q is the trailing particle Δz = z1 – z distance The force acting on the trailing particle is given by the Lorentz force E and B are the electric and magnetic field generated by q1 in the structure R. Bartolini, John Adams Institute, 9 May 2013

12. The variation of the momentum and the energy of the trailing particle is computed assuming that these deviations in energy and momentum are small and the motion of the trailing particle remains uniform, the relative distance does not change as long as the two charged particles travel within the structure C where the wakefields are generated. Integrating the expression for the Lorentz force sensed by the trailing particle as it travels along the structure C, we get the momentum variation Wakefields and wake functions (I) and the energy variation R. Bartolini, John Adams Institute, 9 May 2013

13. The longitudinal wake function is the energy variation per unit charge q and q1 longitudinal wake function [V/C] The wake function w|| describes the total energy lost by the test particle q generated by a single point-like charge q1 and therefore is the Green functionfor the problem: the energy lost by a test particle with charge q due to a collection charge Q = eNp described by a longitudinal charge density ρ is given by the convolution Wakefields and wake functions (II) where Np is the total number of particle of charge e in the bunch. Similarly for the transverse planes we can define transverse wake function and the total momentum transfer imparted by a bunch with charge density ρ is given by

14. charge wake s t head: t = z/c tail: t = (s+z)/c Example: wake of a slot cavity like structure Causality implies that no wake can exist in front of the bunch z s

15. The wakefields can be described conveniently also in the frequency domain. The Fourier transform of the wake function is called coupling impedance longitudinal coupling impedance Coupling impedance transverse coupling impedance Since the induced voltage on a trailing particle is given by the convolution of the wake with the charge distribution, in the frequency domain we have i.e. the effect of the wakefield depends on the product of the coupling impedance times the Fourier transform of the bunch distribution (the bunch spectrum). Analogous formulae hold for the transverse plane. R. Bartolini, John Adams Institute, 9 May 2013

16. The impedance is a function of frequency and its spectrum depends on the accelerator components At low frequencies it is dominated by the resistive wall impedance. High Q resonators (cavities) show up as sharp peaks, and the overall impedance made up of the various components in the ring gives the broadband contribution. At frequencies beyond the cutoff frequency, the wake field propagates freely along the chamber. This is reflected by the roll-off of the broadband contribution. The cut-off frequency is given by ωc = c/b, b is the chamber radius. Overlap of the bunch spectrum with the impedance of the machine generate instabilities Impedance of an accelerator R. Bartolini, John Adams Institute, 9 May 2013

17. There is a strong analogy between wakefields and electronic circuit theories. This can be exploited and wakes can be represented by equivalent circuits. For example, the impedance of a parallel RLC circuit is often associated to the impedance of the so-called high order modes (HOM), single resonance wakes in the vacuum chamber. Resonators and broad band impedance (I) Each mode resembles an RLC - circuit and can, to a good approximation, be treated as such. This circuit has a shunt impedance R, inductance L and capacity C. It can drive multi-bunch instabilities. In a real cavity these parameters cannot easily be separated and we use others which can be measured directly: The resonance frequency ωr, the quality factor Q and the damping rate α: R. Bartolini, John Adams Institute, 9 May 2013

18. Using this RLC model the HOMs can be classified in two main categories. Narrow-band impedances. These modes are characterized by relatively high Q and their spectrum is narrow. The associated wake last for a relatively long time making this modes important for multibunch instabilities. Resonators and broad-band impedance (II) Broad-band impedances. These modes are characterized by a low Q and their spectrum is broader. The associated wake last for a relatively short time making this modes important only for single bunch instabilities. R. Bartolini, John Adams Institute, 9 May 2013

19. Numerical computation of wakefield (a DLS example) For complicated structure the wakefields are computed with numerical codes such as MAFIA (or gdfidl , ABCI, URMEL, ...) MAFIA model of the Primary BPM block and BPM button for Diamond R. Bartolini, John Adams Institute, 9 May 2013

20. Diamond BPM wakefields MAFIA provides the wakefield generated by a short pilot bunch (3mm rms; 1 nC) The impendance of the structure is obtained dividing the spectrum of the field generated by the pilot bunch by the spectrum of the pilot bunch R. Bartolini, John Adams Institute, 9 May 2013

21. Impedance database These calculations should be repeated for the main items in the vacuum pipe. Tapers; Small aperture gap (IDs) Pumping holes (and grilles); Kickers chambers Flanges; Cavity like structures Dipole slot for synchrotron radiation Collimators Resistive wall …. The total wakefield is given by the sum of the wakes of each individual item The total wakefield can be approximated as the sum of the wakefields of many RLC resonators which best fit the total impedance Next improvement is to use directly the MAFIA outputs as input in the tracking code. R. Bartolini, John Adams Institute, 9 May 2013 21/32

22. Equation of motion (I) Particles move in the focussing fields of quads (transverse motion) cavities (longitudinal motion) The one turn map for each single particle reads Transverse Longitudinal Includes RF nonlinear potential, chromaticity (and simplecticity)

23. Equation of motion (II) Radiation damping and diffusion for electrons Transverse Longitudinal These terms guarantee that when tracking a distribution of macroparticles the equilibrium distribution is has the correct equilibrium emittances, beam sizes, divergences, bunch length and energy spread. R. Bartolini, John Adams Institute, 9 May 2013 23/32

24. Equation of motion (II) Wakefields apper in the equations of motion as a lumped kick added to the one turn map (both in longitudinal and transverse) Transverse Longitudinal The kick is computed binning the longitudinal distribution  of the electrons and computing W and W|| from analytical formula or numerical codes R. Bartolini, John Adams Institute, 9 May 2013 24/32

25. Structure of sbtrack - mbtrack generate the initial 6D particles distribution apply the one turn map compute and apply the kick due to the wakefields compute the new 6D particle distribution wakefields model • many BBR simultaneously; RW wakefield, CSR, etc • one turn map extended to consider the full nonlinear motion R. Bartolini, John Adams Institute, 9 May 2013 25/32

26. Good wakefields L’OASIS Gas THz Radiation Laser e- bunch Plasma channel Gas jet nozzle Not all of these wakefields are bad in accelerator applications. In fact, there are few examples were wakes play a positive role Bunches in electron storage rings with longitudinal distribution asymmetrically distorted by wake-fields emit coherent synchrotron radiation at much higher frequencies than bunches with nominal Gaussian distribution. This can be exploited for designing THz and far-infrared synchrotron light sources with revolutionary performances Wakefield-based acceleration schemes. Strong R&D and very promising results. 1 GeV e- beam generated at LBNL in 2006 R. Bartolini, John Adams Institute, 9 May 2013 26/32

27. Laser Plasma Acceleration Tajima & Dawson Phys Rev. Lett. 43 267 (1979) 25 years necessary for developing laser system technology suitable for LPAs R. Bartolini, John Adams Institute, 9 May 2013 27/32

28. LPWA (I) • An ionized electron quivers in the E-field of the laser with a ponderomotive energy: Spatial variation in the ponderomotive energy gives rise to a force, the ponderomotive force: R. Bartolini, John Adams Institute, 9 May 2013 28/32

29. LPWA (II) The ponderomotive force in a laser pulse with intensity of the order of 1018 W/cm2 expels electrons from the region of the pulse to form a trailing plasma wakefield The wakefield is strongest when a resonance condition is met: The electric field within the plasma can reach the wave-breaking limit: Three orders of magnitude larger than the field used in conventional RF accelerators.

30. Bubble regime At very high laser intensities the wake reach the “blow-out” or “bubble” regime. A cavity with strong electric fields is created after the laser pulse. Snap-shot from PIC simulation of bubble acceleration: electron density map, propagation direction z. The typical length scale is the plasma wavelength, thus micrometers. The “bubble” behind the laser can trap nC charge, thus yielding electron beam currents on the scale of 100 kA GeV electron beams have been created with this method R. Bartolini, John Adams Institute, 9 May 2013 30/32

31. Wave analogies R. Bartolini, John Adams Institute, 9 May 2013 31/32

32. Wakefields J. D. Jackson, Classical Electrodynamics, John Wiley & sons. A. Chao, Physics of collective beam instabilities in High energy Accelerators, John Wiley & Sons Laser Plasma Wakefield Accelerators E.Esarey, IEEE Trans. Plasma, 24, 232, (1996) W. Leemans et al., Phil. Trans. R. Soc. A 2006 364, 585-600 Bibliography R. Bartolini, John Adams Institute, 9 May 2013