We generalize a simple Monte Carlo (MC) model for dilute gases to consider the transport behavior of positrons and electrons in Percus-Yevick model liquids under highly nonequilibrium conditions, accounting rigorously for coherent scattering processes. The procedure extends an existing technique [Wojcik and Tachiya, Chem. Phys. Lett. 363, 381 (2002)], using the static structure factor to account for the altered anisotropy of coherent scattering in structured material. We identify the effects of the approximation used in the original method, and we develop a modified method that does not require that approximation. We also present an enhanced MC technique that has been designed to improve the accuracy and flexibility of simulations in spatially varying electric fields. All of the results are found to be in excellent agreement with an independent multiterm Boltzmann equation solution, providing benchmarks for future transport models in liquids and structured systems.
The kinetic theory of non-relativistic positrons in an idealized positron emission tomography PET environment is developed by solving the Boltzmann equation, allowing for coherent and incoherent elastic, inelastic, ionizing and annihilating collisions through positronium formation. An analytic expression is obtained for the positronium formation rate, as a function of distance from a spherical source, in terms of the solutions of the general kinetic eigenvalue problem. Numerical estimates of the positron range - a fundamental limitation on the accuracy of PET, are given for positrons in a model of liquid water, a surrogate for human tissue. Comparisons are made with the 'gas-phase' assumption used in current models in which coherent scattering is suppressed. Our results show that this assumption leads to an error of the order of a factor of approximately 2, emphasizing the need to accurately account for the structure of the medium in PET simulations.