In the Geant4 simulation of the D-Be reaction to calculate neutron yields, the simulation results differ from the reference values by two orders of magnitude

Physics Processes, Models and Cross Sections
2

Hello,

Here are some thoughts on your problem:

  • you are already using the NIST-Manager for defining world_mat;
    you can use it to define Beryllium too - your definition of Beryllium looks right, but NIST manager will defintely do it right! This should do the trick:

G4Material* Be = nist->FindOrBuildMaterial( “G4_Be” );

  • Physics list looks okay too - G4HadronPhysicsQGSP_BIC_AllHP should, be the way to go for Be+d. In newer Geant4 versions, this physics list uses the G4TENDL datasets. Your Geant4 version (10.5.1) is 5 years old and there isn’t even a G4TENDL dataset to download on the download page for your version . I strongly recommend updating since these kind of reactions (deuterons on light ion target) seem to be very sensitive to Geant4 version and Dataset version as you can read about in this post.
  • How do you track/count the neutrons? Maybe it’s a counting problem and the right amount of neutrons are produced but you dont detect them.
  • What experimental values are you referring to? Is there an experiment in a paper you are trying to reproduce?

Regards
Max

1 Like