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

Diffusion de Mie due aux structures tissulaires dans le régime térahertz : Vérification expérimentale et Monte Carlo utilisant l'imagerie polarimétrique diffuse dans des fantômes tissulaires hautement atténuants

Informations fondamentales

  • ID de l'article: 2501.00630
  • Titre: Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
  • Auteurs: Erica Heller, Kuangyi Xu, Zachery Harris, M. Hassan Arbab (Université de Stony Brook)
  • Classification: physics.med-ph physics.bio-ph physics.optics
  • Date de soumission: 31 décembre 2024
  • Lien de l'article: https://arxiv.org/abs/2501.00630

Résumé

Cette étude explore le phénomène de diffusion de Mie induit par les structures tissulaires dans le régime térahertz, en mettant l'accent sur ses applications potentielles dans le diagnostic des maladies. La recherche, par le biais de simulations Monte Carlo et de vérification expérimentale, analyse les caractéristiques de diffusion de la lumière térahertz polarisée dans les tissus hautement atténuants. L'étude révèle que les particules diffusantes de différentes tailles présentent des contrastes marqués en termes de degré de polarisation (DOP) et d'intensité de diffusion diffuse, et valide le potentiel diagnostique de cette technique sur des échantillons de brûlures cutanées porcines ex vivo.

Contexte et motivation de la recherche

Définition du problème

Les changements dans les structures tissulaires constituent une caractéristique importante de nombreux processus pathologiques, tels que le bourgeonnement tumoral aux limites des tumeurs cancéreuses et les lésions de brûlure. Les applications biomédicales traditionnelles du térahertz se concentrent principalement sur les propriétés diélectriques des tissus, en particulier la teneur en eau, et exploitent rarement les caractéristiques de diffusion de la lumière térahertz comme outil de diagnostic.

Importance de la recherche

  1. Besoins diagnostiques: Le bourgeonnement tumoral et les amas cellulaires mal différenciés sont des facteurs pronostiques indépendants de la métastase ganglionnaire et de la survie des patients
  2. Limitations technologiques: Les méthodes existantes telles que la coloration H&E et l'immunohistochimie nécessitent des coupes tissulaires et ne permettent pas un diagnostic en temps réel
  3. Mécanismes de diffusion non élucidés: Les mécanismes de réponse polarimétrique dans le régime térahertz n'ont pas été suffisamment étudiés

Motivation de l'innovation

Avec le développement des ellipsomètres térahertz et des systèmes de mesure polarimétrique, il est urgent de développer des modèles de calcul pour expliquer les phénomènes observés et fournir un soutien théorique aux applications de la biophotonique térahertz.

Contributions principales

  1. Établissement d'un modèle Monte Carlo pour la diffusion de Mie térahertz: Basé sur le modèle développé par Xu et Arbab, avec des améliorations pour les particules diffusantes dans les milieux absorbants
  2. Vérification du potentiel diagnostique des mesures polarimètriques: Démonstration expérimentale du contraste entre particules diffusantes de différentes tailles à l'aide de fantômes tissulaires
  3. Fourniture d'une analyse de décomposition de matrice de Mueller: Utilisation de la méthode de décomposition Lu-Chipman pour analyser les propriétés polarimètriques des tissus
  4. Démonstration d'applications pratiques: Vérification de la faisabilité de la technique sur des échantillons de brûlures cutanées porcines ex vivo
  5. Simplification des exigences de mesure: Démonstration que pour les particules sphériques, une seule polarisation est nécessaire pour construire la matrice de Mueller complète

Détails méthodologiques

Conception du système expérimental

Système spectroscopique temporel térahertz:

  • Émetteur: Laser femtoseconde à 1560 nm excitant une antenne photoconductive (PCA)
  • Récepteur: Bras de détection rotatif 0-160 degrés, incluant un polariseur à grille métallique (WGP)
  • Plage de fréquence: 0,2-1,5 THz

Configuration de mesure de diffusion:

  • Angle sélectionné pour la mesure de diffusion diffuse: 140 degrés (contraste maximal entre puissance cohérente et non-cohérente)
  • Balayage en grille: espacement de 2 mm, zone de 2 cm × 2 cm
  • Enregistrement de 100 mesures par état de polarisation

Préparation des fantômes tissulaires

Sélection des matériaux:

  • Matrice: Solution de gélatine bovine à 15% (simulation du fond tissulaire)
  • Particules diffusantes: Granules de polypropylène (simulation de structures à faible indice de réfraction telles que follicules pileux et glandes sudoripares)
  • Plage de taille de particules: 115-280 μm (correspondant aux paramètres de diffusion de Mie 0,2-20)

Procédé de fabrication:

  1. Chauffage de la solution de gélatine à 40°C
  2. Mélange rapide avec les granules de polypropylène
  3. Versement dans des boîtes de culture et solidification (épaisseur ~5 mm)
  4. Analyse de la distribution granulométrique et de la concentration à l'aide d'un microscope optique et du logiciel ImageJ

Simulation Monte Carlo

Algorithme fondamental: Basé sur le code Monte Carlo pour la lumière polarisée du plan méridien d'après Ramella-Roman et al., modifié pour la diffusion de Mie dans les milieux absorbants.

Calcul des paramètres clés:

  • Coefficient de diffusion: μs=2πρnm2k02j=1(2j+1)(aj2+bj2)\mu_s = \frac{2\pi\rho}{|n_m|^2k_0^2}\sum_{j=1}^{\infty}(2j+1)(|a_j|^2+|b_j|^2)
  • Coefficient d'absorption: μa=2πρ(nm)k02[(j=1(2j+1)(cj2ψj(z)ψj(z)dj2ψj(z)ψj(z)))/np]\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]

Construction de la matrice de Mueller: Calcul des paramètres de Stokes spatiaux à l'aide de quatre états de polarisation d'incidence différents (polarisation linéaire horizontale, polarisation linéaire verticale, polarisation linéaire à 45°, polarisation circulaire droite), permettant la construction de la matrice de Mueller.

Configuration expérimentale

Paramètres des fantômes

Quatre fantômes tissulaires différents ont été préparés:

FantômeTaille de particule (μm)Densité surfacique (particules/mm²)Densité volumique (mm⁻³)
A115310
B13027
C1803,58,5
D2806,510,5

Mesure des propriétés des matériaux

Propriétés diélectriques:

  • Polypropylène: indice de réfraction ~1,51, coefficient d'absorption <2 cm⁻¹
  • Gélatine: indice de réfraction et coefficient d'absorption intermédiaires entre la peau et l'eau

Expérience de brûlure ex vivo

  • Échantillon: Peau porcine ex vivo
  • Préparation de la brûlure: Tige de laiton à 180°C appliquée pendant 5 minutes, créant une brûlure d'épaisseur totale d'un pouce
  • Zone de mesure: 2 × 3,5 cm, taille de pixel de 1 mm
  • Acquisition de données: 10 mesures par pixel, moyenne temporelle de 10 fois

Résultats expérimentaux

Résultats de la simulation Monte Carlo

Efficacité de diffusion et facteur d'asymétrie: Avec l'augmentation de la taille des particules, l'efficacité de diffusion et le facteur d'asymétrie augmentent tous deux dans la bande basse fréquence (<1 THz), indiquant que les particules plus grandes produisent une diffusion avant plus forte.

Décomposition de la matrice de Mueller: Décomposition Lu-Chipman des matrices de Mueller pour les fantômes A et D:

  • Fantôme A: MΔ,Adiag[1,0,921,0,925,0,849]M_{\Delta,A} \approx \text{diag}[1, 0,921, 0,925, 0,849]
  • Fantôme D: MΔ,Ddiag[1,0,586,0,586,0,207]M_{\Delta,D} \approx \text{diag}[1, 0,586, 0,586, 0,207]

Les résultats indiquent que les deux fantômes peuvent être approximés comme une combinaison d'un atténuateur uniforme et d'un dépolariseur non uniforme.

Résultats de vérification expérimentale

Intensité de diffusion diffuse:

  • Les simulations et les expériences montrent toutes deux: plus la taille des particules est grande, plus l'intensité de diffusion est élevée
  • Dépendance fréquentielle: intensité élevée à basse fréquence, décroissance exponentielle avec la fréquence
  • Cohérence des tendances: les résultats de simulation et expérimentaux sont en accord

Contraste du degré de polarisation (DOP):

  • Le fantôme D voit son DOP chuter à environ 0,6 près de 0,6 THz
  • Le DOP du fantôme C diminue d'une manière intermédiaire entre A, B et D
  • Localisation fréquentielle: la position fréquentielle de la chute du DOP est corrélée à la taille des particules

Analyse des tissus brûlés

Contraste spatial: À la fréquence de 0,6 THz, la région brûlée et le tissu sain montrent tous deux des contrastes évidents en termes de DOP et d'intensité de diffusion diffuse.

Caractéristiques spectrales:

  • Le DOP du tissu sain diminue plus rapidement dans la plage 0,2-0,8 THz
  • L'intensité de diffusion diffuse de la région brûlée est significativement plus élevée que celle du tissu sain
  • Cela peut être attribué à la destruction des grandes structures cutanées (follicules pileux, glandes sudoripares) causée par la brûlure et à la perte d'humidité

Travaux connexes

Applications biomédicales du térahertz

Les recherches existantes se concentrent principalement sur:

  1. Diagnostic du cancer: Études spectroscopiques du cancer du sein ex vivo, du carcinome épidermoïde oral, etc.
  2. Classification des brûlures cutanées: Classement des brûlures in vivo et prédiction de la cicatrisation
  3. Analyse des propriétés diélectriques: Principalement axée sur les changements de teneur en eau des tissus

Spectroscopie de diffusion polarisée

Applications largement établies dans d'autres bandes spectrales:

  1. Imagerie endoscopique: Détection des polypes du côlon et des lésions péritonéales
  2. Diagnostic du cancer: Multiples sous-types de cancer incluant le cancer de la peau, le cancer de l'estomac, le cancer du sein, etc.
  3. Modélisation Monte Carlo: Applications des modèles de diffusion lumineuse en biophotonique

Lacunes technologiques

La recherche sur les mécanismes de réponse polarimétrique dans le régime térahertz est relativement en retard, avec un manque de modèles de calcul matures pour expliquer les phénomènes observés.

Conclusions et discussion

Conclusions principales

  1. Vérification du modèle réussie: Les résultats de simulation Monte Carlo sont hautement cohérents avec les données expérimentales en termes de tendances
  2. Potentiel diagnostique confirmé: Le DOP et l'intensité de diffusion diffuse peuvent distinguer les particules diffusantes de différentes tailles
  3. Simplification des mesures: Pour les particules sphériques, une mesure à polarisation unique suffit pour obtenir les informations complètes de la matrice de Mueller
  4. Faisabilité clinique: Démonstration réussie de la capacité de détection des états pathologiques dans les tissus ex vivo

Analyse des limitations

  1. Limitation de la bande passante du système: La bande passante actuelle de 0,2-1,5 THz limite la capacité de résolution des différences de petites particules
  2. Contraintes du rapport signal-bruit: La fiabilité des mesures de DOP au-dessus de 0,8 THz diminue
  3. Effet de la rugosité de surface: Hypothèse d'une surface d'échantillon lisse et plane; les effets de dépolarisation de surface doivent être considérés dans les applications pratiques
  4. Exigences de différence de taille de particule: Une différence de taille de particule suffisamment grande est nécessaire pour produire un contraste de signal évident

Directions de recherche futures

  1. Extension de la bande passante du système: Utilisation de cristaux électro-optiques pour remplacer le PCA, réalisant une bande passante de 8 THz
  2. Vérification sur tissus cancéreux ex vivo: Extension de l'étude des fantômes aux vrais tissus cancéreux
  3. Optimisation du traitement du signal: Développement de techniques de décomposition de signaux d'onde et d'analyse statistique de speckle
  4. Séparation des sources de dépolarisation multiples: Distinction des contributions relatives de la rugosité de surface et des structures internes

Évaluation approfondie

Points forts

  1. Combinaison théorie et expérience: Intégration de la simulation Monte Carlo et de la vérification expérimentale, fournissant un cadre de recherche complet
  2. Innovation méthodologique: Première étude systématique des propriétés polarimètriques de la diffusion de Mie dans le régime térahertz
  3. Valeur pratique élevée: Simplification des exigences de mesure de la matrice de Mueller, améliorant la faisabilité des applications cliniques
  4. Vérification suffisante: Stratégie de vérification progressive des fantômes idéaux aux vrais tissus biologiques, raisonnable
  5. Analyse approfondie: La décomposition Lu-Chipman fournit une compréhension plus profonde des mécanismes physiques

Insuffisances

  1. Diversité limitée des échantillons: Utilisation uniquement de fantômes polypropylène-gélatine; les propriétés des matériaux peuvent différer des vrais tissus
  2. Analyse statistique insuffisante: Manque de tests de signification statistique détaillés et d'analyse d'erreur
  3. Vérification clinique limitée: Vérification sur un seul échantillon de brûlure; des études cliniques à plus grande échelle sont nécessaires
  4. Considération insuffisante de la diffusion complexe: Principalement axée sur les particules sphériques; les structures tissulaires réelles sont plus complexes

Évaluation de l'impact

  1. Valeur académique: Fournit une base théorique importante et des méthodes expérimentales pour la biophotonique térahertz
  2. Avancement technologique: Promotion de l'application de la technologie d'imagerie polarimétrique térahertz dans le diagnostic médical
  3. Signification interdisciplinaire: Connexion de la théorie de la diffusion électromagnétique, de la modélisation Monte Carlo et de l'imagerie biomédicale
  4. Perspectives industrielles: Fourniture d'une voie technologique pour le développement d'équipements médicaux térahertz

Scénarios d'application

  1. Diagnostic des maladies de la peau: Classification des brûlures, détection du cancer de la peau
  2. Identification des limites tumorales: Navigation chirurgicale, analyse pathologique
  3. Caractérisation des structures tissulaires: Analyse des fibres de collagène, des amas cellulaires
  4. Surveillance de l'efficacité thérapeutique: Suivi des changements de structure tissulaire au cours du traitement

Références

L'article cite 41 références pertinentes, couvrant plusieurs domaines incluant la technologie térahertz, l'optique biomédicale, la modélisation Monte Carlo, etc., fournissant une base théorique solide pour la recherche.


Évaluation générale: Cet article est une recherche interdisciplinaire de haute qualité possédant une valeur théorique et pratique importante dans le domaine de l'imagerie biomédicale térahertz. La méthodologie de recherche est rigoureuse, les résultats sont fiables, et elle jette une base solide pour le développement ultérieur de ce domaine.