Reduced Collision Optics Stabilize Four-Point Beam Dynamics in Next-Generation Lepton Collider

A comprehensive electromagnetic modeling and beam dynamics investigation has demonstrated that modifying the interaction point optics of the proposed Future Circular Collider electron-positron ring resolves severe collective beam instabilities that arise when expanding the machine to four collision points. As particle physicists and detector designers evaluate a layout featuring four distinct interaction points instead of the earlier two-point baseline, the complex interplay between electromagnetic wakefields generated by circulating bunches and the localized electromagnetic forces of colliding beams threatened to eliminate the stable operational tune space required to maintain design luminosity. By carrying out full three-dimensional, self-consistent tracking simulations that incorporate refined broadband coupling impedance models of the vacuum chamber, researchers identified that compressing the horizontal betatron function at the collision points from fifteen centimeters down to ten centimeters successfully recovers wide, stable operating regions. This optical mitigation suppresses destructive head-tail dynamics while preserving collision luminosity well above the project baseline threshold of one hundred and eighty times ten to the thirty-fourth inverse square centimeters per second.

The Future Circular Collider project represents a proposed scientific initiative to construct an advanced particle collider complex within an underground circular tunnel spanning roughly one hundred kilometers in the CERN region. Conceived as a major research instrument for the post-Large Hadron Collider era, the circular facility encompasses both a hadron machine and an electron-positron collider, known as FCC-ee, whose technical feasibility study is underway to inform the European Particle Physics Strategy Update scheduled for 2026 to 2027. The lepton collider is planned to operate across four discrete center-of-mass energy configurations tailored to high-precision investigations of key fundamental particles: 45.6 gigaelectronvolts for the Z pole, 80 gigaelectronvolts for W boson pairs, 120 gigaelectronvolts for the Higgs boson, and 182.5 gigaelectronvolts for top quark pairs. Among these operational regimes, the low-energy Z pole configuration is the most vulnerable to collective instability phenomena because its beam parameters combine an exceptionally high single-bunch population of 2.53 times ten to the eleventh particles with ultra-small horizontal and vertical emittances.

To rigorously predict how circulating bunches interact with the accelerator infrastructure, the research team undertook a complete revision of the machine coupling impedance model, matching progressive updates in the mechanical design of the beam vacuum chamber. In the updated model, the primary source of longitudinal and transverse impedance remains the resistive wall effect, which is governed by the multi-layer composition of the vacuum chamber. A non-evaporable getter, or NEG, thin-film coating is deposited inside the copper vacuum pipe to ensure effective vacuum pumping and to suppress electron cloud accumulation in the positron ring. While earlier dynamic evaluations assumed a coating thickness of one hundred nanometers, recent manufacturing guidance from vacuum specialists established that a thickness of one hundred and fifty nanometers is necessary to guarantee uniform deposition across the expansive chamber surface. Numerical evaluations conducted with the specialized code IW2D confirmed that increasing the coating thickness to one hundred and fifty nanometers raises the imaginary reactive component of the resistive wall impedance across both the longitudinal and transverse planes without altering the resistive loss factor.

Beyond the coating thickness adjustment, the realistic cross-sectional profile of the vacuum chamber was thoroughly evaluated to capture the electromagnetic impact of lateral winglet extensions. The actual beam pipe design incorporates these protruding lateral channels on either side of the central circular aperture to house synchrotron radiation absorbers that remain physically shielded from direct exposure to the circulating particle beams. Because standard multi-layer boundary algorithms like IW2D only model idealized circular or planar geometries, the scientists implemented an indirect calculation strategy using the three-dimensional electromagnetic software CST Microwave Studio. By first evaluating a thick, single-layer chamber profile possessing a relatively low conductivity of ten to the fifth Siemens per meter and then rescaling the calculated field response using theoretical two-layer surface impedance formulas, the team compared the winglet model with circular geometry calculations up to thirty gigahertz. The comparative analysis demonstrated that the lateral winglets introduce an approximate ten percent increase in both longitudinal and transverse resistive wall impedance relative to a plain circular chamber. This enhancement was formally incorporated into tracking calculations as a numerical Yokoya form factor of 1.1, while asymmetric quadrupolar wakefield terms were determined to be twenty times weaker than dipolar terms and consequently omitted.

Trapped electromagnetic modes that oscillate within the recessed cavity volume formed between the beam chamber and the corrugated bellows convolutions represent the primary source of broadband impedance. The researchers modeled three successive iterations of the bellows: an initial circular model, a circular chamber incorporating lateral winglets, and a high-fidelity mechanical computer-aided design provided by CERN engineers that incorporated squared bellow convolutions alongside comb-type radiofrequency shielding fingers. Performing finite-difference time-domain wakefield calculations on the realistic model proved computationally demanding, requiring up to one billion mesh cells to properly resolve the microscopic electric contact fingers that maintain electrical continuity across the expansion joint. A major finding emerged from these numerical evaluations: while the simplified circular assembly exhibited a strong resonant impedance peak near eleven gigahertz, the presence of lateral winglets in the realistic model broadened the effective aperture to match the outer bellows cavity diameter, completely suppressing this low-frequency resonance and shifting first-order dipolar modes above twenty gigahertz. Multiplying the single-unit response by a conservative total estimate of twenty thousand bellows across the ring established the broadband wakefield catalog.

The comprehensive impedance database also accounted for other beamline components, including beam position monitors and accelerating radiofrequency cryo-modules. Standard beam position monitors contributed only minor localized impedance perturbations around five gigahertz, while earlier evaluations of photon stoppers confirmed their contributions were negligible. The investigation also examined alternative accelerating systems, contrasting the conceptual design report layout—which features fifty-two single-cell four-hundred-megahertz cavities grouped across thirteen dual-taper cryo-modules—with a newly proposed alternative utilizing two-cell six-hundred-megahertz Slotted Waveguide Elliptical cavities. Although the collimator impedance is still missing and its study remains in progress, the researchers almost doubled the bellows impedance to take into account other important and yet unknown possible impedance sources at this preliminary stage.

Equipped with this complete wake potential database, the authors conducted longitudinal single-bunch beam dynamics simulations using the macroparticle tracking software PyHEADTAIL. To resolve micro-bunch structures accurately, the wake potential of an ultra-short Gaussian bunch with an effective length of 0.4 millimeters was generated and applied as a pseudo-Green function, which is ten times shorter than the minimum nominal bunch length of the accelerator. Tracking the bunch over thousands of turns revealed that doubling the lattice momentum compaction factor compared to initial conceptual designs provides substantial mitigation against longitudinal perturbations. The calculated microwave instability threshold occurred at a bunch population near two times ten to the eleventh particles. At the nominal working intensity of 2.53 times ten to the eleventh particles, the bunch operated in a weak microwave instability regime characterized by an energy spread growth of only ten percent, accompanied by moderate bunch lengthening. The researchers underscored that this modest longitudinal blowing up does not lead to particle losses and is predicted to be entirely suppressed during active collisions due to the stabilizing influence of beamstrahlung.

In contrast to the longitudinal plane, transverse tracking simulations highlighted severe stability constraints governing single-beam operations prior to collision. Evaluating the transverse mode coupling instability, also termed the strong head-tail instability, demonstrated that electromagnetic wakefields induce intensity-dependent shifts in the coherent oscillation eigenmodes of the bunch. As the single-bunch population rises, the zeroth azimuthal mode shifts downward toward the negative-first azimuthal mode until the two modes couple, driving an exponential growth in transverse emittance that culminates in rapid beam loss. When PyHEADTAIL simulations incorporated the transverse impedance alone, mode coupling occurred at bunch populations beneath the nominal design intensity. Crucially, when researchers evaluated both longitudinal and transverse wakefields simultaneously in a self-consistent framework, the longitudinal impedance further lowered the instability threshold. Although longitudinal wakefields lengthen the bunch—which typically relieves transverse wakefield strength—they simultaneously produce a downward synchronous tune shift that accelerates the coalescence of mode zero and mode negative-one. Consequently, uncollided beams cannot reach nominal bunch intensity without dedicated mitigation during single-beam commissioning phases.

The most acute operational challenge identified in the investigation arises from the coherent horizontal-longitudinal, or X-Z, instability under the recently adopted four-interaction-point lattice. The collision geometry in the collider combines very small vertical emittances with a large Piwinski angle and a crab waist crossing scheme. Unlike wakefield-induced collective effects distributed along the vacuum chamber, the X-Z instability is driven by localized, highly non-linear horizontal beam-beam forces exerted at each collision point. According to semi-analytical scaling laws, the intensity threshold of this instability scales proportionally with the momentum compaction factor, bunch length, energy spread, and synchrotron tune, while scaling inversely with the horizontal beam-beam tune shift and the horizontal betatron function at the interaction point. When distributing four interaction points symmetrically along the collider, the effective synchrotron tune per collision segment is reduced by a factor of four while the radiation damping time measured in turns is four times longer. Three-dimensional self-consistent collisions modeled with the simulation code IBB revealed that for the baseline horizontal betatron value of fifteen centimeters, the normalized horizontal beam size blew up catastrophically, leaving no stable betatron tune operating space in the horizontal tune range between 0.53 and 0.58, even when introducing positive chromaticity settings of two or five.

To overcome this performance bottleneck, the researchers evaluated an optical mitigation strategy centered on reducing the horizontal betatron function at the interaction point from fifteen centimeters to ten centimeters. While lowering the betatron function restricts machine dynamic aperture and momentum acceptance, recent lattice developments confirmed that a ten-centimeter configuration remains feasible while reducing the horizontal beam-beam tune shift parameter by roughly a factor of two, to approximately 0.0023. Implementing this optical compression in three-dimensional tracking simulations at single-bunch populations of 1.4 times ten to the eleventh and 2.8 times ten to the eleventh particles successfully established two robust, stable operating bands surrounding a horizontal tune of 0.56. Introducing a positive chromaticity of five further broadened these stable operational windows. Although operating at ten centimeters with positive chromaticity incurs an estimated fifteen percent luminosity penalty per collision point at the higher intensity, the resulting performance remains above the baseline design target of one hundred and eighty times ten to the thirty-fourth inverse square centimeters per second, confirming that refined collision optics can reconcile four-point collider operations with collective wakefield constraints.

Subject of Research: Beam coupling impedance modeling and collective instability mitigation in the Future Circular Collider

Article Title: Impedance modelling and collective effects in the Future Circular e+e− Collider with 4 IPs

Article References: Migliorati, M., Antuono, C., Carideo, E., Zhang, Y., & Zobov, M. (2022). Impedance modelling and collective effects in the Future Circular e+e− Collider with 4 IPs. EPJ Techniques and Instrumentation, 9(1), Article 10. https://doi.org/10.1140/epjti/s40485-022-00084-z

Image Credits: AI Generated

DOI: 10.1140/epjti/s40485-022-00084-z

Keywords: Future Circular Collider, FCC-ee, coupling impedance, beam-beam interaction, collective effects, resistive wall impedance, transverse mode coupling instability, X-Z instability, interaction points, PyHEADTAIL, wakefields, bellows shielding, NEG coating, synchrotron radiation absorbers, microwave instability, betatron function