60 MeV p on Tantalum QGSP_BIC_AllHP errors

Physics Processes, Models and Cross Sections
1

Hello,

I am using the latest Geant4 release. During a simple simulation of 60 MeV protons on a small Tantalum box, using the QGSP_BIC_AllHP physics list there are some most unusual error messages:

// **************************************************************************************************
-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
Process / Model: protonInelastic / ParticleHPInelastic
Primary: proton (2212), E= 980.133, target nucleus (73, 181)
E(initial - final) = -656318 MeV.

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
Process / Model: protonInelastic / ParticleHPInelastic
Primary: proton (2212), E= 964.286, target nucleus (73, 181)
E(initial - final) = -665658 MeV.

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
Process / Model: protonInelastic / ParticleHPInelastic
Primary: proton (2212), E= 983.524, target nucleus (73, 181)
E(initial - final) = -504617 MeV.

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
Process / Model: protonInelastic / ParticleHPInelastic
Primary: proton (2212), E= 993.723, target nucleus (73, 181)
E(initial - final) = -498092 MeV.

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
Process / Model: protonInelastic / ParticleHPInelastic
Primary: proton (2212), E= 995.345, target nucleus (73, 181)
E(initial - final) = -656454 MeV.

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
Process / Model: protonInelastic / ParticleHPInelastic
Primary: proton (2212), E= 984.615, target nucleus (73, 181)
E(initial - final) = -667519 MeV.

*** This is just a warning message. ***
// **************************************************************************************************

I know these are meant to be only warning messages, but keeping in ming that the projectiles are 60 MeV protons, these messages are talking about differences in energy of ~67 GeV and proton energies of 984 (MeV?)…
Also most simulation runs die without a succesfull output file creation.

I have tried also QGSP_INCLXX_HP, where the error messages are more sensible, but also present in a large number.

I would be very grateful for any suggestions.

Thank you.

2

With example Hadr03, QGSP_BIC_AllHP selected, I see these warnings.
But the list of nuclear reactions is strange, even absurd … and can be the cause of the warnings.
I use latest Geant4 development version : ref-05, 31-May-2023.
Unless ‘responsibles’ of this code react quickly, the situation will be unchanged in Geant4 11.2-beta

cbungau.mac.txt (280 Bytes)
cbungau.out.txt (10.7 KB)

3

Hello,

for the time being the best advice would be to use QGSP_BIC_HP or QBBC for low-energy proton physics. QGSP_BIC_AllHP is still experimental and required more work to address issues.

VI

4

Hello,

Thank you - I will switch to BIC_HP then.

Best wishes,
Cristian