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.
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
- 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
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.
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.
- 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
- Limitazioni Tecnologiche: I metodi esistenti come la colorazione H&E e l'immunoistochimica richiedono sezioni tissutali e non consentono una diagnostica in tempo reale
- Meccanismi di Scattering Non Chiariti: I meccanismi della risposta polarimetrica nel regime dei terahertz non sono stati sufficientemente studiati
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.
- 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
- Verifica del Potenziale Diagnostico delle Misurazioni Polarimetriche: Dimostrazione sperimentale del contrasto tra particelle di scattering di diverse dimensioni utilizzando fantomi tissutali
- Analisi di Decomposizione della Matrice di Mueller: Utilizzo del metodo di decomposizione Lu-Chipman per analizzare le proprietà polarimetriche dei tessuti
- Dimostrazione di Applicazioni Pratiche: Verifica della fattibilità della tecnica in campioni di ustioni cutanee suine ex vivo
- 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
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
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:
- Riscaldamento della soluzione di gelatina a 40°C
- Miscelazione rapida con granuli di polipropilene
- Versamento in piastre di coltura per la solidificazione (spessore circa 5 mm)
- Analisi della distribuzione delle dimensioni e della concentrazione utilizzando microscopio ottico e software ImageJ
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=∣nm∣2k022πρ∑j=1∞(2j+1)(∣aj∣2+∣bj∣2)
- Coefficiente di assorbimento: μa=ℜ(nm)k022πρℑ[(∑j=1∞(2j+1)(∣cj∣2ψj(z)ψj′∗(z)−∣dj∣2ψj′(z)ψj∗(z)))/np]
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.
Sono stati preparati quattro fantomi tissutali diversi:
| Fantoma | Dimensione Particella (μm) | Densità Superficiale (particelle/mm²) | Densità Volumetrica (mm⁻³) |
|---|
| A | 115 | 3 | 10 |
| B | 130 | 2 | 7 |
| C | 180 | 3,5 | 8,5 |
| D | 280 | 6,5 | 10,5 |
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
- 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
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Δ,A≈diag[1,0,921,0,925,0,849]
- Fantoma D: MΔ,D≈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.
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
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à
La ricerca esistente si concentra principalmente su:
- Diagnostica del Cancro: Studi spettroscopici su cancro al seno ex vivo, carcinoma squamoso orale e altri
- Classificazione delle Ustioni Cutanee: Classificazione delle ustioni in vivo e previsione della guarigione
- Analisi delle Proprietà Dielettriche: Principalmente focalizzata su variazioni del contenuto di acqua nei tessuti
Ampiamente applicata in altri intervalli di lunghezza d'onda:
- Imaging Endoscopico: Rilevamento di polipi del colon e lesioni peritoneali
- Diagnostica del Cancro: Sottotipi di cancro della pelle, dello stomaco, del seno e altri
- Modellazione Monte Carlo: Applicazione di modelli di scattering luminoso nella biofotonica
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.
- Verifica del Modello Riuscita: I risultati della simulazione Monte Carlo sono altamente coerenti con i dati sperimentali nelle tendenze
- Conferma del Potenziale Diagnostico: Il DOP e l'intensità dello scattering diffuso possono distinguere particelle di scattering di diverse dimensioni
- Semplificazione della Misurazione: Per particelle sferiche, una singola misura di stato di polarizzazione è sufficiente per ottenere informazioni complete della matrice di Mueller
- Fattibilità Clinica: Dimostrazione riuscita della capacità di rilevare stati patologici in tessuti ex vivo
- 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
- Vincoli del Rapporto Segnale-Rumore: L'affidabilità della misurazione del DOP al di sopra di 0,8 THz diminuisce
- 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
- Requisiti di Differenza di Dimensione delle Particelle: È necessaria una differenza di dimensione delle particelle sufficientemente grande per produrre un contrasto di segnale evidente
- 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
- Verifica su Tessuto Canceroso Ex Vivo: Estensione della ricerca su fantomi a veri tessuti cancerosi
- Ottimizzazione dell'Elaborazione dei Segnali: Sviluppo di tecniche di decomposizione dei segnali d'onda e analisi statistica delle speckle
- Separazione della Depolarizzazione da Fonti Multiple: Distinzione dei contributi relativi della rugosità superficiale e delle strutture interne
- Integrazione di Teoria e Esperimento: Combinazione di simulazione Monte Carlo e verifica sperimentale fornisce un quadro di ricerca completo
- Innovazione Metodologica: Primo studio sistematico delle proprietà polarimetriche dello scattering di Mie nel regime dei terahertz
- Alto Valore Pratico: Semplificazione dei requisiti di misurazione della matrice di Mueller, miglioramento della fattibilità dell'applicazione clinica
- Verifica Sufficiente: Strategia di verifica progressiva da fantomi ideali a tessuti biologici reali è razionale
- Analisi Approfondita: La decomposizione Lu-Chipman fornisce una comprensione più profonda dei meccanismi fisici
- Varietà Limitata di Campioni: Utilizzo solo di fantomi polipropilene-gelatina; le proprietà dei materiali potrebbero differire dai veri tessuti
- Analisi Statistica Insufficiente: Mancanza di test di significatività statistica dettagliati e analisi degli errori
- Verifica Clinica Limitata: Verifica solo su un campione di ustione; è necessaria una ricerca clinica su scala più ampia
- Considerazione Insufficiente dello Scattering Complesso: Principalmente focalizzato su particelle sferiche; le strutture tissutali reali sono più complesse
- Valore Accademico: Fornisce una base teorica importante e metodi sperimentali per la biofotonica dei terahertz
- Avanzamento Tecnologico: Promuove l'applicazione della tecnologia di imaging polarimetrico terahertz nella diagnostica medica
- Significato Interdisciplinare: Collega la teoria dello scattering elettromagnetico, la modellazione Monte Carlo e l'imaging biomedico
- Prospettive Industriali: Fornisce un percorso tecnologico per lo sviluppo di dispositivi medici terahertz
- Diagnostica delle Malattie della Pelle: Classificazione delle ustioni, rilevamento del cancro della pelle
- Identificazione dei Margini Tumorali: Navigazione chirurgica, analisi patologica
- Caratterizzazione delle Strutture Tissutali: Analisi delle fibre di collagene, cluster cellulari
- Monitoraggio dell'Efficacia dei Farmaci: Tracciamento dei cambiamenti nelle strutture tissutali durante il trattamento
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.