OSA's Digital Library

Biomedical Optics Express

Biomedical Optics Express

  • Editor: Joseph A. Izatt
  • Vol. 2, Iss. 5 — May. 1, 2011
  • pp: 1169–1183

Modeling hemoglobin at optical frequency using the unconditionally stable fundamental ADI-FDTD method

Ding Yu Heh and Eng Leong Tan  »View Author Affiliations

Biomedical Optics Express, Vol. 2, Issue 5, pp. 1169-1183 (2011)

View Full Text Article

Enhanced HTML    Acrobat PDF (1497 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



This paper presents the modeling of hemoglobin at optical frequency (250 nm – 1000 nm) using the unconditionally stable fundamental alternating-direction-implicit finite-difference time-domain (FADI-FDTD) method. An accurate model based on complex conjugate pole-residue pairs is proposed to model the complex permittivity of hemoglobin at optical frequency. Two hemoglobin concentrations at 15 g/dL and 33 g/dL are considered. The model is then incorporated into the FADI-FDTD method for solving electromagnetic problems involving interaction of light with hemoglobin. The computation of transmission and reflection coefficients of a half space hemoglobin medium using the FADI-FDTD validates the accuracy of our model and method. The specific absorption rate (SAR) distribution of human capillary at optical frequency is also shown. While maintaining accuracy, the unconditionally stable FADI-FDTD method exhibits high efficiency in modeling hemoglobin.

© 2011 OSA

OCIS Codes
(170.3660) Medical optics and biotechnology : Light propagation in tissues
(050.1755) Diffraction and gratings : Computational electromagnetic methods

ToC Category:
Optics of Tissue and Turbid Media

Original Manuscript: January 4, 2011
Revised Manuscript: March 18, 2011
Manuscript Accepted: March 22, 2011
Published: April 12, 2011

Ding Yu Heh and Eng Leong Tan, "Modeling hemoglobin at optical frequency using the unconditionally stable fundamental ADI-FDTD method," Biomed. Opt. Express 2, 1169-1183 (2011)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. W. G. Zijlstra, A. Buursma, and W. P. Meeuwsen-van der Roest, “Absorption spectra of human fetal and adult oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin, and methemoglobin,” Clin. Chem. 37, 1633–1668 (1991).
  2. W. G. Zijlstra, A. Buursma, H. E. Falke, and J. F. Catsburg, “Spectrophotometry of hemoglobin: absorption spectra of rat oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin, and methemoglobin,” Comput. Biochem. Physiol. 107B, 161–166 (1994). [CrossRef]
  3. W. G. Zijlstra, A. Buursma, and O. W. van Assendelft, Visible and Near Infrared Absorption Spectra of Human and Animal Haemoglobin: Determination and Application , (VSP, Zeist, 2000).
  4. M. Cope, “The application of near infrared spectroscopy to non invasive monitoring of cerebral oxygenation in the newborn infant,” PhD Dissertation, (1991).
  5. S. A. Prahl, “Tabulated molar extinction coefficient for hemoglobin in water,” http://omlc.ogi.edu/spectra/hemoglobin/summary.html (1998).
  6. J. G. Kim, M. Xia, and H. Liu, “Extinction coefficients of hemoglobin for near-infrared spectroscopy of tissue,” IEEE Eng. Med. Biol. Mag. 24, 118–121 (2005). [CrossRef] [PubMed]
  7. J. G. Kim and H. Liu, “Variation of haemoglobin extinction coefficients can cause errors in the determination of haemoglobin concentration measured by near-infrared spectroscopy,” Phys. Med. Biol. 52, 6295–6332 (2007). [CrossRef] [PubMed]
  8. M. Meinke and M. Friebel, “Complex refractive index of hemoglobin in the wavelength range from 250 to 1100 nm,” Proc. SPIE 5862, 1–7 (2005).
  9. M. Meinke and M. Friebel, “Determination of the complex refractive index of highly concentrated hemoglobin solutions using transmittance and reflectance measurements,” J. Biomed. Opt. 10, 064019 (2005). [CrossRef]
  10. M. Meinke and M. Friebel, “Model function to calculate the refractive index of native hemoglobin in the wavelength range of 250 nm to 1100 nm dependent on concentration,” Appl. Opt. 45, 2838–2842 (2006). [CrossRef] [PubMed]
  11. K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equation in isotropic media,” IEEE Trans. Antennas Propag. 14, 302–307 (1966). [CrossRef]
  12. A. Taflove and S. C. Hagness, Computational Electrodynamics: the Finite-Difference Time-Domain Method , (Artech House, 2005).
  13. A. Taflove and M. E. Brodwin, “Numerical solution of steady-state electromagnetic scattering problems using the time-dependent Maxwell’s equations,” IEEE Trans. Microw. Theory Tech. 23, 623–630 (1975). [CrossRef]
  14. F. Zheng, Z. Chen, and J. Zhang, “Toward the development of a three-dimensional unconditionally stable finite-difference time-domain method,” IEEE Trans. Microw. Theory Tech. 48, 1550–1558 (2000). [CrossRef]
  15. T. Namiki, “3-D ADI-FDTD method-Unconditionally stable time domain algorithm for solving full vector Maxwell’s equations,” IEEE Trans. Microw. Theory Tech. 48, 1743–1748 (2000). [CrossRef]
  16. E. L. Tan, “Efficient algorithm for the unconditionally stable 3-D ADI-FDTD method,” IEEE Microw. Wireless Comp. Lett. 17, 7–9 (2007). [CrossRef]
  17. E. L. Tan, “Fundamental schemes for efficient unconditionally stable implicit finite-difference time-domain methods,” IEEE Trans. Antennas Propag. 56, 170–177 (2008). [CrossRef]
  18. G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 um wavelength region,” Appl. Opt. 12, 555–563 (1973). [CrossRef] [PubMed]
  19. B. Gustavsen and A. Semlyen, “Rational approximation of frequency domain responses by vector fitting,” IEEE Trans. Power Delivery 14, 1052–1061 (1999). [CrossRef]
  20. D. Y. Heh and E. L. Tan, “Corrected impulse invariance method in z-transform theory for frequency-dependent FDTD methods,” IEEE Trans. Antennas Propag. 57, 2683–2690 (2009). [CrossRef]
  21. E. L. Tan, “Concise current source implementation for efficient 3-D ADI-FDTD method,” IEEE Microw. Wireless Comp. Lett. 17, 748–750 (2007). [CrossRef]
  22. C. M. Furse and O. P. Gahdhi, “A memory efficient method of calculating specific absorption rate in CW FDTD simulations,” IEEE Trans. Biomed. Eng. 43, 558–560 (1996). [CrossRef] [PubMed]
  23. S. Paker and L. Sevgi, “FDTD evaluation of the SAR distribution in a human head near a mobile cellular phone,” Elektrik 6, 227–242 (1998).
  24. I. Laakso, “FDTD method in assessment of human exposure to base station radiation,” Masters Dissertation, (2007).
  25. S. G. Garcez, C. F. Galan, L. H. Bonani, and V. Baranauskas, “Estimating the electromagnetic field effects in biological tissues through the finite-difference time-domain method,” SBMO/IEEE MTT-S Int. Microw. and Optoelectronics Conf. Proc. , 43, 58–62 (2007).
  26. J. D. Jackson, Classical Electrodynamics , 3rd ed. (John Wiley & Sons, 1998).
  27. D. M. Pozar, Microwave Engineering , 3rd ed. (Wiley, 2005).
  28. W. C. Tay, D. Y. Heh, and E. L. Tan, “GPU-accelerated fundamental ADI-FDTD with complex frequency shifted convolutional perfectly matched layer,” PIER M 14, 177–192 (2010). [CrossRef]
  29. K. Boryczko, W. Dzewinel, and D. A. Yuen, “Modeling fibrin aggregation in blood flow with discrete particles,” Comput. Methods Prog. Biomed. 75, 181–194 (2004). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited