Frobenius-Perron Operator Approach to the Beam-Beam Interaction in Circular Colliders

TL;DR

This work reframes beam-beam interactions in circular colliders through the Frobenius-Perron operator for symplectic twist maps, offering a complementary route to the standard Vlasov-Poisson treatment. By applying Renormalization Group reduction to the non-resonant case, it derives a Fokker-Planck-type evolution for the invariant action-density and explicitly computes the incoherent tune shift that governs relaxation. The linearized FP framework yields a discrete matrix formulation for coherent beam-beam resonances, enabling closed-form stability criteria for isolated resonances and detailed stability diagrams, while the analysis reveals how repeated beam-beam encounters can, in principle, increase luminosity per kick by tuning the phase advances. Collectively, the approach provides a robust, scalable tool for analyzing equilibrium densities, resonance structure, and luminosity control, with potential extensions to multi-degree-of-freedom systems and efficient numerical implementations.

Abstract

Unlike most publications devoted to the application of the self-consistent method of the nonlinear Vlasov-Poisson system to the study of beam-beam interaction, in this article an alternative strategy using the elegant approach of the Frobenius-Perron operator for symplectic twist maps has been developed. A detailed analysis of the establishment of an equilibrium density distribution in phase space, as well as the behavior of the perturbed distribution function with respect to the coherent stability of the two beams, has been carried out. Using the Renormalization Group technique for the reduction of the Frobenius-Perron operator, the case where the unperturbed rotation frequency (unperturbed betatron tune) of the map is far from any structural resonance driven by the beam-beam kick perturbation has been analyzed in detail. It has been shown that up to second order in the beam-beam parameter, the renormalized map propagator with nonlinear stabilization describes a random walk of the angle variable, implying that there exists an equilibrium distribution function depending only on the action variable. The linearized Frobenius-Perron operators for each beam imply a discrete form of the linearized Vlasov equations, which essentially comprises a new method for calculating coherent beam-beam instabilities using a matrix mapping technique. In the special case of an isolated coherent beam-beam resonance, a stability criterion for coherent beam-beam resonances has been found in closed form. An intriguing particular concerning the effect of repeated beam-beam collisions on collider luminosity has been derived explicitly. An addition of luminosity per kick (small though, of the order of the beam-beam parameter) in the course of successive beam-beam collisions could be achieved.

Figures (4)

• Figure 1: Dependence of the first-order incoherent tune shift $- \omega_{3-k}^{(u)} \sigma_{3-k}$ as a function of the action variable $J / \sigma_{3-k}^2$ given by Eq. \ref{['IncohTS']}.
• Figure 2: Stability diagram (the shaded region) of a coherent beam-beam coupling resonance of the form given by Eq. \ref{['CohBBReson']}, where $n_1 = n_2 = 1$. The plot is presented in the fractional part of the tune ${\left( \nu_1, \nu_2 \right)}$-space. For demonstrativeness, the beam-beam parameter is taken to be $\lambda_k \sim 4.7712 \times 10^{-5}$.
• Figure 3: Stability diagram (the shaded region) of a coherent nonlinear beam-beam resonance of the form given by Eq. \ref{['CohBBReson']}, where $n_1 = 1$, $n_2 = 3$. The plot is presented in the fractional part of the tune ${\left( \nu_1, \nu_2 \right)}$-space. For demonstrativeness, the beam-beam parameter is taken to be $\lambda_k \sim 4.7712 \times 10^{-5}$.
• Figure 4: Stability diagram (the shaded region) of a coherent beam-beam coupling resonance of the form given by Eq. \ref{['CohBBReson']}, where $n_1 = n_2 = 1$. The plot is presented in the fractional part of the tune ${\left( \nu_1, \nu_2 \right)}$-space. The beam-beam parameter is taken to be $\lambda_k \sim 4.7712 \times 10^{-6}$, which corresponds to the realistic case, where $N_k \sim 4 \times 10^9$.