Mo 17 Walker

Two methods for modelling the propagation of terahertz radiation in a layered structure.: 

Two methods for modelling the propagation of terahertz radiation in a layered structure. GILLIAN C. WALKER1*, ELIZABETH BERRY1, STEPHEN W. SMYE2, NICK N. ZINOV’EV3, ANTHONY J. FITZGERALD1, ROBERT. E. MILES3, MARTYN CHAMBERLAIN3 AND MICHAEL A. SMITH1 1Academic Unit of Medical Physics, University of Leeds, UK 2Department of Medical Physics and Engineering, Leeds Teaching Hospitals NHS Trust, UK 3Institute of Microwaves and Photonics, University of Leeds, UK (*Author for correspondence, email:gcw@medphysics.leeds.ac.uk)

Objectives: 

Objectives Create a modelling tool to simulate the passage of THz radiation through biological tissue. Biological tissue is highly inhomogeneous. The interactions that occur with THz and biological tissue are complex.

Outline: 

Outline Modelling Biological Tissue In Vitro Phantom Thin Film Matrix Model Monte Carlo Model Results Discussion Future Work

Modelling the interaction of THz radiation with biological tissue.: 

Modelling the interaction of THz radiation with biological tissue. THz radiation is being investigated as an imaging tool for skin. It has been shown that TPI can resolve the stratum corneum, epidermis and dermis. (Cole et al. Laser Plasma Generation and Diagnostics, SPIE Proc 2001; 4286.) A three-layer system of parallel sided slabs, each with frequency dependent physical properties could be used to simulate human skin.

The modelling problem: 

The modelling problem Incident THz Spectrum Transmitted THz Spectrum

In vitro phantom : 

In vitro phantom Water/Propanol-1 solution Spacer - 180 m TPX - 2 mm

Physical Properties of Water/Propanol -1: 

Physical Properties of Water/Propanol -1 The absorption coefficient and index of refraction of water and propanol-1 were calculated using the Cole-Cole model. (Kindt et al. Journal of Physical Chemistry 1996;100 :10373-9) These were averaged using volume weighting to give the physical properties of the specific water/propanol-1 solution.

Physical Properties of Water/Propanol-1 : 

Physical Properties of Water/Propanol-1

Thin Film Matrix Model: 

Thin Film Matrix Model A method for calculating the change in electric field as it passes through the layered medium. Implementation of the boundary conditions for the electric and magnetic components of the incident radiation at each boundary result in a matrix formulation of the problem.

Thin Film Matrix Model: 

Thin Film Matrix Model

Thin Film Matrix Model : 

Thin Film Matrix Model

Thin Film Matrix Model : 

Thin Film Matrix Model

Monte Carlo Model: 

Monte Carlo Model Creating a Photon Distribution. A THz pulse is recorded in units proportional to electric field. A photon distribution is created by randomly sampling the spectrum and a Poisson distribution to account for the coherent nature of the radiation. One million photons were included in an incident ensemble.

The Monte Carlo Model : 

The Monte Carlo Model The position of each photon is tracked in the sample. The probability of a photon crossing a boundary within the sample is determined by the Fresnel coefficients. In the water/propanol-1 solution the Beer-Lambert law is sampled as the probability distribution for an absorption event. The number of photons transmitted and reflected is counted to give the output spectra.

Presentation of Results: 

Presentation of Results The results of the Monte Carlo simulation were in expressed as a photon distribution while the experimental results were in arbitrary units proportional to electric field. The experimental results were converted into photon distributions. The number of photons included in each respective ensemble were calculated as a fraction of one million, the fraction given by the experimental area to the incident spectrum area.

Graphical comparison: 

Graphical comparison Transmitted spectra Incident spectra 1 000 000 photons a1 a2 a2/a1*1 000 000 photons

Results - Full Spectrum: 

Results - Full Spectrum

Results - Up to 1 THz: 

Results - Up to 1 THz

Results - Full Spectrum: 

Results - Full Spectrum

Discussion: 

Discussion The model results show good agreement with the experimental results reproducing all major features. Up to 1 THz, where the physical parameters have been verified there is close agreement with model and simulated amplitude. Generally the Thin Film Matrix Model more closely reproduces the experimental results.

Future Work: 

Future Work Investigations into how the absorption coefficient and index of refraction affect the simulated results are being carried out. Implementation of a fitting routine to extract physical parameters from the tissue.

Future Work: 

Future Work 50%-250% of the absorption coefficient Experimental result

Future Work: 

Future Work The Monte Carlo model is to be used to investigate scattering of THz radiation. The models are to be used to simulate results from human skin for both transmission and reflection data, in vitro and in vivo.

Acknowledgements: 

Acknowledgements This work was supported by the Engineering and Physical Sciences Research Council (GR/N39678) We are grateful for the contributions of the members of the EU Teravision project (IST-1999-10154), especially W. Th. Wenckebach, T.U. Delft.