Related papers: Development of a GPU-based Monte Carlo dose calcul…
Particle tracking simulations with space charge effects are very important for high-intensity proton rings. Since they include not only Hamilton mechanics of a single particle but constructing charge densities and solving Poisson equations…
In computed tomography (CT) reconstruction, scattering causes server quality degradation of the reconstructed CT images by introducing streaks and cupping artifacts which reduce the detectability of low contrast objects. Monte Carlo (MC)…
The performance of the Hybrid Monte Carlo algorithm is determined by the speed of sparse matrix-vector multiplication within the context of preconditioned conjugate gradient iteration. We study these operations as implemented for the…
In a cyclotron-based proton therapy facility, the energy changes are performed by means of a degrader of variable thickness. The interaction of the proton beam with the degrader creates energy tails and increases the beam emittance. A…
We discuss the advantages of parallelization by multithreading on graphics processing units (GPUs) for parallel tempering Monte Carlo computer simulations of an exemplified bead-spring model for homopolymers. Since the sampling of a large…
The distribution of absorbed dose in radionuclide therapy with Lu$^{177}$ can be approximated by convolving an image of the time-integrated activity distribution with a dose voxel kernel representing different tissue types. This fast but…
A proton pencil beam is associated with a surrounding low-dose envelope, originating from nuclear interactions. It is important for treatment planning systems to accurately model this envelope when performing dose calculations for pencil…
Targeted Radionuclide Therapy (TRT) is a well-established technique for cancer treatment. In this approach, radionuclides are bound to specific drugs that selectively transport them to the tumor site. Within the ISOLPHARM project, a…
We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core…
We present the methodology of a photon-conserving, spatially-adaptive, ray-tracing radiative transfer algorithm, designed to run on multiple parallel Graphic Processing Units (GPUs). Each GPU has thousands computing cores, making them…
In radiotherapy, a trade-off exists between computational workload/speed and dose calculation accuracy. Calculation methods like pencil-beam convolution can be much faster than Monte-Carlo methods, but less accurate. The dose difference,…
To overcome the problem of inefficient computing time and unreliable results in MCNP5 calculation, a two-step method is adopted to calculate the energy deposition of prompt gamma-rays in detectors for depleted uranium spherical shells under…
We present our experience of porting the code used in the wave-packet convergent-close-coupling (WP-CCC) approach to run on NVIDIA V100 and AMD MI250X GPUs. The WP-CCC approach is a method used in the field of ion-atom collision physics to…
The goal of this manuscript is to estimate Type B uncertainties in absorbed-dose calculations arising from the different implementations in current state-of-the-art Monte Carlo codes of low-energy photon cross-sections (<200 keV). Monte…
The Monte Carlo (MC) simulation method is a powerful tool for radiation physicists, and several general-purpose software packages are commonly applied in a myriad of different radiation physics fields today. In medical physics, charged…
Purpose: Respiratory motion and the interplay effect cause the dose delivered to a patient with spot-scanned proton therapy to differ from the dose planned for during optimization. A new 4D robust optimization methodology was developed…
Simulating water droplets made up of millions of molecules and on timescales as needed in biological and technological applications is challenging due to the difficulty of balancing accuracy with computational capabilities. Most detailed…
We have previously described RapidBrachyMCTPS, a brachytherapy treatment planning toolkit consisting of a graphical user interface (GUI) and a Geant4-based Monte Carlo (MC) dose calculation engine. This work describes the tools that have…
This work presents an updated and extended guide on methods of a proper acceleration of the Monte Carlo integration of stochastic differential equations with the commonly available NVIDIA Graphics Processing Units using the CUDA programming…
We adapt the Quantum Monte Carlo method to the cascaded formalism of quantum optics, allowing us to simulate the emission of photons of known energy. Statistical processing of the photon clicks thus collected agrees with the theory of…