When running neutron-hp with the $G4NEUTRONHP_USE_ONLY_PHOTONEVAPORATION flag on, captures on hydrogen yield two deuterons rather than a single deuteron and a gamma.
Thank you for reporting this problem. It looks like that the problem itself is in the PhotonEvapoation model of Geant4. The model does not produce (2.2MeV) gamma from capturing nucleus. However I made an internal fix within HP, therefore this problem should not appear in final states. Tatsumi