Geant4 Version: 11.2.2
Operating System: Redhat Linux 7.9
Compiler/Version: gnu/12.4.0
CMake Version: 3.25.2
Hi,
I am simulating the neutron yield of deuterons hitting a 3mm thick Beryllium target. I’ve tried a few physics lists, with an energy scan ranging from 1 MeV to 100 MeV. However, I observed an abrupt change in neutron yield when using the QGSP_BIC_AllHP physics list. What could be the reason behind this?
I initially thought it might be due to models or cross-section data being switched. However, upon examining the log messages, it doesn’t seem to support my expectation. By the way, this abrupt change is absent from proton spallation. But they won’t let me post two pictures at one time.
Hadronic Processes for deuteron
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: dInelastic
Model: ParticleHPInelastic: 0 eV /n ---> 200 MeV/n
Model: Binary Light Ion Cascade: 190 MeV/n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: DeuteronHPInelasticXS: 0 eV ---> 200 MeV
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Hadronic Processes for proton
Process: hadElastic
Model: hElasticCHIPS: 0 eV ---> 100 TeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: protonInelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: Binary Cascade: 190 MeV ---> 6 GeV
Model: ParticleHPInelastic: 0 eV ---> 200 MeV
Cr_sctns: ProtonHPInelasticXS: 0 eV ---> 200 MeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Which physics list would be most appropriate for this situation? I wanted to use QGSP_BIC_AllHP because it seem have to adopted finer physics for neutrons and protons at lower energy and there is one paper and another suggest doing so. But the abrupt change behavior got me hesitant.
