Condensed history Monte Carlo methods for photon transport problems

We study methods for accelerating Monte Carlo simulations that retain most of the accuracy of conventional Monte Carlo algorithms. These methods - called Condensed History (CH) methods - have been very successfully used to model the transport of ionizing radiation in turbid systems. Our primary objective is to determine whether or not such methods might apply equally well to the transport of photons in biological tissue. In an attempt to unify the derivations, we invoke results obtained first by Lewis, Goudsmit and Saunderson and later improved by Larsen and Tolar. We outline how two of the most promising of the CH models - one based on satisfying certain similarity relations and the second making use of a scattering phase function that permits only discrete directional changes - can be developed using these approaches. The main idea is to exploit the connection between the space-angle moments of the radiance and the angular moments of the scattering phase function. We compare the results obtained when the two CH models studied are used to simulate an idealized tissue transport problem. The numerical results support our findings based on the theoretical derivations and suggest that CH models should play a useful role in modeling light-tissue interactions.

[1]  C. N. Kelber,et al.  Computing methods in reactor physics , 1970 .

[2]  José M. Fernández-Varea,et al.  On the theory and simulation of multiple elastic scattering of electrons , 1993 .

[3]  Edward W. Larsen A theoretical derivation of the Condensed History Algorithm , 1992 .

[4]  G. I. Bell,et al.  Nuclear Reactor Theory , 1952 .

[5]  J. L. Saunderson,et al.  Multiple Scattering of Electrons , 1940 .

[6]  B. Wilson,et al.  Similarity relations for anisotropic scattering in monte carlo simulations of deeply penetrating neu , 1989 .

[7]  Brian Claude Franke,et al.  Monte Carlo Electron Dose Calculations Using Discrete Scattering Angles and Discrete Energy Losses , 2005 .

[8]  E. Gelbard,et al.  Monte Carlo Principles and Neutron Transport Problems , 2008 .

[9]  Richard W. Hamming,et al.  Numerical Methods for Scientists and Engineers , 1962 .

[10]  H. W. Lewis Multiple Scattering in an Infinite Medium , 1950 .

[11]  W. Nelson,et al.  Monte Carlo Transport of Electrons and Photons , 1988 .

[12]  L. Zou,et al.  Several Theorems for the Trace of Self-conjugate Quaternion Matrix , 2008 .

[13]  T. Mackie,et al.  MMC--a high-performance Monte Carlo code for electron beam treatment planning. , 1995, Physics in medicine and biology.

[14]  D. P. Sloan,et al.  New multigroup Monte Carlo scattering algorithm suitable for neutral- and charged-particle Boltzmann and Fokker-Planck calculations , 1983 .

[15]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[16]  Edward W. Larsen,et al.  Generalized Fokker-Planck Approximations of Particle Transport with Highly Forward-Peaked Scattering , 2001 .

[17]  K. B. Oldham,et al.  An Atlas of Functions. , 1988 .

[18]  Iwan Kawrakow,et al.  On the condensed history technique for electron transport , 1998 .

[19]  Michael S. Patterson,et al.  A discrete method for anisotropic angular sampling in Monte Carlo simulations , 1988 .

[20]  G. C. Pomraning THE FOKKER-PLANCK OPERATOR AS AN ASYMPTOTIC LIMIT , 1992 .

[21]  A Transport Condensed History Algorithm for Electron Monte Carlo Simulations , 2001 .

[22]  I. Kawrakow Accurate condensed history Monte Carlo simulation of electron transport. I. EGSnrc, the new EGS4 version. , 2000, Medical physics.