Рассеяние Ми, вызванное структурами ткани в терагерцовом диапазоне: экспериментальная и Монте-Карло верификация с использованием диффузной поляриметрической визуализации в высокоослабляющих тканевых фантомах

Основная информация

• 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 (Stony Brook University)
• Классификация: physics.med-ph physics.bio-ph physics.optics
• Дата подачи: 31 декабря 2024 г.
• Ссылка на статью: https://arxiv.org/abs/2501.00630

Аннотация

В данном исследовании рассматривается явление рассеяния Ми, вызванное структурами ткани в терагерцовом диапазоне, с акцентом на потенциальные приложения в диагностике заболеваний. Исследование проводилось путём моделирования Монте-Карло и экспериментальной верификации характеристик рассеяния поляризованного терагерцового света в высокоослабляющих тканях. Исследование выявило явный контраст между рассеивающими частицами различных размеров в отношении степени поляризации (DOP) и интенсивности диффузного рассеяния, а также подтвердило диагностический потенциал технологии на образцах ожогов свежей кожи свиньи.

Научная база и мотивация

Определение проблемы

Изменения структуры ткани являются важной характеристикой многих патологических процессов, таких как опухолевые почки (tumor budding) на границах опухолей и ожоговые повреждения. Традиционные приложения терагерцовой биомедицины сосредоточены главным образом на диэлектрических свойствах ткани, в частности на содержании воды, и редко используют рассеивающие свойства терагерцового света в качестве диагностического инструмента.

Научная значимость

1. Потребность в диагностике заболеваний: опухолевые почки и кластеры недифференцированных клеток являются независимыми прогностическими факторами лимфатического метастазирования и выживаемости пациентов
2. Технологические ограничения: существующие методы, такие как окрашивание H&E и иммуногистохимия, требуют срезов ткани и не позволяют проводить диагностику в реальном времени
3. Неясные механизмы рассеяния: механизмы различий в поляризационной реакции в терагерцовом диапазоне ещё недостаточно изучены

Инновационная мотивация

С развитием терагерцовых эллипсометров и систем поляриметрического измерения возникает острая необходимость в вычислительных моделях для объяснения наблюдаемых явлений и обеспечения теоретической поддержки приложений терагерцовой биофотоники.

Основные вклады

1. Разработка модели Монте-Карло для рассеяния Ми в терагерцовом диапазоне: на основе модели, разработанной Xu и Arbab, проведены улучшения для рассеивающих частиц в поглощающей среде
2. Верификация диагностического потенциала поляриметрических измерений: экспериментальное подтверждение контраста для рассеивающих частиц различных размеров на тканевых фантомах
3. Предоставление анализа разложения матрицы Мюллера: использование метода разложения Lu-Chipman для анализа поляризационных свойств ткани
4. Демонстрация практического применения: верификация осуществимости технологии на образцах ожогов свежей кожи свиньи
5. Упрощение требований к измерениям: доказано, что для сферических частиц достаточно измерения одного поляризационного состояния для построения полной матрицы Мюллера

Методологические подробности

Конструкция экспериментальной системы

Система терагерцовой спектроскопии во временной области:

• Передающая часть: использование фемтосекундного лазера 1560 нм для возбуждения фотопроводящей антенны (PCA)
• Приёмная часть: поворотный детекторный рычаг 0-160 градусов, включающий линейный поляризатор (WGP)
• Диапазон частот: 0,2-1,5 ТГц

Конфигурация измерения рассеяния:

• Выбран угол 140 градусов для измерения диффузного рассеяния (максимальный контраст между когерентной и некогерентной мощностью)
• Растровое сканирование: шаг пикселя 2 мм, область 2 см × 2 см
• Запись 100 измерений для каждого поляризационного состояния

Подготовка тканевых фантомов

Выбор материалов:

• Матрица: 15% раствор бычьего желатина (моделирование фона ткани)
• Рассеивающие частицы: полипропиленовые гранулы (моделирование волосяных фолликулов, потовых желёз и других структур с низким показателем преломления)
• Диапазон размеров частиц: 115-280 мкм (соответствует параметрам рассеяния Ми 0,2-20)

Технология подготовки:

1. Нагревание раствора желатина до 40°C
2. Быстрое смешивание с полипропиленовыми гранулами
3. Заливка в чашку Петри для отверждения (толщина ~5 мм)
4. Анализ распределения размеров частиц и концентрации с использованием оптической микроскопии и программного обеспечения ImageJ

Моделирование Монте-Карло

Основной алгоритм: На основе кода Монте-Карло для меридиональной поляризованной световой волны, разработанного Ramella-Roman и соавторами, с модификациями для рассеяния Ми в поглощающей среде.

Расчёт ключевых параметров:

• Коэффициент рассеяния: $\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]$

Построение матрицы Мюллера: Использование четырёх различных состояний падающей поляризации (горизонтальная линейная поляризация, вертикальная линейная поляризация, линейная поляризация под углом 45°, правая круговая поляризация) для расчёта пространственных параметров Стокса и последующего построения матрицы Мюллера.

Экспериментальная установка

Параметры фантомов

Подготовлены четыре различных тканевых фантома:

Измерение характеристик материалов

Диэлектрические свойства:

• Полипропилен: показатель преломления ~1,51, коэффициент поглощения <2 см⁻¹
• Желатин: показатель преломления и коэффициент поглощения находятся между значениями для кожи и воды

Эксперимент с ожогами свежей ткани

• Образец: свежая кожа свиньи
• Подготовка ожога: воздействие латунного стержня при 180°C в течение 5 минут, образование полнослойного ожога размером 1 дюйм
• Область измерения: 2 × 3,5 см, размер пикселя 1 мм
• Сбор данных: 10 измерений на пиксель, 10-кратное усреднение по времени

Результаты экспериментов

Результаты моделирования Монте-Карло

Эффективность рассеяния и фактор асимметрии: С увеличением размера частиц эффективность рассеяния и фактор асимметрии увеличиваются в низкочастотном диапазоне (<1 ТГц), что указывает на более сильное переднее рассеяние для более крупных частиц.

Разложение матрицы Мюллера: Разложение Lu-Chipman матриц Мюллера для фантомов A и D:

• Фантом A: $M_{\Delta,A} \approx \text{diag}[1, 0.921, 0.925, 0.849]$
• Фантом D: $M_{\Delta,D} \approx \text{diag}[1, 0.586, 0.586, 0.207]$

Результаты показывают, что оба фантома могут быть аппроксимированы как комбинация однородного ослабителя и неоднородного деполяризатора.

Результаты экспериментальной верификации

Интенсивность диффузного рассеяния:

• Как моделирование, так и эксперименты показывают: чем больше размер частиц, тем выше интенсивность рассеяния
• Частотная зависимость: высокая интенсивность на низких частотах, экспоненциальное затухание с увеличением частоты
• Согласованность тенденций: результаты моделирования и эксперимента совпадают по тенденции

Контраст степени поляризации (DOP):

• Фантом D показывает снижение DOP до ~0,6 вблизи 0,6 ТГц
• Фантом C демонстрирует промежуточное снижение DOP между A, B и D
• Частотная локализация: частота, при которой происходит снижение DOP, связана с размером частиц

Анализ ожоговой ткани

Пространственный контраст: На частоте 0,6 ТГц область ожога и здоровая ткань демонстрируют явный контраст как в DOP, так и в интенсивности диффузного рассеяния.

Спектральные характеристики:

• DOP здоровой ткани снижается быстрее в диапазоне 0,2-0,8 ТГц
• Интенсивность диффузного рассеяния в области ожога значительно выше, чем в здоровой ткани
• Это может быть обусловлено разрушением крупных структур кожи (волосяные фолликулы, потовые железы) вследствие ожога и потерей влаги

Связанные работы

Приложения терагерцовой биомедицины

Существующие исследования сосредоточены на:

1. Диагностика рака: спектроскопические исследования рака молочной железы, плоскоклеточного рака полости рта и других типов
2. Классификация ожогов кожи: оценка степени ожога in vivo и прогнозирование заживления
3. Анализ диэлектрических свойств: главным образом сосредоточены на изменениях содержания воды в ткани

Поляризационная спектроскопия рассеяния

Широко применяется в других спектральных диапазонах:

1. Эндоскопическая визуализация: обнаружение полипов толстой кишки и поражений брюшины
2. Диагностика рака: меланома, рак желудка, рак молочной железы и другие подтипы рака
3. Моделирование Монте-Карло: применение моделей светорассеяния в биофотонике

Технологический пробел

Исследование механизмов поляризационной реакции в терагерцовом диапазоне отстаёт, отсутствуют зрелые вычислительные модели для объяснения наблюдаемых явлений.

Выводы и обсуждение

Основные выводы

1. Успешная верификация модели: результаты моделирования Монте-Карло высоко согласуются с экспериментальными данными по тенденции
2. Подтверждение диагностического потенциала: DOP и интенсивность диффузного рассеяния способны различать рассеивающие частицы различных размеров
3. Упрощение измерений: для сферических частиц измерение одного поляризационного состояния достаточно для получения полной информации матрицы Мюллера
4. Клиническая осуществимость: успешная демонстрация способности обнаружения патологических состояний на образцах ex vivo

Анализ ограничений

1. Ограничение полосы пропускания системы: текущая полоса пропускания 0,2-1,5 ТГц ограничивает способность разрешения различий в размерах малых частиц
2. Ограничения отношения сигнал-шум: надёжность измерений DOP выше 0,8 ТГц снижается
3. Влияние шероховатости поверхности: предполагается гладкая и ровная поверхность образца, в практических приложениях необходимо учитывать эффект деполяризации поверхности
4. Требование различия размеров частиц: требуется значительное различие в размерах частиц для создания явного контраста сигнала

Направления будущих исследований

1. Расширение полосы пропускания системы: использование электрооптических кристаллов вместо PCA для достижения полосы пропускания 8 ТГц
2. Верификация на образцах раковой ткани ex vivo: расширение исследований фантомов на реальную раковую ткань
3. Оптимизация обработки сигналов: разработка методов разложения волновых сигналов и анализа статистических спеклов
4. Разделение множественных источников деполяризации: различение относительного вклада шероховатости поверхности и внутренних структур

Глубокая оценка

Преимущества

1. Сочетание теории и эксперимента: моделирование Монте-Карло в сочетании с экспериментальной верификацией обеспечивает полную исследовательскую базу
2. Методологические инновации: первое систематическое исследование поляризационных свойств рассеяния Ми в терагерцовом диапазоне
3. Высокая практическая ценность: упрощение требований к измерению матрицы Мюллера повышает осуществимость клинического применения
4. Достаточная верификация: разумная стратегия постепенной верификации от идеальных фантомов к реальным биологическим тканям
5. Глубокий анализ: разложение Lu-Chipman обеспечивает глубокое понимание физических механизмов

Недостатки

1. Ограниченное разнообразие образцов: использование только полипропилен-желатиновых фантомов, характеристики материалов могут отличаться от реальных тканей
2. Недостаточный статистический анализ: отсутствие подробных тестов статистической значимости и анализа ошибок
3. Ограниченная клиническая верификация: верификация проведена только на одном образце ожога, требуются крупномасштабные клинические исследования
4. Недостаточное рассмотрение сложного рассеяния: главным образом рассматриваются сферические частицы, реальные структуры ткани более сложны

Оценка влияния

1. Академическая ценность: обеспечение важной теоретической базы и экспериментальных методов для терагерцовой биофотоники
2. Технологический прогресс: продвижение применения терагерцовой поляриметрической визуализации в медицинской диагностике
3. Междисциплинарное значение: связь теории электромагнитного рассеяния, моделирования Монте-Карло и биомедицинской визуализации
4. Промышленные перспективы: обеспечение технологического пути для разработки терагерцовых медицинских устройств

Применимые сценарии

1. Диагностика кожных заболеваний: классификация ожогов, обнаружение рака кожи
2. Идентификация границ опухоли: навигация при хирургии, патологический анализ
3. Характеризация структуры ткани: анализ коллагеновых волокон, клеточных кластеров
4. Мониторинг эффективности лекарств: отслеживание изменений структуры ткани в процессе лечения

Список литературы

Статья цитирует 41 соответствующую литературу, охватывающую множество областей, включая терагерцовые технологии, биомедицинскую оптику, моделирование Монте-Карло и другие важные работы, обеспечивая прочную теоретическую базу для исследования.

Общая оценка: Это высококачественная междисциплинарная исследовательская статья, имеющая важное теоретическое и практическое значение в области терагерцовой биомедицинской визуализации. Методология исследования строга, результаты надёжны, что создаёт хорошую основу для дальнейшего развития этой области.