We have seen similar effects, using Geant4 10.7.4 with G4ParticleHP, with neutron activation of our germanium detectors, i.e., Ge-70(n,gamma)Ge-71. We thought the problem was with RadioactiveDecay, since the gammas we see are consistent with the Ga-71 X-rays following Ge-71 electron capture. But the problem is not reproducible with rdecay01.
I agree that this is likely related to PhotonEvaporation. Tracing back from my post above, it’s possible that setting
export G4NEUTRONHP_USE_ONLY_PHOTONEVAPORATION=1
or
/run/initialize
/process/had/particle_hp/use_photo_evaporation true
might eliminate the duplicated gammas, at the cost of not getting the evaluated data intensities for the gamma lines. If your simulation only needs energy conservation, and not exact spectra, this may be an option.