严南,周云,孙向上,廖文涛,刘俊雅,蒲越虎.敏感域的形状选取对微剂量动力学模型参数的影响[J].中华放射医学与防护杂志,2023,43(3):209-215
敏感域的形状选取对微剂量动力学模型参数的影响
Influence of selection of domain shape on parameters of microdosimetric kinetic model
投稿时间:2022-09-06  
DOI:10.3760/cma.j.cn112271-20220906-00364
中文关键词:  微剂量动力学模型  剂量加权平均比能  敏感域  传能线密度
英文关键词:Microdosimetric kinetic model  Dose-averaged specific energy  Domain  Linear energy transfer
基金项目:国家重点研发计划(2016YFC0105400)
作者单位E-mail
严南 中国科学院上海应用物理研究所, 上海 201800
中国科学院大学, 北京 100049 
 
周云 中国科学院上海应用物理研究所, 上海 201800
中国科学院大学, 北京 100049 
 
孙向上 中国科学院上海应用物理研究所, 上海 201800
中国科学院大学, 北京 100049 
 
廖文涛 中国科学院上海应用物理研究所, 上海 201800
中国科学院大学, 北京 100049 
 
刘俊雅 中国科学院上海应用物理研究所, 上海 201800
中国科学院大学, 北京 100049 
 
蒲越虎 四川大学华西医院医学装备创新研究中心, 成都 610041 puyuehu@163.com 
摘要点击次数: 1054
全文下载次数: 525
中文摘要:
      目的 基于粒子径迹结构模型,利用三重积分计算粒子在球状敏感域(domain)内的单次事件剂量加权平均比能,并探讨敏感域的形状对微剂量动力学模型(MKM)参数带来的影响及其对应的物理意义。方法 分别假定敏感域为圆柱状和球状。α0、域半径rd和细胞核半径Rn为待定系数,3种带电粒子(3He、12C、20Ne)的核电荷数、动能及其对应的传能线密度(LET)为自变量,D10为因变量。以D10计算值与实验值的残差均方值J2为优化目标,采用稳健最小二乘法分别得到人类唾液腺肿瘤(HSG)细胞和中国仓鼠肺(V79)细胞对应待定系数的最优拟合值即为MKM最优模型参数值。结果 对于HSG细胞,圆柱状敏感域:α0=0.073/Gy,rd=0.29 μm,Rn=4.1 μm,J2=0.039 7 Gy2;球状敏感域:α0=0.023/Gy,rd=0.29 μm,Rn=4.4 μm,J2=0.039 3 Gy2;对于V79细胞,圆柱状敏感域:α0=0.114/Gy,rd=0.25 μm,Rn=3.8 μm,J2=0.097 4 Gy2,球状敏感域:α0=0.095/Gy,rd=0.26 μm,Rn=4.1 μm,J2=0.096 9 Gy2结论 对于同一种细胞,分别选取圆柱状和球状的敏感域,最终计算拟合得到的MKM参数存在明显差异,其中两种形状的敏感域半径拟合值rd相差不大,而球状敏感域拟合得到的α0更小,细胞核半径Rn更大,更接近于荧光显微镜观察的细胞核尺寸。在低LET (<20 keV/μm)区域,根据两种形状敏感域所得参数计算的D10存在明显差异,所以敏感域形状选取会对质子放疗在布拉格峰附近区域的相对生物效应(RBE)计算造成影响。
英文摘要:
      Objective To calculate the single-event dose-averaged specific energy of particles delivered in spherical domains based on the track structure model and using triple integration, and to investigate the influence of the domain shape on the key model parameters of microdosimetric kinetic model (MKM) and its corresponding physical significance.Methods The domains are assumed to be cylinders and spheres, respectively. With α0, domain radius, rd, and nucleus radius, Rn, as undetermined coefficients, the nuclear charge numbers, kinetic energies and their corresponding LETs of three kinds of charged particles (3He, 12C, 20Ne) as independent variables, D10 as dependent variable, the mean value of squared residuals, J2, between the D10 calculated values and D10 experimental values as the optimization objective, the final fitting values of the above undetermined coefficients of human salivary gland (HSG) cells and Chinese hamster lung (V79) cells obtained after iteration by the robust least square method are the optimal model parameter values of MKM.Results For HSG cells, cylindrical domain:α0=0.073/Gy, rd=0.29 μm, Rn=4.1 μm, J2=0.039 7 Gy2; spherical domain:α0=0.023/Gy, rd=0.29 μm, Rn=4.4 μm, J2=0.039 3 Gy2; For V79 cells, cylindrical domain:α0=0.114/Gy, rd=0.25 μm, Rn=3.8 μm, J2=0.097 4 Gy2; spherical domain:α0=0.095/Gy, rd=0.26 μm, Rn=4.1 μm, J2=0.096 9 Gy2.Conclusions For the same type of cells, cylindrical and spherical domains were selected respectively, and there are significant differences in MKM parameters obtained by fitting. The fitting values of the domain radius, rd of the two shapes of domains show no significant difference, while the fitting values of α0 of spherical domains are smaller than those of cylindrical domains, the fitting values of nucleus radius, Rn, of spherical domain are larger than those of cylindrical domains, closer to the nucleus radius observed by fluorescence microscopy. In the low LET (<20 keV/μm) region, D10 calculated according to the parameters of the two different shapes of domains are different, so the selection of the domain shape will cause differences in the relative biological effectiveness(RBE) calculation of proton in the region near Bragg peak.
HTML  查看全文  查看/发表评论  下载PDF阅读器
关闭