BeamBeam Effects and Parameter Optimization for FCCee Dmitry

Beam-Beam Effects and Parameter Optimization for FCC-ee Dmitry Shatilov BINP, Novosibirsk Acknowledgements: K. Ohmi, K. Oide, F. Zimmermann EPS-HEP, Ghent 13 July 2019

What is FCC-ee? FCC-ee is the first stage of the integrated Future Circular Colliders (FCC) programme double ring e+e- collider ~100 km perimeter follows footprint of FCC-hh, except around IPs asymmetric IR layout & optics to limit synchrotron radiation towards the detector presently 2 IPs (alternative layouts with 3 or 4 IPs under study), horizontal crossing angle 30 mrad, crab-waist optics beam energies [Ge. V]: 45. 6 (Z), 80 (WW), 120 (ZH), 175 182. 5 (ttbar) https: //fcc. web. cern. ch/ FCC-ee: Your Questions Answered https: //arxiv. org/abs/1906. 02693 D. Shatilov synchrotron radiation power 50 MW/beam at all beam energies; tapering of arc magnet strengths to match local energy top-up injection scheme; requires booster synchrotron in collider tunnel EPS-HEP 2019, Ghent 1364 contributors from 351 institutes

Introduction To achieve high luminosity it is necessary to fulfill many conditions: Ø Small horizontal emittance: x < 1 nm. Ø Small betatron coupling: y / x 0. 002. Ø Small beta-functions at the IP: y* 1 mm. Ø Large enough dynamic aperture. Ø Large momentum acceptance: > 2% at high energies. Other important issues: Ø RF systems: high beam current & low voltage at Z, low beam current & high voltage at ttbar. Ø Collective and multi-bunch instabilities, electron cloud (esp. at low energy). Ø Misalignments, lattice errors, corrections. Ø MDI, injection, energy calibration, etc. Ø And finally, beam-beam effects. D. Shatilov In this presentation, we discuss beam effects assuming that beams with the required parameters can be obtained. The question is, what limitations does beam-beam impose and how can the other parameters be optimized to increase the luminosity. FCC-ee has unique features (large Piwinski angle and beamstrahlung) that significantly affect the beam dynamics. New types of beam-beam instability were found in simulations, and then mitigated by proper selection of parameters. EPS-HEP 2019, Ghent

Basic Equations Luminosity: x + e Li ∙ x z 2 Piwinski angle: Large Piwinski angle (LPA) e- 2 x z Collision scheme with large Piwinski angle § Li << z => small without hourglass! << z § Crab waist => large y 0. 2 P. Raimondi, 2006 Beam-beam parameters for flat beams, 1 and 1: Proportional to x , does not depend on x Does not depend on x , x D. Shatilov EPS-HEP 2019, Ghent Increase in Np and z in the same proportion: Li , y and L remain unchanged, x drops. Small y is needed to achieve high y. This implies small betatron coupling and small x.

Main Limitations Associated with Beam-Beam § Two new phenomena were recently discovered in simulations: 1) 3 D flip-flop (occurs only in the presence of beamstrahlung). 2) Coherent beam-beam instability. Both instabilities are bound with LPA and horizontal synchro-betatron resonances – satellites of half-integer. Most strongly manifested at low energies. § Beamstrahlung leads to an increase in the energy spread (several times at low energies) and creates long non-Gaussian tails (mainly at high energies). This requires obtaining a large momentum acceptance (especially at high energies) to ensure the necessary beam lifetime. § For high luminosity, an allowable asymmetry in the population of colliding bunches should be small. This imposes strict requirements on the injector and the scheme of its operation. D. Shatilov EPS-HEP 2019, Ghent

3 D Flip-Flop 1) Asymmetry in the bunch currents leads to asymmetry in z due to beamstrahlung (BS). 2) In collision with LPA, asymmetry in z: The threshold depends on the asymmetry of the colliding bunches. But even in symmetrical case the instability arises (with higher Np). a) Enhances synchrotron modulation of the horizontal kick for a longer (weak) bunch, thus amplifying synchro-betatron resonances. b) xw grows quadratically and yw – linearly with decrease of zs, so the footprint expands and can cross more resonances. All this leads to an increase in both emittances of the weak bunch (at the first stage, mainly xw is affected). 3) An increase in xw has two consequences: 1) Weakening of BS for the strong bunch, which makes it shorter and thereby enhances BS for the weak bunch. 2) Growth of yw due to betatron coupling, which leads to asymmetry in the vertical beam sizes. 4) Asymmetry in y enhances BS for the weak bunch and its lengthening, while BS for the opposite bunch weakens and zs shrinks. Thus the asymmetry in z increases even more. 5) Go back to point 2, and the loop is closed. D. Shatilov EPS-HEP 2019, Ghent Density contour plots ( e between successive lines) in the space of normalized betatron amplitudes. All three beam sizes grow slowly, until the footprint touches strong resonance, then the week bunch blows up.

3 D Flip-Flop (continued) Bunch length Vertical beam size y = 0. 64 y = 0. 61 Turn There are two possible scenarios for 3 D flip-flop: 1) Starts from xw growth (e. g. synchro-betatron resonances 2 x - k z = 1), then yw increases due to betatron coupling. 2) Starts from yw growth (e. g. non-optimal y or strength of crab sextupoles). After zs is sufficiently reduced, and xw increased, the resonances 2 x - k z = 1 lead to xw blowup. In all cases, beamstrahlung plays a key role: z dependence on emittances, and x, y dependence on z. D. Shatilov y = 0. 64 y = 0. 61 x / x 0 z / z 0 y / y 0 y = 0. 64 y = 0. 61 Horizontal beam size Turn To avoid 3 D flip-flop: § Mitigation of synchro-betatron resonances, satellites of half-integer. This is also very important for coherent beam-beam instability (see the next slides). § Avoid the vertical blowup: good choice of the working point, CW strength, etc. We need enough room for the footprint. § Minimize asymmetry in the population of colliding bunches. This sets the requirements for the injector. EPS-HEP 2019, Ghent

Coherent Beam-Beam Instability Discovered by K. Ohmi in strong-strong simulations (BBSS). Reproduced in quasi-strong simulations (Lifetrac). Good agreement between the two codes. Recently it was observed at Super. KEKB (K. Ohmi). The effect is 2 D, x increases 5 15 times. Then betatron coupling leads to y growth in the same proportion, and luminosity falls several times. Evolution of the horizontal emittance x / x x (cm) Bunch shape in the horizontal plane at some turns Turn z / z This instability cannot be mitigated by feedback. The only solution: find conditions under which it does not arise. D. Shatilov EPS-HEP 2019, Ghent

Parameter Optimization at Z (45. 6 Ge. V) Coherent instability: x dependence on x and z. Quasi-strong simulations. URF = 250 MV (red) and 100 MV (green, blue). Luminosity vs. betatron tunes, simplified model, weak-strong simulations. Colors from zero (blue) to 2. 3∙ 1036 cm-2 c-1 (red). § 1. 0 This is done by changing FODO arc cell, which also leads to an increase in x. However, y = 1 pm can be achieved. Besides, the threshold of microwave instability is raised. 2 x - 8 z = 1 2 x - 10 z = 1 x (cm) y x + 2 y 0. 5 (0. 57, 0. 61) x § Decrease (and thus x). This leads to a decrease in the energy acceptance. Eventually it can be reduced to 15 cm. =n 1. 0 The range of permissible x for large y is bounded on the right by 0. 57 0. 58. x The distance between resonances is z. The width depends on x and the order of resonances. We need to reduce x / z ratio and increase the order of resonances near the working point. D. Shatilov Increase the momentum compaction factor: z and z grow, x decreases. EPS-HEP 2019, Ghent § Reduce the RF voltage. This decreases z and x in the same proportion, but increases the order of resonances near the w. p. § Neat choice of x between synchro-betatron resonances.

Bootstrapping § When the energy spread is defined mainly by beamstrahlung, the dependence on Np (bunch population) becomes: z 1 / z 0 z 2 / z 0 x = const, E , z , y , L § With the nominal Np = 1. 7 1011 required for high luminosity, z increases 3. 5 times. § If we bring into collision such bunches with the “initial” z (energy spread created only by SR), the beam-beam parameters will be far above the limits. Np = 4. 0 1010 Np = 5. 0 1010 Np = 6. 0 1010 Np = 4. 5 1010 Np = 5. 5 1010 Turn x 1 / x 0 x 2 / x 0 § The beams will be blown up and killed on the transverse aperture, before they are stabilized by the beamstrahlung. § To avoid this, we have to gradually increase the bunch population during collision, so we come to bootstrapping. Turn D. Shatilov EPS-HEP 2019, Ghent

Parameter Optimization at 80 and 120 Ge. V § In order to obtain a resonant depolarization, we need z 0. 05 => momentum compaction factor should be large – same as at low energy. The RF voltage also is determined by z => URF increased to 750 MV. § Another limitation is the HOM power. This sets the upper limit on Np which corresponds to 2000 bunches. § Perform steps 3 – 6 as described below (except that should be 20 cm). WW (80 Ge. V) Here we do not care about polarization and select the parameters as follows: ZH (120 Ge. V) 1) Lattice with small momentum compaction and small emittances. 2) URF is made “small”, but so that RF acceptance still exceeds the energy acceptance, and this determines z. 3) Then x is selected in the range of 0. 56 0. 58, between synchro-betatron resonances. 4) Look for at which the coherent instabilities disappear; in our case 30 cm is enough. 5) With the given x and , the length of interaction area defines the optimal 6) The lattice optimization for the selected . , to maximize the dynamic aperture and energy acceptance. 7) The bunch population is scanned, while the restriction is the lifetime. Thus we determine the maximum Np and luminosity. Single high-energy beamstrahlung photons become important and they impose a limit on Np. D. Shatilov EPS-HEP 2019, Ghent

Parameter Optimization at ttbar (175 182. 5 Ge. V) The major tool for increasing the lifetime is making larger. For flat beams, is inversely proportional to the surface charge density: Lg 10 ( / 0 ) Energy distribution in the logarithmic scale, black line: Gauss with E = 1. 3 E 0 (assuming Length of interaction area E/ E 0 Luminosity is limited by the beamstrahlung lifetime: – fine structure constant – energy acceptance – bending radius of a particle’s trajectory at the IP D. Shatilov ) Luminosity § To reduce beamstrahlung, x should be increased. As a result, Li grows and we have to increase as well. § We also need to keep y small. Thus x is controlled by which was increased to 100 cm. § At such high energies, the coherent instabilities are suppressed by very strong damping and not so large Piwinski angle, so we can allow increase in. § Momentum acceptance is made asymmetrical to match the actual distribution with beamstrahlung. EPS-HEP 2019, Ghent

Parameter Optimization Summary “Low” energies (Z and WW) • • 3 D flip-flop and coherent X-Z instability are dangerous => p URF Resonant depolarization requires large synchrotron tune => p URF Small emittances are required for high luminosity => p Dynamic aperture, momentum acceptance => There are contradictions between the requirements. The optimum was found taking into account the possibility of changing various parameters. “High” energy (ttbar) • • Coherent instabilities are suppressed by strong damping There is no polarization Small emittances are required for high luminosity Lifetime limitation due to beamstrahlung => p => “Medium” energy (ZH) • Coherent instabilities are weaker, but still exist • There is no polarization, small emittances are better Optimal D. Shatilov => p should be comparable with Li => increase with energy. EPS-HEP 2019, Ghent URF is determined by the energy loss per turn. There is no much freedom for optimization.

Basic FCC-ee Parameters from CDR parameter Z W H (ZH) 45. 6 80 120 60 / 60 90 / 90 momentum compaction [10 -5] 1. 48 0. 73 horizontal emittance [nm] 0. 27 0. 84 0. 63 1. 34 1. 46 vertical emittance [pm] 1. 0 1. 7 1. 3 2. 7 2. 9 horizontal beta* [m] 0. 15 0. 2 0. 3 1 vertical beta* [mm] 0. 8 1 1 1. 6 length of interaction area [mm] 0. 42 0. 85 0. 9 RF frequency [MHz] 400 400 (0. 57, 0. 61, 0. 0125) (0. 562, 0. 60, 0. 0253) (0. 565, 0. 60, 0. 0179) (0. 554, 0. 59, 0. 0409) (0. 554, 0. 59, 0. 0436) 415 77 23 7. 5 6. 6 SR energy loss / turn [Ge. V] 0. 036 0. 34 1. 72 7. 8 9. 2 total RF voltage [GV] 0. 10 0. 75 2. 0 4. 0 + 5. 4 = 9. 4 4. 0 + 6. 9 = 10. 9 energy acceptance [%] 1. 3 1. 7 +2. 4 / -2. 8 energy spread (SR / BS) [%] 0. 038 / 0. 132 0. 066 / 0. 131 0. 099 / 0. 165 0. 144 / 0. 186 0. 150 / 0. 192 bunch length (SR / BS) [mm] 3. 5 / 12. 1 3. 0 / 6. 0 3. 15 / 5. 3 2. 01 / 2. 62 1. 97 / 2. 54 Piwinski angle (SR / BS) 8. 2 / 28. 5 3. 5 / 7. 0 3. 4 / 5. 8 0. 8 / 1. 1 0. 8 / 1. 0 crab sextupoles [%] 97 87 80 40 40 bunch intensity [1011] 1. 7 1. 5 1. 8 2. 2 2. 3 number of bunches / beam 16640 2000 328 59 48 beam current [m. A] 1390 147 29 6. 4 5. 4 0. 004 / 0. 133 0. 010 / 0. 113 0. 016 / 0. 118 0. 097 / 0. 128 0. 099 / 0. 126 allowable asymmetry [%] 5 3 3 Luminosity / IP [1034 cm-2 s-1] 230 28 8. 5 1. 8 1. 55 beam energy [Ge. V] arc cell optics tunes, half-ring (x, y, s) longitudinal damping time [ms] beam-beam parameter (x / y) D. Shatilov EPS-HEP 2019, Ghent ttbar 175 182. 5 1. 8 400 + 800

Summary § The main factors limiting the FCC-ee luminosity at high and low energies were recognized and understood. Mitigation techniques have been found. § The parameters have been optimized at each energy separately, taking into account various requirements and limitations. § The injection scheme requirements have been developed. D. Shatilov EPS-HEP 2019, Ghent