Abstract:Objective.This work aims to develop a graphics processing unit (GPU)-accelerated Monte Carlo code for the coupled transport of photon, electron/positron and neutron over a broad range of energies for medical applications.Approach.By separating the MC evolution of radiation into source, transport, and interaction kernels, the branch divergence was alleviated. The memory coalescence was achieved by vectorizing the access pattern in which the secondary particles were archived. To accelerate further particle tracking, ray-tracing hardware acceleration in the Nvidia OptiXTMframework was applied. For photon and electron/positron, the EGSnrc interaction modules were ported as a GPU-optimized configuration. For neutron, a group-wised transport based on NJOY21 preprocessed data was implemented. The developed code was validated against CPU-based FLUKA. Neutron, x-ray and electron beams incident on water and ICRP phantoms were simulated. The neutron energy group and the transport parameters of photon and electron were set to be the same in both codes. A single Nvidia RTX 4090 card was used in this code while all 20 threads of a single Intel Core i9-10900K node were used in FLUKA.Main results.The number of histories was set to ensure that statistical uncertainties lower than 2% for all voxels whose doses were larger than 20% of the maximum. In all cases, the dose differences in the voxels between the codes were within 2.5%. For photons and electrons, the developed code was 150-300 times faster than FLUKA in both geometries. For neutrons, the code was respectively 80 and 135 times faster in the water and ICRP phantoms than FLUKA.Significance.This study offers an appropriate solution for uncoalesced memory access and branch divergence commonly encountered in coupled MC transport on the GPU architecture. The formidable acceleration in computing times and accuracy shown in this study can promise a routine clinical use of MC simulations.

Parallelization of MCATNP MONTE CARLO Particle Transport Code by Using MPI

Study on MPI/OpenMP hybrid parallelism for Monte Carlo neutron transport code

A general-purpose Monte Carlo particle transport code based on inverse transform sampling for radiotherapy dose calculation

Study on Parallel Computing of Radiation Detector Simulation Based on Geant4

Parallel optimization of the Crystal-KMC on Tianhe-2

Study of A Gpu-Based Parallel Computing Method for the Monte Carlo Program

Research on GPU Acceleration for Monte Carlo Criticality Calculation

A GPU-accelerated Monte Carlo code, RT2for coupled transport of photon, electron/positron, and neutron

Scalable and massively parallel Monte Carlo photon transport simulations for heterogeneous computing platforms

Parallelization of 3D Explicit Material Point Method Using OpenMP

Development of MC/DC: a performant, scalable, and portable Python-based Monte Carlo neutron transport code

GPU-Accelerated Monte Carlo Electron Transport Methods: Development and Application for Radiation Dose Calculations Using Six GPU cards

Development of the Integrated Parallelism Strategy for Large Scale Depletion Calculation in the Monte Carlo Code RMC

Parallelization of Kinetic Theory Simulations

A fast GPU-based Monte Carlo simulation of proton transport with detailed modeling of non-elastic interactions

Toward real-time Monte Carlo simulation using a commercial cloud computing infrastructure

Comparison of Two Accelerators for Monte Carlo Radiation Transport Calculations, NVIDIA Tesla M2090 GPU and Intel Xeon Phi 5110p Coprocessor: A Case Study for X-ray CT Imaging Dose Calculation

Performance Portable Monte Carlo Particle Transport on Intel, NVIDIA, and AMD GPUs

Performance Study of Monte Carlo Codes on Xeon Phi Coprocessors — Testing MCNP 6.1 and Profiling ARCHER Geometry Module on the FS7ONNi Problem

Evaluation of Single-Node Performance of Parallel Algorithms for Multigroup Monte Carlo Particle Transport Methods

An OpenCL-based Monte Carlo dose calculation engine (oclMC) for coupled photon-electron transport