Dispersión de Mie debida a estructuras tisulares en el régimen de terahercios: Verificación experimental y de Monte Carlo utilizando imágenes polarimétrica difusa en fantasmas de tejido altamente atenuante

Información Básica

• ID del Artículo: 2501.00630
• Título: Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
• Autores: Erica Heller, Kuangyi Xu, Zachery Harris, M. Hassan Arbab (Universidad de Stony Brook)
• Clasificación: physics.med-ph physics.bio-ph physics.optics
• Fecha de Presentación: 31 de diciembre de 2024
• Enlace del Artículo: https://arxiv.org/abs/2501.00630

Resumen

Este estudio investiga el fenómeno de dispersión de Mie causado por estructuras tisulares en el régimen de terahercios, enfocándose en sus aplicaciones potenciales en el diagnóstico de enfermedades. La investigación analiza las características de dispersión de luz terahercia polarizada en tejidos altamente atenuantes mediante simulación de Monte Carlo y verificación experimental. El estudio revela que las partículas dispersoras de diferentes tamaños muestran contraste notable en el grado de polarización (DOP) e intensidad de dispersión difusa, y verifica el potencial diagnóstico de la técnica en muestras de quemaduras de piel de cerdo ex vivo.

Antecedentes de Investigación y Motivación

Definición del Problema

Los cambios en la estructura tisular son características importantes de muchos procesos patológicos, como el brote tumoral en los márgenes de tumores cancerosos y lesiones por quemaduras. Las aplicaciones biomédicas tradicionales de terahercios se han enfocado principalmente en las propiedades dieléctricas de los tejidos, especialmente el contenido de agua, utilizando raramente las características de dispersión de luz terahercia como herramienta diagnóstica.

Importancia de la Investigación

1. Necesidades de Diagnóstico de Enfermedades: El brote tumoral y los cúmulos de células mal diferenciadas son factores pronósticos independientes de metástasis en ganglios linfáticos y supervivencia del paciente
2. Limitaciones Tecnológicas: Los métodos existentes como la tinción H&E y la inmunohistoquímica requieren secciones de tejido, imposibilitando el diagnóstico en tiempo real
3. Mecanismo de Dispersión Poco Claro: El mecanismo de respuesta de polarización en el régimen de terahercios aún no ha sido suficientemente estudiado

Motivación Innovadora

Con el desarrollo de elipsómetros de terahercios y sistemas de medición de polarización, existe una necesidad urgente de modelos computacionales que expliquen los fenómenos observados, proporcionando apoyo teórico para aplicaciones en biofotónica de terahercios.

Contribuciones Principales

1. Establecimiento de un Modelo de Monte Carlo para Dispersión de Mie en Terahercios: Basado en el modelo desarrollado por Xu y Arbab, se realizaron mejoras para partículas dispersoras en medios absorbentes
2. Verificación del Potencial Diagnóstico de Mediciones de Polarización: Se demostró mediante experimentos con fantasmas de tejido el contraste de partículas dispersoras de diferentes tamaños
3. Análisis de Descomposición de Matriz de Mueller: Se utilizó el método de descomposición Lu-Chipman para analizar las propiedades de polarización del tejido
4. Demostración de Aplicación Práctica: Se verificó la viabilidad de la técnica en muestras de quemaduras de piel de cerdo ex vivo
5. Simplificación de Requisitos de Medición: Se demostró que para partículas esféricas, solo se necesita una medición de estado de polarización para construir la matriz de Mueller completa

Detalles de la Metodología

Diseño del Sistema Experimental

Sistema de Espectroscopia en el Dominio del Tiempo de Terahercios:

• Extremo de Transmisión: Utiliza láser de femtosegundos de 1560 nm para excitar antena conductora fotoinducida (PCA)
• Extremo de Recepción: Brazo de detección rotativo de 0-160 grados, incluye polarizador de rejilla de alambre (WGP)
• Rango de Frecuencia: 0.2-1.5 THz

Configuración de Medición de Dispersión:

• Ángulo seleccionado de 140 grados para medición de dispersión difusa (contraste más evidente entre potencia coherente e incoherente)
• Escaneo de cuadrícula: espaciado de píxeles de 2 mm, área de 2 cm × 2 cm
• Cien mediciones registradas para cada estado de polarización

Preparación de Fantasmas de Tejido

Selección de Materiales:

• Matriz: Solución de gelatina bovina al 15% (simula el fondo del tejido)
• Partículas Dispersoras: Partículas de polipropileno (simulan estructuras de bajo índice de refracción como folículos pilosos y glándulas sudoríparas)
• Rango de Tamaño de Partícula: 115-280 μm (correspondiente a parámetro de dispersión de Mie de 0.2-20)

Proceso de Preparación:

1. Calentar solución de gelatina a 40°C
2. Mezclar rápidamente con partículas de polipropileno
3. Verter en placa de cultivo para solidificar (espesor aproximado de 5 mm)
4. Analizar distribución de tamaño de partícula y concentración utilizando microscopio óptico y software ImageJ

Simulación de Monte Carlo

Algoritmo Principal: Basado en el código de Monte Carlo de luz polarizada en el plano meridional de Ramella-Roman et al., modificado para dispersión de Mie en medios absorbentes.

Cálculo de Parámetros Clave:

• Coeficiente de Dispersión: $\mu_s = \frac{2\pi\rho}{|n_m|^2k_0^2}\sum_{j=1}^{\infty}(2j+1)(|a_j|^2+|b_j|^2)$
• Coeficiente de Absorción: $\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]$

Construcción de Matriz de Mueller: Se calculan parámetros de Stokes espaciales utilizando cuatro estados de polarización incidente diferentes (polarización lineal horizontal, polarización lineal vertical, polarización lineal a 45°, polarización circular derecha), construyendo posteriormente la matriz de Mueller.

Configuración Experimental

Parámetros de Fantasmas

Se prepararon cuatro fantasmas de tejido diferentes:

Medición de Propiedades de Materiales

Propiedades Dieléctricas:

• Polipropileno: Índice de refracción ~1.51, coeficiente de absorción <2 cm⁻¹
• Gelatina: Índice de refracción y coeficiente de absorción intermedios entre piel y agua

Experimento de Quemadura Ex Vivo

• Muestra: Piel de cerdo ex vivo
• Preparación de Quemadura: Varilla de latón a 180°C actuando durante 5 minutos, formando quemadura de espesor total de 1 pulgada
• Área de Medición: 2 × 3.5 cm, tamaño de píxel de 1 mm
• Adquisición de Datos: 10 mediciones por píxel, 10 promedios temporales

Resultados Experimentales

Resultados de Simulación de Monte Carlo

Eficiencia de Dispersión y Factor de Asimetría: Con el aumento del tamaño de partícula, tanto la eficiencia de dispersión como el factor de asimetría aumentan en el segmento de baja frecuencia (<1 THz), indicando que partículas más grandes producen dispersión frontal más fuerte.

Descomposición de Matriz de Mueller: Se realizó descomposición Lu-Chipman de la matriz de Mueller para fantasmas A y D:

• Fantasma A: $M_{\Delta,A} \approx \text{diag}[1, 0.921, 0.925, 0.849]$
• Fantasma D: $M_{\Delta,D} \approx \text{diag}[1, 0.586, 0.586, 0.207]$

Los resultados indican que ambos fantasmas pueden aproximarse como una combinación de atenuador uniforme y despolarizador no uniforme.

Resultados de Verificación Experimental

Intensidad de Dispersión Difusa:

• Tanto la simulación como los experimentos muestran: cuanto mayor es el tamaño de partícula, mayor es la intensidad de dispersión
• Dependencia de Frecuencia: Intensidad alta en baja frecuencia, decayendo exponencialmente con la frecuencia
• Consistencia de Tendencia: Las tendencias de simulación y resultados experimentales coinciden

Contraste de Grado de Polarización (DOP):

• El fantasma D muestra DOP descendente a aproximadamente 0.6 cerca de 0.6 THz
• El DOP del fantasma C disminuye en un grado intermedio entre A, B y D
• Ubicación de Frecuencia: La ubicación de frecuencia de la disminución de DOP está relacionada con el tamaño de partícula

Análisis de Tejido Quemado

Contraste Espacial: A la frecuencia de 0.6 THz, el área quemada versus tejido sano muestra contraste evidente tanto en DOP como en intensidad de dispersión difusa.

Características Espectrales:

• El DOP del tejido sano disminuye más rápidamente en el rango de 0.2-0.8 THz
• La intensidad de dispersión difusa del área quemada es significativamente mayor que la del tejido sano
• Esto puede atribuirse a la destrucción de estructuras grandes de piel (folículos pilosos, glándulas sudoríparas) causada por quemaduras y pérdida de humedad

Trabajos Relacionados

Aplicaciones Biomédicas de Terahercios

La investigación existente se concentra principalmente en:

1. Diagnóstico de Cáncer: Estudios espectroscópicos de cáncer de mama ex vivo, carcinoma de células escamosas oral, etc.
2. Clasificación de Quemaduras de Piel: Clasificación de quemaduras in vivo y predicción de cicatrización
3. Análisis de Propiedades Dieléctricas: Enfoque principal en cambios de contenido de agua en tejidos

Espectroscopia de Dispersión Polarizada

Ampliamente aplicada en otros segmentos de longitud de onda:

1. Imágenes Endoscópicas: Detección de pólipos de colon y lesiones peritoneales
2. Diagnóstico de Cáncer: Múltiples subtipos de cáncer incluyendo cáncer de piel, gástrico y mamario
3. Modelado de Monte Carlo: Aplicación de modelos de dispersión de luz en biofotónica

Brecha Tecnológica

La investigación del mecanismo de respuesta de polarización en el régimen de terahercios está relativamente rezagada, careciendo de modelos computacionales maduros para explicar los fenómenos observados.

Conclusiones y Discusión

Conclusiones Principales

1. Verificación de Modelo Exitosa: Los resultados de simulación de Monte Carlo son altamente consistentes con datos experimentales en tendencia
2. Potencial Diagnóstico Confirmado: DOP e intensidad de dispersión difusa pueden distinguir partículas dispersoras de diferentes tamaños
3. Simplificación de Medición: Para partículas esféricas, la medición de un único estado de polarización es suficiente para obtener información completa de la matriz de Mueller
4. Viabilidad Clínica: Se demostró exitosamente la capacidad de detección de estados patológicos en tejido ex vivo

Análisis de Limitaciones

1. Limitación de Ancho de Banda del Sistema: El ancho de banda actual de 0.2-1.5 THz limita la capacidad de resolución de diferencias de pequeñas partículas
2. Restricción de Relación Señal-Ruido: La confiabilidad de mediciones de DOP por encima de 0.8 THz disminuye
3. Efecto de Rugosidad Superficial: Se asume superficie de muestra lisa y plana; en aplicaciones reales debe considerarse el efecto de despolarización de superficie
4. Requisito de Diferencia de Tamaño de Partícula: Se requiere diferencia de tamaño de partícula relativamente grande para producir contraste de señal evidente

Direcciones de Investigación Futura

1. Expansión de Ancho de Banda del Sistema: Utilizar cristales electroópticos para reemplazar PCA, logrando ancho de banda de 8 THz
2. Verificación de Tejido Canceroso Ex Vivo: Expandir investigación de fantasmas a tejido canceroso real
3. Optimización de Procesamiento de Señal: Desarrollar técnicas de descomposición de señal de onda y análisis estadístico de speckle
4. Separación de Despolarización Multifuente: Distinguir contribuciones relativas de rugosidad superficial versus estructuras internas

Evaluación Profunda

Fortalezas

1. Integración de Teoría y Experimento: La combinación de simulación de Monte Carlo con verificación experimental proporciona un marco de investigación completo
2. Innovación Metodológica: Primer estudio sistemático de propiedades de polarización de dispersión de Mie en el régimen de terahercios
3. Alto Valor Práctico: Simplificación de requisitos de medición de matriz de Mueller, mejorando la viabilidad de aplicación clínica
4. Verificación Suficiente: Estrategia de verificación progresiva desde fantasmas ideales a tejido biológico real es razonable
5. Análisis Profundo: La descomposición Lu-Chipman proporciona comprensión profunda del mecanismo físico

Insuficiencias

1. Diversidad Limitada de Muestras: Solo se utilizaron fantasmas de polipropileno-gelatina; las propiedades de material pueden diferir del tejido real
2. Análisis Estadístico Insuficiente: Falta de pruebas de significancia estadística detalladas y análisis de errores
3. Verificación Clínica Limitada: Solo se realizó verificación en una muestra de quemadura; se requiere investigación clínica a mayor escala
4. Consideración Insuficiente de Dispersión Compleja: Se consideran principalmente partículas esféricas; las estructuras tisulares reales son más complejas

Evaluación de Impacto

1. Valor Académico: Proporciona base teórica importante y métodos experimentales para biofotónica de terahercios
2. Avance Tecnológico: Impulsa la aplicación de tecnología de imágenes polarimétrica de terahercios en diagnóstico médico
3. Significado Interdisciplinario: Conecta teoría de dispersión electromagnética, modelado de Monte Carlo e imágenes biomédicas
4. Perspectiva Industrial: Proporciona ruta tecnológica para desarrollo de dispositivos médicos de terahercios

Escenarios Aplicables

1. Diagnóstico de Enfermedades de Piel: Clasificación de quemaduras, detección de cáncer de piel
2. Identificación de Márgenes Tumorales: Navegación quirúrgica, análisis patológico
3. Caracterización de Estructura Tisular: Análisis de fibras de colágeno, cúmulos celulares
4. Monitoreo de Eficacia de Medicamentos: Seguimiento de cambios de estructura tisular durante tratamiento

Referencias

El artículo cita 41 referencias relacionadas, abarcando trabajos importantes en múltiples campos incluyendo tecnología de terahercios, óptica biomédica y modelado de Monte Carlo, proporcionando base teórica sólida para la investigación.

Evaluación General: Este es un artículo de investigación de alta calidad interdisciplinaria con valor teórico y práctico importante en el campo de imágenes biomédicas de terahercios. La metodología de investigación es rigurosa, los resultados son confiables, y proporciona una base sólida para el desarrollo futuro del campo.