I’m trying to simulate the thermal neutron capture of Gd. This reaction produces Gd metastable and I would like to see the internal conversion electrons released when it returns to the ground state.
The problem I’m currently facing is that I don’t see any metastable states after neutron capture. I’m using the QGSP_BIC_HP physics list, in which I believe the G4RadioActiveDecay process is active, and version 10.6.1.
I’ve tried the example in extended/radioactivedecay/rdecay01. There I can put for example Gd158[233.1] and I see IC electrons. So, the problem seems to be that the nCapture process doesn’t create metastable states of Gd.
Should I add something to my physics list or do you know any other options?
Any advise would be very useful! Thanks!