1 / 24

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.

magnar
Download Presentation

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

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 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. SMITH11Academic Unit of Medical Physics, University of Leeds, UK2Department of Medical Physics and Engineering, Leeds Teaching Hospitals NHS Trust, UK3Institute of Microwaves and Photonics,University of Leeds, UK(*Author for correspondence, email:gcw@medphysics.leeds.ac.uk)

  2. 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.

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

  4. 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.

  5. The modelling problem Incident THz Spectrum Transmitted THz Spectrum

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

  7. 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.

  8. Physical Properties of Water/Propanol-1

  9. 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.

  10. Thin Film Matrix Model

  11. Thin Film Matrix Model

  12. Thin Film Matrix Model

  13. 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.

  14. 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.

  15. 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.

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

  17. Results - Full Spectrum

  18. Results - Up to 1 THz

  19. Results - Full Spectrum

  20. 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.

  21. 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.

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

  23. 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.

  24. 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.

More Related