相对阻止本领(relative stopping power, RSP)

马上去吃汉堡王

已于 2023-07-27 14:50:35 修改

阅读量744

点赞数

分类专栏：数学、物理与算法文章标签：算法

于 2023-07-17 11:42:56 首次发布

本文链接：https://blog.csdn.net/G751888629/article/details/131761436

版权

数学、物理与算法专栏收录该内容

12 篇文章 0 订阅

订阅专栏

接着Bethe-Bloch公式往下写。

Bethe-Bloch给出了入射粒子的平均能量损失率公式 -dE/dx，又称为“组织本领”（stopping power, S），或简称为“能量损失”或“能损”。

The linear stopping power, −dE/dx, is energy loss of a particle per unit distance. The mass stopping power, (1/ρ)dE/dx, is the linear stopping power divided by the density (ρ) of the substance. 下文的阻止本领都默认为质量组织本领。

$S=4\pi N_Ar_e^2m_ec^2z^2\frac{Z}{A}\frac{1}{\beta^2}\cdot (\ln(\frac{2m_ec^2\beta^2}{I(1-\beta^2)})-\beta^2)$

相对组织本领（reletive mass stopping power, RSP）is the ratio of the mass stopping power of a substance to that of a standard substance, usually water in radiology.

物质medium (m)与水water (w)的相对组织本领 RSP 可以表示为

$RSP=\frac{S_m}{S_w}=\frac{N_g^m}{N_g^w}\cdot\frac{\ln(\frac{2m_ec^2\beta^2}{I_m(1-\beta^2)})-\beta^2}{\ln(\frac{2m_ec^2\beta^2}{I_w(1-\beta^2)})-\beta^2}$

$N_g=N_A\cdot\frac{Z}{A}, \rho_e=\rho\cdot N_g$

Ng是单位质量内的电子数量（mass electron density），rho_e是电子密度。

RSP对能量不敏感。

$RSP=\frac{N_g^m}{N_g^w}\cdot\frac{\ln\left (\frac{2m_ec^2\beta^2}{1-\beta^2}\right)-\beta^2-\ln I_m}{\ln\left(\frac{2m_ec^2\beta^2}{1-\beta^2}\right)-\beta^2-\ln I_w}$

$\beta=\sqrt{1-\left ( \frac{1}{E/m_0c^2+1} \right )^2}$

取E=200MeV进行化简，入射粒子为质子

beta = @(E) sqrt(1-(1/(E/938.272+1))^2);
result = @(E) log(2*0.511*1e6*beta(E)^2/(1-beta(E)^2)) - beta(E)^2;

result(200)
% ans = 12.7654

$RSP=\frac{N_g^m}{N_g^w}\cdot\frac{12.77-\ln I_m}{12.77-\ln I_w}=\frac{N_g^m}{N_g^w}\cdot\frac{12.77-\ln I_m}{8.45}$

ps: 取水的平均激发能为75eV。

根据Bragg additivity rule，混合物的平均激发能(mean ionization energy, $I_m$ )为

$\ln I_m=\left(\sum \frac{\omega_iZ_i}{A_i} \ln I_i\right )\left(\sum \frac{\omega_iZ_i}{A_i}\right )^{-1}$

where Zi and Ai are the atomic number and atomic weight of the i th element and ωi is its propotion by weight.

关于双能CT：双能量 CT 成像即利用两种不同能量的 X 射线对物体进行断层扫描，得到不同能量物质衰减的分布。采用双能量 CT 重建算法，可以计算出被扫描物体的等效原子序数和电子密度的分布。

ln(I)与等效原子系数Zeff大体呈现一个线性相关

$\ln(I)=a\cdot Z_{eff}+b$

因此可以用双能CT得到RSP

$RSP=\frac{\rho_{e,m}/\rho_m}{\rho_{e,w}/\rho_w}\cdot\frac{12.77-(a\cdot Z_{eff}+b)}{8.45}$

TPS里电子密度标定曲线得到的是相对电子密度（参考TPS手册和Schneider 1996）

Electron Density: This curve is used by the AAA and MRDC algorithms to convert HU values to relative electron density.

The usual form of the calibration is a bilinear relationship between relative electron density and CT units.

参考：

1. Experimental verification of ion stopping power prediction from dual energy CT data in tissue surrogates

2. The calibration of CT Hounsfield units for radiotherapy treatment planning

3. X-RAY ArIENUATION COEFFICIENTS OF ELEMENTS AND MIXTURES

第1篇论文里有关电子密度的符号写的有错，注意区别电子密度、质量电子密度、相对电子密度。第二篇里Ng的文字定义写错了，但公式没错，Ng is the number of electrons per unit volume（错），应该是 The quantity Ng is the mass electron density,i.e. the number of electrons per unit mass。

2023.7.17