Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!
_Geant4 Version:_Geant4 10.7.2
_Operating System:_Ubuntu 22.04
Compiler/Version:
CMake Version:
Hi all！
I’m currently running a Geant4 simulation for PuO₂, focusing on the spontaneous fission of Pu-240, using the official RadioactiveDecay and G4SFDecay models. My goal is to analyze the energy and angular distributions of neutrons emitted from spontaneous fission.And I use FREYA2.0 to generate the energy and emitted direction of spontaneous fission neutron.
However, I consistently encounter the following fatal error sometimes:
-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had001
issued by : G4CrossSectionDataStore::GetIsoCrossSection
No isotope cross section found for neutron off Element Pu240 in PuO2 Z= 94 A= 240 E(MeV)= 27.3642
*** Fatal Exception *** core dump ***
G4Track (0x626f7a41b0d0) - track ID = 12, parent ID = 1
Particle type : neutron - creator process : RadioactiveDecayBase, creator model : ionIoni_subcut
Kinetic energy : 27.3642 MeV - Momentum direction : (-0.0204457,-0.593415,0.804637)
Step length : 0 fm - total energy deposit : 0 meV
Pre-step point : (1.36457,2.31086,2.95249) - Physical volume : Pu_phys (PuO2)
- defined by : not available
Post-step point : (1.36457,2.31086,2.95249) - Physical volume : Pu_phys (PuO2) - defined by : hadElastic - step status : 7
*** Note: Step information might not be properly updated.
-------- EEEE -------- G4Exception-END --------- EEEE -------
*** G4Exception: Aborting execution ***
Aborted (core dumped)
- Is this issue due to the neutron cross section data for Pu-240 being limited to below a certain energy (e.g., <20 MeV)? If so, is there a recommended way to handle spontaneous fission neutrons with higher energies in Geant4?
- Is there a way to automatically force Geant4 to use the highest available cross section value if the neutron energy is out of range (as some documentation suggests for certain models)?
Thanks in advice!