What physics list is suitable for BNCT

Hello Geant4 experts:
I am simulating the process of Boron Neutron Capture Therapy (BNCT) in Geant4. Are there any suggestion of the physics list in this work?

In this work I need to simulate the process of nCapture in B10(n,a)Li7 channel, track the dose of 2 heavy fragments, and also track the 478keV gamma from Li7* deexcited to simulate the BNCT-SPECT. Beside B10 nCapture, H(n,gamma)D is also a necessary channel to be simulated, and other light nuclear nCapture as well.

Currently my physics list is :

  // Bosons
  G4BosonConstructor bs;

  // Leptons
  G4LeptonConstructor lp;

  // Mesons
  G4MesonConstructor mc;

  // Baryons
  G4BaryonConstructor br;

  // Ions
  G4IonConstructor in;

  // Short Lived
  G4ShortLivedConstructor shortLived;

  // physics process
  RegisterPhysics(new G4EmStandardPhysics);
  RegisterPhysics(new G4HadronElasticPhysicsHP(verb));
  RegisterPhysics(new G4HadronPhysicsQGSP_BIC_HP(verb));

  RegisterPhysics( new G4IonElasticPhysics(verb));
  RegisterPhysics( new G4IonPhysicsXS(verb));

  RegisterPhysics(new G4StoppingPhysics(verb));         
  RegisterPhysics(new G4RadioactiveDecayPhysics(verb));
  RegisterPhysics(new G4IonBinaryCascadePhysics(verb));
  RegisterPhysics(new G4DecayPhysics(verb));
  // RegisterPhysics(new G4EmExtraPhysics);
  // RegisterPhysics(new G4NeutronTrackingCut);

With the above code I simulated the dose distribution of 25 keV monoenergetic neutron induced, we can obviously find that the dose at entrance is higher than that in the target volume (see this figure from top to bottom which is the depth increasing direction, and target is at the middle depth):

For the physics list code above, are there any suggestions, about which one is necessary and which is not, and any other physics is needed to append? Thanks in advance.

_Geant4 Version:_11.0.1
_Operating System:_windows subsystem for linux (WSL2), ubuntu 22.04 LTS
_Compiler/Version:_gcc version 11.4.0
_CMake Version:_3.22.1


I would suggest QGSP_BIC_HP or QGSP_BIC_allHP with thermal cross sections on .


This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.