I am trying to reproduce the experimental total cross section for production of Zn 63 and Zn 65 due to p,n reaction on nat Cu target.
Simulation result:
for proton energy 4.32 MeV, sigma (Zn63): 45.82 mb and sigma (Zn 65): 49.595 mb
for proton energy 4.63 MeV, sigma (Zn63): 75.62 mb and sigma (Zn 65): 70.43 mb
The table II in the attached pdf shows large differences in sigma (Zn63) and sigma (Zn 65) for these above energies.
I am attaching the input file for your reference as well.
The point is that the experimental cross sections are (naturally) given on 63Cu and 65Cu, respectively, and not natural copper (despite the fact that natural copper was actually used). This means that you should more conveniently run two separate simulations, with pure isotopic composition of 63Cu and 65Cu, respectively.
This way, you should rather find about 65 mb (63Zn) and 160 mb (65Zn) at 4.32 MeV and 110 mb (63Zn) and 230 mb (65Zn) at 4.63 MeV. The discrepancy with the reported measurements can be easily understood looking at the huge gradient displayed by the excitation functions in Figs. 2 and 3 up to about 6 MeV.
I gave a run with pure Cu 63 and I got 81.55 mb (63 Zn) for 4.32 MeV; and 95.63 mb for 4.63 MeV. The statistical error is below 1%. Then is there any reason for getting slightly different value than yours ?
The problem is a miscommunication glitch between Flair and FLUKA, triggered by the use of AIR as alternate material for ionisation processes and affecting your material selection in the USRYIELD card, which turns out to be AIR instead of single isotope copper. You can realize this by looking at the output file, where the material index retained in the last field of the USRYIELD card actually corresponds to the air one (as in the material list). As temporary solution before dedicated Flair fixing, please move the AIR definition before the copper definition.
For 65Cu(p,n)65Zn, FLUKA gives 450 mb (vs 422 mb) at 6.16 MeV and 550 mb (vs 492 mb) at 7.23 MeV.
For 63Cu(p,n)63Zn, 92 mb and 142 mb, respectively. Here the energy gradient looks more pronounced and the discrepancy lowers below 30% above 9 MeV.