Radiative transport in the delta-P1 approximation: accuracy of fluence rate and optical penetration depth predictions in turbid semi-infinite media

Abstract

Using the delta-P(1) approximation to the Boltzmann transport equation we develop analytic solutions for the fluence rate produced by planar (1-D) and Gaussian beam (2-D) irradiation of a homogeneous, turbid, semi-infinite medium. To assess the performance of these solutions we compare the predictions for the fluence rate and two metrics of the optical penetration depth with Monte Carlo simulations. We provide results under both refractive-index matched and mismatched conditions for optical properties where the ratio of reduced scattering to absorption lies in the range 0< or =(micro(s')/micro(a))< or =10(4). For planar irradiation, the delta-P(1) approximation provides fluence rate profiles accurate to +/-16% for depths up to six transport mean free paths (l*) over the full range of optical properties. Metrics for optical penetration depth are predicted with an accuracy of +/-4%. For Gaussian irradiation using beam radii r(0) > or =3 l*, the accuracy of the fluence rate predictions is no worse than in the planar irradiation case. For smaller beam radii, the predictions degrade significantly. Specifically for media with (micro(s')/micro(a))=1 irradiated with a beam radius of r(0)=l*, the error in the fluence rate approaches 100%. Nevertheless, the accuracy of the optical penetration depth predictions remains excellent for Gaussian beam irradiation, and degrades to only +/-20% for r(0)=l*. These results show that for a given set of optical properties (micro(s')/micro(a)), the optical penetration depth decreases with a reduction in the beam diameter. Graphs are provided to indicate the optical and geometrical conditions under which one must replace the delta-P(1) results for planar irradiation with those for Gaussian beam irradiation to maintain accurate dosimetry predictions.