_Geant4 Version:_11.3.0
_Operating System:_Ubuntu 22.04.5
Hello everyone,
I’m currently working on a simulation involving a gaseous deuterium (D₂) target bombarded with deuterons, with the goal of modeling a neutron source. Specifically, I’m interested in quantifying the relative yields of two deuteron-deuteron reaction channels:
- d + d → n + ³He
- d + d → d + n + p (deuteron breakup)
My main issue is that so far, I do observe the n + ³He reaction occurring in the simulation, but I do not see any evidence of the breakup reaction channel (d + d → d + n + p) being triggered.
I’m currently using the QGSP_BIC_AllHP physics list. I would like to know:
- Does Geant4 support simulation of this specific breakup channel using any available hadronic or ion models?
- Would using a different physics list (or customizing one) help include this reaction?
- Is there any recommended approach to simulate such three-body breakup processes (possibly through a custom model or external generator) ?
Any insights or suggestions on how to accurately include this reaction channel would be greatly appreciated.
Thank you very much in advance for your help!