|Message: Re: How to add thermal neutron scattering to code?||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)|
I’m a colleague of Florian and I have some new insights into the problem.
The negative-energy errors occurs only for inelastic thermal neutron scattering (G4NeutronHPThermalScattering.cc lines 447-526 for Geant4 10.01, same lines in G4ParticleHPThermalScattering.cc for Geant4 10.04p01) and it is the result of the linear interpolation/extrapolation of the data for sE (line 501). I notice that in our use case the temperature of the used polyethylene is with 50mK much smaller than the lower edge of the interpolation interval which is 264K.
This behavior can be reproduced with the attached macro run1.mac and examples/extended/hardonic/Hadr06 of Geant4 10.04p01 by changing the temperature of ‘water_ts’ from 593K to 50mK on line 92 in DetectorConstruction.cc. With the original temperature of 593K the error does not occur.
My assumption is that the extrapolation breaks down for temperatures far away from the tabulated data. However, I could not find any details about the applicability range of the thermal neutron scattering in the physics reference manual version 10.4.
Can anyone with more knowledge about the thermal neutron scattering comment if my assumption is plausible? Does anyone know the applicability range of the thermal neutron scattering model? Are there any plans to extend them?