New Capabilities of the FLUKA Multi-Purpose Code

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.

Monte Carlo Modeling of Crystal Channeling at High Energies

Philippe Jean Schoofs

Charged particles entering a crystal close to some preferred direction can be trapped in the electromagnetic potential well existing between consecutive planes or strings of atoms. This channeling effect can be used to extract beam particles if the crystal is bent beforehand. Crystal channeling is becoming a reliable and efficient technique for collimating beams and removing halo particles. At the European Organization for Nuclear Research (CERN), the installation of silicon crystals in the Large Hadron Collider (LHC) is under scrutiny by the UA9 collaboration with the goal of investigating if they are a viable option for the collimation system upgrade. This thesis describes a new Monte Carlo model of planar channeling which has been developed from scratch in order to be implemented in the FLUKA code simulating particle transport and interactions. Crystal channels are described through the concept of continuous potential taking into account thermal motion of the lattice atoms and using Moliere screening function. The energy of the particle transverse motion determines whether or not it is trapped between the crystal planes while single Coulomb scattering on lattice atoms can lead to dechanneling. The volume capture and reflection applying to quasi-channeled particles are also modeled. Analogously to dechanneling, single scattering is used to determine the occurrence of volume capture. The parameters of the crystals, such as torsion or miscut, are described as well. For channeled particles, the suppression of electromagnetic and nuclear collisions is implemented. It is stronger for particles oscillating close to the center of the channel and is crucial for a correct evaluation of the rate of dechanneling without having recourse to the use of a macroscopic dechanneling length. The UA9-H8 experiment conducted at CERN aims at investigating new crystal physics as well as characterizing crystals that can be of interest in view of the implementation of crystal collimation at CERN, including in the LHC. This experiment uses silicon strip detectors situated on both sides of the crystal. Putting together upstream and downstream tracks in coincidence and matching at an identical fitted location on the crystal, it yields information about the deflections given to the beam population. Several runs from the UA9-H8 experiment are analyzed and compared to the model results. Channeling and dechanneling rates, as well as angular distributions at crystal exit are shown to be in a very encouraging agreement both for strip and quasi-mosaic crystals.

EPFL2014

Development and validation of data sets and Monte Carlo methods for electron/positron transport at low and medium energies

Cezarina Elena Negreanu

When exposed to ionising radiation, living tissue can potentially suffer somatic and genetic damage - effects depending mainly on the radiation dose or energy absorbed, the type of radiation, and the type and mass of cells affected. It is well known that large doses of radiation lead to high damage of the cell nucleus and additional cell structures, which results in harmful somatic effects, and even rapid death of the individual exposed, while at low doses, cancer is by far the most important possible consequence. Understanding the mechanisms by which low doses of radiation cause cell damage is thus of great significance, not only from this viewpoint but also from that of practical medical physics applications such as radiotherapy treatment planning. Ionising radiation, such as electrons and positrons, begins to cause damage to the genome of a living cell by direct ionisation of atoms, thus depositing energy in the DNA double helix itself. The energy threshold for inducing strand-breaks by electrons, however, is around 7 eV, well below the energy levels required for direct ionization. The low-energy electrons that are set in motion around the tracks of energetic charged particles, for example, are responsible for a multitude of low-energy events (energy transfer of the order of 10 eV), which play a significant role in inducing molecular damage. Assessing the spatial configuration of energy transfer events and the deposited energy spectra, in regions of cellular and sub-cellular dimensions, can be aimed at via the application of appropriate Monte Carlo simulation tools. Such calculations depend primarily on an accurate knowledge of the production and subsequent slowing down of secondary electrons that form the basic structure of the charged particle track. In the above context, an important requirement is the provision of detailed quantitative information concerning the interaction cross sections of electrons over an energy range extending down to low energies, i.e. including the sub-excitation domain. In addition, developing fast particle-transport simulation algorithms to cover the entire slowing down process efficiently is a key aspect. Thus, the present doctoral research has, as global goal, the development and validation of new Monte Carlo calculational tools for electron and positron transport in biological materials, both at high and low energies. More specifically, it aims at providing (i) a comprehensive and accurate set of appropriate cross section data, and (ii) a fast and reliable algorithm for the simulation of charged particle transport. Thus, the first part of the thesis concerns the assessment, further development and validation of standard theoretical models for generating electron and positron cross sections to cover the main interactions of these particles with matter, in particular with the basic atomic components of biomaterials (water, bio-polymers, etc.). This has been done for bremsstrahlung, and both elastic and inelastic scattering, considering a wide range of atomic numbers and high up to thermal incident particle energies. In particular, the excitation cross sections for medium and low energies (down to 1 eV) have been derived by using a new formalism based on many-body field theory. The accuracy of the presently obtained data sets are assessed against other theoretical models, as also a large experimental database for each type of interaction, so that both a comprehensive coverage and adequate accuracies have been ensured for the cross section data sets generated. In the second part of the thesis, an extension of the Monte Carlo code system PENELOPE is first undertaken such that use can be made of elastic scattering differential cross sections which have been made available in numerical form. Thereby, new computational routines (incorporated into the new code PENELAST) prepare the cross sections, needed for a given energy and scattering angle, by applying a fast and accurate sampling technique to a provided data set. The present development will allow various electron and positron cross section data libraries, appropriately formatted, to be used with PENELOPE for benchmarking purposes. The development of a high calculation-speed (Class I) Monte Carlo tool for charged particle transport in biological materials has then been addressed. Thereby, a numerical algorithm for calculating the multiple-scattering angular distributions of high energy electrons and positrons is developed, based on the multiple-scattering theory of Lewis which accounts for energy losses within the continuous slowing down approximation. Partial-wave elastic scattering differential cross sections made available in numerical form, as indicated above, are used for the calculations, the inelastic scattering differential cross sections being obtained from the Sternheimer-Liljequist generalized oscillator strength model implemented in PENELOPE. The new code LEWIS has been used to calculate multiple-scattering angular distributions for given path lengths and can be readily adopted for Class I Monte Carlo simulations. The simultaneous generation of a large number of Legendre expansion coefficients is rendered possible, both rapidly and accurately. Results from LEWIS have been found to be in satisfactory agreement, both with detailed simulations carried out using PENELAST and with various sets of experimental data for high to medium energy electrons. In brief, the present research represents a significant improvement in the quality of Monte Carlo modelling of charged particle slowing down processes, thus contributing to understanding the. role of low-energy secondary electrons in radiation protection studies. It will also allow the further development of a complete Class I Monte Carlo code, which can then be reliably used in practical applications such as radiation treatment planning.

EPFL2005

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