I am simulating a sphere of mainly U233 and H1 in which I generate neutrons in the centre of the sphere. I execute Geant4.10.7 on Ubuntu 20.04.
Thermal scattering is activated with TS_H_of_Water.
The temperature of my material is 250 kelvin. The physics list is QGSP_BIC_HP and I use G4NDL4.6 neutron data.
I tried to set DONOTADJUSTFINALSTATE on and off but both lead to the same results
When I execute my code, I get an error in which a negative energy has been found:
G4WT0 > Final state energy was: E = -1.65487e-10
G4WT0 >
-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had006
issued by : G4HadronElasticProcess::PostStepDoIt
In /usr/local/geant4.10.07/source/processes/hadronic/util/src/G4HadFinalState.cc, line 45:
===> G4HadFinalState: fatal - negative energy
Call for NeutronHPThermalScattering
Target element TS_H_of_Water Z= 1 A= 1
Unrecoverable error in the method ApplyYourself of hadElastic
TrackID= 1 ParentID= 0 neutron
Ekin(GeV)= 4.36253e-13; direction= (0.751898,-0.644773,-0.137543)
Position(mm)= (-30.5531,12.0646,-116.79); material uranium
PhysicalVolume
ApplyYourself failed*** Fatal Exception *** core dump ***
G4WT0 > G4Track (0x14bde227a450) - track ID = 1, parent ID = 0
G4WT0 > Particle type : neutron - creator process : not available
G4WT0 > Kinetic energy : 0.000436253 eV - Momentum direction : (0.751898,-0.644773,-0.137543)
G4WT0 > Step length : 1.54444 mm - total energy deposit : 0 eV
G4WT0 > Pre-step point : (-31.7144,13.0604,-116.578) - Physical volume : USphere2 (uranium)
G4WT0 > - defined by : hadElastic - step status : 4
G4WT0 > Post-step point : (-30.5531,12.0646,-116.79) - Physical volume : USphere2 (uranium)
G4WT0 > - defined by : hadElastic - step status : 4
G4WT0 > *** Note: Step information might not be properly updated.
G4WT0 >
-------- EEEE -------- G4Exception-END --------- EEEE -------G4WT0 >
G4WT0 > *** G4Exception: Aborting execution ***
Is there a reason for TS to lead to this problem at low temperature ?