Azimuthal particle redistribution for the reduction of latent phase-space variance in Monte Carlo simulations

It is well known that the use of a phase space in Monte Carlo simulation introduces a baseline level of variance that cannot be suppressed through the use of standard particle recycling techniques. This variance (termed latent phase-space variance by Sempau et al) can be a significant limiting factor in achieving accurate, low-uncertainty dose scoring results, especially near the surface of a phantom. A BEAMnrc component module (MCTWIST) has been developed to reduce the presence of latent variance in phase-space-based Monte Carlo simulations by implementing azimuthal particle redistribution (APR). For each recycled use of a phase-space particle a random rotation about the beam's central axis is applied, effectively utilizing cylindrical symmetry of the particle fluence and therefore providing a more accurate representation of the source. The MCTWIST module is unique in that no physical component is actually added to the accelerator geometry. Beam modifications are made by directly transforming particle characteristics outside of BEAMnrc/EGSnrc particle transport. Using MCTWIST, we have demonstrated a reduction in latent phase-space variance by more than a factor of 20, for a 10 x 10 cm(2) field, when compared to standard phase-space particle recycling techniques. The reduction in latent variance has enabled the achievement of dramatically smoother in-water dose profiles. This paper outlines the use of MCTWIST in Monte Carlo simulation and quantifies for the first time the latent variance reduction resulting from exploiting cylindrical phase-space symmetry.