Abstract:Purpose: Very fast Monte Carlo (MC) simulations of proton transport have been implemented recently on GPUs. However, these usually use simplified models for non-elastic (NE) proton-nucleus interactions. Our primary goal is to build a GPU-based proton transport MC with detailed modeling of elastic and NE collisions. Methods: Using CUDA, we implemented GPU kernels for these tasks: (1) Simulation of spots from our scanning nozzle configurations, (2) Proton propagation through CT geometry, considering nuclear elastic scattering, multiple scattering, and energy loss straggling, (3) Modeling of the intranuclear cascade stage of NE interactions, (4) Nuclear evaporation simulation, and (5) Statistical error estimates on the dose. To validate our MC, we performed: (1) Secondary particle yield calculations in NE collisions, (2) Dose calculations in homogeneous phantoms, (3) Re-calculations of head and neck plans from a commercial treatment planning system (TPS), and compared with Geant4.9.6p2/TOPAS. Results: Yields, energy and angular distributions of secondaries from NE collisions on various nuclei agree well with the Geant4 Bertini and Binary cascade models. The 3D-gamma pass rate at 2\%-2 mm for treatment plan simulations is typically 98\%. The net calculation time on a NVIDIA GTX680 card, including all data transfers, is $\sim$20 s for $1\times10^7$ proton histories. Conclusions: Our GPU-based MC is the first of its kind to include a detailed nuclear model to handle NE interactions of protons with any nucleus. Dosimetric calculations are in very good agreement with Geant4/TOPAS. Our MC is being used to perform fast routine clinical QA of pencil-beam based treatment plans, and has also been adopted as the dose engine in a clinically-applicable MC-based IMPT TPS. The detailed nuclear modeling will allow us to perform very fast linear energy transfer and neutron dose estimates on the GPU.

Monte Carlo Simulations of Down-Scattered Neutron and Knock-on Deuteron Spectra in Deuterium-Tritium Capsule Implosions

A Monte Carlo Code to Get Response Spectrum of Ions for Neutron Depth Profiling

Thermal-to-fusion Neutron Convertor and Monte Carlo Coupled Simulation of Deuteron/triton Transport and Secondary Products Generation

Use of MCNP to Simulate Nuclear Detection of Inertial Confinement Fusion

Research on Ionization and Secondary Reaction Coupled Monte Carlo Simulation and Its Application

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

Simulation of Neutron Multiplicity Counting Based on Monte Carlo Code RMC

Neutron spectra from photonuclear reactions: Performance testing of Monte-Carlo particle transport simulation codes

Development of a Coupled Depletion Perturbation Theory Methodology in Continuous-Energy Monte Carlo Depletion Simulations

SIMULATION OF PARTICLE INTERACTIONS IN A HIGH INTENSITY RADIO-FREQUENCY QUADRUPOLE FOR MOLECULAR HYDROGEN IONS

Prompt: Probability-Conserved Cross Section Biasing Monte Carlo Particle Transport System

Primary Neutron Spectra in Ion Vlasov-Fokker-Planck Simulations

TMCC: A transient three-dimentional neutron transport code by the direct simulation method

Monte Carlo simulations of muCF processes kinetics in deuterium gas

The TMDICE Monte Carlo shower program and algorithm for jet-fragmentation via coherent medium induced radiations and scatterings

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

neuSIM4: A comprehensive GEANT4 based neutron simulation code

Monte-Carlo Methods for the Neutron Transport Equation

Monte Carlo simulations of charged particle transport in plasmas: A fast single-scattering model

Secondary fusion coupled deuteron/triton transport simulation and thermal-to-fusion neutron convertor measurement

An effective initial particle sampling technique for Monte Carlo reactor transient simulations