Quasi-Monte Carlo method for calculating X-ray scatter in CT

Area of Science:

• Medical Physics
• Computational Physics
• High-Performance Computing

Background:

• Simulating X-ray interactions with matter is crucial for medical imaging and radiation therapy.
• Traditional Monte Carlo methods can be computationally intensive, limiting simulation speed and accuracy.
• Developing efficient algorithms is essential for advancing X-ray-based technologies.

Purpose of the Study:

• To develop and evaluate a novel, efficient algorithm for simulating X-ray photon transport through phantoms.
• To transform X-ray interaction trajectories into a high-dimensional integration problem.
• To provide an integral formula for photon detection probability.

Main Methods:

• Developed the gQMCFRD algorithm, integrating GPU-based quasi-Monte Carlo (gQMC) with forced random detection (FRD).
• gQMC utilizes deterministic low-discrepancy points for enhanced simulation convergence over traditional random sampling.
• Simulated X-ray photon transport and detection probability through a phantom material.

Main Results:

• The gQMCFRD algorithm demonstrated superior efficiency and convergence rates compared to GPU-based Monte Carlo tools (gMCDRR, gMMC, gMCFRD) and MC-GPU.
• Efficiency Improvement Factors ranged from 27 to 37 for gQMCFRD against MC-GPU.
• Results showed excellent agreement across all benchmarked methods, validating the simulation approach.

Conclusions:

• gQMCFRD offers a significant advancement in simulating X-ray photon transport, providing high accuracy and efficiency.
• The algorithm's effectiveness is demonstrated by its substantial efficiency gains over existing methods.
• This work contributes to more accurate and faster simulations in medical physics and related fields.