Significance: Changes in the structure of tissue occur in many disease processes, such as the boundaries of cancerous tumors and burn injuries. Spectroscopic and polarimetric alterations of terahertz light caused by Mie scattering patterns has the potential to be a diagnostic marker.
Aim: We present an analysis of Monte Carlo simulation of Mie scattering of polarized terahertz light from cancerous tumor budding, compare the simulation to experimental results obtained in phantom models, and present an analysis of a polarization-sensitive terahertz scan of an ex vivo porcine burn injury.
Approach: Using a Monte Carlo simulation, we modeled the changes in diffuse intensity and degree of polarization of broadband off-specular terahertz light due to scattering particles in highly attenuating tissue. We extracted the Mueller matrix of the tissue using this model and analyzed the Lu-Chipman product decomposition matrices. We compared this model to experimental data from four phantoms consisting of polypropylene particles of varying sizes embedded in gelatin. Finally, we induced a full-thickness burn injury in ex vivo porcine skin samples and compared experimental data from burned and healthy regions of the tissue.
Results: Simulation revealed contrast in the Stokes vectors and Mueller Matrix elements for varying scattering particle sizes. Experimental phantom results showed contrast between different sizes of scattering particles in degree of polarization and diffuse intensity in agreement with Monte Carlo simulation results. Finally, we demonstrated a similar diffused imaging signal contrast between burned and healthy regions of ex vivo porcine skin.
Conclusion: Polarimetric terahertz imaging has the potential to detect structural changes due to biological disease processes.
Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
- 论文ID: 2501.00630
- 标题: Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
- 作者: Erica Heller, Kuangyi Xu, Zachery Harris, M. Hassan Arbab (Stony Brook University)
- 分类: physics.med-ph physics.bio-ph physics.optics
- 提交时间: 2024年12月31日
- 论文链接: https://arxiv.org/abs/2501.00630
本研究探讨了太赫兹频段下组织结构引起的Mie散射现象,重点关注其在疾病诊断中的潜在应用。研究通过Monte Carlo模拟和实验验证,分析了极化太赫兹光在高衰减组织中的散射特性。研究发现,不同尺寸散射粒子在偏振度(DOP)和漫散射强度方面表现出明显对比,并在离体猪皮肤烧伤样本中验证了该技术的诊断潜力。
组织结构的变化是许多疾病过程的重要特征,如癌症肿瘤边界的肿瘤出芽(tumor budding)和烧伤损伤等。传统的太赫兹生物医学应用主要关注组织的介电性质,特别是水含量,而很少利用太赫兹光的散射特性作为诊断手段。
- 疾病诊断需求:肿瘤出芽和分化不良的细胞簇是淋巴结转移和患者生存的独立预后因子
- 技术局限性:现有方法如H&E染色和免疫组织化学需要组织切片,无法实现实时诊断
- 散射机制未明:太赫兹频段的偏振响应差异机制尚未被充分研究
随着太赫兹椭偏仪和偏振测量系统的发展,迫切需要计算模型来解释观察到的现象,为太赫兹生物光子学应用提供理论支撑。
- 建立了太赫兹Mie散射的Monte Carlo模型:基于Xu和Arbab开发的模型,针对吸收介质中的散射粒子进行了改进
- 验证了偏振测量的诊断潜力:通过组织模体实验证明了不同尺寸散射粒子的对比度
- 提供了Mueller矩阵分解分析:使用Lu-Chipman分解方法分析了组织的偏振特性
- 展示了实际应用:在离体猪皮肤烧伤样本中验证了技术的可行性
- 简化了测量要求:证明对于球形粒子,仅需一种偏振态测量即可构建完整的Mueller矩阵
太赫兹时域光谱系统:
- 发射端:使用1560 nm飞秒激光激发光导天线(PCA)
- 接收端:可旋转0-160度的检测臂,包含线栅偏振器(WGP)
- 频率范围:0.2-1.5 THz
散射测量配置:
- 选择140度角进行漫散射测量(相干功率和非相干功率对比最明显)
- 栅格扫描:2mm像素间距,2cm×2cm区域
- 每个偏振态记录100次测量
材料选择:
- 基质:15%牛皮明胶溶液(模拟组织背景)
- 散射粒子:聚丙烯颗粒(模拟毛囊、汗腺等低折射率结构)
- 粒径范围:115-280 μm(对应Mie散射参数0.2-20)
制备工艺:
- 明胶溶液加热至40°C
- 快速混入聚丙烯颗粒
- 倒入培养皿固化(厚度约5mm)
- 使用光学显微镜和ImageJ软件分析粒径分布和浓度
核心算法:
基于Ramella-Roman等人的子午面偏振光Monte Carlo代码,针对吸收介质中的Mie散射进行了修改。
关键参数计算:
- 散射系数:μs=∣nm∣2k022πρ∑j=1∞(2j+1)(∣aj∣2+∣bj∣2)
- 吸收系数:μa=ℜ(nm)k022πρℑ[(∑j=1∞(2j+1)(∣cj∣2ψj(z)ψj′∗(z)−∣dj∣2ψj′(z)ψj∗(z)))/np]
Mueller矩阵构建:
使用四种不同入射偏振态(水平线偏振、垂直线偏振、45°线偏振、右圆偏振)计算空间Stokes参数,进而构建Mueller矩阵。
制备了四个不同的组织模体:
| 模体 | 粒径(μm) | 表面密度(粒子/mm²) | 体积密度(mm⁻³) |
|---|
| A | 115 | 3 | 10 |
| B | 130 | 2 | 7 |
| C | 180 | 3.5 | 8.5 |
| D | 280 | 6.5 | 10.5 |
介电性质:
- 聚丙烯:折射率~1.51,吸收系数<2 cm⁻¹
- 明胶:折射率和吸收系数介于皮肤和水之间
- 样本:离体猪皮肤
- 烧伤制备:180°C黄铜棒作用5分钟,形成1英寸全层烧伤
- 测量区域:2×3.5 cm,1mm像素尺寸
- 数据采集:每像素10次测量,10次时间平均
散射效率和不对称因子:
随着粒径增大,散射效率和不对称因子在低频段(<1 THz)均增加,表明更大粒子产生更强的前向散射。
Mueller矩阵分解:
对于模体A和D的Mueller矩阵进行Lu-Chipman分解:
- 模体A:MΔ,A≈diag[1,0.921,0.925,0.849]
- 模体D:MΔ,D≈diag[1,0.586,0.586,0.207]
结果表明两个模体都可以近似为均匀衰减器和非均匀去偏振器的组合。
漫散射强度:
- 仿真和实验均显示:粒径越大,散射强度越高
- 频率依赖性:低频强度高,随频率指数衰减
- 趋势一致性:仿真与实验结果趋势相符
偏振度(DOP)对比:
- 模体D在0.6 THz附近DOP降至约0.6
- 模体C的DOP下降程度介于A、B和D之间
- 频率定位:DOP下降的频率位置与粒径相关
空间对比:
在0.6 THz频率下,烧伤区域与健康组织在DOP和漫散射强度方面均显示出明显对比。
频谱特性:
- 健康组织的DOP在0.2-0.8 THz范围内下降更快
- 烧伤区域的漫散射强度显著高于健康组织
- 这可能归因于烧伤导致的大型皮肤结构(毛囊、汗腺)破坏和水分流失
现有研究主要集中在:
- 癌症诊断:离体乳腺癌、口腔鳞状细胞癌等的光谱学研究
- 皮肤烧伤分类:体内烧伤分级和愈合预测
- 介电性质分析:主要关注组织水含量变化
在其他波段已有广泛应用:
- 内窥镜成像:结肠息肉和腹膜病变检测
- 癌症诊断:皮肤癌、胃癌、乳腺癌等多种癌症亚型
- Monte Carlo建模:光散射模型在生物光子学中的应用
太赫兹频段的偏振响应机制研究相对滞后,缺乏成熟的计算模型来解释观察到的现象。
- 模型验证成功:Monte Carlo仿真结果与实验数据在趋势上高度一致
- 诊断潜力确认:DOP和漫散射强度能够区分不同尺寸的散射粒子
- 测量简化:对于球形粒子,单一偏振态测量即可获得完整Mueller矩阵信息
- 临床可行性:在离体组织中成功演示了病理状态检测能力
- 系统带宽限制:当前0.2-1.5 THz带宽限制了小粒径差异的分辨能力
- 信噪比约束:0.8 THz以上的DOP测量可靠性下降
- 表面粗糙度影响:假设样本表面光滑平整,实际应用中需考虑表面去偏振效应
- 粒径差异要求:需要较大的粒径差异才能产生明显的信号对比
- 扩展系统带宽:使用电光晶体替代PCA,实现8 THz带宽
- 离体癌组织验证:将模体研究扩展到真实癌组织
- 信号处理优化:开发波信号分解和统计散斑分析技术
- 多源去偏振分离:区分表面粗糙度和内部结构的相对贡献
- 理论与实验结合:Monte Carlo仿真与实验验证相结合,提供了完整的研究框架
- 方法学创新:首次在太赫兹频段系统研究了Mie散射的偏振特性
- 实用价值高:简化了Mueller矩阵测量要求,提高了临床应用的可行性
- 验证充分:从理想模体到真实生物组织的渐进式验证策略合理
- 分析深入:Lu-Chipman分解提供了物理机制的深层理解
- 样本多样性有限:仅使用了聚丙烯-明胶模体,材料特性可能与真实组织存在差异
- 统计分析不足:缺乏详细的统计显著性检验和误差分析
- 临床验证有限:仅在一个烧伤样本上进行了验证,需要更大规模的临床研究
- 复杂散射考虑不足:主要考虑球形粒子,实际组织结构更加复杂
- 学术价值:为太赫兹生物光子学提供了重要的理论基础和实验方法
- 技术推进:推动了太赫兹偏振成像技术在医学诊断中的应用
- 跨学科意义:连接了电磁散射理论、Monte Carlo建模和生物医学成像
- 产业前景:为太赫兹医疗设备的开发提供了技术路径
- 皮肤病诊断:烧伤分级、皮肤癌检测
- 肿瘤边界识别:手术导航、病理分析
- 组织结构表征:胶原纤维、细胞簇分析
- 药物疗效监测:治疗过程中组织结构变化追踪
论文引用了41篇相关文献,涵盖了太赫兹技术、生物医学光学、Monte Carlo建模等多个领域的重要工作,为研究提供了坚实的理论基础。
总体评价:这是一篇高质量的跨学科研究论文,在太赫兹生物医学成像领域具有重要的理论和实用价值。研究方法严谨,结果可信,为该领域的进一步发展奠定了良好基础。