Abstract:We provide optimal particle split numbers for speeding up TOPAS Monte Carlo simulations of linear accelerator (linac) treatment heads while maintaining accuracy. In addition, we provide a new TOPAS physics module for simulating photoneutron production and transport. TOPAS simulation of a Siemens Oncor linac was used to determine the optimal number of splits for directional bremsstrahlung splitting as a function of the field size for 6MV and 18MV x-ray beams. The linac simulation was validated against published data of lateral dose profiles and percentage depth-dose curves (PDD) for the largest square field (40cm side). In separate simulations, neutron particle split and the custom TOPAS physics module was used to generate and transport photoneutrons, called "TsPhotoNeutron". Verification of accuracy was performed by comparing simulations with published measurements of: 1) neutron yields as a function of beam energy for thick targets of Al, Cu, Ta, W, Pb and concrete; and 2) photoneutron energy spectrum at 40cm laterally from the isocenter of the linac from an 18MV beam with closed jaws and MLC. The optimal number of splits obtained for directional bremsstrahlung splitting enhanced the computational efficiency by two orders of magnitude. The efficiency decreased with increasing beam energy and field size. Calculated lateral profiles in the central region agreed within 1mm/2% from measured data, PDD curves within 1 mm/1%. For the TOPAS physics module, at a split number of 146, the efficiency of computing photoneutron yields was enhanced by a 27.6 factor, whereas it improved the accuracy over existing Geant4 physics modules. This work provides simulation parameters and a new TOPAS physics module to improve the efficiency and accuracy of TOPAS simulations that involve photonuclear processes occurring in high-Z materials found in linac components, patient devices, and treatment rooms, as well as to explore new therapeutic modalities such as very-high-energy electron therapy.

Advancing Nuclear Physics Through TOPS Solvers and Tools

UNEDF: Advanced Scientific Computing Transforms the Low-Energy Nuclear Many-Body Problem

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

UNEDF: Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

UNEDF:Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

PETSc/TAO Developments for GPU-Based Early Exascale Systems

SUNDIALS Time Integrators for Exascale Applications with Many Independent ODE Systems

TOPAS simulation of photoneutrons in radiotherapy: accuracy and speed with variance reduction

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

A Two-Level Parallel Decomposition Approach for Transient Stability Constrained Optimal Power Flow

Scalable Nonlinear Solvers for Fully Implicit Coupled Nuclear Fuel Modeling. Final Report

Enhancement of DFT-calculations at Petascale: Nuclear Magnetic Resonance, Hybrid Density Functional Theory and Car-Parrinello calculations

A robust solver for large-scale heat transfer topology optimization

Compass, the Community Petascale Project for Accelerator Science and Simulation, A Broad Computational Accelerator Physics Initiative

High-Performance Computing for Nuclear Reactor Design and Safety Applications

ytopt: Autotuning Scientific Applications for Energy Efficiency at Large Scales

Eigenvalue Solvers for Modeling Nuclear Reactors on Leadership Class Machines

NSTX-U research advancing the physics of spherical tokamaks

AN INITIATIVE FOR THE DEVELOPMENT AND APPLICATION OF OPEN-SOURCE MULTI-PHYSICS SIMULATION IN SUPPORT OF R&D AND E&T IN NUCLEAR SCIENCE AND TECHNOLOGY

Towards Robust Solvers for Nuclear Fusion Simulations Using JOREK: A Numerical Analysis Perspective

Methods for Simulating TOF-PET in TOPAS Using a Low-Z Medium