Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImpact of photon cross section systematic uncertainties on Monte Carlo-determined depth-dose distributions
87
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This work studies the impact of systematic uncertainties associated to interaction cross sections on depth dose curves determined by Monte Carlo simulations. The corresponding sensitivity factors are quantified by changing cross sections in a given amount and determining the variation in the dose. The influence of total cross sections for all particles, photons and only for Compton scattering is addressed. The PENELOPE code was used in all simulations. It was found that photon cross section sensitivity factors depend on depth. In addition, they are positive and negative for depths below and above an equilibrium depth, respectively. At this depth, sensitivity factors are null. The equilibrium depths found in this work agree very well with the mean free path of the corresponding incident photon energy. Using the sensitivity factors reported here, it is possible to estimate the impact of photon cross section uncertainties on the uncertainty of Monte Carlo-determined depth dose curves.
rate research
Read More
Monte Carlo (MC) method has been recognized the most accurate dose calculation method for radiotherapy. However, its extremely long computation time impedes clinical applications. Recently, a lot of efforts have been made to realize fast MC dose calculation on GPUs. Nonetheless, most of the GPU-based MC dose engines were developed in NVidia CUDA environment. This limits the code portability to other platforms, hindering the introduction of GPU-based MC simulations to clinical practice. The objective of this paper is to develop a fast cross-platform MC dose engine oclMC using OpenCL environment for external beam photon and electron radiotherapy in MeV energy range. Coupled photon-electron MC simulation was implemented with analogue simulations for photon transports and a Class II condensed history scheme for electron transports. To test the accuracy and efficiency of our dose engine oclMC, we compared dose calculation results of oclMC and gDPM, our previously developed GPU-based MC code, for a 15 MeV electron beam and a 6 MV photon beam on a homogenous water phantom, one slab phantom and one half-slab phantom. Satisfactory agreement was observed in all the cases. The average dose differences within 10% isodose line of the maximum dose were 0.48-0.53% for the electron beam cases and 0.15-0.17% for the photon beam cases. In terms of efficiency, our dose engine oclMC was 6-17% slower than gDPM when running both codes on the same NVidia TITAN card due to both different physics particle transport models and different computational environments between CUDA and OpenCL. The cross-platform portability was also validated by successfully running our new dose engine on a set of different compute devices including an Nvidia GPU card, two AMD GPU cards and an Intel CPU card using one or four cores. Computational efficiency among these platforms was compared.
We reconsider the discovery limit of multi-ton direct detection dark matter experiments in the light of recent measurements of the coherent elastic neutrino-nucleus scattering process. Assuming the cross section to be a parameter entirely determined by data, rather than using its Standard Model prediction, we use the COHERENT CsI and LAr data sets to determine WIMP discovery limits. Being based on a data-driven approach, the results are thus free from theoretical assumptions and fall within the WIMP mass regions where XENONnT and DARWIN have best expected sensitivities. We further determine the impact of subleading nuclear form factor and weak mixing angle uncertainties effects on WIMP discovery limits. We point out that these effects, albeit small, should be taken into account. Moreover, to quantify the impact of new physics effects in the neutrino background, we revisit WIMP discovery limits assuming light vector and scalar mediators as well as neutrino magnetic moments/transitions. We stress that the presence of new interactions in the neutrino sector, in general, tend to worsen the WIMP discovery limit.
We study the energy deposition by light and heavy nuclei in tissue-like media as used for cancer therapy. The depth-dose distributions for protons, $^{3}$He, $^{12}$C, $^{20}$Ne, and $^{58}$Ni nuclei are calculated within a Monte Carlo model based on the GEANT4 toolkit. These distributions are compared with each other and with available experimental data. It is demonstrated that nuclear fragmentation reactions essentially reduce the peak-to-plateau ratio of the dose profiles for deeply penetrating energetic ions heavier than $^{3}$He. On the other hand, all projectiles up to $^{20}$Ne were found equally suitable for therapeutic use at low penetration depths.
We recently built an analytical source model for GPU-based MC dose engine. In this paper, we present a sampling strategy to efficiently utilize this source model in GPU-based dose calculation. Our source model was based on a concept of phase-space-ring (PSR). This ring structure makes it effective to account for beam rotational symmetry, but not suitable for dose calculations due to rectangular jaw settings. Hence, we first convert PSR source model to its phase-space let (PSL) representation. Then in dose calculation, different types of sub-sources were separately sampled. Source sampling and particle transport were iterated. So that the particles being sampled and transported simultaneously are of same type and close in energy to alleviate GPU thread divergence. We also present an automatic commissioning approach to adjust the model for a good representation of a clinical linear accelerator . Weighting factors were introduced to adjust relative weights of PSRs, determined by solving a quadratic minimization problem with a non-negativity constraint. We tested the efficiency gain of our model over a previous source model using PSL files. The efficiency was improved by 1.70 ~ 4.41, due to the avoidance of long data reading and transferring. The commissioning problem can be solved in ~20 sec. Its efficacy was tested by comparing the doses computed using the commissioned model and the uncommissioned one, with measurements in different open fields in a water phantom under a clinical Varian Truebeam 6MV beam. For the depth dose curves, the average distance-to-agreement was improved from 0.04~0.28 cm to 0.04~0.12 cm for build-up region and the root-mean-square (RMS) dose difference after build-up region was reduced from 0.32%~0.67% to 0.21%~0.48%. For lateral dose profiles, RMS difference was reduced from 0.31%~2.0% to 0.06%~0.78% at inner beam and from 0.20%~1.25% to 0.10%~0.51% at outer beam.
Monte Carlo (MC) simulation is considered as the most accurate method for radiation dose calculations. Accuracy of a source model for a linear accelerator is critical for the overall dose calculation accuracy. In this paper, we presented an analytical source model that we recently developed for GPU-based MC dose calculations. A key concept called phase-space-ring (PSR) was proposed. It contained a group of particles that are of the same type and close in energy and radial distance to the center of the phase-space plane. The model parameterized probability densities of particle location, direction and energy for each primary photon PSR, scattered photon PSR and electron PSR. For a primary photon PSRs, the particle direction is assumed to be from the beam spot. A finite spot size is modeled with a 2D Gaussian distribution. For a scattered photon PSR, multiple Gaussian components were used to model the particle direction. The direction distribution of an electron PSRs was also modeled as a 2D Gaussian distribution with a large standard deviation. We also developed a method to analyze a phase-space file and derive corresponding model parameters. To test the accuracy of our linac source model, dose distributions of different open fields in a water phantom were calculated using our source model and compared to those directly calculated using the reference phase-space file. The average distance-to-agreement (DTA) was within 1 mm for the depth dose in the build-up region and beam penumbra regions. The root-mean-square (RMS) dose difference was within 1.1% for dose profiles at inner and outer beam regions. The maximal relative difference of output factors was within 0.5%. Good agreements were also found in an IMRT prostate patient case and an IMRT head-and-neck case. These results demonstrated the efficacy of our source model in terms of accurately representing a reference phase-space file.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
