Hardware acceleration of a Monte Carlo simulation for photodynamic therapy [corrected] treatment planning.

Monte Carlo (MC) simulations are being used extensively in the field of medical biophysics, particularly for modeling light propagation in tissues. The high computation time for MC limits its use to solving only the forward solutions for a given source geometry, emission profile, and optical interaction coefficients of the tissue. However, applications such as photodynamic therapy treatment planning or image reconstruction in diffuse optical tomography require solving the inverse problem given a desired dose distribution or absorber distribution, respectively. A faster means for performing MC simulations would enable the use of MC-based models for accomplishing such tasks. To explore this possibility, a digital hardware implementation of a MC simulation based on the Monte Carlo for Multi-Layered media (MCML) software was implemented on a development platform with multiple field-programmable gate arrays (FPGAs). The hardware performed the MC simulation on average 80 times faster and was 45 times more energy efficient than the MCML software executed on a 3-GHz Intel Xeon processor. The resulting isofluence lines closely matched those produced by MCML in software, diverging by only less than 0.1 mm for fluence levels as low as 0.00001 cm(-2) in a skin model.

[1]  J. L. Reuss,et al.  The pulse in reflectance pulse oximetry: Modeling and experimental studies , 2004, Journal of Clinical Monitoring and Computing.

[2]  S. Harding.  Treatment of Age-Related Macular Degeneration with Photodynamic Therapy (TAP) Study Group. Verteporfin therapy for subfoveal choroidal neovascularization in age-related macular degeneration: three-year results of an open-label extension of 2 randomized clinical trials¿TAP Report No. 5 , 2002 .

[3]  Vaughn Betz,et al.  VPR: A new packing, placement and routing tool for FPGA research , 1997, FPL.

[4]  Maya Gokhale,et al.  Monte Carlo Radiative Heat Transfer Simulation on a Reconfigurable Computer , 2004, FPL.

[5]  T. Booth,et al.  MCNP variance reduction overview , 1985 .

[6]  S. A. Prahl,et al.  A Monte Carlo model of light propagation in tissue , 1989, Other Conferences.

[7]  M. Schweiger,et al.  A finite element approach for modeling photon transport in tissue. , 1993, Medical physics.

[8]  M. Blumenkranz,et al.  Verteporfin therapy for subfoveal choroidal neovascularization in age-related macular degeneration: three-year results of an open-label extension of 2 randomized clinical trials--TAP Report no. 5. , 2002 .

[9]  Lothar Lilge,et al.  GPU-accelerated Monte Carlo simulation for photodynamic therapy treatment planning , 2009, European Conference on Biomedical Optics.

[10]  Scott McMillan,et al.  A re-evaluation of the practicality of floating-point operations on FPGAs , 1998, Proceedings. IEEE Symposium on FPGAs for Custom Computing Machines (Cat. No.98TB100251).

[11]  J. Feyh,et al.  MRI-guided laser interstitial thermal therapy (LITT) of head and neck tumors: progress with a new method. , 1996, Journal of clinical laser medicine & surgery.

[12]  R. Ca.  Biological and physical strategies to improve the therapeutic index of photodynamic therapy. , 2008 .

[13]  Quan Liu,et al.  Experimental validation of Monte Carlo modeling of fluorescence in tissues in the UV-visible spectrum. , 2003, Journal of biomedical optics.

[14]  M. Biel.  Advances in photodynamic therapy for the treatment of head and neck cancers , 2006, Lasers in surgery and medicine.

[15]  Philip R. Moorby,et al.  The Verilog Hardware Description Language, 5th Edition , 2002 .

[16]  V. Fanti,et al.  Monte Carlo Computations for Radiotherapy with the use of Dedicated Processors , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[17]  B. Pogue,et al.  Photodynamic Therapy Dosimetry , 2003 .

[18]  L. C. Henyey,et al.  Diffuse radiation in the Galaxy , 1940 .

[19]  James L. Reuss,et al.  Multilayer modeling of reflectance pulse oximetry , 2005, IEEE Transactions on Biomedical Engineering.

[20]  T. Mang.  Dosimetric concepts for PDT. , 2008, Photodiagnosis and photodynamic therapy.

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

[22]  B. Pogue,et al.  Tutorial on diffuse light transport. , 2008, Journal of biomedical optics.

[23]  Lothar Lilge,et al.  FPGA-based Monte Carlo Computation of Light Absorption for Photodynamic Cancer Therapy , 2009, 2009 17th IEEE Symposium on Field Programmable Custom Computing Machines.

[24]  Jonathan Rose,et al.  CALL FOR ARTICLES IEEE Design & Test of Computers Special Issue on Microprocessors , 1996 .

[25]  M. Schweiger,et al.  Image reconstruction in optical tomography. , 1997, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[26]  L. Lilge,et al.  Towards conformal light delivery using tailored cylindrical diffusers: attainable light dose distributions , 2006, Physics in medicine and biology.

[27]  Vaughn Betz,et al.  Architecture and CAD for Deep-Submicron FPGAS , 1999, The Springer International Series in Engineering and Computer Science.

[28]  H. Kahn,et al.  Methods of Reducing Sample Size in Monte Carlo Computations , 1953, Oper. Res..

[29]  G. Lin,et al.  Skin Necrosis due to Photodynamic Action of Benzoporphyrin Depends on Circulating Rather than Tissue Drug Levels: Implications for Control of Photodynamic Therapy , 1998, Photochemistry and photobiology.

[30]  B. Wilson,et al.  The physics, biophysics and technology of photodynamic therapy , 2008, Physics in medicine and biology.

[31]  Thomas M. Keane,et al.  Distributed Monte Carlo simulation of light transportation in tissue , 2006, Proceedings 20th IEEE International Parallel & Distributed Processing Symposium.

[32]  Jonathan Rose,et al.  The Transmogrifier-4: an FPGA-based hardware development system with multi-gigabyte memory capacity and high host and memory bandwidth , 2005, Proceedings. 2005 IEEE International Conference on Field-Programmable Technology, 2005..

[33]  M. Schweiger,et al.  The finite element method for the propagation of light in scattering media: boundary and source conditions. , 1995, Medical physics.

[34]  M. Altschuler,et al.  Optimized interstitial PDT prostate treatment planning with the Cimmino feasibility algorithm. , 2005, Medical physics.

[35]  George A. Constantinides,et al.  A High Throughput FPGA-based Floating Point Conjugate Gradient Implementation , 2008, ARC.

[36]  R. Steiner,et al.  Spatially resolved absolute diffuse reflectance measurements for noninvasive determination of the optical scattering and absorption coefficients of biological tissue. , 1996, Applied optics.

[37]  Michael T. Wilson.  Lethal photosensitisation of oral bacteria and its potential application in the photodynamic therapy of oral infections , 2004, Photochemical & photobiological sciences : Official journal of the European Photochemistry Association and the European Society for Photobiology.

[38]  A. El Gamal,et al.  Architecture of field-programmable gate arrays , 1993, Proc. IEEE.

[39]  S. Stylli,et al.  Photodynamic therapy of cerebral glioma – A review Part II – Clinical studies , 2006, Journal of Clinical Neuroscience.

[40]  L. Grossweiner.  The Science of Phototherapy: An Introduction , 2005 .

[41]  A S Pasciak,et al.  High-speed evaluation of track-structure Monte Carlo electron transport simulations , 2008, Physics in medicine and biology.

[42]  Xiaobo Sharon Hu,et al.  An FPGA Solution for Radiation Dose Calculation , 2006, 2006 14th Annual IEEE Symposium on Field-Programmable Custom Computing Machines.

[43]  N. Ramanujam,et al.  Monte Carlo-based inverse model for calculating tissue optical properties. Part I: Theory and validation on synthetic phantoms. , 2006, Applied optics.

[44]  B. Wilson,et al.  Comparison of the In Vivo Photodynamic Threshold Dose for Photofrin, Mono‐ and Tetrasulfonated Aluminum Phthalocyanine Using a Rat Liver Model , 1998, Photochemistry and photobiology.

[45]  Russell Tessier,et al.  FPGA Architecture: Survey and Challenges , 2008, Found. Trends Electron. Des. Autom..

[46]  C. Ma,et al.  Clinical implementation of a Monte Carlo treatment planning system. , 1999, Medical physics.

[47]  Džena Hidović-Rowe,et al.  Modelling and validation of spectral reflectance for the colon , 2005, Physics in medicine and biology.

[48]  Q. Peng,et al.  5‐Aminolevulinic Acid‐Based Photodynamic Therapy: Principles and Experimental Research , 1997, Photochemistry and photobiology.

[49]  E. Heath,et al.  Dosimetric evaluation of the clinical implementation of the first commercial IMRT Monte Carlo treatment planning system at 6 MV. , 2004, Medical physics.

[50]  U. Nseyo,et al.  Photodynamic therapy (PDT) in the treatment of patients with resistant superficial bladder cancer: a long-term experience. , 1998, Journal of clinical laser medicine & surgery.

[51]  B. Funke,et al.  Transdermal photodynamic therapy—a treatment option for rheumatic destruction of small joints? , 2006, Lasers in surgery and medicine.

[52]  D. Boas,et al.  Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head. , 2002, Optics express.

[53]  Augusto Rendon,et al.  Treatment planning using tailored and standard cylindrical light diffusers for photodynamic therapy of the prostate. , 2008, Physics in medicine and biology.

[54]  Stanley B. Brown,et al.  The present and future role of photodynamic therapy in cancer treatment. , 2004, The Lancet. Oncology.

[55]  T. Dougherty,et al.  HOW DOES PHOTODYNAMIC THERAPY WORK? , 1992, Photochemistry and photobiology.

[56]  Zvonko G. Vranesic,et al.  Fundamentals of Digital Logic with Verilog Design , 1999 .

[57]  B. Overholt,et al.  Five-year efficacy and safety of photodynamic therapy with Photofrin in Barrett's high-grade dysplasia. , 2007, Gastrointestinal endoscopy.

[58]  B. Wilson,et al.  The propagation of optical radiation in tissue I. Models of radiation transport and their application , 1991, Lasers in Medical Science.

[59]  A. Evans,et al.  Vascular‐targeted photodynamic therapy (padoporfin, WST09) for recurrent prostate cancer after failure of external beam radiotherapy: a study of escalating light doses , 2008, BJU international.

[60]  Johannes Swartling,et al.  Realtime light dosimetry software tools for interstitial photodynamic therapy of the human prostate. , 2007, Medical physics.

[61]  I. Tan,et al.  The importance of In situ light dosimetry for photodynamic therapy of oral cavity tumors , 1999, Head & neck.

[62]  Donald E. Thomas,et al.  The Verilog® Hardware Description Language , 1990 .

[63]  M. Dereski,et al.  DEPTH MEASUREMENTS AND HISTOPATHOLOGICAL CHARACTERIZATION OF PHOTODYNAMIC THERAPY GENERATED NORMAL BRAIN NECROSIS AS A FUNCTION OF INCIDENT OPTICAL ENERGY DOSE , 1991, Photochemistry and photobiology.

[64]  Emil Wolf,et al.  Principles of Optics: Contents , 1999 .

[65]  J. Usuda,et al.  Photodynamic therapy (PDT) for lung cancers. , 2006, Journal of thoracic oncology : official publication of the International Association for the Study of Lung Cancer.

[66]  Giovanni De Micheli,et al.  Synthesis and Optimization of Digital Circuits , 1994 .

[67]  A. Johansson,et al.  Influence of treatment-induced changes in tussue absorption on treatment volume during interstitial photodynamic therapy , 2006 .

[68]  S L Jacques,et al.  CONV--convolution for responses to a finite diameter photon beam incident on multi-layered tissues. , 1997, Computer methods and programs in biomedicine.

[69]  Annamaria Kisslinger,et al.  Multiple processor version of a Monte Carlo code for photon transport in turbid media , 2000 .

[70]  M. Nichols,et al.  Design and testing of a white-light, steady-state diffuse reflectance spectrometer for determination of optical properties of highly scattering systems. , 1997, Applied optics.

[71]  Andrew D. Smith,et al.  SIMPROT: Using an empirically determined indel distribution in simulations of protein evolution , 2005, BMC Bioinformatics.

[72]  Pierre L'Ecuyer,et al.  Maximally equidistributed combined Tausworthe generators , 1996, Math. Comput..

[73]  J. F. Briesmeister.  MCNP-A General Monte Carlo N-Particle Transport Code , 1993 .

[74]  V. V. Tuchin.  Light scattering study of tissues , 1997 .

[75]  Sari Fien,et al.  Photodynamic therapy for non-melanoma skin cancer. , 2007, Journal of the National Comprehensive Cancer Network : JNCCN.

[76]  N. Metropolis,et al.  The Monte Carlo method. , 1949 .

[77]  F. Hetzel,et al.  Interstitial photodynamic therapy in the canine prostate. , 1997, British journal of urology.

[78]  T. Dougherty.  An update on photodynamic therapy applications. , 2002, Journal of clinical laser medicine & surgery.

[79]  L. Lilge,et al.  Implicit and explicit dosimetry in photodynamic therapy: a New paradigm , 1997, Lasers in Medical Science.

[80]  Herman H. Goldstine,et al.  The Electronic Numerical Integrator and Computer (ENIAC) , 1996, IEEE Ann. Hist. Comput..

[81]  D. Côté,et al.  Robust concentration determination of optically active molecules in turbid media with validated three-dimensional polarization sensitive Monte Carlo calculations. , 2005, Optics express.

[82]  Tomas Svensson,et al.  Parallel computing with graphics processing units for high-speed Monte Carlo simulation of photon migration. , 2008, Journal of biomedical optics.

[83]  B. Wilson,et al.  A Monte Carlo model for the absorption and flux distributions of light in tissue. , 1983, Medical physics.

[84]  Helen H Liu,et al.  Report of the AAPM Task Group No. 105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning. , 2007, Medical physics.

[85]  C. Gomer,et al.  Cellular targets and molecular responses associated with photodynamic therapy. , 1996, Journal of clinical laser medicine & surgery.