Comparing finite difference time domain and Monte Carlo modeling of human skin interaction with terahertz radiation

Assessing the biological reaction to electromagnetic (EM) radiation of all frequencies and intensities is essential to the understanding of both the potential damage caused by the radiation and the inherent mechanisms within biology that respond, protect, or propagate damage to surrounding tissues. To understand this reaction, one may model the electromagnetic irradiation of tissue phantoms according to empirically measured or intelligently estimated dielectric properties. Of interest in this study is the terahertz region (0.2-2.0 THz), ranging from millimeter to infrared waves, which has been studied only recently due to lack of efficient sources. The specific interaction between this radiation and human tissue is greatly influenced by the significant EM absorption of water across this range, which induces a pronounced heating of the tissue surface. This study compares the Monte Carlo Multi-Layer (MCML) and Finite Difference Time Domain (FDTD) approaches for modeling the terahertz irradiation of human dermal tissue. Two congruent simulations were performed on a one-dimensional tissue model with unit power intensity profile. This works aims to verify the use of either technique for modeling terahertz-tissue interaction for minimally scattering tissues.

[1]  Donald D. Arnone,et al.  Terahertz-pulsed imaging of cancers , 2003, SPIE BiOS.

[2]  J. M. Chamberlain,et al.  The interaction between Terahertz radiation and biological tissue. , 2001, Physics in medicine and biology.

[3]  L Wang,et al.  MCML--Monte Carlo modeling of light transport in multi-layered tissues. , 1995, Computer methods and programs in biomedicine.

[4]  E. Pickwell‐MacPherson,et al.  Simulation of Terahertz Pulse Propagation in Biological Systems , 2004 .

[5]  Robert E. Miles,et al.  Optical properties of tissue measured using terahertz-pulsed imaging , 2003, SPIE Medical Imaging.

[6]  R. E. Miles,et al.  Modelling the propagation of terahertz radiation through a tissue simulating phantom. , 2004, Physics in medicine and biology.

[7]  V. Wallace,et al.  In vivo study of human skin using pulsed terahertz radiation , 2004, Physics in medicine and biology.

[8]  Robert E. Miles,et al.  Two methods for modeling the propagation of terahertz radiation in a layered structure , 2003, SPIE Medical Imaging.

[9]  Ruth M. Woodward,et al.  An analytical model of skin: comparison with experimental results in vivo , 2005, SPIE BiOS.

[10]  DeChang Dai,et al.  Terahertz scattering: comparison of a novel theoretical approach with experiment , 2006, SPIE OPTO.

[11]  P. Taday,et al.  Simulating the response of terahertz radiation to basal cell carcinoma using ex vivo spectroscopy measurements. , 2005, Journal of biomedical optics.

[12]  K. Yee Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[13]  Dwight L. Woolard,et al.  Terahertz spectroscopic characterization of cancer cells , 2005, SPIE BiOS.

[14]  Dustin G. Mixon,et al.  Electromagnetic properties of tissue in the optical region , 2007, SPIE BiOS.