Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms

Basic Information

• Paper ID: 2501.00630
• Title: Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
• Authors: Erica Heller, Kuangyi Xu, Zachery Harris, M. Hassan Arbab (Stony Brook University)
• Classification: physics.med-ph physics.bio-ph physics.optics
• Submission Date: December 31, 2024
• Paper Link: https://arxiv.org/abs/2501.00630

Abstract

This study investigates Mie scattering phenomena induced by tissue structures in the terahertz regime, with emphasis on potential applications in disease diagnosis. Through Monte Carlo simulations and experimental verification, the research analyzes the scattering characteristics of polarized terahertz light in highly attenuating tissue. The findings demonstrate that scattering particles of different sizes exhibit distinct contrasts in degree of polarization (DOP) and diffuse scattering intensity, with diagnostic potential verified in ex vivo porcine skin burn samples.

Research Background and Motivation

Problem Definition

Changes in tissue structure are important characteristics of many disease processes, such as tumor budding at cancer margins and burn injuries. Traditional terahertz biomedical applications have primarily focused on the dielectric properties of tissue, particularly water content, while rarely exploiting the scattering characteristics of terahertz light as a diagnostic tool.

Research Significance

1. Disease Diagnosis Requirements: Tumor budding and poorly differentiated cell clusters are independent prognostic factors for lymph node metastasis and patient survival
2. Technical Limitations: Existing methods such as H&E staining and immunohistochemistry require tissue sections and cannot enable real-time diagnosis
3. Unclear Scattering Mechanisms: The mechanisms underlying polarization response differences in the terahertz regime remain insufficiently studied

Innovation Motivation

With the development of terahertz ellipsometry and polarization measurement systems, there is an urgent need for computational models to explain observed phenomena and provide theoretical support for terahertz biophotonics applications.

Core Contributions

1. Established Monte Carlo model for terahertz Mie scattering: Improved upon the model developed by Xu and Arbab, specifically for scattering particles in absorbing media
2. Verified diagnostic potential of polarization measurements: Demonstrated through tissue phantom experiments the contrast of scattering particles of different sizes
3. Provided Mueller matrix decomposition analysis: Analyzed tissue polarization properties using the Lu-Chipman decomposition method
4. Demonstrated practical applications: Verified technical feasibility in ex vivo porcine skin burn samples
5. Simplified measurement requirements: Demonstrated that for spherical particles, Mueller matrix can be constructed from measurements of a single polarization state

Methodology Details

Experimental System Design

Terahertz Time-Domain Spectroscopy System:

• Transmitter: 1560 nm femtosecond laser exciting photoconductive antenna (PCA)
• Receiver: Rotatable detection arm (0-160 degrees) with wire grid polarizer (WGP)
• Frequency range: 0.2-1.5 THz

Scattering Measurement Configuration:

• Selected 140-degree angle for diffuse scattering measurements (maximum contrast between coherent and incoherent power)
• Grid scanning: 2 mm pixel spacing, 2 cm × 2 cm area
• 100 measurements recorded per polarization state

Tissue Phantom Preparation

Material Selection:

• Matrix: 15% bovine skin gelatin solution (simulating tissue background)
• Scattering particles: polypropylene beads (simulating low refractive index structures such as hair follicles and sweat glands)
• Particle size range: 115-280 μm (corresponding to Mie scattering parameters 0.2-20)

Preparation Protocol:

1. Heat gelatin solution to 40°C
2. Rapidly mix in polypropylene beads
3. Pour into culture dishes for solidification (approximately 5 mm thickness)
4. Analyze particle size distribution and concentration using optical microscopy and ImageJ software

Monte Carlo Simulation

Core Algorithm: Based on the meridional polarized light Monte Carlo code of Ramella-Roman et al., modified for Mie scattering in absorbing media.

Key Parameter Calculations:

• Scattering coefficient: $\mu_s = \frac{2\pi\rho}{|n_m|^2k_0^2}\sum_{j=1}^{\infty}(2j+1)(|a_j|^2+|b_j|^2)$
• Absorption coefficient: $\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]$

Mueller Matrix Construction: Four different incident polarization states (horizontal linear polarization, vertical linear polarization, 45° linear polarization, right circular polarization) were used to calculate spatial Stokes parameters, from which the Mueller matrix was constructed.

Experimental Setup

Phantom Parameters

Four tissue phantoms with different characteristics were prepared:

Material Characterization

Dielectric Properties:

• Polypropylene: refractive index ~1.51, absorption coefficient <2 cm⁻¹
• Gelatin: refractive index and absorption coefficient intermediate between skin and water

Ex Vivo Burn Experiments

• Sample: ex vivo porcine skin
• Burn preparation: 180°C brass rod applied for 5 minutes, creating 1-inch full-thickness burn
• Measurement area: 2×3.5 cm, 1 mm pixel size
• Data acquisition: 10 measurements per pixel, 10-time temporal averaging

Experimental Results

Monte Carlo Simulation Results

Scattering Efficiency and Asymmetry Factor: With increasing particle size, both scattering efficiency and asymmetry factor increase at low frequencies (<1 THz), indicating that larger particles produce stronger forward scattering.

Mueller Matrix Decomposition: Lu-Chipman decomposition of Mueller matrices for phantoms A and D:

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

Results indicate both phantoms can be approximated as combinations of uniform attenuators and non-uniform depolarizers.

Experimental Verification Results

Diffuse Scattering Intensity:

• Both simulation and experiment show: larger particle size yields higher scattering intensity
• Frequency dependence: high intensity at low frequencies, exponential decay with increasing frequency
• Trend consistency: simulation and experimental results show consistent trends

Degree of Polarization (DOP) Contrast:

• Phantom D shows DOP decreasing to approximately 0.6 near 0.6 THz
• Phantom C's DOP decrease is intermediate between A, B, and D
• Frequency localization: the frequency position of DOP decrease is related to particle size

Burn Tissue Analysis

Spatial Contrast: At 0.6 THz frequency, burn regions and healthy tissue show distinct contrasts in both DOP and diffuse scattering intensity.

Spectral Characteristics:

• Healthy tissue DOP decreases more rapidly in the 0.2-0.8 THz range
• Burn region diffuse scattering intensity is significantly higher than healthy tissue
• This may be attributed to destruction of large skin structures (hair follicles, sweat glands) and water loss caused by burns

Related Work

Terahertz Biomedical Applications

Existing research primarily focuses on:

1. Cancer Diagnosis: Spectroscopic studies of ex vivo breast cancer, oral squamous cell carcinoma, etc.
2. Skin Burn Classification: In vivo burn grading and healing prediction
3. Dielectric Property Analysis: Primarily concerned with tissue water content changes

Polarized Scattering Spectroscopy

Widely applied in other spectral regions:

1. Endoscopic Imaging: Colonic polyp and peritoneal lesion detection
2. Cancer Diagnosis: Multiple cancer subtypes including skin cancer, gastric cancer, and breast cancer
3. Monte Carlo Modeling: Light scattering models in biophotonics applications

Technical Gaps

Research on polarization response mechanisms in the terahertz regime lags behind, lacking mature computational models to explain observed phenomena.

Conclusions and Discussion

Main Conclusions

1. Successful Model Verification: Monte Carlo simulation results show high consistency with experimental data in trend
2. Diagnostic Potential Confirmed: DOP and diffuse scattering intensity can distinguish scattering particles of different sizes
3. Measurement Simplification: For spherical particles, complete Mueller matrix information can be obtained from single polarization state measurements
4. Clinical Feasibility: Successfully demonstrated pathological state detection capability in ex vivo tissue

Limitations Analysis

1. System Bandwidth Limitation: Current 0.2-1.5 THz bandwidth limits resolution of small particle size differences
2. Signal-to-Noise Ratio Constraints: DOP measurement reliability decreases above 0.8 THz
3. Surface Roughness Effects: Assumes smooth sample surfaces; actual applications must consider surface depolarization effects
4. Particle Size Difference Requirements: Requires substantial particle size differences to produce obvious signal contrast

Future Research Directions

1. Extended System Bandwidth: Use electro-optic crystals to replace PCA, achieving 8 THz bandwidth
2. Ex Vivo Cancer Tissue Verification: Extend phantom studies to real cancer tissue
3. Signal Processing Optimization: Develop wavelet signal decomposition and statistical speckle analysis techniques
4. Multi-Source Depolarization Separation: Distinguish relative contributions of surface roughness and internal structures

In-Depth Evaluation

Strengths

1. Theory-Experiment Integration: Combines Monte Carlo simulation with experimental verification, providing a complete research framework
2. Methodological Innovation: First systematic study of polarization properties of Mie scattering in the terahertz regime
3. High Practical Value: Simplifies Mueller matrix measurement requirements, improving clinical application feasibility
4. Sufficient Verification: Reasonable progressive verification strategy from ideal phantoms to real biological tissue
5. In-Depth Analysis: Lu-Chipman decomposition provides deeper understanding of physical mechanisms

Limitations

1. Limited Sample Diversity: Uses only polypropylene-gelatin phantoms; material properties may differ from real tissue
2. Insufficient Statistical Analysis: Lacks detailed statistical significance testing and error analysis
3. Limited Clinical Verification: Verification on only one burn sample; requires larger-scale clinical studies
4. Insufficient Complex Scattering Consideration: Primarily considers spherical particles; actual tissue structures are more complex

Impact Assessment

1. Academic Value: Provides important theoretical foundation and experimental methods for terahertz biophotonics
2. Technical Advancement: Promotes application of terahertz polarization imaging technology in medical diagnosis
3. Interdisciplinary Significance: Connects electromagnetic scattering theory, Monte Carlo modeling, and biomedical imaging
4. Industrial Prospects: Provides technical pathway for development of terahertz medical devices

Applicable Scenarios

1. Dermatological Diagnosis: Burn grading, skin cancer detection
2. Tumor Margin Identification: Surgical navigation, pathological analysis
3. Tissue Structure Characterization: Collagen fiber, cell cluster analysis
4. Drug Efficacy Monitoring: Tracking tissue structure changes during treatment

References

The paper cites 41 relevant references covering important works in terahertz technology, biomedical optics, Monte Carlo modeling, and other fields, providing a solid theoretical foundation for the research.

Overall Assessment: This is a high-quality interdisciplinary research paper with significant theoretical and practical value in the field of terahertz biomedical imaging. The research methodology is rigorous, results are credible, and it establishes a solid foundation for further development in this field.