Comparison of a Semi-Analytic Variance Reduction Technique to Classical Monte Carlo Variance Reduction Techniques for High Aspect Ratio Pencil Beam Collimators for Emission Tomography Applications
Abstract
A semi-analytic variance reduction technique developed for collimated gamma emission tomography problems was compared to classic Monte Carlo variance reduction techniques within the Monte Carlo N Particle Transport (MCNP) code. In the semi-analytic technique, a computationally efficient, non-analog, monodirectional source biased Monte Carlo simulation is first performed. Analytical expressions or empirical values are then used to correct for solid angle and field-of-view effects introduced by the non-analog source definition. This variance reduction technique was compared with deterministic transport sphere (DXTRAN) and geometry splitting variance reduction schemes to determine the accuracy and computational savings of each technique relative to an analog pulse height tally (F8 tally) at 1, 3, 5, and 10 mm collimator aperture radii. For radii 0.5 mm and 0.3 mm a DXTRAN sphere was used in place of an analog F8 tally, due to large particle history demands, to analyze the accuracy of the semi-analytic variance reduction technique. The computational savings and accuracy were evaluated for six to seven photopeaks depending on the method used. The monodirectional source biasing technique overestimated the count rates by approximately 9%-19% when the radius is less than 3 mm, but the technique overestimated by a factor of 2 to 7, when the radius is greater than or equal to 3 mm. The monodirectionally source biased technique offered computational saving factors on the order of 108-1013 over 1-10 mm collimator radii studied. DXTRAN and geometry splitting methods yielded higher accuracy, but computational savings range from approximately 0.13 to 2.2 and 0.07 to 2.9, respectively indicating marginal improvement at best.
Recommended Citation
S. Kilby et al., "Comparison of a Semi-Analytic Variance Reduction Technique to Classical Monte Carlo Variance Reduction Techniques for High Aspect Ratio Pencil Beam Collimators for Emission Tomography Applications," Nuclear Instruments and Methods in Physics Research, Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 1001, Elsevier, Jun 2021.
The definitive version is available at https://doi.org/10.1016/j.nima.2021.165236
Department(s)
Nuclear Engineering and Radiation Science
Keywords and Phrases
Emission Tomography Transport; Monte Carlo N Particle Transport (MCNP); Variance Reduction
International Standard Serial Number (ISSN)
0168-9002
Document Type
Article - Journal
Document Version
Citation
File Type
text
Language(s)
English
Rights
© 2021 North-Holland, All rights reserved.
Publication Date
11 Jun 2021
Comments
U.S. Department of Energy, Grant NRC-HQ-13-G-38-0026