Hello everyone, while studying chem6, I came across the following code in the macro file beam.in:
It seems that the author intended to use the electronic solvation subtype model. However, in the literature I have read, it was mentioned that this model is not suitable for the chemical reaction stage.The title of the reference paper is “Geant4-DNA example applications for track structure simulations in liquid water: A report from the Geant4-DNA Project”.
Since the addition of this model has little impact on the output results, I would like to understand the effect of adding the “e-SolvationSubType” model on the program results and whether there is any information available about this model that can be referred to.
I wrote this paragraph to add some pictures
Sorry for the misunderstanding! The paper says that this process does not apply for the physics constructors and not included in the table 1. This process is included only when chemistry simulation is activated (in the chemistry constructors)
Hello, I couldn’t find any reference to the e-SolvationSubType, neither in the BookForApplicationDevelopers nor online. It appears that people pay less attention to this option in their models. Do you happen to know how to find more information about this particular model?
Thank you so much, this paper help me a lot!
I really hope that the paper can be cited in the readme.txt of the example code in future release
but I still have some questions about the setting of parameters that haven’t discussed in the paper. For instance, in the primaryKiller setting, T2 seems to be 1% larger than T1. What is the reason for specifying this value? Also, when setting the specified energy for incident particles, why is the first one 2 keV and the next one 3.5 keV without following a consistent pattern?
Thank you for your assistance!
here are the table in the paper,and the code of particle number
I speculate that the setting of this value is related to Linear Energy Transfer (LET), but how exactly are they correlated?
the G-value we calculate is called track segment G-value. This means we just count G-value for a piece of the particle track. At different LETs, this piece is different. The reason of this way to calculate G-value is time consuming. Please do some simulation to see the difference.