|Message: Re: Does Geant4 model proton boron fusion reaction?||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)|
Many thanks for the reply.
I incorporated the physics list from Hadr06 and activated QGSP_BIC_HP based on the information you provided. However, still absorbed dose in a boron layer in my geometry is very low compared with published data obtained from MCNPX.
I then tried running example Hadr03 for a proton beam of 80 MeV in a phantom with B10 or B11 material: /control/verbose 2 /run/verbose 2 #
/testhadr/det/setIsotopeMat B10 5 10 2.08 g/cm3 /testhadr/det/setSize 10 cm # /run/initialize # /gun/particle proton /gun/energy 80 MeV # /process/list /process/inactivate hadElastic
# #/run/printProgress 10000 # /run/beamOn 100000
The output is attached. While the number of alpha particles produced are much larger when water material is used, I am pretty sure I am not picking up proton boron fusion reaction and that's why I see discrepancy with MCNPX results. After the proton reacts with the boron (11B), the boron changes to carbon (12C) in an excited state. The excited carbon nucleus is split into alpha particle of ~3.76 MeV and beryllium (8Be). 8Be is then divided into the two alpha particles of ~2.74 MeV each.
Could you please give me some feedback, how I could get fusion reaction?
Many thanks! Leila