Penetration of THz Radiation in Biological Tissue

Terahertz technology has revolutionized many research fields. Much of the research in biomedical applications focuses on radiation-tissue interaction using submillimeter wavelengths in the THz gap range. This technology has a promising future, but good equipment is expensive.
Therefore, this Demonstration investigates a modified bio-heat model equation, which includes the calculation of factors that affect the rate of heat generation per unit volume in tumor cell tissue (). This model analyzes the power effect on the THz radiation penetration in biological tissue by analyzing the scattering, which depends on the angle and wavelength. We study and test our model using THz frequencies and submillimeter wavelengths in the range from 0.1 mm to 30 mm. Because the interaction between THz radiation and tissue involves photon radiation absorption and scattering, the most important parameter is power density (), which is a function of the power delivered by the source and the spot size. In fact, the product of fluence and absorption coefficient equals the heat from the source, which is the amount of energy deposited in a unit volume of tissue.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


When tissue is penetrated, the dominant effects are scattering and absorption by polarizable molecules. The most important parameter in the interaction is the power density, (). This is a function of power source delivered from a THz radiation source to the tissue, where input power is (W), and the spot has a radius (mm). The propagation of the scattered radiation is described by the photon transport equation
where the rate of radiance change ) at a point is indicated by in the direction as a result of light incident in direction . This equation is solved using the simplest isotropic phase function, , where is the attenuation coefficient (), is the scattering coefficient, and is the absorption coefficient. At larger distances from the dipole, the field is perpendicular to the direction of propagation and falls off as . The scattering amplitude is given by
where is the potential energy. The amplitude depends on the potential, the scattering angle , and the target area . We can expand the plane waves in terms of Legendre polynomials to find an expression from the scattering amplitude and the differential cross section
The total cross section is . We verify directly that , which is called the optical theorem.
Snapshots 1, 2, and 3: The graphs show that the change from a shorter to a longer wavelength increases the heat production rate in a tumor cell volume slightly (measured in ). This means that shorter wavelengths contribute to the heat production more than the longer wavelengths. The longer wavelengths remain constant, with variation of tissue depth due to penetration of THz photon radiation , and also due to properties of tissue optical depth, . Other factors include the variation of the THz radiation power source and the variation of the wavelength , as well as the radiation photon incidence angle . Because the heat production is caused by the THz radiation beam propagating through the tissue with strong water absorption, the intensity is attenuated exponentially, due to the low scattering factor.
The following assumptions are made in this model: First, the frequency of the scattering radiation remains exactly the same as that of the incident radiation. Therefore, the scattering process is elastic. Second, any interaction between the scattered radiation particles themselves is neglected. Third, possible multiple scattering processes are likewise neglected. Fourth, the incident beam width is much larger than a typical range of the scattering potential, so that the particle will have a well-defined momentum.
[1] W. Hoppe, W. Lohmann, H. Markl, and H. Ziegler, Biophysics, New York: Springer-Verlag, 1983.
[2] G. Videen and D. Ngo, "Light Scattering from a Cell," in Optics of Biological Particles, NATO Science Series, Series II: Mathematics, Physics and Chemistry, Vol. 238, (A. Hoekstra, V. Maltsev, and G. Videen, eds.), New York: Springer, 2007.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+