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

Scattering di Mie dovuto a strutture tissutali nel regime dei terahertz: Verifica sperimentale e Monte Carlo utilizzando imaging polarimetrico diffuso in fantomi tissutali altamente attenuanti

Informazioni Fondamentali

  • ID Articolo: 2501.00630
  • Titolo: Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
  • Autori: Erica Heller, Kuangyi Xu, Zachery Harris, M. Hassan Arbab (Stony Brook University)
  • Classificazione: physics.med-ph physics.bio-ph physics.optics
  • Data di Presentazione: 31 dicembre 2024
  • Link Articolo: https://arxiv.org/abs/2501.00630

Riassunto

Questo studio esamina il fenomeno dello scattering di Mie indotto da strutture tissutali nel regime dei terahertz, con particolare attenzione alle potenziali applicazioni nella diagnostica delle malattie. La ricerca analizza le caratteristiche di scattering della luce terahertz polarizzata in tessuti altamente attenuanti attraverso simulazioni Monte Carlo e verifica sperimentale. Lo studio rivela contrasti significativi nelle particelle di scattering di diverse dimensioni in termini di grado di polarizzazione (DOP) e intensità di scattering diffuso, e verifica il potenziale diagnostico della tecnica in campioni di ustioni cutanee suine ex vivo.

Contesto e Motivazione della Ricerca

Definizione del Problema

Le variazioni nelle strutture tissutali rappresentano caratteristiche importanti di molti processi patologici, come il tumor budding ai margini dei tumori cancerosi e i danni da ustione. Le applicazioni biomediche tradizionali dei terahertz si concentrano principalmente sulle proprietà dielettriche dei tessuti, in particolare sul contenuto di acqua, e raramente sfruttano le caratteristiche di scattering della luce terahertz come mezzo diagnostico.

Importanza della Ricerca

  1. Esigenze di Diagnostica delle Malattie: Il tumor budding e i cluster cellulari scarsamente differenziati sono fattori prognostici indipendenti per la metastasi linfonodale e la sopravvivenza del paziente
  2. Limitazioni Tecnologiche: I metodi esistenti come la colorazione H&E e l'immunoistochimica richiedono sezioni tissutali e non consentono una diagnostica in tempo reale
  3. Meccanismi di Scattering Non Chiariti: I meccanismi della risposta polarimetrica nel regime dei terahertz non sono stati sufficientemente studiati

Motivazione dell'Innovazione

Con lo sviluppo degli ellipsometri terahertz e dei sistemi di misurazione polarimetrica, è urgente sviluppare modelli computazionali per spiegare i fenomeni osservati e fornire supporto teorico alle applicazioni della biofotonica terahertz.

Contributi Principali

  1. Sviluppo di un Modello Monte Carlo per lo Scattering di Mie nei Terahertz: Miglioramento del modello sviluppato da Xu e Arbab per particelle di scattering in mezzi assorbenti
  2. Verifica del Potenziale Diagnostico delle Misurazioni Polarimetriche: Dimostrazione sperimentale del contrasto tra particelle di scattering di diverse dimensioni utilizzando fantomi tissutali
  3. Analisi di Decomposizione della Matrice di Mueller: Utilizzo del metodo di decomposizione Lu-Chipman per analizzare le proprietà polarimetriche dei tessuti
  4. Dimostrazione di Applicazioni Pratiche: Verifica della fattibilità della tecnica in campioni di ustioni cutanee suine ex vivo
  5. Semplificazione dei Requisiti di Misurazione: Dimostrazione che per particelle sferiche è sufficiente una singola misura di stato di polarizzazione per costruire la matrice di Mueller completa

Dettagli Metodologici

Progettazione del Sistema Sperimentale

Sistema Spettroscopico nel Dominio del Tempo dei Terahertz:

  • Estremità di trasmissione: Laser a femtosecondi da 1560 nm che eccita un'antenna fotoconducente (PCA)
  • Estremità di ricezione: Braccio di rilevamento ruotabile da 0-160 gradi, con polarizzatore a griglia di filo (WGP)
  • Intervallo di frequenza: 0,2-1,5 THz

Configurazione di Misurazione dello Scattering:

  • Selezione di un angolo di 140 gradi per la misurazione dello scattering diffuso (contrasto massimo tra potenza coerente e incoerente)
  • Scansione a griglia: spaziatura di pixel di 2 mm, area di 2 cm × 2 cm
  • 100 misurazioni registrate per ogni stato di polarizzazione

Preparazione dei Fantomi Tissutali

Scelta dei Materiali:

  • Matrice: Soluzione di gelatina bovina al 15% (simula lo sfondo tissutale)
  • Particelle di scattering: Granuli di polipropilene (simulano strutture a basso indice di rifrazione come follicoli piliferi e ghiandole sudoripare)
  • Intervallo di dimensioni: 115-280 μm (corrispondente a parametri di scattering di Mie di 0,2-20)

Processo di Preparazione:

  1. Riscaldamento della soluzione di gelatina a 40°C
  2. Miscelazione rapida con granuli di polipropilene
  3. Versamento in piastre di coltura per la solidificazione (spessore circa 5 mm)
  4. Analisi della distribuzione delle dimensioni e della concentrazione utilizzando microscopio ottico e software ImageJ

Simulazione Monte Carlo

Algoritmo Principale: Basato sul codice Monte Carlo per la luce polarizzata nel piano meridiano di Ramella-Roman et al., modificato per lo scattering di Mie in mezzi assorbenti.

Calcolo dei Parametri Chiave:

  • Coefficiente di scattering: μ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)
  • Coefficiente di assorbimento: μ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]

Costruzione della Matrice di Mueller: Calcolo dei parametri di Stokes spaziali utilizzando quattro diversi stati di polarizzazione incidente (polarizzazione lineare orizzontale, polarizzazione lineare verticale, polarizzazione lineare a 45°, polarizzazione circolare destra), da cui si costruisce la matrice di Mueller.

Configurazione Sperimentale

Parametri dei Fantomi

Sono stati preparati quattro fantomi tissutali diversi:

FantomaDimensione Particella (μm)Densità Superficiale (particelle/mm²)Densità Volumetrica (mm⁻³)
A115310
B13027
C1803,58,5
D2806,510,5

Misurazione delle Proprietà dei Materiali

Proprietà Dielettriche:

  • Polipropilene: indice di rifrazione ~1,51, coefficiente di assorbimento <2 cm⁻¹
  • Gelatina: indice di rifrazione e coefficiente di assorbimento intermedi tra la pelle e l'acqua

Esperimento di Ustione Ex Vivo

  • Campione: Pelle suina ex vivo
  • Preparazione dell'ustione: Asta di ottone a 180°C applicata per 5 minuti, creando un'ustione a tutto spessore di 1 pollice
  • Area di misurazione: 2 × 3,5 cm, dimensione di pixel di 1 mm
  • Acquisizione dati: 10 misurazioni per pixel, media temporale di 10 volte

Risultati Sperimentali

Risultati della Simulazione Monte Carlo

Efficienza di Scattering e Fattore di Asimmetria: Con l'aumento delle dimensioni delle particelle, sia l'efficienza di scattering che il fattore di asimmetria aumentano nel segmento a bassa frequenza (<1 THz), indicando che particelle più grandi producono uno scattering più forte in avanti.

Decomposizione della Matrice di Mueller: Decomposizione Lu-Chipman della matrice di Mueller per i fantomi A e D:

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

I risultati indicano che entrambi i fantomi possono essere approssimati come una combinazione di un attenuatore uniforme e un depolarizzatore non uniforme.

Risultati della Verifica Sperimentale

Intensità dello Scattering Diffuso:

  • Sia la simulazione che l'esperimento mostrano: maggiore è la dimensione delle particelle, maggiore è l'intensità di scattering
  • Dipendenza dalla frequenza: intensità elevata a bassa frequenza, decadimento esponenziale con la frequenza
  • Coerenza delle tendenze: le tendenze della simulazione e dei risultati sperimentali sono coerenti

Contrasto del Grado di Polarizzazione (DOP):

  • Il fantoma D mostra un DOP che scende a circa 0,6 vicino a 0,6 THz
  • Il DOP del fantoma C mostra una diminuzione intermedia tra A, B e D
  • Posizionamento in frequenza: la posizione in frequenza della diminuzione del DOP è correlata alle dimensioni delle particelle

Analisi dei Tessuti Ustionati

Contrasto Spaziale: Alla frequenza di 0,6 THz, l'area ustionata e il tessuto sano mostrano contrasti evidenti sia nel DOP che nell'intensità dello scattering diffuso.

Caratteristiche Spettrali:

  • Il DOP del tessuto sano diminuisce più rapidamente nell'intervallo 0,2-0,8 THz
  • L'intensità dello scattering diffuso nell'area ustionata è significativamente superiore al tessuto sano
  • Ciò potrebbe essere attribuito alla distruzione delle grandi strutture cutanee (follicoli piliferi, ghiandole sudoripare) causata dall'ustione e alla perdita di umidità

Lavori Correlati

Applicazioni Biomediche dei Terahertz

La ricerca esistente si concentra principalmente su:

  1. Diagnostica del Cancro: Studi spettroscopici su cancro al seno ex vivo, carcinoma squamoso orale e altri
  2. Classificazione delle Ustioni Cutanee: Classificazione delle ustioni in vivo e previsione della guarigione
  3. Analisi delle Proprietà Dielettriche: Principalmente focalizzata su variazioni del contenuto di acqua nei tessuti

Spettroscopia di Scattering Polarimetrico

Ampiamente applicata in altri intervalli di lunghezza d'onda:

  1. Imaging Endoscopico: Rilevamento di polipi del colon e lesioni peritoneali
  2. Diagnostica del Cancro: Sottotipi di cancro della pelle, dello stomaco, del seno e altri
  3. Modellazione Monte Carlo: Applicazione di modelli di scattering luminoso nella biofotonica

Lacune Tecnologiche

La ricerca dei meccanismi di risposta polarimetrica nel regime dei terahertz è relativamente arretrata, con una mancanza di modelli computazionali maturi per spiegare i fenomeni osservati.

Conclusioni e Discussione

Conclusioni Principali

  1. Verifica del Modello Riuscita: I risultati della simulazione Monte Carlo sono altamente coerenti con i dati sperimentali nelle tendenze
  2. Conferma del Potenziale Diagnostico: Il DOP e l'intensità dello scattering diffuso possono distinguere particelle di scattering di diverse dimensioni
  3. Semplificazione della Misurazione: Per particelle sferiche, una singola misura di stato di polarizzazione è sufficiente per ottenere informazioni complete della matrice di Mueller
  4. Fattibilità Clinica: Dimostrazione riuscita della capacità di rilevare stati patologici in tessuti ex vivo

Analisi delle Limitazioni

  1. Limitazioni della Larghezza di Banda del Sistema: L'attuale larghezza di banda di 0,2-1,5 THz limita la capacità di risoluzione delle differenze di piccole particelle
  2. Vincoli del Rapporto Segnale-Rumore: L'affidabilità della misurazione del DOP al di sopra di 0,8 THz diminuisce
  3. Effetto della Rugosità Superficiale: Si assume che la superficie del campione sia liscia e piana; gli effetti di depolarizzazione della superficie devono essere considerati nelle applicazioni pratiche
  4. Requisiti di Differenza di Dimensione delle Particelle: È necessaria una differenza di dimensione delle particelle sufficientemente grande per produrre un contrasto di segnale evidente

Direzioni di Ricerca Futura

  1. Espansione della Larghezza di Banda del Sistema: Utilizzo di cristalli elettro-ottici al posto di PCA per realizzare una larghezza di banda di 8 THz
  2. Verifica su Tessuto Canceroso Ex Vivo: Estensione della ricerca su fantomi a veri tessuti cancerosi
  3. Ottimizzazione dell'Elaborazione dei Segnali: Sviluppo di tecniche di decomposizione dei segnali d'onda e analisi statistica delle speckle
  4. Separazione della Depolarizzazione da Fonti Multiple: Distinzione dei contributi relativi della rugosità superficiale e delle strutture interne

Valutazione Approfondita

Punti di Forza

  1. Integrazione di Teoria e Esperimento: Combinazione di simulazione Monte Carlo e verifica sperimentale fornisce un quadro di ricerca completo
  2. Innovazione Metodologica: Primo studio sistematico delle proprietà polarimetriche dello scattering di Mie nel regime dei terahertz
  3. Alto Valore Pratico: Semplificazione dei requisiti di misurazione della matrice di Mueller, miglioramento della fattibilità dell'applicazione clinica
  4. Verifica Sufficiente: Strategia di verifica progressiva da fantomi ideali a tessuti biologici reali è razionale
  5. Analisi Approfondita: La decomposizione Lu-Chipman fornisce una comprensione più profonda dei meccanismi fisici

Insufficienze

  1. Varietà Limitata di Campioni: Utilizzo solo di fantomi polipropilene-gelatina; le proprietà dei materiali potrebbero differire dai veri tessuti
  2. Analisi Statistica Insufficiente: Mancanza di test di significatività statistica dettagliati e analisi degli errori
  3. Verifica Clinica Limitata: Verifica solo su un campione di ustione; è necessaria una ricerca clinica su scala più ampia
  4. Considerazione Insufficiente dello Scattering Complesso: Principalmente focalizzato su particelle sferiche; le strutture tissutali reali sono più complesse

Valutazione dell'Impatto

  1. Valore Accademico: Fornisce una base teorica importante e metodi sperimentali per la biofotonica dei terahertz
  2. Avanzamento Tecnologico: Promuove l'applicazione della tecnologia di imaging polarimetrico terahertz nella diagnostica medica
  3. Significato Interdisciplinare: Collega la teoria dello scattering elettromagnetico, la modellazione Monte Carlo e l'imaging biomedico
  4. Prospettive Industriali: Fornisce un percorso tecnologico per lo sviluppo di dispositivi medici terahertz

Scenari Applicabili

  1. Diagnostica delle Malattie della Pelle: Classificazione delle ustioni, rilevamento del cancro della pelle
  2. Identificazione dei Margini Tumorali: Navigazione chirurgica, analisi patologica
  3. Caratterizzazione delle Strutture Tissutali: Analisi delle fibre di collagene, cluster cellulari
  4. Monitoraggio dell'Efficacia dei Farmaci: Tracciamento dei cambiamenti nelle strutture tissutali durante il trattamento

Bibliografia

L'articolo cita 41 riferimenti correlati, coprendo importanti lavori in più campi tra cui tecnologia terahertz, ottica biomedica e modellazione Monte Carlo, fornendo una base teorica solida per la ricerca.


Valutazione Complessiva: Questo è un articolo di ricerca di alta qualità interdisciplinare con importante valore teorico e pratico nel campo dell'imaging biomedico terahertz. Il metodo di ricerca è rigoroso, i risultati sono affidabili e fornisce una base solida per il ulteriore sviluppo del campo.