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: 11.1.1
Operating System: Windows 7 Professional
Compiler/Version: Visual Studio 2022
CMake Version: 3.25.3
On a semi-related occurrence to another one of my posts (DNAPhysics and RadioactiveDecayPhysics not simulating relaxation of excited daughter isomer), I have noticed a problem with the atomic relaxation of the daughter product of Tb-161, which decays by beta minus emission.
In the below screenshot you’ll see the particles produced by G4RadioActiveDecay (simulated in vacuum, tested with both the dna_opt2 and std_opt3 physics lists). In 9 out of 10 runs, the Dy-161 remains in an excited state and does not relax to the ground state - this is crucial for the simulation of Auger electron emissions of this radionuclide.
I have even tried to force its instantaneous decay (which is why the mean life time displays as 0 ps) by increasing G4NuclideTable::GetInstance()->SetThresholdOfHalfLife to above is actual mean life of 4.53 ns, but still only 1 out of 10 times is the relaxation modeled.
Am I missing something? I use the exact same code for other radionuclides (e.g. I-123, I-125, I-131, etc) and there’s not problem.
