QGSP_BIC_AllH and QGSP_BIC_HP in Hadr03

Physics Processes, Models and Cross Sections
1

I have a question!

I would like to use “QGSP_BIC_AllHP” for cross-section data of the C-12(p,X)C-11 and C-12(p,X)C-10 reactions induced by proton irradiation.

When I compared the results of “QGSP_BIC_HP” and “QGSP_BIC_AllHP” in Hadr03, the number of produced isotopes was significantly lower with “QGSP_BIC_AllHP”, and important isotopes such as C-11 and C-10 were not displayed.

What could be the reason for this?

I’m using Geant4-v11.3.0

/testhadr/det/setMat G4_POLYETHYLENE
/run/initialize
/gun/particle proton
/gun/energy 100 MeV
/process/inactivate hadElastic

QGSP_BIC_AllHP
List of generated particles:
C12: 35875 Emean = 558.97 keV ( 10.005 eV → 22.322 MeV)
C13: 392 Emean = 554.29 keV ( 2.38 keV → 10.248 MeV)
alpha: 32303 Emean = 6.7452 MeV ( 11.048 keV → 80.348 MeV)
deuteron: 1 Emean = 7.9866 MeV ( 7.9866 MeV → 7.9866 MeV)
neutron: 32777 Emean = 21.197 MeV ( 1.7471 keV → 88.566 MeV)
proton: 40640 Emean = 31.676 MeV ( 16.662 keV → 99.995 MeV)

QGSP_BIC_HP
List of generated particles:
B10: 778 Emean = 5.8711 MeV ( 7.7599 keV → 27.556 MeV)
B10[2154.270]: 242 Emean = 4.7453 MeV ( 142.52 keV → 27.66 MeV)
B10[718.380]: 795 Emean = 5.1114 MeV ( 41.836 keV → 26.542 MeV)
B11: 1882 Emean = 4.4434 MeV ( 36.204 keV → 23.955 MeV)
B12: 15 Emean = 3.2843 MeV ( 164.19 keV → 11.064 MeV)
B8: 78 Emean = 4.2532 MeV ( 99.819 keV → 21.054 MeV)
Be10: 128 Emean = 7.6735 MeV ( 155.1 keV → 20.563 MeV)
Be6: 196 Emean = 6.6564 MeV ( 152.65 keV → 43.673 MeV)
Be7: 1567 Emean = 4.6486 MeV ( 3.2651 keV → 26.888 MeV)
Be9: 82 Emean = 7.4995 MeV ( 202.2 keV → 28.256 MeV)
C10: 425 Emean = 7.8394 MeV ( 240.19 keV → 25.194 MeV)
C11: 4868 Emean = 4.039 MeV ( 16.469 keV → 28.503 MeV)
C12: 35884 Emean = 631.91 keV ( 26.39 eV → 20.4 MeV)
C13: 375 Emean = 522.43 keV ( 792.72 eV → 5.0173 MeV)
C9: 9 Emean = 7.9307 MeV ( 925.99 keV → 19.976 MeV)
He3: 3864 Emean = 7.7023 MeV ( 15.252 keV → 48.565 MeV)
He6: 45 Emean = 6.023 MeV ( 655.01 keV → 21.474 MeV)
Li6: 1253 Emean = 5.984 MeV ( 8.6493 keV → 41.941 MeV)
Li7: 775 Emean = 4.4629 MeV ( 23.539 keV → 25.532 MeV)
Li8: 8 Emean = 5.9729 MeV ( 649.56 keV → 12.901 MeV)
N12: 33 Emean = 6.0394 MeV ( 1.2463 MeV → 14.869 MeV)
N13: 2 Emean = 5.8985 MeV ( 5.1363 MeV → 6.6607 MeV)
alpha: 40623 Emean = 6.4897 MeV ( 341.03 eV → 66.895 MeV)
deuteron: 5820 Emean = 7.5987 MeV ( 6.087 keV → 63.085 MeV)
gamma: 13193 Emean = 2.4967 MeV ( 45.238 eV → 20.75 MeV)
neutron: 25739 Emean = 21.058 MeV ( 4.157 keV → 91.86 MeV)
proton: 56622 Emean = 30.306 MeV ( 8.1482 keV → 99.992 MeV)
triton: 1420 Emean = 6.6221 MeV ( 19.67 keV → 38.407 MeV)

3

AllHP uses the ENDF and TENDL evaluated databases. There is no ENDF data for C10, C11 (in fact only a few isotopes are present in ENDF) neither in TENDL. C10 and C11 are not present in natural C, so it is not surprising that there is no experimental data and their cross sections are 0.
HP uses its own model, not based on evaluated databases.

4

Thank you for your response!

I believe the data contained in /path/to/geant4-v11.3.0-install/share/Geant4/data/G4TENDL1.4/Proton/Inelastic/F02 corresponds to (p,X) reactions.

When I checked inside 6_12_Carbon.z, it seemed to contain the same data as shown in the following link (TENDL-2021):
https://www-nds.iaea.org/exfor/servlet/E4sGetIntSection?SectID=16394704&req=75674&e4up=0&PenSectID=19781785&pen=0

I believe that the “AWP : 1.091641+1” corresponds to C-11.

Therefore, if C-12 is present, I would expect C-11 to be produced…

5

Right, this line marks the start of the C11 XS data
6.011000e+03 1.091641e+01 0 1
and comparing the numbers to other products, the XS are negligible. Same for C10

6

Thank you!

Indeed, with 100 MeV proton irradiation, the neutron yield is 7.194589e−1, and C-11 yield is 4.972447e−2, which is about 1/15th of the neutron yield. While small, It doesn’t seem negligible.

In fact, although the model is different, when using “QGSP_BIC_HP”, particles with yields differing by about 100 are still displayed.

Is “QGSP_BIC_AllHP” set to exclude products when the cross sectional area falls below a certain threshold?

Additionally, is it possible to use other ENDF in Geant4?