I’m trying to produce normalized differential fluence spectra within a simulated aluminum-PNTD stack and a 975 MeV/n 28-Si beam. I haven’t done this since 2010. Going back to it in 2024, I’m not finding the heavy ions fragmenting as they should:
The above is behind 20 g/cm2. I’ve attempted to modify my input file to align it with modern FLUKA standards with no success (found some suggestions here). Please find attached my old and new input files.
How would one implement heavy ion fragmentation in 2024? Many thanks in advance.
You just forgot to use the standard executable embedding the RQMD library (flukadpm), which allows FLUKA to perform ion fragmentation in the concerned energy range:
The spectrum you produced is correct. I don’t understand why I’m not producing the same. The documentation states that FLUKA no longer uses environment variables, so that cannot be the issue.
This is mysterious. I’m using the fluka_4-4.1.x86-Linux-gfor9_amd64.deb build. Any further suggestions will be greatly appreciated.
UPDATE: Running from the command line, I can confirm that the regular fluka exec is being used:
If you launch the simulation from Flair, you have to select the <flukadpm> executable in the Exe field of the Run tab.
If you launch it from the command line, you have to add as argument of rfluka
-e /usr/loca/fluka/bin/flukadpm
The output file (*001.out) will tell you if RQMD (together with DPMJET) has been actually linked:
—> requested particle production model RQMD
is the crucial record.
I ended up running my input file from the command line last night via /usr/local/fluka/bin/rfluka -M 1 -d BB5_Si_Al.inp, forcing it to use flukadpm with the -d flag. Success:
