|Message: Alpha energies in Radioactive Decay are wrong!||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)|
It appears that all of the alpha energies in the "radioactive decay" are calculated incorrectly.
This was noticed for a handful of isotopes.
As an example, I've been trying to model the decay of 221Th. The kinetic energy of an alpha particle in g.s.->g.s. alpha decay should be 8470 keV.
To check this number, we took the atomic masses from AME2003:
M_a(221Th) = 205877.112 MeV
M_a(217Ra) = 202140.085 MeV
M_(alpha) = 3727.379 MeV
And derived the Q-value:
Q = 205877.112e3 - 202140.085e3 - 3727.379e3 - 2*511 = 8626 keV =>
E_alpha = 8626/(1+4./217) = 8470 keV
This agrees to the measured alpha energy (see, for example, https://www-nds.iaea.org/relnsd/vcharthtml/VChartHTML.html - 8470(5) keV)
Still, Geant4 generates alpha particles of 8509 keV for this transition.
We've tried to reverse-injeneer what's happening in Geant4 in order to find what's wrong.
Judging by the fact that G4 always returns the masses of bare nuclei when asked for the total energy of a nucleus at rest, we supposed that the NUCLEAR masses (atomic mass - electrons + binding energy) may have been used for the calculation of the alpha energy INSTEAD OF ATOMIC.
We have used the following approximation formula to deduce the M_n
M_n(217Ra) = 202095.797 MeV
M_n(221Th) = 205831.843 MeV
So, if indeed the nuclear masses were by mistake taken instead of the atomic ones, and the whole lot of the electrons was ignored:
"Q" = 205831.843e3 - 202095.797e3 - 3727.379e3 = 8667 keV =>
"E_alpha" = 8510 keV
As Geant4 proposes alpha-energy of 8509 keV instead of 8470, I assume that indeed there is a mistake and an atomic mass is replaced by the nuclear mass. Same calculations can be carried out for various nuclides.
It is also wrong to use the nuclear masses, because practically none of them were measured directly, but are derived from the estimates similar to the one I've attached. It's the atomic masses that are measured.
Also, it is wrong that an object "ion" has a mass such as a bare nucleus. That's not what's happening in real life, they are atoms…! Fully stripping a nucleus like this is a challenge.
I failed to find where the masses are tabulated in Geant4. Could you tell me where is this? I guess, one could modify those files?
Are the tabulated masses the nuclear ones or the atomic ones (such as from AME)? If latter, where is the function that converts them?
|Inline Depth:||Outline Depth:||Add message:|