4. METHODS
FLUKA treats an arbitrary three-dimensional configuration of materials in geometric regions bounded by first- and second-degree surfaces. Derived from the Combinatorial Geometry package, it has been entirely rewritten. A completely new, fast tracking strategy has been developed, with special attention to charged particle transport, especially in magnetic fields. New bodies have been introduced, resulting in increased rounding accuracy and speed. Input preparation has been made much easier by the possibility to use names instead of numbers, free format and nested parentheses. The distance to nearest boundary is taken into account for improved performance.
Repetitive structures (lattices) and voxel geometries can also be handled.
The FLUKA hadron-nucleon interaction models are based on resonance production and decay below a few GeV, and on the Dual Parton model above. Two models are used also in hadron-nucleus interactions. At momenta below 3-5 GeV/c the PEANUT package includes a very detailed Generalised Intra-Nuclear Cascade (GINC) and a preequilibrium stage, while at high energies the Gribov-Glauber multiple collision mechanism is included in a less refined GINC. Both modules are followed by equilibrium processes: evaporation, fission, Fermi break-up, gamma deexcitation. FLUKA can also simulate photonuclear interactions (described by Vector Meson Dominance, Delta Resonance, Quasi-Deuteron and Giant Dipole Resonance). The PEANUT model is set to become the default at all energies: please read the release notes for further details about this possibility.
Hadron elastic scattering is described by means of parameterised nucleon-nucleon cross sections, tabulated nucleon-nucleus cross sections and tabulated phase shift data for pion-proton and phase-shift analysis for kaon-proton scattering. Detailed kinematics of elastic scattering is performed on hydrogen nuclei and transport of proton recoils.
Nuclear interactions ge
nerated by ions are treated through interfaces to external event generators: DPMJET-2.5 or DPMJET-3, with a special initialisation procedure, above 5 GeV per nucleon, modified Rqmd-2.4 between 0.1 and 5 GeV per nucleon, and BME (Boltzmann Master Equation) below 0.1 GeV per nucleon. Pre-compiled libraries for these event generators are included in the distributed packages: the source code is not yet included, pending finalization of proper licensing.
Particle transport includes time-dependence.
Transport of charged particles is based on an original treatment of multiple Coulomb scattering and of ionisation fluctuations which allows the code to handle accurately some challenging problems such as electron backscattering and energy deposition in thin layers even in the few keV energy range.
Energy loss of charged particles is based on the Bethe-Bloch theory, with optional delta-ray production and transport with account for spin effects and ionisation fluctuations. Shell and other low-energy corrections are derived from Ziegler, the density effect is according to Sternheimer.
For all charged particles (hadrons and muons as well as electrons and positrons) a special transport algorithm, based on Moliere's theory of multiple Coulomb scattering improved by Bethe, accounts for correlations between lateral and longitudinal displacement and the deflection angle, between projected angles, and between projected step length and total deflection.
The algorithm includes an accurate treatment of boundaries and curved trajectories in magnetic fields, an automatic control of the step, a path length correction, spin-relativistic effects at the level of the second Born approximation, nuclear size effects (scattering suppression) on option, and a correction for cross section variation with energy over the step.
Bremsstrahlung and electron pair production at high energy by heavy charged particles are treated as a continuous energy loss and deposition or
as discrete processes depending on user choice. Muon photonuclear interactions are simulated with or without transport of the produced secondaries.
Differences between positrons and electrons are taken into account concerning both stopping power and bremsstrahlung. The bremsstrahlung differential cross sections of Seltzer and Berger have been extended to include the finite value at "tip" energy, and the angular distribution of bremsstrahlung photons is sampled accurately. The Landau-Pomeranchuk-Migdal suppression effect and the Ter-Mikaelyan polarisation effect in the soft part of the bremsstrahlung spectrum are also implemented.
Positron annihilation is simulated both in flight and at rest. Delta-ray production by positrons and electrons is described via Bhabha and Moller scattering. The lowest transport limit for electrons is 1 keV. Although in high-Z materials the Moliere multiple scattering model becomes unreliable below 20-30 keV, a single-scattering option is available which allows to obtain satisfactory results in any material also in this low energy range.
Photon interactions include pair production with actual angular distribution of electrons and positrons, Compton effect with account for atomic bonds through use of inelastic Hartree-Fock form factors, photoelectric effect with actual photoelectron angular distribution, detailed interaction on six K and L single sub-shells, optional emission of fluorescence photons and an approximate treatment of Auger electrons, and Rayleigh effect. Photon polarisation can be taken into account for Compton, Rayleigh and photoelectric effects. Photohadron production is modelled according to the Vector Meson Dominance Model, modified and improved using PEANUT below 770 MeV, Quasideuteron interactions and Giant Dipole Resonance. Photomuon production is described according to Tsai.
For neutrons with energy lower than 20 MeV, FLUKA uses its own neutron cross section libraries (P5 Legendre
angular expansion, 260 or 72 neutron energy groups), containing more than 200 different materials, selected for their interest in physics, dosimetry and accelerator engineering and derived from the most recently evaluated data. Gamma-ray generation and different temperatures are available. Doppler broadening is applied for temperatures above 0 K.
The neutron transport is based on standard multigroup transport with photon and fission neutron generation, detailed kinematics of elastic scattering on hydrogen nuclei, transport of proton recoils and protons from 14-N(n,p)14-C reaction. Capture photons are generated according to the multigroup treatment, but transported with the more accurate electromagnetic package of FLUKA which performs continuous transport in energy and allows for secondary electron generation. The 2.226 MeV gamma line from capture in hydrogen is generated as an actual precise energy, and the whole gamma cascade is available for Xenon and Cadmium isotopes. For nuclei other than hydrogen, kerma factors are used to calculate energy deposition (including from low-energy fission). Pointwise cross section transport is available for a few nuclei and reactions.
Electron, muon, and tau (anti)neutrinos are produced and tracked on option, without interactions, but neutrino interactions are implemented, independently from tracking.
Generation and transport are available (on user's request) of Cherenkov and scintillation radiation. Transport of light of given wavelength in materials can be simulated with user-defined optical properties.
FLUKA has extended scoring capabilities, requiring in most cases no user-written code. Quantities which can be scored include star (hadron inelastic interaction) density by producing particle and region, energy density by region, total or from electrons/photons only, and energy and momentum transfer density in a geometry-independent binning structure (Cartesian or cylindrical), averaged ove
r the run or event by event. The step size is independent of bin size. Energy deposition can be weighted by a quenching factor (Birks law). Scoring can be done in a time window. It is possible to simulate coincidences and anti-coincidences.
Fluence and current can be scored as a function of energy and angle, via boundary-crossing, collision and track-length estimators coincident with regions or region boundaries.
Track-length fluence can be scored in a binning structure (Cartesian or cylindrical) independent of geometry.
Particle yield from a target is available, or differential cross section with respect to several different kinematic variables.
Other scoring possibilities include residual nuclei, fission density, momentum transfer density, neutron balance, unweighted energy deposition.
All quantities from radioactive decay of residual nuclei can be scored according to user-defined irradiation and cooling time profiles (decay radiation transport is provided on request).
FLUKA can be used in analog mode or with a variety of variance reduction options. These include: Leading particle biasing for electrons and photons: region dependent, below user-defined energy threshold and for selected physical effects; Russian Roulette and splitting at boundary crossing based on region relative importance; region-dependent multiplicity tuning in high energy nuclear interactions; region-dependent biased downscattering and non-analog absorption of low-energy neutrons; biased decay length for increased daughter production, biased inelastic nuclear interaction length; biased interaction lengths for electron and photon electromagnetic interactions; biased angular distribution of decay secondary particles; region-dependent weight window in three energy ranges (and energy group dependent for low energy neutrons).