Impedance & Instabilities The concept of wakefields and impedance Wakefield effects and their relation to important beam parameters Beam-pipe geometry and materials and their impact on impedance An introduction to beam instabilities (including ion effects) Rainer Wanzenberg CAS Vacuum for Particle Accelerator June 7, 2017 Hotel Orenas Slott, Glumslov, Sweden
The concept of Wakefields and Impedance Wake Field = the track left by a moving body (as a ship) in a fluid (as water); broadly : a track or path left Impedance = Fourier Transform (Wake Field) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 2
Electric Field of a Bunch Point charge in a beam pipe Gaussian Bunch, Step out transition Opening angle Wakefields Rainer Wanzenberg Impedance & Instabilities June 2017 Page 3
Geometric wakefields and resistive wall wakes Geometric Wake Resistive Wall Wake Material roughness, Oxide layers σ = 57 10 6 σ = 37 10 ( Ωm) σ = 1.5 10 (1 6 Ω m = 10 ( Ωm) 6 ( Ωm) 6 Ω 1 1 1 mm m Cu Al Steel 2 ) M. Dohlus, M.I. Ivanyan, V.M. Tsakanov Surface Roughness Study for the TESLA-FEL, DESY-TESLA-FEL-2000-26 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 4
Short and Long Range Wakefields Example: small cavity Transient effect: within one bunch, head interacts with the tail Long range wakefields: One or many modes are excited in the cavity bunch interacts with other bunches Rainer Wanzenberg Impedance & Instabilities June 2017 Page 5
Effects on a test charge Change of energy Bunch charge q 1 Test charge q 2 depends on the total bunch charge q 1 q 2 s q 1 Equation of motion q 1 : q 2 : z=c t z = c t - s Rainer Wanzenberg Impedance & Instabilities June 2017 Page 6
Effects on a test charge (cont.) Lorentz force: r F = q 2 r r ( E + c u B) z r q 2 The electric and magnetic fields are generated by the bunch charge q 1 test charge q 2 change of energy transverse kick Bunch: q 1 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 7
Approximations Numerical calculation of a wakefield of a Gaussian bunch traversing a cavity The concept of a Wakefield assumes the Rigid Beam Approximation The wakefield does not affect the motion of the beam The wakefield does not affect the motion of the test charge (only the energy or momentum change is calculated) The interaction of the beam with the environment is not self consistent. Nevertheless, one can use results from wake field calculations for a turn by turn tracking code - sometimes cutting a beam into slices. Rainer Wanzenberg Impedance & Instabilities June 2017 Page 8
Wakefield effects and their relation to important beam parameters Beam Parameter: Total bunch charge q 1 Bunch length (and shape) q 2 s q 1 Transverse dimensions (Emittance, Beta-function) Number of bunches Total beam current Synchrotron tune Betatron tune Damping time (Synchrotron Rad.) Chromaticity Transient effects: Wakefield depends on bunch charge and charge distribution geometry change or material properties Rainer Wanzenberg Impedance & Instabilities June 2017 Page 9
Example: Beam parameters of PETRA III Design Parameter PETRA III Energy / GeV 6 Circumference /m 2304 RF Frequency / MHz 500 RF harmonic number 3840 RF Voltage / MV 20 Experimental Hall Momentum compaction 1.22 10-3 Synchrotron tune 0.049 Total current / ma 100 Number of bunches 960 40 Bunch population / 10 10 0.5 12 Bunch separation / ns 8 192 Undulators Emittance (horz. / vert.) /nm 1.2 / 0.01 Bunch length / mm 12 Damping time H/V/L / ms 16 / 16 / 8 RF system Coupled bunch instabilities threshold current ~ 10 ma powerful broadband feedback neccessary Rainer Wanzenberg Impedance & Instabilities June 2017 Page 10
Wakepotential Longitudinal Wakepotential: t (Unit V/C) Electric field of a Gaussian bunch time of the head: t = z/c test charge tail: t = (s+z)/c (or test charge) energy loss/gain: E = e q1 Wz( s) s z Rainer Wanzenberg Impedance & Instabilities June 2017 Page 11
Transverse Wakepotential Lorentz Force Longitudinal Wakepotential Transverse Wakepotential Change of momentum of a test charge q 2 Kick on an electron (q 2 ) due to the bunch charge q 1 beam energy Rainer Wanzenberg Impedance & Instabilities June 2017 Page 12
Example: TESLA cavity cavity length ~ 1 m rms bunch length σ = 50 µm σ = 250 µm σ = 1 mm E= e 1 nc (-20 kv/nc) = 20 kev Energy spread: 20 kev / 20 MeV = 0.1 % s/σ assuming an accelerating gradient of 20 MV/m Rainer Wanzenberg Impedance & Instabilities June 2017 Page 13
Wake: point charge versus bunch Longitudinal Wakepotential The fields E and B are generated by the charge distribution ρ and the current density j. They are solutions of the Maxwell equations and have to obey several boundary conditions. Line charge density Rainer Wanzenberg Impedance & Instabilities June 2017 Page 14
Longitudinal Impedance Longitudinal Impedance = Fourier transform of the point charge wake Wakepotential of a Gaussian bunch = convolution of point charge wake and charge distribution Fourier transform of the wakepotential of a Gaussian bunch Impedance Fourier transform of the charge distribution Rainer Wanzenberg Impedance & Instabilities June 2017 Page 15
Transverse Impedance Transverse Impedance = Fourier transform of the point charge transverse wake Panofsky-Wenzel-Theorem Relation between the transverse and longitudinal Wakepotential Maxwell equation Rainer Wanzenberg Impedance & Instabilities June 2017 Page 16
Multipole expansion of the wake Longitudinal Wakepotential Multipole expansion in Cartesian coordinates: Using the Panofsky-Wenzel-Theorem: Transverse dipole wake Rainer Wanzenberg Impedance & Instabilities June 2017 Page 17
Wakepotential of a cavity Numerical calculation with the CST Studio suite (commercial 3D code) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 18
Numerical calculations There exist several numerical codes to calculate wakefields. Examples are: Non commercial codes (2D, r-z geometry) ABCI, Yong Ho Chin, KEK http://abci.kek.jp/abci.htm Echo 2D, Igor Zagorodnov, DESY http://www.desy.de/~zagor/wakefieldcode_echoz/ Recently development: Echo3D (version 1.0) Commercial codes (3D) GdfidL CST (Particle Studio, Microwave Studio) The Maxwell equations are solved on a grid Rainer Wanzenberg Impedance & Instabilities June 2017 Page 19
Impedance of a cavity Fourier transform of the time domain calculation Impedance of a cavity modes! The results for the impedance can strongly depend on the mesh and the range of the wakepotential! Rainer Wanzenberg Impedance & Instabilities June 2017 Page 20
Rainer Wanzenberg Impedance & Instabilities June 2017 Page 21 ( ) ( ) cos 2 0 t k t V m ω = > i(t) V(t) ( ) = = R LC k C m 2 1 2 ω no losses: ( ) ( ) t q t i δ = bunch current: with voltage: with losses: m m m Q k L Q C Q R ω ω ω 2 = = = longitudinal impedance (one mode): ( ) ( ) ( ) 1 1 1 2 1 1 + = + + = = ω ω ω ω ω ω ω ω ω ω m m m Q i k Q R C i L i I V Z shunt impedance k q E 2 = loss parameter k Equivalent circuit model for the longitudinal wake
Frequency domain analysis i(t) V(t) Q -value C = ( 2k ) LC = ω 1 2 m R = ω m Q C 2 GHz frequency loss parameter quality R/Q R ω = 2π m 2 k = V 4 Q Q = 2k ω m f m ( U ) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 22
Frequency domain post processing ( ) How to obtain the 2 k = V 4 U loss parameter of a mode 2 E = q k Voltage = integrated electric field with phase factor Point cahrge wake potential of one mode (s > 0) Stored energy Rainer Wanzenberg Impedance & Instabilities June 2017 Page 23
Frequency domain post processing: Dipole modes Numerical Solution of the Maxwell Equations Mode frequency Electric and magnetic field Basic post processing Q-value stored energy in the E, B field of the mode voltage with phase factor Loss parameter at beam offset r, R/Q. R For Dipole modes Tesla cavity, 1.62 GHz Dipole mode Rainer Wanzenberg Impedance & Instabilities June 2017 Page 24
Wakefield effects and their relation to important beam parameters Beam Parameter: Total bunch charge q 1 Bunch length (and shape) Transverse dimensions (Emittance, Beta-function) Number of bunches Total beam current Kick on the test charge: Synchrotron tune Betatron tune Damping time (Synchrotron Rad.) Chromaticity Ω = ω β + i τ 1 Betatron frequency growth / damping rate Rainer Wanzenberg Impedance & Instabilities June 2017 Page 25
Equation of motion (rigid beam approximation) sum: n turns C = Circumference of the storage ring dipole wake ansatz for the solution: Equation of motion is transformed: Impedance is related to betatron frequency shift and instability growth rate Rainer Wanzenberg Impedance & Instabilities June 2017 Page 26
Beam-pipe geometry and materials and their impact on impedance Geometric Wake Resistive Wall Wake Loss parameter and Kick parameter Line charge density: Loss parameter of a mode: Kick parameter and dipole wake Rainer Wanzenberg Impedance & Instabilities June 2017 Page 27
Example: CMS experimental chamber Transient heating revolution frequency 2.4 V / nc Number of bunches single bunch charge rms bunch length HL-LHC parameters: P ~ 100 W (2800 bunches, 35 nc) CERN-ATS-Note-2013-018 TECH Rainer Wanzenberg Impedance & Instabilities June 2017 Page 28
Example: PETRA III Vacuum chamber Resistive wall impedance analytic formula (valid for long bunches) 6 σ = 57 10 ( Ωm) 6 σ = 37 10 ( Ωm) 6 σ = 1.5 10 ( Ωm) (1 Ω m = 10 6 Ω 1 1 1 mm m Cu Al Steel b = pipe radius L = pipe length 2 ) Arc: Al, 80 mm x 40 mm Loss and Kick parameter/length: b = 2 cm Al chamber σ z = 10 mm Loss parameter 1.1 V / (nc m) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 29
Example: European XFEL http://www.xfel.eu/ First laser light May 4, 2017 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 30
Example: European XFEL 3 rd harmonic RF dogleg 4 accelerator modules 12 accelerator modules main linac collimator Gun laser heater BC0 Before the Undulator: Impedance BC1 bunch compressors BC2 In the Undulator: Impedance SASE1 total energy loss 35.3 MeV total energy spread 15.4 MeV collimators 80% is related to material properties energy spread 14 MeV ( 35 sections) elliptical pipe cavities 274/412 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 31
An Introduction to beam instabilities Longitudinal Energy spread TESLA cavity Transverse Head Tail Kick on tail due to the wake of the head Rainer Wanzenberg Impedance & Instabilities June 2017 Page 32
Head tail instability Storage ring Circumference C Wakefield: Tail Head Synchrotron oscillations: Positions of Head and Tail are exchanged after one synchrotron oscillation period Equation of motion 0 T s /2 (time for a synchrotron period) N/2 = bunch population of the head Wake potential of the storage ring (also the total kick parameter could be used here) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 33
Head Tail Instability (cont.) phasor notation: 0 T s /2 0 T s Stability requires pure imaginary eigen values of the matrix which can be translated into a criteria for: Rainer Wanzenberg Impedance & Instabilities June 2017 Page 34
Single bunch instability: TMCI TMCI = Transverse Mode Coupling Instability mode coupling Macro particle model Tail Head Modes of the charge distribution (Vlasov equation with wakefields) Measurement (2009) PETRA III betatron tune synchrotron side band I B T0 Q = β k Q < β 4π E / e β Q s 1 2 Tune shift with intensity Kick parameter Q = 37.13, x Q = 30.31, z Q = 0.049, s f f x x s = 16.9 khz = 39.5 khz f = 6.1 khz Rainer Wanzenberg Impedance & Instabilities June 2017 Page 35
Multi bunch instabilities Impedance: HOMs in PETRA Cavities Petra III, h = 3840, rf- buckets N b = 960, 8 ns (125 MHz) N b = 40, 192 ns Multi bunch modes f = 500 MHz, R/Q = 830 Ohm at PETRA III feedback systems are needed often it helps to detune HOMs f = 728 MHz, R/Q = 89 Ohm Rainer Wanzenberg Impedance & Instabilities June 2017 Page 36
Ions Equation for an ideal gas: Boltzmann constant Residual gas density at room temperature (300 K) Ion density bunch population Typical cross section for ionization: 2 Mbarn = 2 x 10-18 cm 2 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 37
Example PETRA III: ion density Residual gas density Ion density compare with 20.8 x10 6 electrons /cm average electron density bunch to bunch distance 8 ns oder 2.4 m = 0.24 ions/cm one bunch with 5 x 10 9 electrons (100 ma, 960 Bunche) = 230 ions/cm 960 Bunchen (after one turn) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 38
Ion Optics: linear model Drift Quad (beam as a lense) Trapped ions in the beam: + critical ion mass number N b = bunch population σ x, σ y beam dimensions r p = 1.535 10-18 m Ions: A = 2 H 2 A = 16 CH 4 A = 18 H 2 O A = 28 CO, N 2 A = 32 O 2 A = 44 CO 2 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 39
Ion effects Effects of trapped ions on the beam: increased emittance, betatron tune shifts, reduced beam lifetime The effect of the ion cloud on the beam can be modelled as a broad band resonator wake field *) (linear approximation of the force) Non linear interaction between beam ion cloud Multi turn effects (trapped ions) + + + + + + + + Single pass effects (fast ion instability) *) L. Wang et al. Phys. Rev. STAB 14, 084401 (2011) Suppression of beam-ion instability in electron rings with multibunch train beam fillings Rainer Wanzenberg Impedance & Instabilities June 2017 Page 40
Example Ion effects at PETRA III PETRA III vertical emittance increase (data from June 2015) larger vert. emittance (~ factor 2) strange tune spectrum + additional lines in the multibunch spectrum 2014: Installation of new vacuum chambers The effects disappeared with improved vacuum conditions (conditioning) Rainer Wanzenberg Impedance & Instabilities June 2017 Page 41
Conclusion: Less is More less ions less impedance more beam current more luminosity more brilliance vacuum chamber design small loss and kick parameters I B T0 Q = β k < β 4π E / e Q s 1 2 Rainer Wanzenberg Impedance & Instabilities June 2017 Page 42
Thank you for your attention! Rainer Wanzenberg Impedance & Instabilities June 2017 Page 43