Photon migration in turbid media using a cumulant approximation to radiative transfer

A photon transport model for light migration in turbid media based on a cumulant approximation to radiative transfer is presented for image reconstruction inside an infinite medium or a bounded medium with a planar geometry. This model treats weak inhomogeneities through a Born approximation of the Boltzmann radiative transfer equation and uses the second-order cumulant solution of photon density to the Boltzmann equation as the Green's function for the uniform background. It provides the correct behavior of photon migration at early times and reduces at long times to the center-moved diffusion approximation. At early times, it agrees much better with the result from the Monte Carlo simulation than the diffusion approximation. Both approximations agree well with the Monte Carlo simulation at later times. The weight function for image reconstruction under this proposed model is shown to have a strong dependence at both early and later times on absorption and/or scattering inhomogeneities located in the propagation direction of and close to the source, or in the field of view of and close to the detector. This effect originates from the initial ballistic motion of incident photons, which is substantially underestimated by the diffusion approximation.