[Paper Review] The physics models of FLUKA: status and recent development

A. Fassò, A. Ferrari+10 more|ArXiv.org|Jun 27, 2003

High-Energy Particle Collisions ResearchPhysics and Astronomy126 citations

TL;DR

This paper presents the physics models underlying the FLUKA Monte Carlo code, emphasizing its microscopic, non-tuned approach to modeling hadron-nucleon and nucleus-nucleus interactions from low to ultra-high energies. It details the integration of advanced models like RQMD and QMD for ion collisions, validates performance against experimental data, and outlines future coupling with preequilibrium and Boltzmann-based codes for improved fragmentation and low-energy ion transport.

ABSTRACT

A description of the intermediate and high energy hadronic interaction models used in the FLUKA code is given. Benchmarking against experimental data is also reported in order to validate the model performances. Finally the most recent developments and perspectives for nucleus-nucleus interactions are described together with some comparisons with experimental data.

Motivation & Objective

• To provide a comprehensive overview of the physics models implemented in the FLUKA Monte Carlo code for particle transport up to 10,000 TeV.
• To validate the accuracy of FLUKA’s intermediate- and high-energy hadronic interaction models through benchmarking against experimental data.
• To describe the recent integration of ion interaction models, replacing the outdated superposition model approximation for nucleus-nucleus collisions.
• To present the development and validation of a new QMD-based model for simulating A–A collisions, including initial nuclear configurations and dynamical evolution.
• To outline future integration paths with preequilibrium and Boltzmann Master Equation codes for improved low- and intermediate-energy ion transport.

Proposed method

• FLUKA employs microscopic, self-consistent models for hadron-nucleon interactions, avoiding tuning to integral observables and relying on fundamental conservation laws.
• Intermediate-energy hadron-nucleon interactions are modeled via the isobar model, with resonances such as Δ(1232) described through Breit-Wigner formalism and phase-shift fits.
• High-energy inelastic processes are simulated using quark-parton string models and the Dual Parton Model (DPM), with leading two-chain diagrams used to describe multi-particle final states.
• For nucleus-nucleus collisions, the superposition model has been replaced by a new ion interaction model based on the RQMD-2.4 code, which treats full A–A dynamics.
• A new QMD code is developed from scratch, using Gaussian wave packets for nucleons, a Skyrme-type non-relativistic Hamiltonian with surface and symmetry terms, and Coulomb repulsion.
• Initial nuclear configurations are selected based on stability, density, and momentum distribution, with only those avoiding spurious nucleon emission being stored for simulation.

Experimental results

Research questions

• RQ1 How accurately do FLUKA’s microscopic models reproduce experimental data for hadron-nucleon interactions across a broad energy range?
• RQ2 What improvements does the new RQMD-based ion interaction model bring over the previous superposition model approximation in nucleus-nucleus collisions?
• RQ3 How well does the newly developed QMD code reproduce ground state nuclear properties and dynamical evolution in A–A collisions?
• RQ4 What are the key physical features—such as rms radius and density oscillations—observed during the time evolution of nuclear configurations in the QMD model?
• RQ5 How can the coupling of FLUKA with preequilibrium and Boltzmann Master Equation codes enhance the description of fragmentation and low-energy ion interactions?

Key findings

• The FLUKA code successfully reproduces experimental data for double differential neutron production at 400 MeV/n and fragment production at 1.05 GeV/n using the modified RQMD-2.4 model.
• The root mean square radii of $^{20}$Ne, $^{40}$Ca, and $^{90}$Zr nuclei exhibit small oscillations over time due to finite Fermi momentum and dynamic motion, indicating non-static initial configurations.
• The radial density profiles of $^{20}$Ne nuclei show oscillations around mean values during time evolution, reflecting dynamic nuclear responses under internal forces.
• The new QMD model achieves reasonable reproduction of nuclear ground state properties using a few parameters, with Gaussian widths increasing with mass number to improve accuracy across the nuclear chart.
• Initial configurations are carefully selected to avoid spurious nucleon emission and to maintain correct density and momentum distributions over hundreds of fm/c.
• The integration of the new ion interaction model renders the superposition model obsolete, enabling more accurate simulation of A–A collisions at high energies.