Monte Carlo simulations of ultra high vacuum and synchrotron radiation for particle accelerators

With preparation of Hi-Lumi LHC fully underway, and the FCC machines under study, accelerators will reach unprecedented energies and along with it very large amount of synchrotron radiation (SR). This will desorb photoelectrons and molecules from accelerator walls, which contribute to electron cloud buildup and increase the residual pressure - both effects reducing the beam lifetime. In current accelerators these two effects are among the principal limiting factors, therefore precise calculation of synchrotron radiation and pressure properties are very important, desirably in the early design phase. This PhD project shows the modernization and a major upgrade of two codes, Molflow and Synrad, originally written by R. Kersevan in the 1990s, which are based on the test-particle Monte Carlo method and allow ultra-high vacuum and synchrotron radiation calculations. The new versions contain new physics, and are built as an all-in-one package - available to the public. Existing vacuum calculation methods are overviewed, then the steady-state and time-dependent algorithms behind the ultra-high vacuum simulator Molflow are presented. Some practices to tackle the most common problems that arise when simulating large systems are also discussed. Results are compared to theory, and validated through two experiments. Next the the main steps of synchrotron radiation simulations are presented. Properties of SR are summarized, along with optimizations that allow simulating the rather complex underlying physics at a higher speed. The resulting software's photon generation algorithm is benchmarked against published data. The phenomenon of photon stimulated desorption and its literature is overviewed, then two dedicated photodesorption experiments carried out in KEK (Tsukuba, Japan) are presented: one with six room-temperature samples and an other at liquid nitrogen temperature. A simple synchrotron radiation calculation is performed for the LHeC interaction region, allowing to compare Synrad+ results with published analytic calculations. Then the calculations are repeated for a more precise geometry description. The pressure profile of a crotch absorber of the recently started Max IV light source is calculated using Molflow+ and Synrad+ together. Finally the pressure analysis of the SuperKEKB interaction region is presented, consisting of modeling the vacuum chamber and the optics, calculating synchrotron radiation, then performing vacuum simulations. It is confirmed that pressure is expected to meet the design requirements during operation of the machine.

Electron cloud studies for the LHC and future proton colliders

César Octavio Domínguez Sánchez De La Blanca

The Large Hadron Collider (LHC) is the world’s largest and most powerful particle collider. Its main objectives are to explore the validity of the standard model of particle physics and to look for new physics beyond it, at unprecedented collision energies and rates. A good luminosity performance is imperative to attain these goals. In the last stage of the LHC commissioning (2011-2012), the limiting factor to achieving the design bunch spacing of 25 ns has been the electron cloud effects. The electron cloud is also expected to be the most important luminosity limitation after the first Long Shut-Down of the LHC (LS1), when the machine should be operated at higher energy and with 25-ns spacing, as well as for the planned luminosity upgrade (HL-LHC) and future high energy proton colliders (HE-LHC and VHE-LHC). This thesis contributes to the understanding of the electron cloud observations during the first run of the LHC (2010-2012), presents the first beam dynamics analysis for the next generation of high energy hadron colliders, and assists in the prediction of how electron clouds will impact the performance of the future high-luminosity and high-energy machines. In particular, the thesis discusses a method to benchmark pressure measurements at the LHC against electron cloud build-up simulations for identifying the most relevant surface parameters. This method allowed monitoring the effectiveness of LHC “scrubbing runs”, revealing that in the warm regions, the maximum Secondary Electron Yield, δmax, decreased from an initial value of about 1.9 down to about 1.2 (with a low-energy electron reflectivity R ≈ 0.3), thanks to surface conditioning. In addition, the “map formalism”, a good approximation to quickly explain and predict electron cloud effects, has been further developed and applied, for the first time, to optimize the scrubbing process at the LHC. For the HL-LHC, several novel filling schemes have been analyzed in terms of luminosity performance and electron cloud activity. Only a few of them are compatible with an electron cloud activity lower than for the baseline scenario. We highlight a promising option which could be a good fallback scenario in case the electron cloud effects prevent the injection of the baseline beam. This option could also be considered for the nominal LHC after the LS1 if electron cloud turns out to be a serious obstacle. Regarding the future high-energy proton colliders (HE- and VHE-LHC), in the frame of this thesis a performance model was developed to predict the luminosity as a function of time and to optimize the beam parameters, carrying out the first ever performance analysis for these machines. Several scenarios have been considered, including round and flat beams as well as different bunch spacings. The parameters presented in this thesis have been submitted as an input for the most recent update of the European Strategy for Particle Physics. Finally, we also report the electron cloud studies performed for both high energy machines. The large amount of primary photoelectrons generated by synchrotron radiation at these high energies motivates the consideration of high efficiency photon stops as well as other mitigation techniques (e.g. a-C coatings and clearing electrodes). Although for both machines (HE- and VHE-LHC) a tentative bunch spacing of 25 ns has been considered as baseline assumption, the results of this thesis suggest the possibility of going down to 5 ns, since such a beam would present several advantages.

EPFL2013

New Capabilities of the FLUKA Multi-Purpose Code

Francesco Cerutti, Nikolaos Charitonidis, André Donadon Servelle, Philippe Jean Schoofs

FLUKA is a general purpose Monte Carlo code able to describe the transport and interaction of any particle and nucleus type in complex geometries over an energy range extending from thermal neutrons to ultrarelativistic hadron collisions. It has many different applications in accelerator design, detector studies, dosimetry, radiation protection, medical physics, and space research. In 2019, CERN and INFN, as FLUKA copyright holders, together decided to end their formal collaboration framework, allowing them henceforth to pursue different pathways aimed at meeting the evolving requirements of the FLUKA user community, and at ensuring the long term sustainability of the code. To this end, CERN set up the FLUKA.CERN Collaboration (1) . This paper illustrates the physics processes that have been newly released or are currently implemented in the code distributed by the FLUKA.CERN Collaboration (2) under new licensing conditions that are meant to further facilitate access to the code, as well as intercomparisons. The description of coherent effects experienced by high energy hadron beams in crystal devices, relevant to promising beam manipulation techniques, and the charged particle tracking in vacuum regions subject to an electric field, overcoming a former lack, have already been made available to the users. Other features, namely the different kinds of low energy deuteron interactions as well as the synchrotron radiation emission in the course of charged particle transport in vacuum regions subject to magnetic fields, are currently undergoing systematic testing and benchmarking prior to release. FLUKA is widely used to evaluate radiobiological effects, with the powerful support of the Flair graphical interface, whose new generation (Available at http://flair.cem) offers now additional capabilities, e.g., advanced 3D visualization with photorealistic rendering and support for industry-standard volume visualization of medical phantoms. FLUKA has also been playing an extensive role in the characterization of radiation environments in which electronics operate. In parallel, it has been used to evaluate the response of electronics to a variety of conditions not included in radiation testing guidelines and standards for space and accelerators, and not accessible through conventional ground level testing. Instructive results have been obtained from Single Event Effects (SEE) simulations and benchmarks, when possible, for various radiation types and energies. The code has reached a high level of maturity, from which the FLUKA.CERN Collaboration is planning a substantial evolution of its present architecture. Moving towards a modern programming language allows to overcome fundamental constraints that limited development options. Our long term goal, in addition to improving and extending its physics performances with even more rigorous scientific oversight, is to modernize its structure to integrate independent contributions more easily and to formalize quality assurance through state-of-the-art software deployment techniques. This includes a continuous integration pipeline to automatically validate the codebase as well as automatic processing and analysis of a tailored physics-case test suite. With regard to the aforementioned objectives, several paths are currently envisaged, like finding synergies with Geant4, both at the core structure and interface level, this way offering the user the possibility to run with the same input different Monte Carlo codes and crosscheck the results.

FRONTIERS MEDIA SA2022

High Intensity Beam Issues in the CERN Proton Synchrotron

Sandra Aumon

This PhD work is about limitations of high intensity proton beams observed in the CERN Proton Synchrotron (PS) and, in particular, about issues at injection and transition energies. With its 53 years, the CERN PS would have to operate beyond the limit of its performance to match the future requirements. Beam instabilities driven by transverse impedance and aperture restrictions are important issues for the operation and for the High-Luminosity LHC upgrade which foresees an intensity increase delivered by the injectors. The main subject of the thesis concerns the study of a fast transverse instability occurring at transition energy. The proton beams crossing this energy range are particularly sensitive to wake forces because of the slow synchrotron motion. This instability can cause a strong vertical emittance blow-up and severe losses in less than a synchrotron period. Experimental observations show that the particles at the peak density of the beam longitudinal distribution oscillate in the vertical plane due to a short range wake field and following a travelling wave of about 700 MHz. In order to perform measurements, a dedicated single bunch beam was set up with a zero chromaticity plateau around transition energy. Extensive measurements were performed of the dynamics of instability in order to compute rise times and intensity thresholds. These measurements were done for several peak densities and the results show that the longer the bunch length, the higher is the threshold in intensity. Other measurements performed with a small negative chromaticity and another working point show that the intensity threshold can be pushed at higher values—in this case, the threshold was increased by 20%. The particularity of this work is that the instability is triggered during the acceleration. At transition energy, the momentum compaction factor is zero and the exchange of particles between the head and the tail of the beam from synchrotron motion, which is a natural way to damp instabilities, vanishes for few turns. Therefore the measurements at which η the instability is triggered are fundamental to know in which longitudinal regime the instability develops. Macro particle simulations were performed in order reproduce the dynamics of the instability with the HEADTAIL code, that simulates the beam interaction with the impedance of the machine. In the case of the PS, a very detailed impedance model does not exist, therefore a very simple resonator impedance was considered. The parameters of the code were adapted in order to simulate the beam as close as possible to the experimental conditions. The simulations showed that the travelling wave is well reproduced and therefore rise times were extracted for different beam intensities and compared to the measurements. A good agreement is found for a such simple impedance model. The intensity thresholds are reproduced for zero chromaticity within 30% and it is clear that some damping mechanisms occurring in the measurements are not reproduced by the code. In particular, octupolar components of the field, non-linear coupling and transverse space charge are not taken into account. These limitations could cause some discrepancies between model and measurements. This study allows to establish a broadband impedance model which can be use to predict the dynamics of the instability. A few experiments were done also with the use a gamma transition jump. In this case, the instability appears in the adiabatic regime and the rising of the instability is following the peak density of the beam. Both the intensity threshold measurements and the simulations show that the gamma jump is by far the most efficient way to increase significantly the intensity threshold. The choice of adequate working point appears also to be a cure of the instability. This study with the support of measurements provides a rough impedance model confirmed by simulations and an understanding of the different mechanisms triggering the instability. The last part of the thesis is dedicated to proton beam losses studies at injection energy. An eventual upgrade of the high intensity beam is limited by large beam losses at injection in the PS. Losses were measured with the current beam loss monitor system of the PS that they occur while the incoming beam enters in the machine and then during several hundred turns while the beam is circulating. Optics measurements were performed to check the mismatch of the transfer line between the PSBooster and the PS. A significant horizontal dispersion mismatch was found and in particular a difference in optics between the four injection lines of the PS. Aperture bottlenecks were found at the injection septum in both transverse planes and at the maximum of the injection bump. With the help of Monte Carlo simulations, the high radiation outside the ring induced by the losses was understood. However, the mismatch does not explain the continuous losses after the beam is injected. A possible explanation is that the space charge forces contribute to make particles cross resonances and are lost while they are oscillating at high amplitudes. The transverse phase space is refilled due to the combination of the synchrotron motion and space charge. Several solutions were proposed to improve the situation. A new optics of the transfer line and at injection could make the beam size smaller. Some special quadrupoles already installed in the PS could be used to adapt the optics of the incoming beam. Increase the injection energy from 1.4 GeV to 2 GeV kinetic would cause a gain of 65% in space charge forces, and these studies are currently ongoing for the upgrade of the PS in the framework of high-intensity LHC upgrade project. High intensity beams are also perturbed by space charge forces and transverse instabilities. Coherent tune shift measurements were performed at injection and extraction in order to evaluate the transverse imaginary part of the impedance and the contribution of the space charge to betatron frequency shift with the beam intensity.

EPFL2012