Research Papers

Computational Simulation of Temperature Elevations in Tumors Using Monte Carlo Method and Comparison to Experimental Measurements in Laser Photothermal Therapy

[+] Author and Article Information
Liang Zhu

e-mail: zliang@umbc.edu
Department of Mechanical Engineering,
University of Maryland Baltimore County,
1000 Hilltop Circle,
Baltimore, MD 21250

1Corresponding author.

Contributed by the Bioengineering Division of ASME for publication in the Journal of Biomechanical Engineering. Manuscript received May 31, 2013; final manuscript received August 27, 2013; accepted manuscript posted September 12, 2013; published online October 10, 2013. Assoc. Editor: Ram Devireddy.

J Biomech Eng 135(12), 121007 (Oct 10, 2013) (11 pages) Paper No: BIO-13-1251; doi: 10.1115/1.4025388 History: Received May 31, 2013; Revised August 27, 2013; Accepted September 12, 2013

Accurate simulation of temperature distribution in tumors induced by gold nanorods during laser photothermal therapy relies on precise measurements of thermal, optical, and physiological properties of the tumor with or without nanorods present. In this study, a computational Monte Carlo simulation algorithm is developed to simulate photon propagation in a spherical tumor to calculate laser energy absorption in the tumor and examine the effects of the absorption (μa) and scattering (μs) coefficients of tumors on the generated heating pattern in the tumor. The laser-generated energy deposition distribution is then incorporated into a 3D finite-element model of prostatic tumors embedded in a mouse body to simulate temperature elevations during laser photothermal therapy using gold nanorods. The simulated temperature elevations are compared with measured temperatures in PC3 prostatic tumors in our previous in vivo experimental studies to extract the optical properties of PC3 tumors containing different concentrations of gold nanorods. It has been shown that the total laser energy deposited in the tumor is dominated by μa, while both μa and μs shift the distribution of the energy deposition in the tumor. Three sets of μa and μs are extracted, representing the corresponding optical properties of PC3 tumors containing different concentrations of nanorods to laser irradiance at 808 nm wavelength. With the injection of 0.1 cc of a 250 optical density (OD) nanorod solution, the total laser energy absorption rate is increased by 30% from the case of injecting 0.1 cc of a 50 OD nanorod solution, and by 125% from the control case without nanorod injection. Based on the simulated temperature elevations in the tumor, it is likely that after heating for 15 min, permanent thermal damage occurs in the tumor injected with the 250 OD nanorod solution, while thermal damage to the control tumor and the one injected with the 50 OD nanorod solution may be incomplete.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Bernardi, R. J., Lowery, A. R., Thompson, P. A., Blaney, S. M., and West, J. L., 2008, “Immunonanoshells for Targeted Photothermal Ablation in Medullo-Blastoma and Glioma: An In Vitro Evaluation Using Human Cell Lines,” J. Neuro-Oncol., 86(2), pp. 165–172. [CrossRef]
Gobin, A. M., Moon, J. J., and West, J. L., 2008, “Ephrin Al-Targeted Nano-Shells for Photothermal Ablation of Prostate Cancer Cells,” Int. J. Nanomed., 3(3), pp. 351–358.
Melancon, M. P., Lu, W., Yang, Z., Zhang, R., Cheng, Z., Elliot, A. M., Stafford, J., Olson, T., Zhang, J. Z., and Li, C., 2008, “In Vitro and In Vivo Targeting of Hollow Gold Nanoshells Directed at Epidermal Growth Factor Receptor for Photothermal Ablation Therapy,” Mol. Cancer Ther., 7(6), pp. 1730–1739. [CrossRef] [PubMed]
Loo, C. H., Hirsch, L., Lee, M. H., Chang, E., West, J., Halas, N., and Drezek, R., 2005, “Gold Nanoshell Bioconjugates for Molecular Imaging in Living Cells,” Opt. Lett., 30, pp. 1012–1014. [CrossRef] [PubMed]
Loo, C. L., Lowery, A., Halas, N., West, J., and Drezek, R., 2005, “Immunotargeted Nanoshells for Integrated Cancer Imaging and Therapy,” Nano Lett., 5, pp. 709–711. [CrossRef] [PubMed]
Lal, S., Clare, S. E., and Halas, N. J., 2008, “Nanoshell-Enabled Photothermal Cancer Therapy: Impending Clinical Impact,” Acc. Chem. Res., 41(12), pp. 1842–1851. [CrossRef] [PubMed]
Skrabalak, S. E., Chen, J., Lu, X., Li, X., and Xia, Y., 2007, “Gold Nanocages for Biomedical Applications,” Adv. Mater., 19, pp. 3177–3184. [CrossRef] [PubMed]
Reidenbach, H.-D., 2007, “Laser Safety,” Springer Handbook of Laser and Optics, F.Traeger, ed., Springer, New York.
Mohammed, Y., and Verhey, J. F., 2005, “A Finite Element Method Model to Simulate Laser Interstitial Thermo Therapy in Anatomical Inhomogenous Regions,” Biomed. Eng. Online, 4(2), pp. 10–45.
Welch, A. J., and van Gemert, M. J. C., 1995, “Optical-Thermal Response of Laser-Irradiated Tissue,” A. J.Welch and M. J. C.van Gemert, eds., Plenum Press, New York.
Engin, K., 1994, “Biological Rationale for Hyperthermia in Cancer Treatment (II),” Neoplasma, 41(5), pp. 277–283. [PubMed]
Dewhirst, M. V., Vujaskovic, Z., Jones, E., and Thrall, D., 2005, “Resetting the Biologic Rationale for Thermal Therapy,” Int. J. Hypotherm., 21(8), pp. 779–790. [CrossRef]
El-Sayed, I. H., Huang, X., and El-Sayed, M. A., 2006, “Selective Laser Photo Thermal Therapy of Epithelial Carcinoma Using Anti-EGFR Antibody Conjugated Gold Nanoparticles,” Cancer Lett., 239, pp. 129–135. [CrossRef] [PubMed]
O'Neal, D. P., Hirsch, L. R., Halas, N. J., Payne, J. D., and West, J. L., 2004, “Photo Thermal Tumor Ablation in Mice Using Near Infrared-Absorbing Nanoparticles,” Cancer Lett., 209, pp. 171–176. [CrossRef] [PubMed]
Xie, H., Gill-Sharp, K. L., and O'Neal, D. P., 2007, “Quantitative Estimation of Gold Nanoshell Concentrations in Whole Blood Using Dynamic Light Scattering,” Nanomed Nanotech. Biol. Med., 3, pp. 89–94. [CrossRef]
Weissleder, R., 2001, “A Clearer Vision for In Vivo Imaging,” Nat. Biotechnol., 19(4), pp. 316–317. [CrossRef] [PubMed]
Huang, X. W., Jain, P. K., El-Sayed, I. H., and El-Sayed, M. A., 2008, “Plasmonic Photothermal Therapy (PPTT) Using Gold Nanoparticles,” Lasers Med. Sci., 23, pp. 217–228. [CrossRef] [PubMed]
Ballou, B. L., Lagerholm, C., Ernst, L. A., Bruchez, M. P., and Waggoner, A. S., 2004, “Non-Invasive Imaging of Quantum Dots in Mice,” Bioconjug. Chem., 15, pp. 79–86. [CrossRef] [PubMed]
Choi, M. R., Stanton-Maxey, K. J., Stanley, J. K., Levin, C. S., Bardhan, R., Akin, D., Badve, S., Sturgis, J., Robinson, J. P., Bashir, R., Halas, N. J., and Clare, S. E., 2007, “A Cellular Trojan Horse for Delivery of Therapeutic Nanoparticles into Tumors,” Nano Lett., 7, pp. 3759–3765. [CrossRef] [PubMed]
Khlebtsov, B. N., Zharov, V., Melnikov, A., Tuchin, V., and Khlebtsov, N., 2006, “Optical Amplification of Photothermal Therapy With Gold Nanoparticles and Nanoclusters,” Nanotechnology, 17, pp. 5167–5179. [CrossRef]
Stern, J. M., Stanfield, J., Kabbani, W., Hsieh, J. T., and Cadeddu, J. A., 2008, “Selective Prostate Cancer Thermal Ablation With Laser Activated Gold Nano-Shells,” J. Urol., 179, pp. 748–753. [CrossRef] [PubMed]
Krag, D. N., Fuller, S. P., Oligino, L., Pero, S. C., Weaver, D. L., Soden, A. L., Hebert, C., Mills, S., Liu, C., and Peterson, D., 2002, “Phage-Displayed Random Peptide Libraries in Mice: Toxicity After Serial Panning,” Cancer Chemother. Pharmacol., 50, pp. 325–332. [CrossRef] [PubMed]
Hirsch, L. R., Stafford, R. J., Bankson, J. A., Sershen, S. R., Rivera, B., Price, R. E., Hazle, J. D., Halas, N. J., and West, J. L., 2003, “Nanoshell-Mediated Near-Infrared Thermal Therapy of Tumors Under Magnetic Resonance Guidance,” Proc. Natl. Acad. Sci. USA, 100(23), pp. 13549–13554. [CrossRef]
Qin, Z., and Bischof, J. C., 2010, “One-Dimensional Experimental Setup to Study the Heating of Nanoparticle Laden Systems,” Proceedings of the ASME Summer Bioengineering Engineering Conference, Naples, FL, June 16–19, 2010, ASME Paper No. SBC2010-19676.
Elliott, A. M., Stafford, R. J., Schwartz, J., Wang, J., Shetty, A. M., Bourgoyne, C., O'Neal, P., and Hazle, J. D., 2007, “Laser-Induced Thermal Response and Characterization of Nanoparticles for Cancer Treatment Using Magnetic Resonance Thermal Imaging,” Med. Phys., 34, pp. 3102–3108. [CrossRef] [PubMed]
Anvari, B., Rastegar, S., and Motamedi, M., 1994, “Modeling of Intraluminal Heating of Biological Tissue: Implications for Treatment of Benign Prostatic Hyperplasia,” IEEE Trans. Biomed. Eng., 41(9), pp. 854–864. [CrossRef] [PubMed]
Wang, L., Jacques, S. L., and Zheng, L., 1995, “MCML—Monte Carlo Modeling of Light Transport in Multi-Layered Tissues,” Comput. Meth. Prog. Biomed., 47(2), pp. 131–146. [CrossRef]
Wilson, B. C., and Adam, G., 1983, “A Monte Carlo Model for the Absorption and Flux Distributions of Light in Tissue,” Med. Phys., 10(6), pp. 824–830. [CrossRef] [PubMed]
Prahl, S. A., Keijzer, M., Jacques, S. L., and Welch, A. J., 1989, “A Monte Carlo Model of Light Propagation in Tissue,” SPIE Proc. Dosim. Laser Rad. Med. Biol., 5, pp. 102–111.
Flock, S. W., Wilson, B. C., and Patterson, M. S., 1989, “Monte Carlo Modeling of Light Propagation in Highly Scattering Tissues-II: Comparison With Measurements in Phantoms,” IEEE Trans. Biomed. Eng., 36(12), pp. 1169–1173. [CrossRef] [PubMed]
Keijzer, M., Pickering, J. W., and van Gemert, M. J. C., 1991, “Laser Beam Diameter for Port Wine Stain Treatment Lasers,” Laser. Surg. Med., 11, pp. 601–605. [CrossRef]
Keijzer, M., Jacques, S. L., Prahl, S. A., and Welch, A. J., 1989, “Light Distributions in Artery Tissue: Monte Carlo Simulations for Finite-Diameter Laser Beams,” Laser. Surg. Med., 9, pp. 148–154. [CrossRef]
Jacques, S., 1989, “Time-Resolved Propagation of Ultra-Short Laser Pulses Within Turbid Tissues,” Appl. Opt., 28, pp. 2223–2229. [CrossRef] [PubMed]
Jacques, S., 1989, “Time-Resolved Reflectance Spectroscopy in Turbid Tissues,” IEEE Trans. Biomed. Eng., 36, pp. 1155–1161. [CrossRef] [PubMed]
Wang, L., and Jacques, S. L., 1993, “Hybrid Model of Monte Carlo Simulation and Diffusion Theory for Light Reflectance by Turbid Media,” J. Opt. Soc., 10, pp. 1746–1752. [CrossRef]
Wang, L., and Jacques, S. L., 1994, “Optimized Radial and Angular Positions in Monte Carlo Modeling,” Med. Phys., 21, pp. 1081–1083. [CrossRef] [PubMed]
Wang, L., Jacques, S. L., and Zheng, L. Q., 1997, “CONV—Convolution for Responses to a Finite Diameter Photon Beam Incident on Multi-Layered Tissues,” Comput. Meth. Prog. Biomed., 54, pp. 141–150. [CrossRef]
Gardner, C. W., and Welch, A. J., 1994, “Monte Carlo Simulation of Light Transport in Tissue: Unscattered Absorption Events,” Appl. Opt., 33(13), pp. 2743–2745. [CrossRef] [PubMed]
Feng, Y. F., Fuentes, D., Hawkins, A., Bass, J., Rylander, M. N., Elliott, A., Shetty, A., Stafford, R. J., and Oden, J. T., 2010, “Nanoshell-Mediated Laser Surgery Simulation for Prostate Cancer Treatment,” Eng. Comp., 25(1), pp. 3–13. [CrossRef]
Gardner, C. M., Jacques, S. L., and WelchA. J., 1996, “Light Transport in Tissue: Accurate Expressions for One-Dimensional Fluence Rate and Escape Function Based Upon Monte Carlo Simulation,” Laser. Surg. Med., 18(2), pp. 129–138. [CrossRef]
Chen, Y., 2012, “Monte Carlo Simulation of Enhanced Laser Absorption in Tumors With Gold Nano-rod Inclusions: A Study of Laser Thermal Therapy Applied in Tumor Ablation,” M.S. thesis, University of Maryland Baltimore County, Baltimore, MD.
Mourant, J. R., Fuselier, T., Boyer, J., Johnson, T. M., and Bigio, I. J., 1997, “Predictions and Measurements of Scattering and Absorption Over Broad Wavelength Ranges in Tissue Phantoms,” Appl. Opt., 36(4), pp. 949–957. [CrossRef] [PubMed]
Piao, D. B., Bartels, K. E., Jiang, Z., Holyoak, G. R., Ritchey, J. W., Xu, G., Bunting, C. F., and Slobodov, G., 2010, “Alternative Transrectal Prostate Imaging: A Diffuse Optical Tomography Method,” IEEE J. Sel. Topics Quant. Electron., 16(4), pp. 715–729. [CrossRef]
Piao, D. B., Jiang, Z., Bartels, K. E., Holyoak, G. R., Ritchey, J. W., Xu, G., Bunting, C. F., and Slobodov, G., 2009, “In Vivo Trans-Rectal Ultrasound-Coupled Near-Infrared Optical Tomography of Intact Normal Canine Prostate,” J. Innovat. Opt. Health Sci., 2(3), pp. 215–225. [CrossRef]
Manuchehrabadi, N. A., Attaluri, A., Cai, H., Edziah, R., Lalanne, E., Bieberich, C., Ma, R., Johnson, A. M., and Zhu, L., 2012, “MicroCT Imaging and In Vivo Temperature Elevations in Implanted Prostatic Tumors in Laser Photothermal Therapy Using Gold Nanorods,” ASME J. Nanotech. Eng. Med., 3(2), p. 021003. [CrossRef]
Pennes, H. H., 1948, “Analysis of Tissue and Arterial Blood Temperature in the Resting Human Forearm,” J. Appl. Physiol., 1, pp. 93–122. [PubMed]
Fujita, S. T., Tamazawa, M., and Kuroda, K., 1998, “Effects of Blood Perfusion Rate on the Optimization of RF-Capacitive Hyperthermia,” IEEE Trans. Biomed. Eng., 45(9), pp. 1182–1186. [CrossRef] [PubMed]
Gore, J. P., and Xu, L. X., 2003, “Thermal Imaging for Biological and Medical Diagnostics,” Biomedical Photonics Handbook, T.Vo-Dinh, ed., CRC Press, Boca Raton, FL.
von Maltzahn, G., Park, J. H., Agrawal, A., and Bandaru, N. K., 2009, “Computationally Guided Photothermal Tumor Therapy Using Long-Circulating Gold Nanorod Antennas,” Cancer Res. Suppl. Data, 69(9), pp. 3892–3900. [CrossRef]
Manuchehrabadi, N., Attaluri, A., Cai, H., Edziah, R., Lalanne, E., Bieberich, C., Ma, R., Johnson, A. M., and Zhu, L., 2013, “Tumor Shrinkage Studies and Histological Analyses after Laser Photothermal Therapy Using Gold Nanorods,” International Journal of Biomedical Engineering and Technology, 12(2), pp. 157–174. [CrossRef]
Jain, P. K., Lee, K. S., El-Sayed, I. H., and El-SayedM. A., 2006, “Calculated Absorption and Scattering Properties of Gold Nanoparticles of Different Size, Shape, and Composition: Applications in Biological Imaging and Biomedicine,” J. Phys. Chem., 110, pp. 7238–7248.
Arnifield, M. R., Mathew, R. P., Tulip, J., and Mcphee, M. S., 1992, “Analysis of Tissue Optical Coefficients Using an Approximate Equation Valid for Comparable Absorption and Scattering,” Phys. Med. Biol., 37, pp. 1219–1230. [CrossRef] [PubMed]
Song, C. W., 1984, “Effect of Local Hyperthermia on Blood Flow and Microenvironment: A Review,” Cancer Res., 44, pp. 4721s–4730s. [PubMed]
Zhu, L., 2009, “Heat Transfer Applications in Biological Systems.,” Biomedical Engineering & Design Handbook, Volume 1: Bioengineering Fundamentals, M.Kutz, editor-in-chief, McGraw-Hill, New York, pp. 2.33–2.67.
Moritz, A. R., and Henriques, F. C.1947, “The Relative Importance of Time and Surface Temperature in the Causation of Cutaneous Burns,” Am. J. Pathol., 23, pp. 695–720. [PubMed]


Grahic Jump Location
Fig. 2

The left panel shows the subregions of the top surface of a spherical tumor irradiated by a laser spot of 7 mm in diameter. The right panel illustrates individual grid elements which are 0.2 × 0.2 × 0.2 mm3 in volume in the tumor region.

Grahic Jump Location
Fig. 1

Possible path trajectories for incident photons on the tumor top surface

Grahic Jump Location
Fig. 3

The generated mouse model and the two embedded tumors. The thermocouple paths are illustrated in the spherical tumor as solid and dash lines.

Grahic Jump Location
Fig. 8

Panel (a) gives the SAR contour map along the centerline in the mouse model, including the mouse body and the tumor. An enlarged SAR distribution in the tumor in the 250 OD injection case is shown in panel (b). The top view of the SAR distribution in the 3D structure of the mouse model is illustrated in panel (c).

Grahic Jump Location
Fig. 4

The effect of the absorption coefficient on the SAR distribution in the spherical tumor while the scattering coefficient μs is kept unchanged as 3 cm−1

Grahic Jump Location
Fig. 5

The effect of the scattering coefficient on the SAR distribution in the spherical tumor while the absorption coefficient μa is kept unchanged as 1.5 cm−1

Grahic Jump Location
Fig. 6

The obtained objective function values for various combinations of the absorption and scattering coefficients in the 250 OD injection case. The objective function is the smallest when μa = 1.1 cm−1 and μs = 7 cm−1.

Grahic Jump Location
Fig. 7

Experimental (symbols) and theoretical (lines) temperature distribution profiles in the tumors during laser photothermal therapy in the 250 OD injection case. Temperature mappings are along two tumor paths shown on the right bottom of the figure.

Grahic Jump Location
Fig. 9

The temperature contours inside the tumor based on the obtained optical properties of the tumor. (a) No nanorod injection, (b) the tumor is injected with 0.1 cc of the 50 OD nanorod solution, and (c) the tumor is injected with 0.1 cc of the 250 OD nanorod solution.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In