Appropriate Physicslist for Neutron production

Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!

_Geant4 Version:_11.01.p03
_Operating System:_Sonoma
Compiler/Version:
_CMake Version:_3.28.3

Dear Geant4 Users,
I hope you are all doing well.

I have a question regarding the appropriate PhysicsList for simulating high-energy gamma interactions with metals (e.g., lead), particularly for studying the production of secondary neutrons.
I have searched the forum, but I am still uncertain about the most suitable PhysicsList for this purpose. Below is a simple simulation where a 20 MeV gamma particle is incident on lead (cyan), and I am checking whether secondary neutrons reach the scoring box (red).

When using the commonly recommended “G4HadronPhysicsQGSP_BIC”, no neutrons are produced. However, when I use “GammaNuclearPhysics(“gamma”)”, I can observe neutron production, as shown in the attached image.

Energy, Particle name, Process name

So, my main question is:

  • What is the most appropriate Physics List for accurately simulating secondary neutron production from high-energy gamma interactions with metals?
  • Am I missing any important considerations regarding the Physics List selection?

I would appreciate any insights or recommendations from experienced users.

Best,
yHuh