|Message: Reaction kinematics for neutrons and Lithium6: questions and a hard patch||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
This post is an attempt to tie together the unresolved issues surrounding the incorrect G4NeutronHPInelastic nuclear reaction kinematics for the 6Li(n,t)4He reaction (Q=4.78 MeV). This issue has appeared previously in posts #1032(hadronic), #461(hadronic), #530(physics). As a reminder, G4NeutronHPInelastic is unable to correctly partition the kinetic energy to the triton and alpha in the 2 body final state. While it seems that the momentum is conserved (triton and alpha are emitted with equal and opposite momenta direction), Geant4 attempts to conserve energy by producing non-physical particles, usually gamma rays. Using geant4.9.3 + QGSP_HP_BERT (Fedora12, gcc 4.4.2) to model the 6Li(n,t)4He reaction, G4NeutronHPInelastic produces:
1) triton with 1.56 MeV; should be 2.730 MeV
2) alpha with 1.17 MeV; should be 2.050 MeV
3) 3 non-physical gammas (2 with 1.03 MeV, 1 with 0.00212 MeV )
It should be noted that similar problems exist for the 10B(n,a)7Li reaction, although the kinetic energies of the charged products are off by only ~30 keV in these cases.
Questions: 1) It has been suggested by Dr. Koi that this is due to ENDF.VI classification of these reactions, and Geant4's incorrect ability to simulate them correctly. Is there a patch forthcoming in the near future? 2) Will this issue be resolved when Geant4 moves to ENDF.VII or ENDL based libraries? Is there any progress on either of these fronts?
For those interested, I have developed a "hard" patch in the form of a new G4VDiscreteProcess that can be attached to the triton/alpha/gamma in the user's physics list that resolves the issue. The patch:
1) Detects triton/alpha/gamma from the 6Li(n,t)4He reaction
2) Correctly sets the KE energy of triton/alpha after a forced 0.1 nm step
3) Kills the gamma rays
4) Does not interfere with particle tracking after the first step
I have done some initial work to validate the patch for the case of thin-film based thermal neutron semiconductor detectors. Geant4 simulations have been validated against analytic calculations presented by McGregor et al, NIM A 500 (2003) 272-308. The initial results are excellent, matching calculated detector efficiencies quite well with no adverse effects in Geant4.
I am sure that this is not the best or "correct" way to patch this problem, but it seems to be a good working solution. For those interested, I have attached the header and implementation files, along with initial validation results. This is certainly a work in progress, so please use caution. Any comments on improving the patch are quite welcome.
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/12/21/22.02-4415-ronInelasticCorrection.cc http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/12/21/22.02-95131-ronInelasticCorrection.hh http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2009/12/21/22.02-76286-Li6CorrectionSummary.pdf