Semianalytical algorithms to study longitudinal beam instabilities in double rf systems

Double rf systems are critical for achieving the parameters of fourth-generation light sources. These systems, comprising both main and harmonic rf cavities, relax statistical collective effects but also introduce instabilities, such as Robinson and periodic transient beam loading (PTBL) instabilities. In this paper, we provide semianalytical algorithms designed to predict and analyze these instabilities with improved accuracy and robustness. The algorithms leverage recent advancements in the field, offering a computationally efficient and accurate complement to multibunch tracking simulations. Using the SOLEIL II project as a case study, we demonstrate how these algorithms can optimize rf cavity parameters in high-dimensional parameter spaces, thereby maximizing the Touschek lifetime. An open-source Python package, (Algorithms for Longitudinal Multibunch Beam Stability), is provided as an accessible tool for double rf system stability analysis.