FLUKA has been able to track charged particles in magnetic fields since the late 1980s, thus allowing users from a relatively early stage on to model charged particle transport e.g. in complex particle accelerator components, such as dipole or quadrupole magnets and beam kickers. The availability of user routines provides full customization capabilities to tackle the transport of charged particles in arbitrary user-defined magnetic fields.

Instead, the transport of charged particles in electric fields [2] was a longstanding pending ingredient in FLUKA’s transport algorithm. This lack prevented users from modelling electrostatic septa in particle accelerator extraction studies, for example. Users were forced to resort to a proxy magnetic field effectively providing the desired force (permissible whenever the change of the charged particle energy is negligible). To overcome this shortcoming, a dedicated module for charged particle transport in electric fields in vacuum has been recently included in FLUKA.

The transport of charged particles in both magnetic and electric (static) fields relies on the same adaptive stepping algorithm to solve the resulting equation of motion, offering users an accurate control of both the particle step size and the tolerance with which boundaries to nearby regions can be intercepted by the (effectively curved) particle steps. FLUKA’s module for charged particle tracking in electric fields accounts for the variation (decrease or increase) of the particle’s kinetic energy, and it includes an algorithm which treats situations where the electrostatic force opposes the direction of motion, accurately resolving the trajectory turning point where relevant.

Following the example proposed by A. Bielajew [3], Figure 1 displays the excellent agreement between the analytical (solid curves) and FLUKA (symbols) electron (red) and positron (blue) trajectories starting at the origin of coordinates along the direction indicated by the respective arrows in the yz plane, in an electric field of 511 kV/cm along the y axis.

Coherent interactions of charged particles in bent crystals can be used for beam collimation or extraction purposes. Crystals have long been under study at CERN and in other research centers [4, 5] leading to very promising results for the Large Hadron Collider (LHC) collimation [6, 7]. In this context, crystals are planned to be part of the collimation system for High Luminosity LHC Ion runs [8]. In parallel with this activity, a module able to reproduce the various coherent effects such as channeling, volume capture and volume reflection was added to FLUKA [9]. This module also accounts for the change in both interaction and energy loss rate in these regimes. In the following, relevant aspects of the model are introduced and the most recent developments are detailed.

In a FLUKA geometry, a crystal is a region for which several characteristics such as torsion, temperature, radius of curvature and channel orientation are specified. These are set up independently and it is therefore the user’s responsibility to ensure that microscopic features of the device are consistent with its geometry. However, this decoupling offers an additional degree of freedom, for instance the possibility to introduce a deliberate mismatch between crystal edges and plane orientation to reproduce miscut ⁴ .

As the crystal channel direction is defined, the incoming angle of each entering particle can be calculated. When the angle is small enough, the particle becomes trapped in between the transverse Molière potential barriers sitting at each atomic center. The particle enters the channeling regime and is transported along the channel. For a bent crystal, this implies tracking along an arc which is implemented through the use of a proxy magnetic field. In addition, the channeled particle oscillates transversally between the channeling planes with an amplitude depending on its transverse momentum.

This oscillation is important for the computation of reduced interaction rate in channeling. Firstly, the electronic stopping power for channeled projectiles is reduced compared with the one of an amorphous medium. The FLUKA algorithm uses a prescription similar to that of Esbensen et al. [10], maintaining the contribution of distant collisions to the stopping power unaltered and scaling that of close collisions ⁵ by the fraction of the channel electron density subtended by the oscillating trajectory of the particle. The reduction factor is typically of the order of 2, as shown in Figure 2 for the case of 180 GeV protons through a 4 mm-long Si crystal.

Secondly, Coulomb scattering of channeled particles on the screened potential of target atoms is computed on a single event basis rather than in multiple scattering. Contributions from small impact parameter collisions are suppressed, with a cutoff corresponding to the minimum distance between the particle oscillatory path and the atomic centers.

Thirdly, nuclear interactions are scaled to the average nuclear density along the particle oscillatory path. This leads to a substantial reduction of the nuclear interaction rate [11] for the channeled population at large.

Coherent effects are not limited to channeling. Some particles can experience volume effects inside the crystal: a reflection on the crystal planes or a capture in their current channel after interaction. Both take place when the particle trajectory is almost tangent with the planes, which can only happen when their incoming angle has the same sign as the crystal bending. In order to reproduce these mechanisms, the equations of motion of the particle are solved in the potential of interest, in this case an electrostatic Molière potential to which a centrifugal term is added due to bending.

The volume reflection naturally results from this tracking: on Figure 3 the trajectory of a quasi-channeled particle is displayed close to its tangency point with the planes, as a function of the longitudinal coordinate. After crossing several planes, the particle acquires a small angular deflection in the direction opposite to bending. While progressing through this region, quasi-channeled particles can undergo a disrupting scattering event, which can lead to the capture of the particle if its transverse momentum drops below the Molière potential energy. A captured particle becomes channeled from that point on until it is dechanneled or reaches the end of the crystal.

Microscopic treatment of the coherent effects of quasi-channeled particles is an important recent development of the model. The nuclear interaction rate is altered for quasi-channeled particles, akin to what is implemented for channeled ones. However, quasi-channeled particles cross atomic planes and the resulting interaction rate is higher than in amorphous mode. This increase is apparent on Figure 4 displaying the nuclear interaction rate as a function of the incoming angle of the particle for a crystal bent towards positive angles. The main dip reflects the suppression of nuclear interactions for channeled particles. Next to it, the interaction rate comes back to amorphous regime values, with the excess on the right reflecting the increased interaction probability of quasi-channeled particles, in accordance with experimental data. The extent of the excess, however, remains slightly lower than the data shows and will be the focus of upcoming developments.

Another longstanding shortcoming of FLUKA’s transport algorithm concerned the lack of a dedicated model for nuclear reactions of low energy deuterons (below 150 MeV/n). Up to version 4.1, FLUKA included a coarse tentative model to effectively split the deuteron into its constituent neutron and proton, each taking ∼50% of the kinetic energy of the deuteron. FLUKA was therefore not recommended for deuteron-beam applications.

Recently, deuteron transport has gained increased attention, for applications as diverse as the production of radioisotopes for medical applications using deuteron beams [12, 13], their radiobiological effects and associated radiation protection aspects [14], as well as for nuclear fusion studies [15]. Furthermore, deuteron production from hadron- and nucleus-nucleus inelastic interactions is often far from negligible; an accurate description of their subsequent nuclear interactions is necessary. Consequently, a dedicated model for low energy deuteron nuclear reactions has been developed and will be included in FLUKA 4.2.

The model consists of the following main components: deuteron stripping to the continuum is modelled on the basis of the Serber model [16] for the (uncorrelated) energy and angle distribution of the stripped nucleon, coupled to FLUKA’s hadron-nucleus interaction model PEANUT (Pre-Equilibrium Approach to NUclear Thermalization) [17] in order to properly describe the nuclear inelastic interaction of the remaining nucleon with the target nucleus; deuteron elastic break-up is described within the zero-range distorted wave Born approximation (DWBA) for nuclei with atomic number Z ≥ 3 [18]; deuteron stripping to discrete levels is also accounted for within the DWBA [19] for the (d,p) and (d,n) reaction channels on a selection of light isotopes for which this interaction mechanism is relevant, including ^6,7Li, ⁹Be, and ¹²C. More details and benchmarking of this model will be provided in an upcoming dedicated publication.

Furthermore, a complete-fusion channel is included via the Boltzmann Master Equation (BME) event generator, already employed in FLUKA to describe nucleus-nucleus inelastic interactions up to 150 MeV/n [20]. Whereas the three model components described above imply at most one nucleon penetrating the nucleus, the complete-fusion channel effectively permits both nucleons to penetrate the target nucleus, thus providing a larger excitation energy and accounting for the experimentally observed emission of light nuclear fragments beyond a neutron or proton (d, t, ³He, ⁴He, etc.).

As an example showcasing the performance of the new low-energy deuteron nuclear interaction model in FLUKA, Figure 5 displays the double-differential (energy and angle) cross section for the emission of protons from ⁵⁸Ni under 99.6 MeV deuteron bombardment at various emission angles, resolved as a function of the proton energy, comparing in absolute units the FLUKA simulation results (solid curves) with available experimental data [21] (in symbols). The simulated proton energy spectra exhibit a remarkably good agreement with the experimental data, reproducing their essential features: a peak at half the deuteron kinetic energy, particularly intense for forward emission angles and due to deuteron stripping to the continuum, as well as a peak at low proton energies from the evaporation stage of the conjugate neutron-nucleus inelastic interaction. Figure 5 shows that good agreement is preserved throughout the angular domain covered by the experimental data.

The emission of synchrotron radiation (SR) photons is possible in FLUKA by means of a special source term. This allows the user to define up to two arcs of a charged particle trajectory in a dipole field, along which SR photons are sampled as source particles for the simulation. The radiating charged particle itself is not explicitly tracked. Such a scheme is useful for applications where SR emission is dominated by beam particles, for example when modelling arcs of electron storage rings. This approach is however not suitable for cases where secondary and higher-generation particles lose a significant fraction of their energy through SR emission. This is for instance the case for high energy muon storage rings and colliders, where SR emission by decay electrons and positrons cannot be neglected when assessing the heat load distribution in magnets or the beam-induced background in detectors.

To overcome this drawback, the FLUKA tracking algorithm for charged particles in a magnetic field in vacuum has been enhanced by including the emission of SR photons as a supplementary discrete process during tracking.

SR photons are explicitly generated and the corresponding energy loss of the radiating charged particle is taken into account. The user can control the SR emission mechanism via dedicated production thresholds.

The new feature is illustrated in Figure 6, showing the photon flux in a bending magnet of a 1.5 TeV muon storage ring. The two plots compare different FLUKA simulations with and without SR emission by decay electrons. In the latter case, without the new emission model, the photon flux in the vacuum chamber, shielding and magnet is exclusively due to the impact of decay electrons on the machine aperture. The dipole field bends all decay electrons towards the center of the ring, which leads to a localized energy deposition hot spot on one side of the magnet. In the other case, SR photons are generated along the electron trajectories and can impact either on the inner or outer side of the vacuum chamber, resulting in a larger spread of the energy released in the decay of beam particles. At such or even higher energies, neglecting this phenomenon may lead to a non-faithful representation of the secondary particle distribution in the magnets, energy deposition and background to experiments.

The FLUKA advanced graphical interface Flair ⁶ [22] made its first appearance in 2007. It is a fully integrated working environment that assists the user in all stages of a Monte Carlo simulation, from constructing the input and the geometry to running the simulation, processing the output files and plotting the results. Thanks to its advanced capabilities and intuitive use, it quickly became the most popular way of working with FLUKA. Flair is in constant evolution, continuously expanding its features.

Up to version 2, Flair was based on Python v2, which the Python community decided to phase out in 2020 ⁷ . Profiting from this opportunity, it was decided not only to migrate to v3, but also to restructure the code so as to facilitate its further cooperative development. The most notable improvement concerns the decoupling of the graphical interface from the underlying functionality, which was achieved by the introduction of an intermediate abstract layer allowing Flair to support other simulation packages besides FLUKA. While support for MCNP inputs [23] is partial, in the case of Geant4 it is in a quite advanced state, thanks to a dedicated Geant4 application called Moira (see Section 6).

Particular attention has been paid to computational performance, to further optimize Flair both in terms of CPU-time and numerical precision. A considerable speedup (up to a factor 10) has been attained during the loading and displaying of large input files, thanks to the use of a tk:Text module replacing the previous use of a tk:Canvas ⁸ . Likewise, the performance and the numerical precision of the geometry editor has been greatly optimized. This has been achieved through the use of z-buffers for 3D rendering, the implementation of an adaptive anti-aliasing algorithm, and the storage of quadratic surfaces also in reduced form, in order to minimize numerical rounding errors introduced by geometric transformations, to name but a few.

Among the new features of Flair v3, we highlight the graphical tool for planning radiation protection interventions: starting from a FLUKA Cartesian 3D dose map, the user can define a path which an operator will follow in space and time by means of a spline to obtain a report detailing the dose acquired along the path. In addition, short animated movies can be created displaying the 3D geometry (and the simulation results) by means of a camera object moving along a user-defined 3D spline. To reinforce the consistency of the input interpretation by FLUKA and Flair, the Flair input file now persistently embeds all the FLUKA input information, which can still be directly manipulated in the FLUKA input file.

As high-throughput computational power is getting more accessible, the use of Monte Carlo-based treatment planning systems (TPS) becomes a viable option for clinical applications of radiation therapy. Codes like FLUKA have important advantages over classical treatment planning systems as they can take a wide range of physical effects (e.g., full secondary particle generation due to implants) accurately into account. Yet, the best physics results cannot be fully exploited if they are not put into an anatomic context which speaks to the medical practitioner.

A dedicated effort has been invested in making Flair v3 a self-contained working environment for the use of FLUKA in medical physics applications. Recent developments include an improved DICOM handling ⁹ , facilitating the generation of voxel geometries with built-in automatic conversion from Hounsfield units to material composition. For instance, boolean operations can now be applied on predefined regions of interest (ROI) to customize the material assignment in the simulation. Moreover, Flair now includes a DICOM editor to perform common tasks like anonymizing patient information and, in general, the editing of DICOM tags. The embedded RTViewer helps the user to compare the planned dose calculated by a treatment planning system with the dose simulated by FLUKA. An automatic dose-volume-histogram generator (DVH) is included. RTPLAN-modality DICOM files can be imported to generate the corresponding primary-source cards for FLUKA. Furthermore, RTDOSE-modality DICOM files can be converted to FLUKA’s standard 3D mesh scoring (USRBIN), which can not only be readily visualized in 2D and 3D in Flair, but also be used as a volumetric primary source for further simulations. The opposite conversion from USRBIN to RTDOSE is also implemented.

While support for 3D visualization is becoming increasingly available in the domain of Monte Carlo transport codes, visual quality often trails behind current industry standards. In addition, common surface rendering techniques are not suitable for medical visualization due to the loss of contextual information in one dimension. This stems from the fact that there is no well-defined geometry but the input is a discrete cloud of density values obtained from CT or MRI scans. Consequently, there is no defined geometric surface that one can work with nor are the data continuous in space due to the discretization originating from the data acquisition. The most common way is to assign each point a proxy geometry, e.g. a small cube, and render these with standard techniques. However, this leads to the loss of crucial information as intricate 3D anatomic relationships are not presented anymore to the attending physician. To overcome this drawback, the completely new visualization module in Flair v3 now includes a number of volume rendering techniques for medical applications. These methods allow for proper reconstruction of a 3D function from the discrete dataset via ray-marching based algorithms. Texture-based slicing methods would offer performance advantages due to the possibility of GPU implementation. Nonetheless, for large data sets one quickly reaches the limits of standard consumer hardware and, in addition, ray-casting methods seamlessly integrate with Flair’s new photorealistic, interactive ray tracer for standard geometries. Figures 7 and 8 illustrate the use of standard surface visualization methods versus Flair v3’s new medical visualization module using the same dataset. As can be seen, volumetric rendering retains 3D anatomic structures from the scans even on a very fine level and permits a superior evaluation of the dose distribution in a voxel geometry by the medical physicist or physician. Visualization of the geometry is immediate with progressive refinement based on so-called space filling curves opening up real-time user feedback while profiting from optimized use of CPU cache lines.

A key task carried out with FLUKA is the simulation of the radiation field in the CERN accelerator complex, as well as in experimental and test facilities. To this end, FLUKA is used to calculate 3D space distributions of the key radiation level quantities that are needed to assess radiation impact on electronic equipment, both in terms of lifetime degradation and stochastic SEEs.

Lifetime effects are quantified with Total Ionising Dose (TID), i.e., the energy loss due to ionising radiation per unit mass, and the equivalent fluence of 1 MeV neutrons in Silicon, which is used to express Displacement Damage effects in material lattices. On the other hand, FLUKA provides the High Energy Hadron-equivalent (HEH-eq) fluence, i.e., the fluence of all hadrons with E > 20 MeV plus the neutrons in the 0.2–20 MeV range with energy-dependent Weibull weights, and the thermal neutron equivalent fluence, where neutrons are weighted proportionally to the inverse of their velocity. These quantities are directly linked to the SEE evaluation.

Figure 9 shows a representative example of radiation environment characterization. It displays the HEH-eq fluence in a portion of LHC tunnel located between 216 and 260 m from Interaction Point (IP) 1, where beam collisions take place. In this case, FLUKA is employed to generate and transport products of inelastic proton-proton interactions from IP1, inducing radiation showers through machine elements, the LHC tunnel infrastructure and nearby underground areas. This way the space distribution of HEH-eq fluence can be obtained and 1D and 2D projections can be extracted for an area of interest (as shown in Figure 9), which in this case consists of the highlighted tunnel section and the nearby RR17 alcove hosting electronic equipment. The simulation is scaled to an integrated luminosity of 360 fb⁻¹, so as to represent annual radiation levels in the final years of LHC operation after the High-Luminosity LHC upgrade. These results, along with similar calculations in other LHC areas, serve as basis for the definition of radiation level specifications for electronics. They are essential to define qualification requirements in accordance with a standard Radiation Hardness Assurance procedure [26].

As introduced above, radiation effects can be classified as either cumulative or stochastic. Monte Carlo codes are particularly suited to the prediction and better understanding of the latter which are of random and single-event-like nature. This capability is especially relevant for high energy accelerator applications in which the related response to the passage of a broad range of particles and energies [27] needs to be derived from a relatively limited data set, typically consisting of experimental results from 200 MeV protons.

In particular, codes like FLUKA can accurately simulate energy deposition distributions from single-particle events in sensitive volumes (SV) of microelectronic components. After the energy deposition step, additional dedicated simulations are performed to determine charge collection and related circuit response, which in turn are used to predict the SEE cross section of the device. In first approximation however, FLUKA can be used for this purpose as well, modelling the SV of the device as a Rectangular Parallelepiped (RPP) [28] or Integral RPP (IRPP)[29], assuming the energy deposited in this volume directly translates to ionization charge collected by the circuit with 100% efficiency. Under this assumption, an SEE will occur if the collected charge is above a certain critical charge value (Q _crit ) in the RPP case or typically determined experimentally through heavy ion SEE testing, in the IRPP case.

This method has been applied to demonstrate that all hadrons have an approximately equal SEE cross section above 20 MeV [30] in a mixed-field environment. For highly sensitive devices (i.e., low Q c r i t with small SVs, this finding has been validated within a factor 2 up to hundreds of GeV [31]. The resulting SEE rate prediction is typically estimated within ±50%, which is considered acceptable for mixed radiation field applications. Moreover, the related predictions were successfully benchmarked against experimental measurements, [30].

Various exceptions to the above formalism were studied with FLUKA. This includes for example differences between pion- and nucleon-induced SEE cross section [32], the enhanced SEE cross section of low-sensitivity devices (i.e. high Q c r i t with high-Z material near the SV in the hadron energy range from 200 MeV to 3 GeV due to fission reactions [33] and the impact of directly ionizing light charged particles such as protons [34] and muons [35]. Additional studies investigated SEEs induced by electro-and photonuclear reactions [36] and, via a Geant4-based extension, by neutrons in the 0.1–10 MeV neutron range [37].

Beyond accelerator applications, similar FLUKA models were used to predict the dominance of protons versus heavy ions for nuclear induced events in space [38], as well as to use experimental SEE results obtained in hadron environments to predict or bound the error rate induced by heavy ions [39]. Likewise, FLUKA was used to study fragmentation and indirect energy deposition from high-energy heavy ions [40].

As a final illustration of the FLUKA capabilities to gauge radiation effects on electronics, the following lines present results of the simulation of SEEs induced by mono-energetic photons in the 1 MeV-1 GeV range in an RPP model of a 40 nm commercial SRAM ¹⁰ . The SEE cross section is obtained from the integral of the distribution of the energy deposited by the photons in the SV above a fixed threshold referred to as Q _crit , assuming Q [fC] = E_dep [GeV] × 10⁶/22.5 in Silicon. Its simulated value is shown in Figure 10 as a function of Q _crit . Considering a critical charge of 1 fC, Figure  10 shows that photons below 10 MeV are not capable of inducing SEUs. Indeed, they interact mostly through Compton scattering which deposits low energy in the SV.

Over 10 MeV, photo-nuclear reactions lead to higher energy deposition values which can trigger SEUs. The photon-induced SEU cross section is represented in Figure 11 as a function of the photon energy and compared to the neutron-induced one, in both cases for Q _crit = 1 fC. This shows, as expected, that the latter is higher by roughly two orders of magnitude. Still, for specific radiation environments such as those generated by high energy electrons (e.g., downstream the CLEAR beam dump [41], absorbing a beam of 200 MeV electrons) the relative proportion of high energy photons can be prevalent, hence potentially having a significant or even dominant contribution to the overall SEU rate.

The FLUKA source code lifecycle is managed using the CERN GitLab instance ¹² . The latter offers a state-of-the-art infrastructure for full version control based on the git version control system [42] as well as Continuous Integration (CI). A continually growing number of test cases is included in the CI framework, aimed both at probing new development features and ensuring the consistency with respect to the established FLUKA version. These test cases run automatically for every pushed commit. A further in-depth validation at the model level and for full-scale test cases is performed within the FLUKAVAL framework (see next subsection). As an additional quality assurance check prior to public releases, the code is deployed and intensively tested in the respective computing clusters of the collaborating institutes, monitoring its performance and stability in real-world applications.

The FLUKAVAL framework is a Python-based command line application that enables semi-automatic batch submission, processing, validation and reporting of test cases. Any FLUKA input and even model-level based tests can be adapted and integrated in a few steps into the framework. FLUKAVAL makes extensive use of the git version control system to store simulation and reference data sets on the CERN GitLab instance. In its current version, FLUKAVAL is optimised for the submission of a large number of test cases to the CERN Batch Service (LXBATCH), which provides distributed computing power with over 220,000 CPU cores.

FLUKA output data are stored in standardised JSON format, which is compatible with the Geant Validation Portal ¹³ , to simplify future direct comparisons of identical test cases with Geant4.

The results are processed, validated and compiled in a comprehensive evaluated report, which allows the direct comparison of results from the FLUKA version under validation with results obtained with previous FLUKA versions, other radiation transport codes and/or experimental data. A quantitative analysis is performed in addition to a qualitative comparison whenever possible. The resulting report is stored in the CERN document management system once the proper approval process is positively completed.

One of the ways to validate a FLUKA version and achieve quality assurance is to directly and quantitatively compare results from the version under development to results obtained with previous released versions of the code. A growing number of test cases is available within the FLUKAVAL framework, each of them aiming to investigate selected features of the code. Specific acceptance criteria are chosen for each of the test cases based on the observable/physics quantity to validate.

Two kinds of test cases are available at present in order to validate the different features of the code, namely on the model/interaction level and full-scale test cases. The first kind of test cases aims to validate a FLUKA version at the model/interaction level by running specific drivers that have access to internal functionalities of the code. Among others, these drivers allow to estimate reaction cross sections, double differential cross sections and yields as a function of energy for any kind of projectile and known target element either as a single isotope or in its natural isotopic composition. These tailored test cases provide a fine control over the changes and improvements of the implemented physics model. Full-scale test cases, instead, aim to validate the code with FLUKA simulations based on real case applications. These test cases provide a coverage of different fields of application and FLUKA functionalities. They can test several observables such as energy deposition, energy and/or spatial particle fluence distributions, radionuclide inventory and residual dose in a more elaborate case. That can include, for example, complex geometries and material definitions, tracking in magnetic fields, special sources and custom user routines.

Comparisons to other codes also provide useful input for the development process. The FLUKAVAL framework can include simulation results of other Monte Carlo radiation transport simulations codes and compare them quantitatively to FLUKA results. In addition, the FLUKA code is regularly included in international Monte Carlo code inter-comparisons, for example the particle production inter-comparison at the SATIF workshop series [43, 44].

Comparisons to experimental data are vital for the quality assurance of FLUKA. The FLUKAVAL framework can embed experimental data and allows for a quantitative comparison with FLUKA results. The experimental data in the FLUKAVAL framework ranges from the single interaction level data, e.g., radionuclide production cross-sections from the EXFOR database, up to measurements in large experiments.

All the above is complemented by a long-standing initiative to generate new experimental data for benchmarks of radiation transport simulation codes with a focus on the relevance to real-world applications. Examples are neutron spectra at the CERN-EU High-Energy Fields (CERF) facility [45–47], neutron spectra and activation at the CERN Shielding Benchmark Facility (CSBF) [48, 49], activation at the CERN High-energy AcceleRator (CHARM) test facility [50] and activation due to stray radiation from the annihilation of low-energy antiprotons in steel [51].

The first stage of the development project, which is named Moira, is presently in an exploratory phase. To foster synergies with other Monte Carlo codes and benefit from existing development efforts, the new code will be based on the Geant4 framework [52–54]. The technical functionalities of FLUKA are translated into a Geant4-based application, incorporating a wide spectrum of existing FLUKA capabilities (e.g., scoring, biasing, radioactive decay calculations). Where available, existing features of Geant4 will be utilized.

The new code will continue to support combinatorial geometries based on the present FLUKA input format, using Flair as an external geometry navigator. At the same time, it will offer access to Geant4-compatible geometry descriptions, for example in the GDML format. A versatile set of built-in scoring options will be provided, which can be used with any type of geometry description. Moira will be closely integrated with Flair in order to preserve an intuitive visual user experience through a graphical user interface and an advanced geometry viewer.

In a second stage, the FLUKA physics models will be migrated. When the Moira program reaches a sufficient level of maturity, it will be distributed as a new FLUKA generation. The goal is to provide a fully integrated package to the user with a fixed set of physics models that have undergone strict quality assurance processes and thus will deliver well identifiable physics results, being consistent with the ones of FLUKA v4 at that time. The development effort until the first release is expected to expand over several years. Nevertheless, in the interest of the Monte Carlo user community, it is planned to distribute intermediate Moira versions at earlier stages, enabling the use of the FLUKA input files also for Geant4 calculations. The modernization and modularization of the program, profiting in particular from existing interfaces in Geant4, will facilitate contributions from a wider community.