PKA spectrum calculation

Versions

FLUKA: 4.5.1
Flair:
Operating system:
Compiler:

Description

I am trying to create a PKA spectra for the whole region (simple thin slab) and was curious if running the model several times with arcdpa parameters, c=0, b= -1 for different threshold displacement energies can give me dpa values that I can process to get a spectra?

Are there better ways of doing this such as extracting from a track file or using a subroutine ?

Input files

Please upload all relevant files. (FLUKA input file, Flair project file, user routines, and data files)

Dear Sanjeeb Pal,

You can score the number of displaced atoms using the generalized particles DPA-*. You can use a USRBIN card with region scoring to get the average number per region, after region volume normalization. You can change the threshold with a MAT-PROP card and run simulation with different thresholds to perform a threshold scan. I suppose you would then obtain what you wish?

As a complementary clarification, DPA scoring will not give you at all a PKA energy spectrum. To obtain the latter, you need to make use of the mgdraw user routine and identify the residual nucleus out of any nuclear reaction and the respective kinetic energy, as discussed for instance here.