2025-11-13T23:19:15.294850

Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms

Heller, Xu, Harris et al.

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.

academic

テラヘルツ領域における組織構造に起因するMie散乱：高減衰組織ファントムを用いた拡散偏光測定による実験的およびMonte Carlo検証

基本情報

論文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（ストーニーブルック大学）
分類: 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度角を選択（コヒーレント電力と非コヒーレント電力のコントラストが最も明確）
グリッドスキャン：2 mmピクセル間隔、2 cm × 2 cm領域
各偏光状態について100回の測定を記録

組織ファントムの製造

材料の選択：

マトリックス：15%ウシ皮膚ゼラチン溶液（組織背景をシミュレート）
散乱粒子：ポリプロピレン粒子（毛嚢、汗腺などの低屈折率構造をシミュレート）
粒径範囲：115～280 μm（Mie散乱パラメータ0.2～20に対応）

製造工程：

ゼラチン溶液を40°Cに加熱
ポリプロピレン粒子を急速に混入
培養皿に注ぎ込み固化（厚さ約5 mm）
光学顕微鏡およびImageJソフトウェアを使用して粒径分布と濃度を分析

Monte Carlo シミュレーション

コアアルゴリズム： Ramella-Roman等が開発した子午面偏光光Monte Carloコードに基づき、吸収媒質中のMie散乱に対して修正を加えた。

主要パラメータの計算：

散乱係数： $\mu_s = \frac{2\pi\rho}{|n_m|^2k_0^2}\sum_{j=1}^{\infty}(2j+1)(|a_j|^2+|b_j|^2)$
吸収係数： $\mu_a = \frac{2\pi\rho}{\Re(n_m)k_0^2}\Im[(\sum_{j=1}^{\infty}(2j+1)(|c_j|^2\psi_j(z)\psi_j'^*(z)-|d_j|^2\psi_j'(z)\psi_j^*(z)))/n_p]$