Proton Upset Monte Carlo Simulation
O'Neill, Patrick M.; Kouba, Coy K.; Foster, Charles C.
2009-01-01
The Proton Upset Monte Carlo Simulation (PROPSET) program calculates the frequency of on-orbit upsets in computer chips (for given orbits such as Low Earth Orbit, Lunar Orbit, and the like) from proton bombardment based on the results of heavy ion testing alone. The software simulates the bombardment of modern microelectronic components (computer chips) with high-energy (.200 MeV) protons. The nuclear interaction of the proton with the silicon of the chip is modeled and nuclear fragments from this interaction are tracked using Monte Carlo techniques to produce statistically accurate predictions.
Tseung, H. Wan Chan; J. Ma; Beltran, C.
2014-01-01
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, consid...
GPU-based fast Monte Carlo dose calculation for proton therapy
Jia, Xun; Schümann, Jan; Paganetti, Harald; Jiang, Steve B.
2012-12-01
Accurate radiation dose calculation is essential for successful proton radiotherapy. Monte Carlo (MC) simulation is considered to be the most accurate method. However, the long computation time limits it from routine clinical applications. Recently, graphics processing units (GPUs) have been widely used to accelerate computationally intensive tasks in radiotherapy. We have developed a fast MC dose calculation package, gPMC, for proton dose calculation on a GPU. In gPMC, proton transport is modeled by the class II condensed history simulation scheme with a continuous slowing down approximation. Ionization, elastic and inelastic proton nucleus interactions are considered. Energy straggling and multiple scattering are modeled. Secondary electrons are not transported and their energies are locally deposited. After an inelastic nuclear interaction event, a variety of products are generated using an empirical model. Among them, charged nuclear fragments are terminated with energy locally deposited. Secondary protons are stored in a stack and transported after finishing transport of the primary protons, while secondary neutral particles are neglected. gPMC is implemented on the GPU under the CUDA platform. We have validated gPMC using the TOPAS/Geant4 MC code as the gold standard. For various cases including homogeneous and inhomogeneous phantoms as well as a patient case, good agreements between gPMC and TOPAS/Geant4 are observed. The gamma passing rate for the 2%/2 mm criterion is over 98.7% in the region with dose greater than 10% maximum dose in all cases, excluding low-density air regions. With gPMC it takes only 6-22 s to simulate 10 million source protons to achieve ˜1% relative statistical uncertainty, depending on the phantoms and energy. This is an extremely high efficiency compared to the computational time of tens of CPU hours for TOPAS/Geant4. Our fast GPU-based code can thus facilitate the routine use of MC dose calculation in proton therapy.
A GPU-accelerated and Monte Carlo-based intensity modulated proton therapy optimization system
Purpose: Conventional spot scanning intensity modulated proton therapy (IMPT) treatment planning systems (TPSs) optimize proton spot weights based on analytical dose calculations. These analytical dose calculations have been shown to have severe limitations in heterogeneous materials. Monte Carlo (MC) methods do not have these limitations; however, MC-based systems have been of limited clinical use due to the large number of beam spots in IMPT and the extremely long calculation time of traditional MC techniques. In this work, the authors present a clinically applicable IMPT TPS that utilizes a very fast MC calculation. Methods: An in-house graphics processing unit (GPU)-based MC dose calculation engine was employed to generate the dose influence map for each proton spot. With the MC generated influence map, a modified least-squares optimization method was used to achieve the desired dose volume histograms (DVHs). The intrinsic CT image resolution was adopted for voxelization in simulation and optimization to preserve spatial resolution. The optimizations were computed on a multi-GPU framework to mitigate the memory limitation issues for the large dose influence maps that resulted from maintaining the intrinsic CT resolution. The effects of tail cutoff and starting condition were studied and minimized in this work. Results: For relatively large and complex three-field head and neck cases, i.e., >100 000 spots with a target volume of ∼1000 cm3 and multiple surrounding critical structures, the optimization together with the initial MC dose influence map calculation was done in a clinically viable time frame (less than 30 min) on a GPU cluster consisting of 24 Nvidia GeForce GTX Titan cards. The in-house MC TPS plans were comparable to a commercial TPS plans based on DVH comparisons. Conclusions: A MC-based treatment planning system was developed. The treatment planning can be performed in a clinically viable time frame on a hardware system costing around 45 000
A GPU-accelerated and Monte Carlo-based intensity modulated proton therapy optimization system
Ma, Jiasen, E-mail: ma.jiasen@mayo.edu; Beltran, Chris; Seum Wan Chan Tseung, Hok; Herman, Michael G. [Department of Radiation Oncology, Division of Medical Physics, Mayo Clinic, 200 First Street Southwest, Rochester, Minnesota 55905 (United States)
2014-12-15
Purpose: Conventional spot scanning intensity modulated proton therapy (IMPT) treatment planning systems (TPSs) optimize proton spot weights based on analytical dose calculations. These analytical dose calculations have been shown to have severe limitations in heterogeneous materials. Monte Carlo (MC) methods do not have these limitations; however, MC-based systems have been of limited clinical use due to the large number of beam spots in IMPT and the extremely long calculation time of traditional MC techniques. In this work, the authors present a clinically applicable IMPT TPS that utilizes a very fast MC calculation. Methods: An in-house graphics processing unit (GPU)-based MC dose calculation engine was employed to generate the dose influence map for each proton spot. With the MC generated influence map, a modified least-squares optimization method was used to achieve the desired dose volume histograms (DVHs). The intrinsic CT image resolution was adopted for voxelization in simulation and optimization to preserve spatial resolution. The optimizations were computed on a multi-GPU framework to mitigate the memory limitation issues for the large dose influence maps that resulted from maintaining the intrinsic CT resolution. The effects of tail cutoff and starting condition were studied and minimized in this work. Results: For relatively large and complex three-field head and neck cases, i.e., >100 000 spots with a target volume of ∼1000 cm{sup 3} and multiple surrounding critical structures, the optimization together with the initial MC dose influence map calculation was done in a clinically viable time frame (less than 30 min) on a GPU cluster consisting of 24 Nvidia GeForce GTX Titan cards. The in-house MC TPS plans were comparable to a commercial TPS plans based on DVH comparisons. Conclusions: A MC-based treatment planning system was developed. The treatment planning can be performed in a clinically viable time frame on a hardware system costing around 45
Dosimetric investigation of proton therapy on CT-based patient data using Monte Carlo simulation
Chongsan, T.; Liamsuwan, T.; Tangboonduangjit, P.
2016-03-01
The aim of radiotherapy is to deliver high radiation dose to the tumor with low radiation dose to healthy tissues. Protons have Bragg peaks that give high radiation dose to the tumor but low exit dose or dose tail. Therefore, proton therapy is promising for treating deep- seated tumors and tumors locating close to organs at risk. Moreover, the physical characteristic of protons is suitable for treating cancer in pediatric patients. This work developed a computational platform for calculating proton dose distribution using the Monte Carlo (MC) technique and patient's anatomical data. The studied case is a pediatric patient with a primary brain tumor. PHITS will be used for MC simulation. Therefore, patient-specific CT-DICOM files were converted to the PHITS input. A MATLAB optimization program was developed to create a beam delivery control file for this study. The optimization program requires the proton beam data. All these data were calculated in this work using analytical formulas and the calculation accuracy was tested, before the beam delivery control file is used for MC simulation. This study will be useful for researchers aiming to investigate proton dose distribution in patients but do not have access to proton therapy machines.
Tseung, H Wan Chan; Beltran, C
2014-01-01
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, en...
Parodi, K.; Ferrari, A.; Sommerer, F.; Paganetti, H.
2007-07-01
Clinical investigations on post-irradiation PET/CT (positron emission tomography/computed tomography) imaging for in vivo verification of treatment delivery and, in particular, beam range in proton therapy are underway at Massachusetts General Hospital (MGH). Within this project, we have developed a Monte Carlo framework for CT-based calculation of dose and irradiation-induced positron emitter distributions. Initial proton beam information is provided by a separate Geant4 Monte Carlo simulation modelling the treatment head. Particle transport in the patient is performed in the CT voxel geometry using the FLUKA Monte Carlo code. The implementation uses a discrete number of different tissue types with composition and mean density deduced from the CT scan. Scaling factors are introduced to account for the continuous Hounsfield unit dependence of the mass density and of the relative stopping power ratio to water used by the treatment planning system (XiO (Computerized Medical Systems Inc.)). Resulting Monte Carlo dose distributions are generally found in good correspondence with calculations of the treatment planning program, except a few cases (e.g. in the presence of air/tissue interfaces). Whereas dose is computed using standard FLUKA utilities, positron emitter distributions are calculated by internally combining proton fluence with experimental and evaluated cross-sections yielding 11C, 15O, 14O, 13N, 38K and 30P. Simulated positron emitter distributions yield PET images in good agreement with measurements. In this paper, we describe in detail the specific implementation of the FLUKA calculation framework, which may be easily adapted to handle arbitrary phase spaces of proton beams delivered by other facilities or include more reaction channels based on additional cross-section data. Further, we demonstrate the effects of different acquisition time regimes (e.g., PET imaging during or after irradiation) on the intensity and spatial distribution of the irradiation
Wuerl, Matthias
2016-08-01
Matthias Wuerl presents two essential steps to implement offline PET monitoring of proton dose delivery at a clinical facility, namely the setting up of an accurate Monte Carlo model of the clinical beamline and the experimental validation of positron emitter production cross-sections. In the first part, the field size dependence of the dose output is described for scanned proton beams. Both the Monte Carlo and an analytical computational beam model were able to accurately predict target dose, while the latter tends to overestimate dose in normal tissue. In the second part, the author presents PET measurements of different phantom materials, which were activated by the proton beam. The results indicate that for an irradiation with a high number of protons for the sake of good statistics, dead time losses of the PET scanner may become important and lead to an underestimation of positron-emitter production yields.
We implemented the simplified Monte Carlo (SMC) method on graphics processing unit (GPU) architecture under the computer-unified device architecture platform developed by NVIDIA. The GPU-based SMC was clinically applied for four patients with head and neck, lung, or prostate cancer. The results were compared to those obtained by a traditional CPU-based SMC with respect to the computation time and discrepancy. In the CPU- and GPU-based SMC calculations, the estimated mean statistical errors of the calculated doses in the planning target volume region were within 0.5% rms. The dose distributions calculated by the GPU- and CPU-based SMCs were similar, within statistical errors. The GPU-based SMC showed 12.30–16.00 times faster performance than the CPU-based SMC. The computation time per beam arrangement using the GPU-based SMC for the clinical cases ranged 9–67 s. The results demonstrate the successful application of the GPU-based SMC to a clinical proton treatment planning. (note)
Wan Chan Tseung, H; Ma, J; Beltran, C [Mayo Clinic, Rochester, MN (United States)
2014-06-15
Purpose: To build a GPU-based Monte Carlo (MC) simulation of proton transport with detailed modeling of elastic and non-elastic (NE) protonnucleus interactions, for use in a very fast and cost-effective proton therapy treatment plan verification system. Methods: Using the CUDA framework, we implemented kernels for the following tasks: (1) Simulation of beam spots from our possible scanning nozzle configurations, (2) Proton propagation through CT geometry, taking into account nuclear elastic and multiple scattering, as well as energy straggling, (3) Bertini-style modeling of the intranuclear cascade stage of NE interactions, and (4) Simulation of nuclear evaporation. To validate our MC, we performed: (1) Secondary particle yield calculations in NE collisions with therapeutically-relevant nuclei, (2) Pencil-beam dose calculations in homogeneous phantoms, (3) A large number of treatment plan dose recalculations, and compared with Geant4.9.6p2/TOPAS. A workflow was devised for calculating plans from a commercially available treatment planning system, with scripts for reading DICOM files and generating inputs for our MC. Results: Yields, energy and angular distributions of secondaries from NE collisions on various nuclei are in good agreement with the Geant4.9.6p2 Bertini and Binary cascade models. The 3D-gamma pass rate at 2%–2mm for 70–230 MeV pencil-beam dose distributions in water, soft tissue, bone and Ti phantoms is 100%. The pass rate at 2%–2mm for treatment plan calculations is typically above 98%. The net computational time on a NVIDIA GTX680 card, including all CPU-GPU data transfers, is around 20s for 1×10{sup 7} proton histories. Conclusion: Our GPU-based proton transport MC is the first of its kind to include a detailed nuclear model to handle NE interactions on any nucleus. Dosimetric calculations demonstrate very good agreement with Geant4.9.6p2/TOPAS. Our MC is being integrated into a framework to perform fast routine clinical QA of pencil
Purpose: To build a GPU-based Monte Carlo (MC) simulation of proton transport with detailed modeling of elastic and non-elastic (NE) protonnucleus interactions, for use in a very fast and cost-effective proton therapy treatment plan verification system. Methods: Using the CUDA framework, we implemented kernels for the following tasks: (1) Simulation of beam spots from our possible scanning nozzle configurations, (2) Proton propagation through CT geometry, taking into account nuclear elastic and multiple scattering, as well as energy straggling, (3) Bertini-style modeling of the intranuclear cascade stage of NE interactions, and (4) Simulation of nuclear evaporation. To validate our MC, we performed: (1) Secondary particle yield calculations in NE collisions with therapeutically-relevant nuclei, (2) Pencil-beam dose calculations in homogeneous phantoms, (3) A large number of treatment plan dose recalculations, and compared with Geant4.9.6p2/TOPAS. A workflow was devised for calculating plans from a commercially available treatment planning system, with scripts for reading DICOM files and generating inputs for our MC. Results: Yields, energy and angular distributions of secondaries from NE collisions on various nuclei are in good agreement with the Geant4.9.6p2 Bertini and Binary cascade models. The 3D-gamma pass rate at 2%–2mm for 70–230 MeV pencil-beam dose distributions in water, soft tissue, bone and Ti phantoms is 100%. The pass rate at 2%–2mm for treatment plan calculations is typically above 98%. The net computational time on a NVIDIA GTX680 card, including all CPU-GPU data transfers, is around 20s for 1×107 proton histories. Conclusion: Our GPU-based proton transport MC is the first of its kind to include a detailed nuclear model to handle NE interactions on any nucleus. Dosimetric calculations demonstrate very good agreement with Geant4.9.6p2/TOPAS. Our MC is being integrated into a framework to perform fast routine clinical QA of pencil
Kurosu, Keita [Department of Medical Physics and Engineering, Osaka University Graduate School of Medicine, Suita, Osaka 565-0871 (Japan); Department of Radiation Oncology, Osaka University Graduate School of Medicine, Suita, Osaka 565-0871 (Japan); Takashina, Masaaki; Koizumi, Masahiko [Department of Medical Physics and Engineering, Osaka University Graduate School of Medicine, Suita, Osaka 565-0871 (Japan); Das, Indra J. [Department of Radiation Oncology, Indiana University School of Medicine, Indianapolis, IN 46202 (United States); Moskvin, Vadim P., E-mail: vadim.p.moskvin@gmail.com [Department of Radiation Oncology, Indiana University School of Medicine, Indianapolis, IN 46202 (United States)
2014-10-01
Although three general-purpose Monte Carlo (MC) simulation tools: Geant4, FLUKA and PHITS have been used extensively, differences in calculation results have been reported. The major causes are the implementation of the physical model, preset value of the ionization potential or definition of the maximum step size. In order to achieve artifact free MC simulation, an optimized parameters list for each simulation system is required. Several authors have already proposed the optimized lists, but those studies were performed with a simple system such as only a water phantom. Since particle beams have a transport, interaction and electromagnetic processes during beam delivery, establishment of an optimized parameters-list for whole beam delivery system is therefore of major importance. The purpose of this study was to determine the optimized parameters list for GATE and PHITS using proton treatment nozzle computational model. The simulation was performed with the broad scanning proton beam. The influences of the customizing parameters on the percentage depth dose (PDD) profile and the proton range were investigated by comparison with the result of FLUKA, and then the optimal parameters were determined. The PDD profile and the proton range obtained from our optimized parameters list showed different characteristics from the results obtained with simple system. This led to the conclusion that the physical model, particle transport mechanics and different geometry-based descriptions need accurate customization in planning computational experiments for artifact-free MC simulation.
Although three general-purpose Monte Carlo (MC) simulation tools: Geant4, FLUKA and PHITS have been used extensively, differences in calculation results have been reported. The major causes are the implementation of the physical model, preset value of the ionization potential or definition of the maximum step size. In order to achieve artifact free MC simulation, an optimized parameters list for each simulation system is required. Several authors have already proposed the optimized lists, but those studies were performed with a simple system such as only a water phantom. Since particle beams have a transport, interaction and electromagnetic processes during beam delivery, establishment of an optimized parameters-list for whole beam delivery system is therefore of major importance. The purpose of this study was to determine the optimized parameters list for GATE and PHITS using proton treatment nozzle computational model. The simulation was performed with the broad scanning proton beam. The influences of the customizing parameters on the percentage depth dose (PDD) profile and the proton range were investigated by comparison with the result of FLUKA, and then the optimal parameters were determined. The PDD profile and the proton range obtained from our optimized parameters list showed different characteristics from the results obtained with simple system. This led to the conclusion that the physical model, particle transport mechanics and different geometry-based descriptions need accurate customization in planning computational experiments for artifact-free MC simulation
Purpose: For proton radiation therapy, Monte Carlo simulation (MCS) methods are recognized as the gold-standard dose calculation approach. Although previously unrealistic due to limitations in available computing power, GPU-based applications allow MCS of proton treatment fields to be performed in routine clinical use, on time scales comparable to that of conventional pencil-beam algorithms. This study focuses on validating the results of our GPU-based code (gPMC) versus fully implemented proton therapy based MCS code (TOPAS) for clinical patient cases. Methods: Two treatment sites were selected to provide clinical cases for this study: head-and-neck cases due to anatomical geometrical complexity (air cavities and density heterogeneities), making dose calculation very challenging, and prostate cases due to higher proton energies used and close proximity of the treatment target to sensitive organs at risk. Both gPMC and TOPAS methods were used to calculate 3-dimensional dose distributions for all patients in this study. Comparisons were performed based on target coverage indices (mean dose, V90 and D90) and gamma index distributions for 2% of the prescription dose and 2mm. Results: For seven out of eight studied cases, mean target dose, V90 and D90 differed less than 2% between TOPAS and gPMC dose distributions. Gamma index analysis for all prostate patients resulted in passing rate of more than 99% of voxels in the target. Four out of five head-neck-cases showed passing rate of gamma index for the target of more than 99%, the fifth having a gamma index passing rate of 93%. Conclusion: Our current work showed excellent agreement between our GPU-based MCS code and fully implemented proton therapy based MC code for a group of dosimetrically challenging patient cases
A Monte Carlo dose calculation algorithm for proton therapy
A Monte Carlo (MC) code (VMCpro) for treatment planning in proton beam therapy of cancer is introduced. It is based on ideas of the Voxel Monte Carlo algorithm for photons and electrons and is applicable to human tissue for clinical proton energies. In the present paper the implementation of electromagnetic and nuclear interactions is described. They are modeled by a Class II condensed history algorithm with continuous energy loss, ionization, multiple scattering, range straggling, δ-electron transport, nuclear elastic proton nucleus scattering and inelastic proton nucleus reactions. VMCpro is faster than the general purpose MC codes FLUKA by a factor of 13 and GEANT4 by a factor of 35 for simulations in a phantom with inhomogeneities. For dose calculations in patients the speed improvement is larger, because VMCpro has only a weak dependency on the heterogeneity of the calculation grid. Dose distributions produced with VMCpro are in agreement with GEANT4 results. Integrated or broad beam depth dose curves show maximum deviations not larger than 1% or 0.5 mm in regions with large dose gradients for the examples presented here
Currently, many proton therapy facilities are used for radiotherapy for treating cancer. The main advantage of proton therapy is the absence of exit dose, which offers a highly conformal dose to treatment target as well as better normal organ sparing. The most of treatment planning system (TPS) in proton therapy calculates dose distribution using a pencil beam algorithm (PBA). PBA is suitable for clinical proton therapy because of the fast computation time. However PBA shows accuracy limitations mainly because of the one-dimensional density scaling of proton pencil beams in water. Recently, we developed Monte Carlo simulation tools for the design of proton therapy facility at National Cancer Center (NCC) using GEANT4 toolkit (version GEANT4.9.2p02). Monte Carlo simulation is expected to reproduce precise influences of complex geometry and material varieties which are difficult to introduce to the PBA. The data format of Monte Carlo simulation result has different from DICOM-RT. Consequently we need we analysis software for comparing between TPS and Monte Carlo simulation. The main objective of this research is to develop an analysis toolkit for verifying precision and accuracy of the proton treatment planning system and to analyze dose calculating algorithm of the proton therapy using Monte Carlo simulation. In this work, we conclude that we developed an analysis software for GEANT4-based medical application. This toolkit is capable of evaluating the accuracy of calculated dose by TPS with Monte Carlo simulation.
Kim, Dae Hyun; Suh, Tae Suk [Dept. of Biomedical Engineering, College of Medicine, The Catholic University of Korea, Seoul (Korea, Republic of); Park, Sey Joon; Yoo, Seung Hoon; Lee, Se Byeong [Proton Therapy Center, National Cancer Center, Goyang (Korea, Republic of); Shin, Jung Wook [Dept. of Radiation Oncology, University of California, SanFrancisco (United States)
2011-11-15
Currently, many proton therapy facilities are used for radiotherapy for treating cancer. The main advantage of proton therapy is the absence of exit dose, which offers a highly conformal dose to treatment target as well as better normal organ sparing. The most of treatment planning system (TPS) in proton therapy calculates dose distribution using a pencil beam algorithm (PBA). PBA is suitable for clinical proton therapy because of the fast computation time. However PBA shows accuracy limitations mainly because of the one-dimensional density scaling of proton pencil beams in water. Recently, we developed Monte Carlo simulation tools for the design of proton therapy facility at National Cancer Center (NCC) using GEANT4 toolkit (version GEANT4.9.2p02). Monte Carlo simulation is expected to reproduce precise influences of complex geometry and material varieties which are difficult to introduce to the PBA. The data format of Monte Carlo simulation result has different from DICOM-RT. Consequently we need we analysis software for comparing between TPS and Monte Carlo simulation. The main objective of this research is to develop an analysis toolkit for verifying precision and accuracy of the proton treatment planning system and to analyze dose calculating algorithm of the proton therapy using Monte Carlo simulation. In this work, we conclude that we developed an analysis software for GEANT4-based medical application. This toolkit is capable of evaluating the accuracy of calculated dose by TPS with Monte Carlo simulation.
Biniam Yohannes Tesfamicael
2014-03-01
Full Text Available Purpose: To construct a dose monitoring system based on an endorectal balloon coupled to thin scintillating fibers to study the dose to the rectum in proton therapy of prostate cancer.Method: A Geant4 Monte Carlo toolkit was used to simulate the proton therapy of prostate cancer, with an endorectal balloon and a set of scintillating fibers for immobilization and dosimetry measurements, respectively.Results: A linear response of the fibers to the dose delivered was observed to within less than 2%. Results obtained show that fibers close to the prostate recorded higher dose, with the closest fiber recording about one-third of the dose to the target. A 1/r2 (r is defined as center-to-center distance between the prostate and the fibers decrease was observed as one goes toward the frontal and distal regions. A very low dose was recorded by the fibers beneath the balloon which is a clear indication that the overall volume of the rectal wall that is exposed to a higher dose is relatively minimized. Further analysis showed a relatively linear relationship between the dose to the target and the dose to the top fibers (total 17, with a slope of (-0.07 ± 0.07 at large number of events per degree of rotation of the modulator wheel (i.e., dose.Conclusion: Thin (1 mm × 1 mm, long (1 m scintillating fibers were found to be ideal for real time in-vivo dose measurement to the rectum during proton therapy of prostate cancer. The linear response of the fibers to the dose delivered makes them good candidates as dosimeters. With thorough calibration and the ability to define a good correlation between the dose to the target and the dose to the fibers, such dosimeters can be used for real time dose verification to the target.-----------------------------------Cite this article as: Tesfamicael BY, Avery S, Gueye P, Lyons D, Mahesh M. Scintillating fiber based in-vivo dose monitoring system to the rectum in proton therapy of prostate cancer: A Geant4 Monte Carlo
Benchmarking of proton transport in Super Monte Carlo simulation program
Full text of the publication follows. The Monte Carlo (MC) method has been traditionally applied in nuclear design and analysis due to its capability of dealing with complicated geometries and multi-dimensional physics problems as well as obtaining accurate results. The Super Monte Carlo Simulation Program (SuperMC) is developed by FDS Team in China for fusion, fission, and other nuclear applications. The simulations of radiation transport, isotope burn-up, material activation, radiation dose, and biology damage could be performed using SuperMC. Complicated geometries and the whole physical process of various types of particles in broad energy scale can be well handled. Bi-directional automatic conversion between general CAD models and full-formed input files of SuperMC is supported by MCAM, which is a CAD/image-based automatic modeling program for neutronics and radiation transport simulation. Mixed visualization of dynamical 3D dataset and geometry model is supported by RVIS, which is a nuclear radiation virtual simulation and assessment system. Continuous-energy cross section data from hybrid evaluated nuclear data library HENDL are utilized to support simulation. Neutronic fixed source and critical design parameters calculates for reactors of complex geometry and material distribution based on the transport of neutron and photon have been achieved in our former version of SuperMC. Recently, the proton transport has also been integrated in SuperMC in the energy region up to 10 GeV. The physical processes considered for proton transport include electromagnetic processes and hadronic processes. The electromagnetic processes include ionization, multiple scattering, Bremsstrahlung, and pair production processes. Public evaluated data from HENDL are used in some electromagnetic processes. In hadronic physics, the Bertini intra-nuclear cascade model with excitons, preequilibrium model, nucleus explosion model, fission model, and evaporation model are incorporated to
Tseung, H Wan Chan; Kreofsky, C R; Ma, D; Beltran, C
2016-01-01
Purpose: To demonstrate the feasibility of fast Monte Carlo (MC) based inverse biological planning for the treatment of head and neck tumors in spot-scanning proton therapy. Methods: Recently, a fast and accurate Graphics Processor Unit (GPU)-based MC simulation of proton transport was developed and used as the dose calculation engine in a GPU-accelerated IMPT optimizer. Besides dose, the dose-averaged linear energy transfer (LETd) can be simultaneously scored, which makes biological dose (BD) optimization possible. To convert from LETd to BD, a linear relation was assumed. Using this novel optimizer, inverse biological planning was applied to 4 patients: 2 small and 1 large thyroid tumor targets, and 1 glioma case. To create these plans, constraints were placed to maintain the physical dose (PD) within 1.25 times the prescription while maximizing target BD. For comparison, conventional IMRT and IMPT plans were created for each case in Eclipse (Varian, Inc). The same critical structure PD constraints were use...
Purpose: To develop a clinically applicable intensity modulated proton therapy (IMPT) optimization system that utilizes more accurate Monte Carlo (MC) dose calculation, rather than analytical dose calculation. Methods: A very fast in-house graphics processing unit (GPU) based MC dose calculation engine was employed to generate the dose influence map for each proton spot. With the MC generated influence map, a modified gradient based optimization method was used to achieve the desired dose volume histograms (DVH). The intrinsic CT image resolution was adopted for voxelization in simulation and optimization to preserve the spatial resolution. The optimizations were computed on a multi-GPU framework to mitigate the memory limitation issues for the large dose influence maps that Result from maintaining the intrinsic CT resolution and large number of proton spots. The dose effects were studied particularly in cases with heterogeneous materials in comparison with the commercial treatment planning system (TPS). Results: For a relatively large and complex three-field bi-lateral head and neck case (i.e. >100K spots with a target volume of ∼1000 cc and multiple surrounding critical structures), the optimization together with the initial MC dose influence map calculation can be done in a clinically viable time frame (i.e. less than 15 minutes) on a GPU cluster consisting of 24 Nvidia GeForce GTX Titan cards. The DVHs of the MC TPS plan compare favorably with those of a commercial treatment planning system. Conclusion: A GPU accelerated and MC-based IMPT optimization system was developed. The dose calculation and plan optimization can be performed in less than 15 minutes on a hardware system costing less than 45,000 dollars. The fast calculation and optimization makes the system easily expandable to robust and multi-criteria optimization. This work was funded in part by a grant from Varian Medical Systems, Inc
Peterson, S W; Polf, J; Archambault, L; Beddar, S [Department of Radiation Physics, Unit 94, University of Texas M D Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX 77030 (United States); Bues, M; Ciangaru, G; Smith, A [Proton Therapy Center, Unit 130, University of Texas M D Anderson Cancer Center, 1840 Old Spanish Trail, Houston, TX 77030 (United States)], E-mail: swpeters@mdanderson.org
2009-05-21
The purpose of this study is to validate the accuracy of a Monte Carlo calculation model of a proton magnetic beam scanning delivery nozzle developed using the Geant4 toolkit. The Monte Carlo model was used to produce depth dose and lateral profiles, which were compared to data measured in the clinical scanning treatment nozzle at several energies. Comparisons were also made between measured and simulated off-axis profiles to test the accuracy of the model's magnetic steering. Comparison of the 80% distal dose fall-off values for the measured and simulated depth dose profiles agreed to within 1 mm for the beam energies evaluated. Agreement of the full width at half maximum values for the measured and simulated lateral fluence profiles was within 1.3 mm for all energies. The position of measured and simulated spot positions for the magnetically steered beams agreed to within 0.7 mm of each other. Based on these results, we found that the Geant4 Monte Carlo model of the beam scanning nozzle has the ability to accurately predict depth dose profiles, lateral profiles perpendicular to the beam axis and magnetic steering of a proton beam during beam scanning proton therapy.
Proton therapy facilities are shielded to limit the amount of secondary radiation to which patients, occupational workers and members of the general public are exposed. The most commonly applied shielding design methods for proton therapy facilities comprise semi-empirical and analytical methods to estimate the neutron dose equivalent. This study compares the results of these methods with a detailed simulation of a proton therapy facility by using the Monte Carlo technique. A comparison of neutron dose equivalent values predicted by the various methods reveals the superior accuracy of the Monte Carlo predictions in locations where the calculations converge. However, the reliability of the overall shielding design increases if simulation results, for which solutions have not converged, e.g. owing to too few particle histories, can be excluded, and deterministic models are being used at these locations. Criteria to accept or reject Monte Carlo calculations in such complex structures are not well understood. An optimum rejection criterion would allow all converging solutions of Monte Carlo simulation to be taken into account, and reject all solutions with uncertainties larger than the design safety margins. In this study, the optimum rejection criterion of 10% was found. The mean ratio was 26, 62% of all receptor locations showed a ratio between 0.9 and 10, and 92% were between 1 and 100. (authors)
Reconstruction for proton computed tomography by tracing proton trajectories: A Monte Carlo study
Proton computed tomography (pCT) has been explored in the past decades because of its unique imaging characteristics, low radiation dose, and its possible use for treatment planning and on-line target localization in proton therapy. However, reconstruction of pCT images is challenging because the proton path within the object to be imaged is statistically affected by multiple Coulomb scattering. In this paper, we employ GEANT4-based Monte Carlo simulations of the two-dimensional pCT reconstruction of an elliptical phantom to investigate the possible use of the algebraic reconstruction technique (ART) with three different path-estimation methods for pCT reconstruction. The first method assumes a straight-line path (SLP) connecting the proton entry and exit positions, the second method adapts the most-likely path (MLP) theoretically determined for a uniform medium, and the third method employs a cubic spline path (CSP). The ART reconstructions showed progressive improvement of spatial resolution when going from the SLP [2 line pairs (lp) cm-1] to the curved CSP and MLP path estimates (5 lp cm-1). The MLP-based ART algorithm had the fastest convergence and smallest residual error of all three estimates. This work demonstrates the advantage of tracking curved proton paths in conjunction with the ART algorithm and curved path estimates
Proton therapy Monte Carlo SRNA-VOX code
Ilić Radovan D.
2012-01-01
The most powerful feature of the Monte Carlo method is the possibility of simulating all individual particle interactions in three dimensions and performing numerical experiments with a preset error. These facts were the motivation behind the development of a general-purpose Monte Carlo SRNA program for proton transport simulation in technical systems described by standard geometrical forms (plane, sphere, cone, cylinder, cube). Some of the possible applications of the SRNA program are:...
Monte Carlo calculations of positron emitter yields in proton radiotherapy
Positron emission tomography (PET) is a promising tool for monitoring the three-dimensional dose distribution in charged particle radiotherapy. PET imaging during or shortly after proton treatment is based on the detection of annihilation photons following the β+-decay of radionuclides resulting from nuclear reactions in the irradiated tissue. Therapy monitoring is achieved by comparing the measured spatial distribution of irradiation-induced β+-activity with the predicted distribution based on the treatment plan. The accuracy of the calculated distribution depends on the correctness of the computational models, implemented in the employed Monte Carlo (MC) codes that describe the interactions of the charged particle beam with matter and the production of radionuclides and secondary particles. However, no well-established theoretical models exist for predicting the nuclear interactions and so phenomenological models are typically used based on parameters derived from experimental data. Unfortunately, the experimental data presently available are insufficient to validate such phenomenological hadronic interaction models. Hence, a comparison among the models used by the different MC packages is desirable. In this work, starting from a common geometry, we compare the performances of MCNPX, GATE and PHITS MC codes in predicting the amount and spatial distribution of proton-induced activity, at therapeutic energies, to the already experimentally validated PET modelling based on the FLUKA MC code. In particular, we show how the amount of β+-emitters produced in tissue-like media depends on the physics model and cross-sectional data used to describe the proton nuclear interactions, thus calling for future experimental campaigns aiming at supporting improvements of MC modelling for clinical application of PET monitoring. (authors)
Monte Carlo comparison of x-ray and proton CT for range calculations of proton therapy beams
Proton computed tomography (CT) has been described as a solution for imaging the proton stopping power of patient tissues, therefore reducing the uncertainty of the conversion of x-ray CT images to relative stopping power (RSP) maps and its associated margins. This study aimed to investigate this assertion under the assumption of ideal detection systems. We have developed a Monte Carlo framework to assess proton CT performances for the main steps of a proton therapy treatment planning, i.e. proton or x-ray CT imaging, conversion to RSP maps based on the calibration of a tissue phantom, and proton dose simulations. Irradiations of a computational phantom with pencil beams were simulated on various anatomical sites and the proton range was assessed on the reference, the proton CT-based and the x-ray CT-based material maps. Errors on the tissue’s RSP reconstructed from proton CT were found to be significantly smaller and less dependent on the tissue distribution. The imaging dose was also found to be much more uniform and conformal to the primary beam. The mean absolute deviation for range calculations based on x-ray CT varies from 0.18 to 2.01 mm depending on the localization, while it is smaller than 0.1 mm for proton CT. Under the assumption of a perfect detection system, proton range predictions based on proton CT are therefore both more accurate and more uniform than those based on x-ray CT. (paper)
Proton therapy Monte Carlo SRNA-VOX code
Ilić Radovan D.
2012-01-01
Full Text Available The most powerful feature of the Monte Carlo method is the possibility of simulating all individual particle interactions in three dimensions and performing numerical experiments with a preset error. These facts were the motivation behind the development of a general-purpose Monte Carlo SRNA program for proton transport simulation in technical systems described by standard geometrical forms (plane, sphere, cone, cylinder, cube. Some of the possible applications of the SRNA program are: (a a general code for proton transport modeling, (b design of accelerator-driven systems, (c simulation of proton scattering and degrading shapes and composition, (d research on proton detectors; and (e radiation protection at accelerator installations. This wide range of possible applications of the program demands the development of various versions of SRNA-VOX codes for proton transport modeling in voxelized geometries and has, finally, resulted in the ISTAR package for the calculation of deposited energy distribution in patients on the basis of CT data in radiotherapy. All of the said codes are capable of using 3-D proton sources with an arbitrary energy spectrum in an interval of 100 keV to 250 MeV.
Kurosu, Keita; Das, Indra J.; Moskvin, Vadim P.
2016-01-01
Spot scanning, owing to its superior dose-shaping capability, provides unsurpassed dose conformity, in particular for complex targets. However, the robustness of the delivered dose distribution and prescription has to be verified. Monte Carlo (MC) simulation has the potential to generate significant advantages for high-precise particle therapy, especially for medium containing inhomogeneities. However, the inherent choice of computational parameters in MC simulation codes of GATE, PHITS and FLUKA that is observed for uniform scanning proton beam needs to be evaluated. This means that the relationship between the effect of input parameters and the calculation results should be carefully scrutinized. The objective of this study was, therefore, to determine the optimal parameters for the spot scanning proton beam for both GATE and PHITS codes by using data from FLUKA simulation as a reference. The proton beam scanning system of the Indiana University Health Proton Therapy Center was modeled in FLUKA, and the geometry was subsequently and identically transferred to GATE and PHITS. Although the beam transport is managed by spot scanning system, the spot location is always set at the center of a water phantom of 600 × 600 × 300 mm3, which is placed after the treatment nozzle. The percentage depth dose (PDD) is computed along the central axis using 0.5 × 0.5 × 0.5 mm3 voxels in the water phantom. The PDDs and the proton ranges obtained with several computational parameters are then compared to those of FLUKA, and optimal parameters are determined from the accuracy of the proton range, suppressed dose deviation, and computational time minimization. Our results indicate that the optimized parameters are different from those for uniform scanning, suggesting that the gold standard for setting computational parameters for any proton therapy application cannot be determined consistently since the impact of setting parameters depends on the proton irradiation technique. We
Microdosimetry of the full slowing down of protons using Monte Carlo track structure simulations.
Liamsuwan, T; Uehara, S; Nikjoo, H
2015-09-01
The article investigates two approaches in microdosimetric calculations based on Monte Carlo track structure (MCTS) simulations of a 160-MeV proton beam. In the first approach, microdosimetric parameters of the proton beam were obtained using the weighted sum of proton energy distributions and microdosimetric parameters of proton track segments (TSMs). In the second approach, phase spaces of energy depositions obtained using MCTS simulations in the full slowing down (FSD) mode were used for the microdosimetric calculations. Targets of interest were water cylinders of 2.3-100 nm in diameters and heights. Frequency-averaged lineal energies ([Formula: see text]) obtained using both approaches agreed within the statistical uncertainties. Discrepancies beyond this level were observed for dose-averaged lineal energies ([Formula: see text]) towards the Bragg peak region due to the small number of proton energies used in the TSM approach and different energy deposition patterns in the TSM and FSD of protons. PMID:25904698
Quantum Monte Carlo study of the protonated water dimer
Dagrada, Mario; Saitta, Antonino M; Sorella, Sandro; Mauri, Francesco
2013-01-01
We report an extensive theoretical study of the protonated water dimer (Zundel ion) by means of the highly correlated variational Monte Carlo and lattice regularized Monte Carlo approaches. This system represents the simplest model for proton transfer (PT) and a correct description of its properties is essential in order to understand the PT mechanism in more complex acqueous systems. Our Jastrow correlated AGP wave function ensures an accurate treatment of electron correlations. Exploiting the advantages of contracting the primitive basis set over atomic hybrid orbitals, we are able to limit dramatically the number of variational parameters with a systematic control on the numerical precision, crucial in order to simulate larger systems. We investigate energetics and geometrical properties of the Zundel ion as a function of the oxygen-oxygen distance, taken as reaction coordinate. In both cases, our QMC results are found in excellent agreement with coupled cluster CCSD(T) technique, the quantum chemistry "go...
Clinical implementation of full Monte Carlo dose calculation in proton beam therapy
Paganetti, Harald; Jiang, Hongyu; Parodi, Katia; Slopsema, Roelf; Engelsman, Martijn [Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114 (United States)
2008-09-07
The goal of this work was to facilitate the clinical use of Monte Carlo proton dose calculation to support routine treatment planning and delivery. The Monte Carlo code Geant4 was used to simulate the treatment head setup, including a time-dependent simulation of modulator wheels (for broad beam modulation) and magnetic field settings (for beam scanning). Any patient-field-specific setup can be modeled according to the treatment control system of the facility. The code was benchmarked against phantom measurements. Using a simulation of the ionization chamber reading in the treatment head allows the Monte Carlo dose to be specified in absolute units (Gy per ionization chamber reading). Next, the capability of reading CT data information was implemented into the Monte Carlo code to model patient anatomy. To allow time-efficient dose calculation, the standard Geant4 tracking algorithm was modified. Finally, a software link of the Monte Carlo dose engine to the patient database and the commercial planning system was established to allow data exchange, thus completing the implementation of the proton Monte Carlo dose calculation engine ('DoC++'). Monte Carlo re-calculated plans are a valuable tool to revisit decisions in the planning process. Identification of clinically significant differences between Monte Carlo and pencil-beam-based dose calculations may also drive improvements of current pencil-beam methods. As an example, four patients (29 fields in total) with tumors in the head and neck regions were analyzed. Differences between the pencil-beam algorithm and Monte Carlo were identified in particular near the end of range, both due to dose degradation and overall differences in range prediction due to bony anatomy in the beam path. Further, the Monte Carlo reports dose-to-tissue as compared to dose-to-water by the planning system. Our implementation is tailored to a specific Monte Carlo code and the treatment planning system XiO (Computerized Medical
Krueger, Rachel A. [Department of Chemistry, California Institute of Technology, Pasadena, California 91125 (United States); Haibach, Frederick G. [Confluent Science, Wilbraham, Massachusetts 01095 (United States); Fry, Dana L.; Gomez, Maria A., E-mail: magomez@mtholyoke.edu [Department of Chemistry, Mount Holyoke College, South Hadley, Massachusetts 01075 (United States)
2015-04-21
A centrality measure based on the time of first returns rather than the number of steps is developed and applied to finding proton traps and access points to proton highways in the doped perovskite oxides: AZr{sub 0.875}D{sub 0.125}O{sub 3}, where A is Ba or Sr and the dopant D is Y or Al. The high centrality region near the dopant is wider in the SrZrO{sub 3} systems than the BaZrO{sub 3} systems. In the aluminum-doped systems, a region of intermediate centrality (secondary region) is found in a plane away from the dopant. Kinetic Monte Carlo (kMC) trajectories show that this secondary region is an entry to fast conduction planes in the aluminum-doped systems in contrast to the highest centrality area near the dopant trap. The yttrium-doped systems do not show this secondary region because the fast conduction routes are in the same plane as the dopant and hence already in the high centrality trapped area. This centrality measure complements kMC by highlighting key areas in trajectories. The limiting activation barriers found via kMC are in very good agreement with experiments and related to the barriers to escape dopant traps.
Monte Carlo simulations of a novel Micromegas 2D array for proton dosimetry
Dolney, D.; Ainsley, C.; Hollebeek, R.; Maughan, R.
2016-02-01
Modern proton therapy affords control of the delivery of radiotherapeutic dose on fine length and temporal scales. The authors have developed a novel detector technology based on Micromesh Gaseous Structure (Micromegas) that is uniquely tailored for applications using therapeutic proton beams. An implementation of a prototype Micromegas detector for Monte Carlo using Geant4 is presented here. Comparison of simulation results with measurements demonstrates agreement in relative dose along the proton longitudinal dose profile to be 1%. The effect of a radioactive calibration source embedded in the chamber gas is demonstrated by measurements and reproduced by simulations, also at the 1% level. Our Monte Carlo simulations are shown to reproduce the time structure of ionization pulses produced by a double-scattering delivery system.
A Monte Carlo track structure code for low energy protons
Endo, S; Nikjoo, H; Uehara, S; Hoshi, M; Ishikawa, M; Shizuma, K
2002-01-01
A code is described for simulation of protons (100 eV to 10 MeV) track structure in water vapor. The code simulates molecular interaction by interaction for the transport of primary ions and secondary electrons in the form of ionizations and excitations. When a low velocity ion collides with the atoms or molecules of a target, the ion may also capture or lose electrons. The probabilities for these processes are described by the quantity cross-section. Although proton track simulation at energies above Bragg peak (>0.3 MeV) has been achieved to a high degree of precision, simulations at energies near or below the Bragg peak have only been attempted recently because of the lack of relevant cross-section data. As the hydrogen atom has a different ionization cross-section from that of a proton, charge exchange processes need to be considered in order to calculate stopping power for low energy protons. In this paper, we have used state-of-the-art Monte Carlo track simulation techniques, in conjunction with the pub...
Srna - Monte Carlo codes for proton transport simulation in combined and voxelized geometries
Ilić Radovan D.
2002-01-01
Full Text Available This paper describes new Monte Carlo codes for proton transport simulations in complex geometrical forms and in materials of different composition. The SRNA codes were developed for three dimensional (3D dose distribution calculation in proton therapy and dosimetry. The model of these codes is based on the theory of proton multiple scattering and a simple model of compound nucleus decay. The developed package consists of two codes: SRNA-2KG and SRNA-VOX. The first code simulates proton transport in combined geometry that can be described by planes and second order surfaces. The second one uses the voxelized geometry of material zones and is specifically adopted for the application of patient computer tomography data. Transition probabilities for both codes are given by the SRNADAT program. In this paper, we will present the models and algorithms of our programs, as well as the results of the numerical experiments we have carried out applying them, along with the results of proton transport simulation obtained through the PETRA and GEANT programs. The simulation of the proton beam characterization by means of the Multi-Layer Faraday Cup and spatial distribution of positron emitters obtained by our program indicate the imminent application of Monte Carlo techniques in clinical practice.
Srna-Monte Carlo codes for proton transport simulation in combined and voxelized geometries
This paper describes new Monte Carlo codes for proton transport simulations in complex geometrical forms and in materials of different composition. The SRNA codes were developed for three dimensional (3D) dose distribution calculation in proton therapy and dosimetry. The model of these codes is based on the theory of proton multiple scattering and a simple model of compound nucleus decay. The developed package consists of two codes: SRNA-2KG and SRNA-VOX. The first code simulates proton transport in combined geometry that can be described by planes and second order surfaces. The second one uses the voxelized geometry of material zones and is specifically adopted for the application of patient computer tomography data. Transition probabilities for both codes are given by the SRNADAT program. In this paper, we will present the models and algorithms of our programs, as well as the results of the numerical experiments we have carried out applying them, along with the results of proton transport simulation obtained through the PETRA and GEANT programs. The simulation of the proton beam characterization by means of the Multi-Layer Faraday Cup and spatial distribution of positron emitters obtained by our program indicate the imminent application of Monte Carlo techniques in clinical practice. (author)
TOPAS: An innovative proton Monte Carlo platform for research and clinical applications
Perl, J.; Shin, J.; Schuemann, J.; Faddegon, B.; Paganetti, H. [SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025 (United States); University of California San Francisco Comprehensive Cancer Center, 1600 Divisadero Street, San Francisco, California 94143-1708 (United States); Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114 (United States); University of California San Francisco Comprehensive Cancer Center, 1600 Divisadero Street, San Francisco, California 94143-1708 (United States); Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114 (United States)
2012-11-15
Purpose: While Monte Carlo particle transport has proven useful in many areas (treatment head design, dose calculation, shielding design, and imaging studies) and has been particularly important for proton therapy (due to the conformal dose distributions and a finite beam range in the patient), the available general purpose Monte Carlo codes in proton therapy have been overly complex for most clinical medical physicists. The learning process has large costs not only in time but also in reliability. To address this issue, we developed an innovative proton Monte Carlo platform and tested the tool in a variety of proton therapy applications. Methods: Our approach was to take one of the already-established general purpose Monte Carlo codes and wrap and extend it to create a specialized user-friendly tool for proton therapy. The resulting tool, TOol for PArticle Simulation (TOPAS), should make Monte Carlo simulation more readily available for research and clinical physicists. TOPAS can model a passive scattering or scanning beam treatment head, model a patient geometry based on computed tomography (CT) images, score dose, fluence, etc., save and restart a phase space, provides advanced graphics, and is fully four-dimensional (4D) to handle variations in beam delivery and patient geometry during treatment. A custom-designed TOPAS parameter control system was placed at the heart of the code to meet requirements for ease of use, reliability, and repeatability without sacrificing flexibility. Results: We built and tested the TOPAS code. We have shown that the TOPAS parameter system provides easy yet flexible control over all key simulation areas such as geometry setup, particle source setup, scoring setup, etc. Through design consistency, we have insured that user experience gained in configuring one component, scorer or filter applies equally well to configuring any other component, scorer or filter. We have incorporated key lessons from safety management, proactively
Nanodosimetric verification in proton therapy: Monte Carlo Codes Comparison
Full text: Nanodosimetry strives to develop a novel dosimetry concept suitable for advanced modalities of cancer radiotherapy, such as proton therapy. This project aims to evaluate the plausibility of the physical models implemented in the Geant4 Very Low Energy (Geant4-DNA) extensions by comparing nanodosimetric quantities calculated with Geant4-DNA and the PTB Monte Carlo track structure code. Nanodosimetric track structure parameters were calculated for cylindrical targets representing DNA and nucleosome segments and converted into the probability of producing a DSB using the model proposed by Garty et al. [1]. Monoenergetic protons and electrons of energies typical for 6-electron spectra were considered as primary particles. Good agreement was found between the two codes for electrons of energies above 200 eV. Below this energy Geant4-DNA produced slightly higher numbers of ionisations in the sensitive volumes and higher probabilities for DSB formation. For protons, Geant4-DNA also gave higher numbers of ionisations and DSB probabilities, particularly in the low energy range, while a satisfactory agreement was found for energies higher than I MeV. Comparing two codes can be useful as any observed divergence in results between the two codes provides valuable information as to where further consideration of the underlying physical models used in each code may be required. Consistently it was seen that the largest difference between the codes was in the low energy ranges for each particle type. (author)
Monte Carlo calculated stopping-power ratios, water/air, for clinical proton dosimetry (50-250 MeV)
Calculations of stopping power ratios, water to air, for the determination of absorbed dose to water in clinical proton beams using ionization chamber measurements have been undertaken using the Monte Carlo method. A computer code to simulate the transport of protons in water (PETRA) has been used to calculate Sw,air-data under different degrees of complexity, ranging from values based on primary protons only to data including secondary electrons and high-energy secondary protons produced in nonelastic nuclear collisions. All numerical data are based on ICRU 49 proton stopping powers. Calculations using primary protons have been compared to the simple continuous slowing-down approximation (c.s.d.a.) analytical technique used in proton dosimetry protocols, not finding significant differences that justify elaborate Monte Carlo simulations except beyond the mean range of the protons (the far side of the Bragg peak). The influence of nuclear nonelastic processes, through the detailed generation and transport of secondary protons, on the calculated stopping-power ratios has been found to be negligible. The effect of alpha particles has also been analysed, finding differences smaller than 0.1% from the results excluding them. Discrepancies of up to 0.6% in the plateau region have been found, however, when the production and transport of secondary electrons are taken into account. The large influence of nonelastic nuclear interactions on proton depth-dose distributions shows that the removal of primary protons from the incident beam decreases the peak-to-plateau ratio by a large factor, up to 40% at 250 MeV. It is therefore emphasized that nonelastic nuclear reactions should be included in Monte Carlo simulations of proton beam depth-dose distributions. (author)
Purpose: To construct a dose monitoring system based on an endorectal balloon coupled to thin scintillating fibers to study the dose delivered to the rectum during prostate cancer proton therapy Methods: The Geant4 Monte Carlo toolkit version 9.6p02 was used to simulate prostate cancer proton therapy treatments of an endorectal balloon (for immobilization of a 2.9 cm diameter prostate gland) and a set of 34 scintillating fibers symmetrically placed around the balloon and perpendicular to the proton beam direction (for dosimetry measurements) Results: A linear response of the fibers to the dose delivered was observed within <2%, a property that makes them good candidates for real time dosimetry. Results obtained show that the closest fiber recorded about 1/3 of the dose to the target with a 1/r2 decrease in the dose distribution as one goes toward the frontal and distal top fibers. Very low dose was recorded by the bottom fibers (about 45 times comparatively), which is a clear indication that the overall volume of the rectal wall that is exposed to a higher dose is relatively minimized. Further analysis indicated a simple scaling relationship between the dose to the prostate and the dose to the top fibers (a linear fit gave a slope of −0.07±0.07 MeV per treatment Gy) Conclusion: Thin (1 mm × 1 mm × 100 cm) long scintillating fibers were found to be ideal for real time in-vivo dose measurement to the rectum for prostate cancer proton therapy. The linear response of the fibers to the dose delivered makes them good candidates of dosimeters. With thorough calibration and the ability to define a good correlation between the dose to the target and the dose to the fibers, such dosimeters can be used for real time dose verification to the target
Monte Carlo Calculations Supporting Patient Plan Verification in Proton Therapy.
Lima, Thiago V M; Dosanjh, Manjit; Ferrari, Alfredo; Molineli, Silvia; Ciocca, Mario; Mairani, Andrea
2016-01-01
Patient's treatment plan verification covers substantial amount of the quality assurance (QA) resources; this is especially true for Intensity-Modulated Proton Therapy (IMPT). The use of Monte Carlo (MC) simulations in supporting QA has been widely discussed, and several methods have been proposed. In this paper, we studied an alternative approach from the one being currently applied clinically at Centro Nazionale di Adroterapia Oncologica (CNAO). We reanalyzed the previously published data (Molinelli et al. (1)), where 9 patient plans were investigated in which the warning QA threshold of 3% mean dose deviation was crossed. The possibility that these differences between measurement and calculated dose were related to dose modeling (Treatment Planning Systems (TPS) vs. MC), limitations on dose delivery system, or detectors mispositioning was originally explored, but other factors, such as the geometric description of the detectors, were not ruled out. For the purpose of this work, we compared ionization chambers' measurements with different MC simulation results. It was also studied that some physical effects were introduced by this new approach, for example, inter-detector interference and the delta ray thresholds. The simulations accounting for a detailed geometry typically are superior (statistical difference - p-value around 0.01) to most of the MC simulations used at CNAO (only inferior to the shift approach used). No real improvement was observed in reducing the current delta ray threshold used (100 keV), and no significant interference between ion chambers in the phantom were detected (p-value 0.81). In conclusion, it was observed that the detailed geometrical description improves the agreement between measurement and MC calculations in some cases. But in other cases, position uncertainty represents the dominant uncertainty. The inter-chamber disturbance was not detected for the therapeutic protons energies, and the results from the current delta threshold
Monte Carlo calculations supporting patient plan verification in proton therapy
Thiago Viana Miranda Lima
2016-03-01
Full Text Available Patient’s treatment plan verification covers substantial amount of the quality assurance (QA resources, this is especially true for Intensity Modulated Proton Therapy (IMPT. The use of Monte Carlo (MC simulations in supporting QA has been widely discussed and several methods have been proposed. In this paper we studied an alternative approach from the one being currently applied clinically at Centro Nazionale di Adroterapia Oncologica (CNAO. We reanalysed the previously published data (Molinelli et al. 2013, where 9 patient plans were investigated in which the warning QA threshold of 3% mean dose deviation was crossed. The possibility that these differences between measurement and calculated dose were related to dose modelling (Treatment Planning Systems (TPS vs MC, limitations on dose delivery system or detectors mispositioning was originally explored but other factors such as the geometric description of the detectors were not ruled out. For the purpose of this work we compared ionisation-chambers measurements with different MC simulations results. It was also studied some physical effects introduced by this new approach for example inter detector interference and the delta ray thresholds. The simulations accounting for a detailed geometry typically are superior (statistical difference - p-value around 0.01 to most of the MC simulations used at CNAO (only inferior to the shift approach used. No real improvement were observed in reducing the current delta-ray threshold used (100 keV and no significant interference between ion chambers in the phantom were detected (p-value 0.81. In conclusion, it was observed that the detailed geometrical description improves the agreement between measurement and MC calculations in some cases. But in other cases position uncertainty represents the dominant uncertainty. The inter chamber disturbance was not detected for the therapeutic protons energies and the results from the current delta threshold are
Monte Carlo Calculations Supporting Patient Plan Verification in Proton Therapy
Lima, Thiago V. M.; Dosanjh, Manjit; Ferrari, Alfredo; Molineli, Silvia; Ciocca, Mario; Mairani, Andrea
2016-01-01
Patient’s treatment plan verification covers substantial amount of the quality assurance (QA) resources; this is especially true for Intensity-Modulated Proton Therapy (IMPT). The use of Monte Carlo (MC) simulations in supporting QA has been widely discussed, and several methods have been proposed. In this paper, we studied an alternative approach from the one being currently applied clinically at Centro Nazionale di Adroterapia Oncologica (CNAO). We reanalyzed the previously published data (Molinelli et al. (1)), where 9 patient plans were investigated in which the warning QA threshold of 3% mean dose deviation was crossed. The possibility that these differences between measurement and calculated dose were related to dose modeling (Treatment Planning Systems (TPS) vs. MC), limitations on dose delivery system, or detectors mispositioning was originally explored, but other factors, such as the geometric description of the detectors, were not ruled out. For the purpose of this work, we compared ionization chambers’ measurements with different MC simulation results. It was also studied that some physical effects were introduced by this new approach, for example, inter-detector interference and the delta ray thresholds. The simulations accounting for a detailed geometry typically are superior (statistical difference – p-value around 0.01) to most of the MC simulations used at CNAO (only inferior to the shift approach used). No real improvement was observed in reducing the current delta ray threshold used (100 keV), and no significant interference between ion chambers in the phantom were detected (p-value 0.81). In conclusion, it was observed that the detailed geometrical description improves the agreement between measurement and MC calculations in some cases. But in other cases, position uncertainty represents the dominant uncertainty. The inter-chamber disturbance was not detected for the therapeutic protons energies, and the results from the current delta
SU-E-T-239: Monte Carlo Modelling of SMC Proton Nozzles Using TOPAS
Purpose: To expedite and cross-check the commissioning of the proton therapy nozzles at Samsung Medical Center using TOPAS. Methods: We have two different types of nozzles at Samsung Medical Center (SMC), a multi-purpose nozzle and a pencil beam scanning dedicated nozzle. Both nozzles have been modelled in Monte Carlo simulation by using TOPAS based on the vendor-provided geometry. The multi-purpose nozzle is mainly composed of wobbling magnets, scatterers, ridge filters and multi-leaf collimators (MLC). Including patient specific apertures and compensators, all the parts of the nozzle have been implemented in TOPAS following the geometry information from the vendor.The dedicated scanning nozzle has a simpler structure than the multi-purpose nozzle with a vacuum pipe at the down stream of the nozzle.A simple water tank volume has been implemented to measure the dosimetric characteristics of proton beams from the nozzles. Results: We have simulated the two proton beam nozzles at SMC. Two different ridge filters have been tested for the spread-out Bragg peak (SOBP) generation of wobbling mode in the multi-purpose nozzle. The spot sizes and lateral penumbra in two nozzles have been simulated and analyzed using a double Gaussian model. Using parallel geometry, both the depth dose curve and dose profile have been measured simultaneously. Conclusion: The proton therapy nozzles at SMC have been successfully modelled in Monte Carlo simulation using TOPAS. We will perform a validation with measured base data and then use the MC simulation to interpolate/extrapolate the measured data. We believe it will expedite the commissioning process of the proton therapy nozzles at SMC
Srna-Monte Carlo codes for proton transport simulation in combined and voxelized geometries
Ilic, R D; Stankovic, S J
2002-01-01
This paper describes new Monte Carlo codes for proton transport simulations in complex geometrical forms and in materials of different composition. The SRNA codes were developed for three dimensional (3D) dose distribution calculation in proton therapy and dosimetry. The model of these codes is based on the theory of proton multiple scattering and a simple model of compound nucleus decay. The developed package consists of two codes: SRNA-2KG and SRNA-VOX. The first code simulates proton transport in combined geometry that can be described by planes and second order surfaces. The second one uses the voxelized geometry of material zones and is specifically adopted for the application of patient computer tomography data. Transition probabilities for both codes are given by the SRNADAT program. In this paper, we will present the models and algorithms of our programs, as well as the results of the numerical experiments we have carried out applying them, along with the results of proton transport simulation obtaine...
Monte Carlo simulation in proton computed tomography: a study of image reconstruction technique
Inocente, Guilherme Franco; Stenico, Gabriela V.; Hormaza, Joel Mesa [Universidade Estadual Paulista Julio de Mesquita Filho (UNESP), Botucatu, SP (Brazil). Inst. de Biociencias. Dept. de Fisica e Biofisica
2012-07-01
Full text: The radiation method is one of the most used for cancer treatment. In this context arises therapy with proton beams in front of conventional radiotherapy. It is known that with proton therapy there are more advantages to the patient treated when compared with more conventional methods. The dose distributed along the path, especially in healthy tissues - neighbor the tumor, is smaller and the accuracy of treatment is much better. To carry out the treatment, the patient undergoes a plan through images for visualization and location of the target volume. The main method for obtaining these images is computed tomography X-ray (XCT). For treatment with proton beam this imaging technique can to generate some uncertainties. The purpose of this project is to study the feasibility of reconstructing images generated from the irradiation with proton beams, thereby reducing some inaccuracies, as it will be the same type of radiation as treatment planning, and also to drastically reduce some errors location, since the planning can be done at the same place and just before where the patient is treated. This study aims to obtain a relationship between the intrinsic property of the interaction of photons and protons with matter. For this we use computational simulation based on Monte Carlo method with the code SRIM 2008 and MCNPX v.2.5.0, to reconstruct images using the technique used in conventional computed tomography. (author)
Monte Carlo and Analytical Calculation of Lateral Deflection of Proton Beams in Homogeneous Targets
Proton radiation therapy is a precise form of radiation therapy, but the avoidance of damage to critical normal tissues and the prevention of geographical tumor misses require accurate knowledge of the dose delivered to the patient and the verification of his position demand a precise imaging technique. In proton therapy facilities, the X-ray Computed Tomography (xCT) is the preferred technique for the planning treatment of patients. This situation has been changing nowadays with the development of proton accelerators for health care and the increase in the number of treated patients. In fact, protons could be more efficient than xCT for this task. One essential difficulty in pCT image reconstruction systems came from the scattering of the protons inside the target due to the numerous small-angle deflections by nuclear Coulomb fields. The purpose of this study is the comparison of an analytical formulation for the determination of beam lateral deflection, based on Moliere's theory and Rutherford scattering with Monte Carlo calculations by SRIM 2008 and MCNPX codes.
Comparison of linear energy transfer scoring techniques in Monte Carlo simulations of proton beams
Monte Carlo (MC) simulations are commonly used to study linear energy transfer (LET) distributions in therapeutic proton beams. Various techniques have been used to score LET in MC simulations. The goal of this work was to compare LET distributions obtained using different LET scoring techniques and examine the sensitivity of these distributions to changes in commonly adjusted simulation parameters. We used three different techniques to score average proton LET in TOPAS, which is a MC platform based on the Geant4 simulation toolkit. We determined the sensitivity of each scoring technique to variations in the range production thresholds for secondary electrons and protons. We also compared the depth-LET distributions that we acquired using each technique in a simple monoenergetic proton beam and in a more clinically relevant modulated proton therapy beam. Distributions of both fluence-averaged LET (LETΦ) and dose-averaged LET (LETD) were studied. We found that LETD values varied more between different scoring techniques than the LETΦ values did, and different LET scoring techniques showed different sensitivities to changes in simulation parameters. (note)
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data
Ilic, Radovan D [Laboratory of Physics (010), Vinca Institute of Nuclear Sciences, PO Box 522, 11001 Belgrade (Serbia and Montenegro); Spasic-Jokic, Vesna [Laboratory of Physics (010), Vinca Institute of Nuclear Sciences, PO Box 522, 11001 Belgrade (Serbia and Montenegro); Belicev, Petar [Laboratory of Physics (010), Vinca Institute of Nuclear Sciences, PO Box 522, 11001 Belgrade (Serbia and Montenegro); Dragovic, Milos [Center for Nuclear Medicine MEDICA NUCLEARE, Bulevar Despota Stefana 69, 11000 Belgrade (Serbia and Montenegro)
2005-03-07
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour.
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data
Ilic, Radovan D.; Spasic-Jokic, Vesna; Belicev, Petar; Dragovic, Milos
2005-03-01
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour.
The Monte Carlo SRNA-VOX code for 3D proton dose distribution in voxelized geometry using CT data
This paper describes the application of the SRNA Monte Carlo package for proton transport simulations in complex geometry and different material compositions. The SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The decay of proton induced compound nuclei was simulated by the Russian MSDM model and our own using ICRU 63 data. The developed package consists of two codes: the SRNA-2KG, which simulates proton transport in combinatorial geometry and the SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of the proton beam characterization by multi-layer Faraday cup, spatial distribution of positron emitters obtained by the SRNA-2KG code and intercomparison of computational codes in radiation dosimetry, indicate immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in the SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumour
Ilić Radovan D.
2004-01-01
Full Text Available This paper describes the application of SRNA Monte Carlo package for proton transport simulations in complex geometry and different material composition. SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The compound nuclei decay was simulated by our own and the Russian MSDM models using ICRU 63 data. The developed package consists of two codes SRNA-2KG, which simulates proton transport in the combinatorial geometry and SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield’s data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of proton beam characterization by Multi-Layer Faraday Cup, spatial distribution of positron emitters obtained by SRNA-2KG code, and intercomparison of computational codes in radiation dosimetry, indicate the immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in SRNA pack age, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumor.
This paper describes the application of SRNA Monte Carlo package for proton transport simulations in complex geometry and different material composition. SRNA package was developed for 3D dose distribution calculation in proton therapy and dosimetry and it was based on the theory of multiple scattering. The compound nuclei decay was simulated by our own and the Russian MSDM models using ICRU 63 data. The developed package consists of two codes: SRNA-2KG, which simulates proton transport in the combinatorial geometry and SRNA-VOX, which uses the voxelized geometry using the CT data and conversion of the Hounsfield's data to tissue elemental composition. Transition probabilities for both codes are prepared by the SRNADAT code. The simulation of proton beam characterization by Multi-Layer Faraday Cup, spatial distribution of positron emitters obtained by SRNA-2KG code, and intercomparison of computational codes in radiation dosimetry, indicate the immediate application of the Monte Carlo techniques in clinical practice. In this paper, we briefly present the physical model implemented in SRNA package, the ISTAR proton dose planning software, as well as the results of the numerical experiments with proton beams to obtain 3D dose distribution in the eye and breast tumor. (author)
Zheng, Na; Xu, Hai-Bo
2015-10-01
An empirical numerical model that includes nuclear absorption, multiple Coulomb scattering and energy loss is presented for the calculation of transmission through thick objects in high energy proton radiography. In this numerical model the angular distributions are treated as Gaussians in the laboratory frame. A Monte Carlo program based on the Geant4 toolkit was developed and used for high energy proton radiography experiment simulations and verification of the empirical numerical model. The two models are used to calculate the transmission fraction of carbon and lead step-wedges in proton radiography at 24 GeV/c, and to calculate radial transmission of the French Test Object in proton radiography at 24 GeV/c with different angular cuts. It is shown that the results of the two models agree with each other, and an analysis of the slight differences is given. Supported by NSAF (11176001) and Science and Technology Developing Foundation of China Academy of Engineering Physics (2012A0202006)
Al-Sulaiti, Leena; Shipley, David; Thomas, Russell; Kacperek, Andrzej; Regan, Patrick; Palmans, Hugo
2010-07-01
The accurate conversion of dose to various materials used in clinical proton dosimetry to dose-to-water is based on fluence correction factors, accounting for attenuation of primary protons and production of secondary particles due to non-elastic nuclear interactions. This work aims to investigate the depth dose distribution and the fluence correction with respect to water or graphite at water equivalent depths (WED) in different target materials relevant for dosimetry such as polymethyl methacrylate (PMMA), graphite, A-150, aluminium and copper at 60 and 200 MeV. This was done through a comparison between Monte Carlo simulation using MCNPX 2.5.0, analytical model calculations and experimental measurements at Clatterbridge Centre of Oncology (CCO) in a 60 MeV modulated and un-modulated proton beam. MCNPX simulations indicated small fluence corrections for all materials with respect to graphite and water in 60 and 200 MeV except for aluminium. The analytical calculations showed an increase in the fluence correction factor to a few percent for all materials with respect to water at 200 MeV. The experimental measurements for 60 MeV un-modulated beam indicated a good agreement with MCNPX. For the modulated beam the fluence correction factor was found to be decreasing below unity by up to few percent with depth for aluminium and copper but almost constant and unity for A-150.
Monte Carlo calculation of lateral deflection of proton beams in homogeneous targets
Proton radiation therapy is a precise form of radiation therapy, but the avoidance of damage to critical normal tissues and the prevention of geographical tumor misses require accurate knowledge of the dose delivered to the patient and the verification of his position demand a precise imaging technique . In proton therapy facilities, the x-ray computed tomography (xCT) is the preferred technique for the planning treatment of patients. This situation has been changing nowadays with the development of proton accelerators for health care and the increase in the number of treated patients. In fact, protons could be more efficient than xCT for this task . One essential difficulty in pCT image reconstruction systems came from the scattering of the protons inside the target due to the numerous small-angle deflections by nuclear Coulomb fields, because high energy proton are slowed down when they lose their kinetic energy, mainly through ionization of the medium during multiple Coulomb scatterings (MCS). The purpose of this study is to determine the lateral deflection of proton beam (pencil beam), in the range of energy between 100 MeV and 200 MeV with Monte Carlo calculations due MCNPX (Monte Carlo N-Particle eXtended) v2.50. (author)
Lai, Bo-Lun; Sheu, Rong-Jiun; Lin, Uei-Tyng
2015-05-01
Monte Carlo simulations are generally considered the most accurate method for complex accelerator shielding analysis. Simplified models based on point-source line-of-sight approximation are often preferable in practice because they are intuitive and easy to use. A set of shielding data, including source terms and attenuation lengths for several common targets (iron, graphite, tissue, and copper) and shielding materials (concrete, iron, and lead) were generated by performing Monte Carlo simulations for 100-300 MeV protons. Possible applications and a proper use of the data set were demonstrated through a practical case study, in which shielding analysis on a typical proton treatment room was conducted. A thorough and consistent comparison between the predictions of our point-source line-of-sight model and those obtained by Monte Carlo simulations for a 360° dose distribution around the room perimeter showed that the data set can yield fairly accurate or conservative estimates for the transmitted doses, except for those near the maze exit. In addition, this study demonstrated that appropriate coupling between the generated source term and empirical formulae for radiation streaming can be used to predict a reasonable dose distribution along the maze. This case study proved the effectiveness and advantage of applying the data set to a quick shielding design and dose evaluation for proton therapy accelerators. PMID:25811254
Computed tomography with a low-intensity proton flux: results of a Monte Carlo simulation study
Schulte, Reinhard W.; Klock, Margio C. L.; Bashkirov, Vladimir; Evseev, Ivan G.; de Assis, Joaquim T.; Yevseyeva, Olga; Lopes, Ricardo T.; Li, Tianfang; Williams, David C.; Wroe, Andrew J.; Schelin, Hugo R.
2004-10-01
Conformal proton radiation therapy requires accurate prediction of the Bragg peak position. This problem may be solved by using protons rather than conventional x-rays to determine the relative electron density distribution via proton computed tomography (proton CT). However, proton CT has its own limitations, which need to be carefully studied before this technique can be introduced into routine clinical practice. In this work, we have used analytical relationships as well as the Monte Carlo simulation tool GEANT4 to study the principal resolution limits of proton CT. The GEANT4 simulations were validated by comparing them to predictions of the Bethe Bloch theory and Tschalar's theory of energy loss straggling, and were found to be in good agreement. The relationship between phantom thickness, initial energy, and the relative electron density uncertainty was systematically investigated to estimate the number of protons and dose needed to obtain a given density resolution. The predictions of this study were verified by simulating the performance of a hypothetical proton CT scanner when imaging a cylindrical water phantom with embedded density inhomogeneities. We show that a reasonable density resolution can be achieved with a relatively small number of protons, thus providing a possible dose advantage over x-ray CT.
Monte Carlo study of secondary electron production from gold nanoparticle in proton beam irradiation
Jeff Gao
2014-03-01
Full Text Available Purpose: In this study, we examined some characteristics of secondary electrons produced by gold nanoparticle (NP during proton beam irradiation.Method: By using the Geant4 Monte Carlo simulation toolkit, we simulated the NP at the range from radius (r of 17.5 nm, 25 nm, 35 nm to r = 50 nm. The proton beam energies used were 20MeV, 50MeV, and 100MeV. Findings on secondary electron production and their average kinetic energy are presented in this paper. Results: Firstly, for NP with a finite size, the secondary electron production increase with decreasing incident proton beam energy and secondary buildup existed outside NP. Secondly, the average kinetic energy of secondary electrons produced by a gold NP increased with incident proton beam energy. Thirdly, the larger the NP size, the more the secondary electron production.Conclusion: Collectively, our results suggest that apart from biological uptake efficiency, we should take the secondary electron production effect into account when considering the potential use of NPs in proton beam irradiation.-----------------------------------------------Cite this article as: Gao J, Zheng Y. Monte Carlo study of secondary electron production from gold nanoparticle in proton beam irradiation. Int J Cancer Ther Oncol 2014; 2(2:02025.DOI: http://dx.doi.org/10.14319/ijcto.0202.5
The proton therapy nozzles at Samsung Medical Center: A Monte Carlo simulation study using TOPAS
Chung, Kwangzoo; Kim, Jinsung; Kim, Dae-Hyun; Ahn, Sunghwan; Han, Youngyih
2015-07-01
To expedite the commissioning process of the proton therapy system at Samsung Medical Center (SMC), we have developed a Monte Carlo simulation model of the proton therapy nozzles by using TOol for PArticle Simulation (TOPAS). At SMC proton therapy center, we have two gantry rooms with different types of nozzles: a multi-purpose nozzle and a dedicated scanning nozzle. Each nozzle has been modeled in detail following the geometry information provided by the manufacturer, Sumitomo Heavy Industries, Ltd. For this purpose, the novel features of TOPAS, such as the time feature or the ridge filter class, have been used, and the appropriate physics models for proton nozzle simulation have been defined. Dosimetric properties, like percent depth dose curve, spreadout Bragg peak (SOBP), and beam spot size, have been simulated and verified against measured beam data. Beyond the Monte Carlo nozzle modeling, we have developed an interface between TOPAS and the treatment planning system (TPS), RayStation. An exported radiotherapy (RT) plan from the TPS is interpreted by using an interface and is then translated into the TOPAS input text. The developed Monte Carlo nozzle model can be used to estimate the non-beam performance, such as the neutron background, of the nozzles. Furthermore, the nozzle model can be used to study the mechanical optimization of the design of the nozzle.
Comparison of some popular Monte Carlo solution for proton transportation within pCT problem
Evseev, Ivan; Assis, Joaquim T. de; Yevseyeva, Olga [Universidade do Estado do Rio de Janeiro (UERJ), Nova Friburgo, RJ (Brazil). Inst. Politecnico], E-mail: evseev@iprj.uerj.br, E-mail: joaquim@iprj.uerj.br, E-mail: yevseyeva@iprj.uerj.br; Lopes, Ricardo T.; Cardoso, Jose J.B.; Silva, Ademir X. da [Universidade Federal do Rio de Janeiro (UFRJ), RJ (Brazil). Coordenacao dos Programas de Pos-graduacao de Engenharia (COPPE). Lab. de Instrumentacao Nuclear], E-mail: ricardo@lin.ufrj.br, E-mail: jjbrum@oi.com.br, E-mail: ademir@con.ufrj.br; Vinagre Filho, Ubirajara M. [Instituto de Engenharia Nuclear IEN/CNEN-RJ, Rio de Janeiro, RJ (Brazil)], E-mail: bira@ien.gov.br; Hormaza, Joel M. [UNESP, Botucatu, SP (Brazil). Inst. de Biociencias], E-mail: jmesa@ibb.unesp.br; Schelin, Hugo R.; Paschuk, Sergei A.; Setti, Joao A.P.; Milhoretto, Edney [Universidade Tecnologica Federal do Parana, Curitiba, PR (Brazil)], E-mail: schelin@cpgei.cefetpr.br, E-mail: sergei@utfpr.edu.br, E-mail: jsetti@gmail.com, E-mail: edneymilhoretto@yahoo.com
2007-07-01
The proton transport in matter is described by the Boltzmann kinetic equation for the proton flux density. This equation, however, does not have a general analytical solution. Some approximate analytical solutions have been developed within a number of significant simplifications. Alternatively, the Monte Carlo simulations are widely used. Current work is devoted to the discussion of the proton energy spectra obtained by simulation with SRIM2006, GEANT4 and MCNPX packages. The simulations have been performed considering some further applications of the obtained results in computed tomography with proton beam (pCT). Thus the initial and outgoing proton energies (3 / 300 MeV) as well as the thickness of irradiated target (water and aluminum phantoms within 90% of the full range for a given proton beam energy) were considered in the interval of values typical for pCT applications. One from the most interesting results of this comparison is that while the MCNPX spectra are in a good agreement with analytical description within Fokker-Plank approximation and the GEANT4 simulated spectra are slightly shifted from them the SRIM2006 simulations predict a notably higher mean energy loss for protons. (author)
Cirrone, G.A.P., E-mail: cirrone@lns.infn.it [Laboratori Nazionali del Sud - National Instiute for Nuclear Physics INFN (INFN-LNS), Via S.Sofia 64, 95100 Catania (Italy); Bucciolini, M. [Department of ' Fisiopatologia Clinica' , University of Florence, V.le Morgagni 85, I-50134 Florence (Italy); Bruzzi, M. [Energetic Department, University of Florence, Via S. Marta 3, I-50139 Florence (Italy); Candiano, G. [Laboratorio di Tecnologie Oncologiche HSR, Giglio Contrada, Pietrapollastra-Pisciotto, 90015 Cefalu, Palermo (Italy); Civinini, C. [National Institute for Nuclear Physics INFN, Section of Florence, Via G. Sansone 1, Sesto Fiorentino, I-50019 Florence (Italy); Cuttone, G. [Laboratori Nazionali del Sud - National Instiute for Nuclear Physics INFN (INFN-LNS), Via S.Sofia 64, 95100 Catania (Italy); Guarino, P. [Nuclear Engineering Department, University of Palermo, Via... Palermo (Italy); Laboratori Nazionali del Sud - National Instiute for Nuclear Physics INFN (INFN-LNS), Via S.Sofia 64, 95100 Catania (Italy); Lo Presti, D. [Physics Department, University of Catania, Via S. Sofia 64, I-95123, Catania (Italy); Mazzaglia, S.E. [Laboratori Nazionali del Sud - National Instiute for Nuclear Physics INFN (INFN-LNS), Via S.Sofia 64, 95100 Catania (Italy); Pallotta, S. [Department of ' Fisiopatologia Clinica' , University of Florence, V.le Morgagni 85, I-50134 Florence (Italy); Randazzo, N. [National Institute for Nuclear Physics INFN, Section of Catania, Via S.Sofia 64, 95123 Catania (Italy); Sipala, V. [National Institute for Nuclear Physics INFN, Section of Catania, Via S.Sofia 64, 95123 Catania (Italy); Physics Department, University of Catania, Via S. Sofia 64, I-95123, Catania (Italy); Stancampiano, C. [National Institute for Nuclear Physics INFN, Section of Catania, Via S.Sofia 64, 95123 Catania (Italy); and others
2011-12-01
In this paper the use of the Filtered Back Projection (FBP) Algorithm, in order to reconstruct tomographic images using the high energy (200-250 MeV) proton beams, is investigated. The algorithm has been studied in detail with a Monte Carlo approach and image quality has been analysed and compared with the total absorbed dose. A proton Computed Tomography (pCT) apparatus, developed by our group, has been fully simulated to exploit the power of the Geant4 Monte Carlo toolkit. From the simulation of the apparatus, a set of tomographic images of a test phantom has been reconstructed using the FBP at different absorbed dose values. The images have been evaluated in terms of homogeneity, noise, contrast, spatial and density resolution.
In this paper the use of the Filtered Back Projection (FBP) Algorithm, in order to reconstruct tomographic images using the high energy (200-250 MeV) proton beams, is investigated. The algorithm has been studied in detail with a Monte Carlo approach and image quality has been analysed and compared with the total absorbed dose. A proton Computed Tomography (pCT) apparatus, developed by our group, has been fully simulated to exploit the power of the Geant4 Monte Carlo toolkit. From the simulation of the apparatus, a set of tomographic images of a test phantom has been reconstructed using the FBP at different absorbed dose values. The images have been evaluated in terms of homogeneity, noise, contrast, spatial and density resolution.
Maria Grazia PiaINFN Sezione di Genova; Marcia BegalliState University Rio de Janeiro; Anton LechnerVienna University of Technology; Lina QuintieriINFN Laboratori Nazionali di Frascati; Paolo SaraccoINFN Sezione di Genova
2014-01-01
The issue of how epistemic uncertainties affect the outcome of Monte Carlo simulation is discussed by means of a concrete use case: the simulation of the longitudinal energy deposition profile of low energy protons. A variety of electromagnetic and hadronic physics models is investigated, and their effects are analyzed. Possible systematic effects are highlighted. The results identify requirements for experimental measurements capable of reducing epistemic uncertainties in the physics models.
Pia, Maria Grazia; Lechner, Anton; Quintieri, Lina; Saracco, Paolo
2010-01-01
The issue of how epistemic uncertainties affect the outcome of Monte Carlo simulation is discussed by means of a concrete use case: the simulation of the longitudinal energy deposition profile of low energy protons. A variety of electromagnetic and hadronic physics models is investigated, and their effects are analyzed. Possible systematic effects are highlighted. The results identify requirements for experimental measurements capable of reducing epistemic uncertainties in the physics models.
Feasibility Study of Neutron Dose for Real Time Image Guided Proton Therapy: A Monte Carlo Study
Kim, Jin Sung; Shin, Jung Suk; Kim, Daehyun; Shin, EunHyuk; Chung, Kwangzoo; Cho, Sungkoo; Ahn, Sung Hwan; Ju, Sanggyu; Chung, Yoonsun; Jung, Sang Hoon; Han, Youngyih
2015-01-01
Two full rotating gantry with different nozzles (Multipurpose nozzle with MLC, Scanning Dedicated nozzle) with conventional cyclotron system is installed and under commissioning for various proton treatment options at Samsung Medical Center in Korea. The purpose of this study is to investigate neutron dose equivalent per therapeutic dose, H/D, to x-ray imaging equipment under various treatment conditions with monte carlo simulation. At first, we investigated H/D with the various modifications...
Wu, D; X. T. He; Yu, W.; Fritzsche, S.
2016-01-01
A Monte-Carlo approach to proton stopping in warm dense matter is implemented into an existing particle-in-cell code. The model is based on multiple binary-collisions among electron-electron, electron-ion and ion-ion, taking into account contributions from both free and bound electrons, and allows to calculate particle stopping in much more natural manner. At low temperature limit, when ``all'' electron are bounded at the nucleus, the stopping power converges to the predictions of Bethe-Bloch...
The Proton Therapy Nozzles at Samsung Medical Center: A Monte Carlo Simulation Study using TOPAS
Chung, Kwangzoo; Kim, Dae-Hyun; Ahn, Sunghwan; Han, Youngyih
2015-01-01
To expedite the commissioning process of the proton therapy system at Samsung Medical Center (SMC), we have developed a Monte Carlo simulation model of the proton therapy nozzles using TOPAS. At SMC proton therapy center, we have two gantry rooms with different types of nozzles; a multi-purpose nozzle and a dedicated scanning nozzle. Each nozzle has been modeled in detail following the geometry information provided by the manufacturer, Sumitomo Heavy Industries, Ltd. For this purpose, novel features of TOPAS, such as the time feature or the ridge filter class, have been used. And the appropriate physics models for proton nozzle simulation were defined. Dosimetric properties, like percent depth dose curve, spread-out Bragg peak (SOBP), beam spot size, have been simulated and verified against measured beam data. Beyond the Monte Carlo nozzle modeling, we have developed an interface between TOPAS and the treatment planning system (TPS), RayStation. An exported RT plan data from the TPS has been interpreted by th...
Monte Carlo calculated stopping power ratio water/air for clinical proton therapy
In order to compute stopping-power ratios water/air for use in clinical proton dosimetry a Monte Carlo code has been developed. The main difference between the present code and other codes for proton transport is the inclusion of the detailed production of secondary electrons along the proton track. For this purpose the code is a Class-II type, where single proton-electron collisions yielding energy losses larger than a specific cut-off are considered individually. Proton multiple scattering is sampled from the complete Moliere distribution. To take into account in an approximate way the effect of inelastic nuclear collisions the fraction of the incident energy that is converted to kinetic energy of charged particles in the interaction is deposited on the spot. The energy that goes to neutral particles is assumed to leave the scoring geometry without any energy deposition. Stopping-power ratios are calculated in-line, i.e. during the transport, thereby reducing the uncertainty of the calculated value. The production and transport of the secondary electrons is used to determine an additional contribution to the stopping-power ratios obtained using the proton spectra alone
In this paper we present Monte Carlo studies of intensity modulated radiation therapy using laser-accelerated proton beams. Laser-accelerated protons coming out of a solid high-density target have broad energy and angular spectra leading to dose distributions that cannot be directly used for therapeutic applications. Through the introduction of a spectrometer-like particle selection system that delivers small pencil beams of protons with desired energy spectra it is feasible to use laser-accelerated protons for intensity modulated radiotherapy. The method presented in this paper is a three-dimensional modulation in which the proton energy spectrum and intensity of each individual beamlet are modulated to yield a homogeneous dose in both the longitudinal and lateral directions. As an evaluation of the efficacy of this method, it has been applied to two prostate cases using a variety of beam arrangements. We have performed a comparison study between intensity modulated photon plans and those for laser-accelerated protons. For identical beam arrangements and the same optimization parameters, proton plans exhibit superior coverage of the target and sparing of neighbouring critical structures. Dose-volume histogram analysis of the resulting dose distributions shows up to 50% reduction of dose to the critical structures. As the number of fields is decreased, the proton modality exhibits a better preservation of the optimization requirements on the target and critical structures. It is shown that for a two-beam arrangement (parallel-opposed) it is possible to achieve both superior target coverage with 5% dose inhomogeneity within the target and excellent sparing of surrounding tissue
Coevolution Based Adaptive Monte Carlo Localization (CEAMCL
Luo Ronghua
2008-11-01
Full Text Available An adaptive Monte Carlo localization algorithm based on coevolution mechanism of ecological species is proposed. Samples are clustered into species, each of which represents a hypothesis of the robot's pose. Since the coevolution between the species ensures that the multiple distinct hypotheses can be tracked stably, the problem of premature convergence when using MCL in highly symmetric environments can be solved. And the sample size can be adjusted adaptively over time according to the uncertainty of the robot's pose by using the population growth model. In addition, by using the crossover and mutation operators in evolutionary computation, intra-species evolution can drive the samples move towards the regions where the desired posterior density is large. So a small size of samples can represent the desired density well enough to make precise localization. The new algorithm is termed coevolution based adaptive Monte Carlo localization (CEAMCL. Experiments have been carried out to prove the efficiency of the new localization algorithm.
Coevolution Based Adaptive Monte Carlo Localization (CEAMCL)
Luo Ronghua; Hong Bingrong
2004-01-01
An adaptive Monte Carlo localization algorithm based on coevolution mechanism of ecological species is proposed. Samples are clustered into species, each of which represents a hypothesis of the robot's pose. Since the coevolution between the species ensures that the multiple distinct hypotheses can be tracked stably, the problem of premature convergence when using MCL in highly symmetric environments can be solved. And the sample size can be adjusted adaptively over time according to the unce...
Purpose: Output dependence on field size for uniform scanning beams, and the accuracy of treatment planning system (TPS) calculation are not well studied. The purpose of this work is to investigate the dependence of output on field size for uniform scanning beams and compare it among TPS calculation, measurements and Monte Carlo simulations. Methods: Field size dependence was studied using various field sizes between 2.5 cm diameter to 10 cm diameter. The field size factor was studied for a number of proton range and modulation combinations based on output at the center of spread out Bragg peak normalized to a 10 cm diameter field. Three methods were used and compared in this study: 1) TPS calculation, 2) ionization chamber measurement, and 3) Monte Carlos simulation. The XiO TPS (Electa, St. Louis) was used to calculate the output factor using a pencil beam algorithm; a pinpoint ionization chamber was used for measurements; and the Fluka code was used for Monte Carlo simulations. Results: The field size factor varied with proton beam parameters, such as range, modulation, and calibration depth, and could decrease over 10% from a 10 cm to 3 cm diameter field for a large range proton beam. The XiO TPS predicted the field size factor relatively well at large field size, but could differ from measurements by 5% or more for small field and large range beams. Monte Carlo simulations predicted the field size factor within 1.5% of measurements. Conclusion: Output factor can vary largely with field size, and needs to be accounted for accurate proton beam delivery. This is especially important for small field beams such as in stereotactic proton therapy, where the field size dependence is large and TPS calculation is inaccurate. Measurements or Monte Carlo simulations are recommended for output determination for such cases
Shielding properties of iron at high energy proton accelerators studied by a Monte Carlo code
Shielding properties of a lateral iron shield and of iron and concrete shields at angles between 5deg and 30deg are studied by means of the Monte Carlo program FLUNEV (DESY-D3 version of the FLUKA code extended for emission and transport of low energy neutrons). The following quantities were calculated for a high energy proton beam hitting an extended iron target: total and partial dose equivalents, attenuation coefficients, neutron spectra, star densities (compared also with the CASIM code) and quality factors. The dependence of the dose equivalent on the energy of primary protons, the effect of a concrete layer behind a lateral iron shielding and the total number of neutrons produced in the target were also estimated. (orig.)
TH-A-19A-10: Fast Four Dimensional Monte Carlo Dose Computations for Proton Therapy of Lung Cancer
Purpose: To develop and validate a fast and accurate four dimensional (4D) Monte Carlo (MC) dose computation system for proton therapy of lung cancer and other thoracic and abdominal malignancies in which the delivered dose distributions can be affected by respiratory motion of the patient. Methods: A 4D computer tomography (CT) scan for a lung cancer patient treated with protons in our clinic was used to create a time dependent patient model using our in-house, MCNPX-based Monte Carlo system (“MC2”). The beam line configurations for two passively scattered proton beams used in the actual treatment were extracted from the clinical treatment plan and a set of input files was created automatically using MC2. A full MC simulation of the beam line was computed using MCNPX and a set of phase space files for each beam was collected at the distal surface of the range compensator. The particles from these phase space files were transported through the 10 voxelized patient models corresponding to the 10 phases of the breathing cycle in the 4DCT, using MCNPX and an accelerated (fast) MC code called “FDC”, developed by us and which is based on the track repeating algorithm. The accuracy of the fast algorithm was assessed by comparing the two time dependent dose distributions. Results: The error of less than 1% in 100% of the voxels in all phases of the breathing cycle was achieved using this method with a speedup of more than 1000 times. Conclusion: The proposed method, which uses full MC to simulate the beam line and the accelerated MC code FDC for the time consuming particle transport inside the complex, time dependent, geometry of the patient shows excellent accuracy together with an extraordinary speed
TH-A-19A-10: Fast Four Dimensional Monte Carlo Dose Computations for Proton Therapy of Lung Cancer
Mirkovic, D; Titt, U; Mohan, R [U.T M.D. Anderson Cancer Center, Houston, TX (United States); Yepes, P [Rice University, Houston, TX (United States)
2014-06-15
Purpose: To develop and validate a fast and accurate four dimensional (4D) Monte Carlo (MC) dose computation system for proton therapy of lung cancer and other thoracic and abdominal malignancies in which the delivered dose distributions can be affected by respiratory motion of the patient. Methods: A 4D computer tomography (CT) scan for a lung cancer patient treated with protons in our clinic was used to create a time dependent patient model using our in-house, MCNPX-based Monte Carlo system (“MC{sup 2}”). The beam line configurations for two passively scattered proton beams used in the actual treatment were extracted from the clinical treatment plan and a set of input files was created automatically using MC{sup 2}. A full MC simulation of the beam line was computed using MCNPX and a set of phase space files for each beam was collected at the distal surface of the range compensator. The particles from these phase space files were transported through the 10 voxelized patient models corresponding to the 10 phases of the breathing cycle in the 4DCT, using MCNPX and an accelerated (fast) MC code called “FDC”, developed by us and which is based on the track repeating algorithm. The accuracy of the fast algorithm was assessed by comparing the two time dependent dose distributions. Results: The error of less than 1% in 100% of the voxels in all phases of the breathing cycle was achieved using this method with a speedup of more than 1000 times. Conclusion: The proposed method, which uses full MC to simulate the beam line and the accelerated MC code FDC for the time consuming particle transport inside the complex, time dependent, geometry of the patient shows excellent accuracy together with an extraordinary speed.
von Wittenau, A; Aufderheide, M B; Henderson, G L
2010-05-07
Given the cost and lead-times involved in high-energy proton radiography, it is prudent to model proposed radiographic experiments to see if the images predicted would return useful information. We recently modified our raytracing transmission radiography modeling code HADES to perform simplified Monte Carlo simulations of the transport of protons in a proton radiography beamline. Beamline objects include the initial diffuser, vacuum magnetic fields, windows, angle-selecting collimators, and objects described as distorted 2D (planar or cylindrical) meshes or as distorted 3D hexahedral meshes. We present an overview of the algorithms used for the modeling and code timings for simulations through typical 2D and 3D meshes. We next calculate expected changes in image blur as scattering materials are placed upstream and downstream of a resolution test object (a 3 mm thick sheet of tantalum, into which 0.4 mm wide slits have been cut), and as the current supplied to the focusing magnets is varied. We compare and contrast the resulting simulations with the results of measurements obtained at the 800 MeV Los Alamos LANSCE Line-C proton radiography facility.
Schach von Wittenau, Alexis E., E-mail: schachvonwittenau1@llnl.gov [Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States); Aufderheide, Maurice; Henderson, Gary [Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States)
2011-10-01
Given the cost and lead-times involved in high-energy proton radiography, it is prudent to model proposed radiographic experiments to see if the images predicted would return useful information. We recently modified our raytracing transmission radiography modeling code HADES to perform simplified Monte Carlo simulations of the transport of protons in a proton radiography beamline. Beamline objects include the initial diffuser, vacuum magnetic fields, windows, angle-selecting collimators, and objects described as distorted 2D (planar or cylindrical) meshes or as distorted 3D hexahedral meshes. We describe the algorithms used for simulations through typical 2D and 3D meshes. We calculate expected changes in image blur as scattering materials are placed upstream and downstream of a resolution test object (a 3 mm thick sheet of tantalum, into which 0.4 mm wide slits have been cut), and as the current supplied to the focusing magnets is varied. We compare and contrast the resulting simulations with the results of measurements obtained at the 800 MeV Los Alamos LANSCE Line-C proton radiography facility.
Monte Carlo approach for hadron azimuthal correlations in high energy proton and nuclear collisions
Ayala, Alejandro; Jalilian-Marian, Jamal; Magnin, J; Tejeda-Yeomans, Maria Elena
2012-01-01
We use a Monte Carlo approach to study hadron azimuthal angular correlations in high energy proton-proton and central nucleus-nucleus collisions at the BNL Relativistic Heavy Ion Collider (RHIC) energies at mid-rapidity. We build a hadron event generator that incorporates the production of $2\\to 2$ and $2\\to 3$ parton processes and their evolution into hadron states. For nucleus-nucleus collisions we include the effect of parton energy loss in the Quark-Gluon Plasma using a modified fragmentation function approach. In the presence of the medium, for the case when three partons are produced in the hard scattering, we analyze the Monte Carlo sample in parton and hadron momentum bins to reconstruct the angular correlations. We characterize this sample by the number of partons that are able to hadronize by fragmentation within the selected bins. In the nuclear environment the model allows hadronization by fragmentation only for partons with momentum above a threshold $p_T^{{\\tiny{thresh}}}=2.4$ GeV. We argue that...
Monte Carlo approach for hadron azimuthal correlations in high energy proton and nuclear collisions
Ayala, Alejandro; Dominguez, Isabel; Jalilian-Marian, Jamal; Magnin, J.; Tejeda-Yeomans, Maria Elena
2012-09-01
We use a Monte Carlo approach to study hadron azimuthal angular correlations in high-energy proton-proton and central nucleus-nucleus collisions at the BNL Relativistic Heavy Ion Collider energies at midrapidity. We build a hadron event generator that incorporates the production of 2→2 and 2→3 parton processes and their evolution into hadron states. For nucleus-nucleus collisions we include the effect of parton energy loss in the quark-gluon plasma using a modified fragmentation function approach. In the presence of the medium, for the case when three partons are produced in the hard scattering, we analyze the Monte Carlo sample in parton and hadron momentum bins to reconstruct the angular correlations. We characterize this sample by the number of partons that are able to hadronize by fragmentation within the selected bins. In the nuclear environment the model allows hadronization by fragmentation only for partons with momentum above a threshold pTthresh=2.4 GeV. We argue that one should treat properly the effect of those partons with momentum below the threshold, because their interaction with the medium may lead to showers of low-momentum hadrons along the direction of motion of the original partons as the medium becomes diluted.
A pencil beam algorithm for intensity modulated proton therapy derived from Monte Carlo simulations.
Soukup, Martin; Fippel, Matthias; Alber, Markus
2005-11-01
A pencil beam algorithm as a component of an optimization algorithm for intensity modulated proton therapy (IMPT) is presented. The pencil beam algorithm is tuned to the special accuracy requirements of IMPT, where in heterogeneous geometries both the position and distortion of the Bragg peak and the lateral scatter pose problems which are amplified by the spot weight optimization. Heterogeneity corrections are implemented by a multiple raytracing approach using fluence-weighted sub-spots. In order to derive nuclear interaction corrections, Monte Carlo simulations were performed. The contribution of long ranged products of nuclear interactions is taken into account by a fit to the Monte Carlo results. Energy-dependent stopping power ratios are also implemented. Scatter in optional beam line accessories such as range shifters or ripple filters is taken into account. The collimator can also be included, but without additional scattering. Finally, dose distributions are benchmarked against Monte Carlo simulations, showing 3%/1 mm agreement for simple heterogeneous phantoms. In the case of more complicated phantoms, principal shortcomings of pencil beam algorithms are evident. The influence of these effects on IMPT dose distributions is shown in clinical examples. PMID:16237243
Proton beam radiotherapy is an effective and non-invasive treatment for uveal melanoma. Recent research efforts have focused on improving the dosimetric accuracy of treatment planning and overcoming the present limitation of relative analytical dose calculations. Monte Carlo algorithms have been shown to accurately predict dose per monitor unit (D/MU) values, but this has yet to be shown for analytical algorithms dedicated to ocular proton therapy, which are typically less computationally expensive than Monte Carlo algorithms. The objective of this study was to determine if an analytical method could predict absolute dose distributions and D/MU values for a variety of treatment fields like those used in ocular proton therapy. To accomplish this objective, we used a previously validated Monte Carlo model of an ocular nozzle to develop an analytical algorithm to predict three-dimensional distributions of D/MU values from pristine Bragg peaks and therapeutically useful spread-out Bragg peaks (SOBPs). Results demonstrated generally good agreement between the analytical and Monte Carlo absolute dose calculations. While agreement in the proximal region decreased for beams with less penetrating Bragg peaks compared with the open-beam condition, the difference was shown to be largely attributable to edge-scattered protons. A method for including this effect in any future analytical algorithm was proposed. Comparisons of D/MU values showed typical agreement to within 0.5%. We conclude that analytical algorithms can be employed to accurately predict absolute proton dose distributions delivered by an ocular nozzle.
Effect of elemental compositions on Monte Carlo dose calculations in proton therapy of eye tumors
Rasouli, Fatemeh S.; Farhad Masoudi, S.; Keshazare, Shiva; Jette, David
2015-12-01
Recent studies in eye plaque brachytherapy have found considerable differences between the dosimetric results by using a water phantom, and a complete human eye model. Since the eye continues to be simulated as water-equivalent tissue in the proton therapy literature, a similar study for investigating such a difference in treating eye tumors by protons is indispensable. The present study inquires into this effect in proton therapy utilizing Monte Carlo simulations. A three-dimensional eye model with elemental compositions is simulated and used to examine the dose deposition to the phantom. The beam is planned to pass through a designed beam line to moderate the protons to the desired energies for ocular treatments. The results are compared with similar irradiation to a water phantom, as well as to a material with uniform density throughout the whole volume. Spread-out Bragg peaks (SOBPs) are created by adding pristine peaks to cover a typical tumor volume. Moreover, the corresponding beam parameters recommended by the ICRU are calculated, and the isodose curves are computed. The results show that the maximum dose deposited in ocular media is approximately 5-7% more than in the water phantom, and about 1-1.5% less than in the homogenized material of density 1.05 g cm-3. Furthermore, there is about a 0.2 mm shift in the Bragg peak due to the tissue composition difference between the models. It is found that using the weighted dose profiles optimized in a water phantom for the realistic eye model leads to a small disturbance of the SOBP plateau dose. In spite of the plaque brachytherapy results for treatment of eye tumors, it is found that the differences between the simplified models presented in this work, especially the phantom containing the homogenized material, are not clinically significant in proton therapy. Taking into account the intrinsic uncertainty of the patient dose calculation for protons, and practical problems corresponding to applying patient
A Monte Carlo pencil beam scanning model for proton treatment plan simulation using GATE/GEANT4
Grevillot, L; Freud, N; Sarrut, D [Universite de Lyon, CREATIS, CNRS UMR5220, Inserm U1044, INSA-Lyon, Universite Lyon 1, Centre Leon Berard, Lyon (France); Bertrand, D; Dessy, F, E-mail: loic.grevillot@creatis.insa-lyon.fr [IBA, B-1348, Louvain-la Neuve (Belgium)
2011-08-21
This work proposes a generic method for modeling scanned ion beam delivery systems, without simulation of the treatment nozzle and based exclusively on beam data library (BDL) measurements required for treatment planning systems (TPS). To this aim, new tools dedicated to treatment plan simulation were implemented in the Gate Monte Carlo platform. The method was applied to a dedicated nozzle from IBA for proton pencil beam scanning delivery. Optical and energy parameters of the system were modeled using a set of proton depth-dose profiles and spot sizes measured at 27 therapeutic energies. For further validation of the beam model, specific 2D and 3D plans were produced and then measured with appropriate dosimetric tools. Dose contributions from secondary particles produced by nuclear interactions were also investigated using field size factor experiments. Pristine Bragg peaks were reproduced with 0.7 mm range and 0.2 mm spot size accuracy. A 32 cm range spread-out Bragg peak with 10 cm modulation was reproduced with 0.8 mm range accuracy and a maximum point-to-point dose difference of less than 2%. A 2D test pattern consisting of a combination of homogeneous and high-gradient dose regions passed a 2%/2 mm gamma index comparison for 97% of the points. In conclusion, the generic modeling method proposed for scanned ion beam delivery systems was applicable to an IBA proton therapy system. The key advantage of the method is that it only requires BDL measurements of the system. The validation tests performed so far demonstrated that the beam model achieves clinical performance, paving the way for further studies toward TPS benchmarking. The method involves new sources that are available in the new Gate release V6.1 and could be further applied to other particle therapy systems delivering protons or other types of ions like carbon.
Purpose: Monte Carlo simulation plays an important role for proton Pencil Beam Scanning (PBS) technique. However, MC simulation demands high computing power and is limited to few large proton centers that can afford a computer cluster. We study the feasibility of utilizing cloud computing in the MC simulation of PBS beams. Methods: A GATE/GEANT4 based MC simulation software was installed on a commercial cloud computing virtual machine (Linux 64-bits, Amazon EC2). Single spot Integral Depth Dose (IDD) curves and in-air transverse profiles were used to tune the source parameters to simulate an IBA machine. With the use of StarCluster software developed at MIT, a Linux cluster with 2–100 nodes can be conveniently launched in the cloud. A proton PBS plan was then exported to the cloud where the MC simulation was run. Results: The simulated PBS plan has a field size of 10×10cm2, 20cm range, 10cm modulation, and contains over 10,000 beam spots. EC2 instance type m1.medium was selected considering the CPU/memory requirement and 40 instances were used to form a Linux cluster. To minimize cost, master node was created with on-demand instance and worker nodes were created with spot-instance. The hourly cost for the 40-node cluster was $0.63 and the projected cost for a 100-node cluster was $1.41. Ten million events were simulated to plot PDD and profile, with each job containing 500k events. The simulation completed within 1 hour and an overall statistical uncertainty of < 2% was achieved. Good agreement between MC simulation and measurement was observed. Conclusion: Cloud computing is a cost-effective and easy to maintain platform to run proton PBS MC simulation. When proton MC packages such as GATE and TOPAS are combined with cloud computing, it will greatly facilitate the pursuing of PBS MC studies, especially for newly established proton centers or individual researchers
Accelerated GPU based SPECT Monte Carlo simulations.
Garcia, Marie-Paule; Bert, Julien; Benoit, Didier; Bardiès, Manuel; Visvikis, Dimitris
2016-06-01
Monte Carlo (MC) modelling is widely used in the field of single photon emission computed tomography (SPECT) as it is a reliable technique to simulate very high quality scans. This technique provides very accurate modelling of the radiation transport and particle interactions in a heterogeneous medium. Various MC codes exist for nuclear medicine imaging simulations. Recently, new strategies exploiting the computing capabilities of graphical processing units (GPU) have been proposed. This work aims at evaluating the accuracy of such GPU implementation strategies in comparison to standard MC codes in the context of SPECT imaging. GATE was considered the reference MC toolkit and used to evaluate the performance of newly developed GPU Geant4-based Monte Carlo simulation (GGEMS) modules for SPECT imaging. Radioisotopes with different photon energies were used with these various CPU and GPU Geant4-based MC codes in order to assess the best strategy for each configuration. Three different isotopes were considered: (99m) Tc, (111)In and (131)I, using a low energy high resolution (LEHR) collimator, a medium energy general purpose (MEGP) collimator and a high energy general purpose (HEGP) collimator respectively. Point source, uniform source, cylindrical phantom and anthropomorphic phantom acquisitions were simulated using a model of the GE infinia II 3/8" gamma camera. Both simulation platforms yielded a similar system sensitivity and image statistical quality for the various combinations. The overall acceleration factor between GATE and GGEMS platform derived from the same cylindrical phantom acquisition was between 18 and 27 for the different radioisotopes. Besides, a full MC simulation using an anthropomorphic phantom showed the full potential of the GGEMS platform, with a resulting acceleration factor up to 71. The good agreement with reference codes and the acceleration factors obtained support the use of GPU implementation strategies for improving computational
Accelerated GPU based SPECT Monte Carlo simulations
Garcia, Marie-Paule; Bert, Julien; Benoit, Didier; Bardiès, Manuel; Visvikis, Dimitris
2016-06-01
Monte Carlo (MC) modelling is widely used in the field of single photon emission computed tomography (SPECT) as it is a reliable technique to simulate very high quality scans. This technique provides very accurate modelling of the radiation transport and particle interactions in a heterogeneous medium. Various MC codes exist for nuclear medicine imaging simulations. Recently, new strategies exploiting the computing capabilities of graphical processing units (GPU) have been proposed. This work aims at evaluating the accuracy of such GPU implementation strategies in comparison to standard MC codes in the context of SPECT imaging. GATE was considered the reference MC toolkit and used to evaluate the performance of newly developed GPU Geant4-based Monte Carlo simulation (GGEMS) modules for SPECT imaging. Radioisotopes with different photon energies were used with these various CPU and GPU Geant4-based MC codes in order to assess the best strategy for each configuration. Three different isotopes were considered: 99m Tc, 111In and 131I, using a low energy high resolution (LEHR) collimator, a medium energy general purpose (MEGP) collimator and a high energy general purpose (HEGP) collimator respectively. Point source, uniform source, cylindrical phantom and anthropomorphic phantom acquisitions were simulated using a model of the GE infinia II 3/8" gamma camera. Both simulation platforms yielded a similar system sensitivity and image statistical quality for the various combinations. The overall acceleration factor between GATE and GGEMS platform derived from the same cylindrical phantom acquisition was between 18 and 27 for the different radioisotopes. Besides, a full MC simulation using an anthropomorphic phantom showed the full potential of the GGEMS platform, with a resulting acceleration factor up to 71. The good agreement with reference codes and the acceleration factors obtained support the use of GPU implementation strategies for improving computational efficiency
Krause, Claudius
2012-04-15
High energy proton-proton collisions lead to a large amount of secondary particles to be measured in a detector. A final state containing top quarks is of particular interest. But top quarks are only produced in a small fraction of the collisions. Hence, criteria must be defined to separate events containing top quarks from the background. From detectors, we record signals, for example hits in the tracker system or deposits in the calorimeters. In order to obtain the momentum of the particles, we apply algorithms to reconstruct tracks in space. More sophisticated algorithms are needed to identify the flavour of quarks, such as b-tagging. Several steps are needed to test these algorithms. Collision products of proton-proton events are generated using Monte Carlo techniques and their passage through the detector is simulated. After that, the algorithms are applied and the signal efficiency and the mistagging rate can be obtained. There are, however, many different approaches and algorithms realized in programs, so the question arises if the choice of the Monte Carlo generator influences the measured quantities. In this thesis, two commonly used Monte Carlo generators, SHERPA and MadGraph/MadEvent, are compared and the differences in the selection efficiency of semimuonic tt events are estimated. In addition, the distributions of kinematic variables are shown. A special chapter about the matching of matrix elements with parton showers is included. The main algorithms, CKKW for SHERPA and MLM for MadGraph/MadEvent, are introduced.
Use of Monte Carlo software to aid design of a proton therapy nozzle
Swanepoel, M.W. [Medical Radiation Group, iThemba LABS, P.O. Box 22, Somerset West 7129 (South Africa)], E-mail: mark@tlabs.ac.za; Jones, D.T.L. [Medical Radiation Group, iThemba LABS, P.O. Box 22, Somerset West 7129 (South Africa)
2007-09-21
A second proton therapy nozzle is being developed at iThemba LABS to irradiate lesions in the body, thus complementing an existing facility for head and neck treatments. A passive scattering system is being developed, the complexity of which necessitates Monte Carlo simulations. We have used MCNPX to set the apertures and spacing of collimators, to model dose distributions in water, to check and modify beam scattering and energy modulating components, and to check radiation shields. The comprehensive shielding model was adapted for other problems by reducing the types of particles transported, limiting the extent and complexity of the geometry, and where possible killing particles by setting their importance to zero. Our results appear to indicate that the Rossi and Greisen description of multiple Coulomb scattering as used in MCNPX predicts high-Z, large angle scattering acceptably well for modeling proton therapy nozzles. MCNPX is easy to learn and implement, but has disadvantages when used to model therapy nozzles: (1) it does not yet offer a true capability to model electromagnetic interactions, (2) it cannot model moving components, and (3) it uses energy rather than range cut-offs for particles. Hence a GEANT4 model of the new nozzle is also being implemented.
Improved efficiency in Monte Carlo simulation for passive-scattering proton therapy
The aim of this work was to improve the computational efficiency of Monte Carlo simulations when tracking protons through a proton therapy treatment head. Two proton therapy facilities were considered, the Francis H Burr Proton Therapy Center (FHBPTC) at the Massachusetts General Hospital and the Crocker Lab eye treatment facility used by University of California at San Francisco (UCSFETF). The computational efficiency was evaluated for phase space files scored at the exit of the treatment head to determine optimal parameters to improve efficiency while maintaining accuracy in the dose calculation.For FHBPTC, particles were split by a factor of 8 upstream of the second scatterer and upstream of the aperture. The radius of the region for Russian roulette was set to 2.5 or 1.5 times the radius of the aperture and a secondary particle production cut (PC) of 50 mm was applied. For UCSFETF, particles were split a factor of 16 upstream of a water absorber column and upstream of the aperture. Here, the radius of the region for Russian roulette was set to 4 times the radius of the aperture and a PC of 0.05 mm was applied. In both setups, the cylindrical symmetry of the proton beam was exploited to position the split particles randomly spaced around the beam axis.When simulating a phase space for subsequent water phantom simulations, efficiency gains between a factor of 19.9 ± 0.1 and 52.21 ± 0.04 for the FHTPC setups and 57.3 ± 0.5 for the UCSFETF setups were obtained. For a phase space used as input for simulations in a patient geometry, the gain was a factor of 78.6 ± 7.5. Lateral-dose curves in water were within the accepted clinical tolerance of 2%, with statistical uncertainties of 0.5% for the two facilities. For the patient geometry and by considering the 2% and 2mm criteria, 98.4% of the voxels showed a gamma index lower than unity. An analysis of the dose distribution resulted in systematic deviations below of 0.88% for 20
Chatzikyriakou, Eleni; Smyrnis, Chris; Chatwin, Chris
2014-01-01
Carbon-based nanotechnology electronics can provide high performance, low-power and low-weight solutions, which are very suitable for innovative aerospace applications. However, its application in the space environment where there is a radiation hazard, requires an assessment of the response of such electronic products to the background irradiance. To explore the potential of carbon-based nanotechnology, Monte Carlo simulations of radiation interacting with a gate-all-around carbo...
Wu, D; Yu, W; Fritzsche, S
2016-01-01
A Monte-Carlo approach to proton stopping in warm dense matter is implemented into an existing particle-in-cell code. The model is based on multiple binary-collisions among electron-electron, electron-ion and ion-ion, taking into account contributions from both free and bound electrons, and allows to calculate particle stopping in much more natural manner. At low temperature limit, when ``all'' electron are bounded at the nucleus, the stopping power converges to the predictions of Bethe-Bloch theory, which shows good consistency with data provided by the NIST. With the rising of temperatures, more and more bound electron are ionized, thus giving rise to an increased stopping power to cold matter, which is consistent with the report of a recently experimental measurement [Phys. Rev. Lett. 114, 215002 (2015)]. When temperature is further increased, with ionizations reaching the maximum, lowered stopping power is observed, which is due to the suppression of collision frequency between projected proton beam and h...
8 particle histories were run. Results: Range measurements of the Monte-Carlo simulations matched the measured data within 1mm. Distal fall-off of the simulated fields matched within <1mm. Lateral penumbra and field size measurements of the standard-sized square and half-beam blocked fields matched within 1mm at all three planes compared. A small difference was seen in the in-air profiles at doses <0%. The suspected cause of the difference was the aperture shape. The measured data utilized a divergent aperture. The Monte-Carlo calculation used a non-divergent aperture. Conclusion: The validation measurements indicate that we were able to accurately model the MEVION s250 Proton therapy system using Monte-Carlo Calculations. This may reduce the commissioning time for future users. Purpose: Monte-Carlo modeling is an important tool for understanding the behavior of therapeutic proton beams in a heterogeneous media such as the patient. To gain confidence that a Monte-Carlo model is accurate in complex geometries and media, it must first be compared with measurement in simple situations. This study documents the validation of our Monte-Carlo Model. Methods: A model of the MEVION s250 Proton therapy system was created in the TOPAS Monte-Carlo environment using machine geometry and field shaping system information provided by the vendor. For each of 24 options, validation of the TOPAS model was performed by comparing the dose scored by TOPAS to the dose measurements obtained during the commissioning of the treatment planning system. The measurements compared consisted of: pristine peak depth-dose profiles, in-air profiles for a standard-sized square field (20cm×20cm or 10cm×10cm depending on the maximum field size for each option) at isocenter and at 20cm upstream and downstream of isocenter, and in-air profiles with a half-beam blocked aperture at isocenter and at 20cm upstream and downstream of isocenter. For all Monte-Carlo simulations,
CERN Summer Student Report 2016 Monte Carlo Data Base Improvement
Caciulescu, Alexandru Razvan
2016-01-01
During my Summer Student project I worked on improving the Monte Carlo Data Base and MonALISA services for the ALICE Collaboration. The project included learning the infrastructure for tracking and monitoring of the Monte Carlo productions as well as developing a new RESTful API for seamless integration with the JIRA issue tracking framework.
Feasibility Study of Neutron Dose for Real Time Image Guided Proton Therapy: A Monte Carlo Study
Kim, Jin Sung; Kim, Daehyun; Shin, EunHyuk; Chung, Kwangzoo; Cho, Sungkoo; Ahn, Sung Hwan; Ju, Sanggyu; Chung, Yoonsun; Jung, Sang Hoon; Han, Youngyih
2015-01-01
Two full rotating gantry with different nozzles (Multipurpose nozzle with MLC, Scanning Dedicated nozzle) with conventional cyclotron system is installed and under commissioning for various proton treatment options at Samsung Medical Center in Korea. The purpose of this study is to investigate neutron dose equivalent per therapeutic dose, H/D, to x-ray imaging equipment under various treatment conditions with monte carlo simulation. At first, we investigated H/D with the various modifications of the beam line devices (Scattering, Scanning, Multi-leaf collimator, Aperture, Compensator) at isocenter, 20, 40, 60 cm distance from isocenter and compared with other research groups. Next, we investigated the neutron dose at x-ray equipments used for real time imaging with various treatment conditions. Our investigation showed the 0.07 ~ 0.19 mSv/Gy at x-ray imaging equipments according to various treatment options and intestingly 50% neutron dose reduction effect of flat panel detector was observed due to multi- lea...
Formal quality control for a proton Monte Carlo system in radiation therapy
TOPAS (TOol for PArticle Simulation) is a Monte Carlo particle transport tool being released to a wide variety of proton therapy users worldwide. Because TOPAS provides unprecedented ease in 4D placement of geometry components, beam sources and scoring, including options to place geometry components, beam sources or scorers within each other, Quality Control (QC) for TOPAS is both critical and challenging. All simulation details (geometry, particle sources, scoring, physics settings, time-dependent motions, gating, etc.) are specified in the TOPAS Parameter Control System (which catches many user errors). QC includes Unit and End-to-End Testing. Each code unit is tested (each geometry component, particle source option, scoring option, etc.) and these unit testing procedures are shared with end users so they can reproduce tests. End-to-End testing of several full clinical setups is routinely performed. End-to-End testing presents a challenge since one cannot anticipate all the ways users will combine TOPAS flexible units for their specific project. Automated checking catches geometry overlaps and some other problematic setups, but one can never rule out the potential for problems when users combine units in new setups. QC is ultimately a partnership between the tool developer and the user. Key is that the developer be clear to the end user about what has been tested and what has not.
Monte Carlo simulations play a crucial role for in-vivo treatment monitoring based on PET and prompt gamma imaging in proton and carbon-ion therapies. The accuracy of the nuclear fragmentation models implemented in these codes might affect the quality of the treatment verification. In this paper, we investigate the nuclear models implemented in GATE/Geant4 and FLUKA by comparing the angular and energy distributions of secondary particles exiting a homogeneous target of PMMA. Comparison results were restricted to fragmentation of 16O and 12C. Despite the very simple target and set-up, substantial discrepancies were observed between the two codes. For instance, the number of high energy (>1 MeV) prompt gammas exiting the target was about twice as large with GATE/Geant4 than with FLUKA both for proton and carbon ion beams. Such differences were not observed for the predicted annihilation photon production yields, for which ratios of 1.09 and 1.20 were obtained between GATE and FLUKA for the proton beam and the carbon ion beam, respectively. For neutrons and protons, discrepancies from 14% (exiting protons–carbon ion beam) to 57% (exiting neutrons–proton beam) have been identified in production yields as well as in the energy spectra for neutrons. (paper)
Cho, S; Shin, E H; Kim, J; Ahn, S H; Chung, K; Kim, D-H; Han, Y; Choi, D H [Samsung Medical Center, Seoul (Korea, Republic of)
2015-06-15
Purpose: To evaluate the shielding wall design to protect patients, staff and member of the general public for secondary neutron using a simply analytic solution, multi-Monte Carlo code MCNPX, ANISN and FLUKA. Methods: An analytical and multi-Monte Carlo method were calculated for proton facility (Sumitomo Heavy Industry Ltd.) at Samsung Medical Center in Korea. The NCRP-144 analytical evaluation methods, which produced conservative estimates on the dose equivalent values for the shielding, were used for analytical evaluations. Then, the radiation transport was simulated with the multi-Monte Carlo code. The neutron dose at evaluation point is got by the value using the production of the simulation value and the neutron dose coefficient introduced in ICRP-74. Results: The evaluation points of accelerator control room and control room entrance are mainly influenced by the point of the proton beam loss. So the neutron dose equivalent of accelerator control room for evaluation point is 0.651, 1.530, 0.912, 0.943 mSv/yr and the entrance of cyclotron room is 0.465, 0.790, 0.522, 0.453 mSv/yr with calculation by the method of NCRP-144 formalism, ANISN, FLUKA and MCNP, respectively. The most of Result of MCNPX and FLUKA using the complicated geometry showed smaller values than Result of ANISN. Conclusion: The neutron shielding for a proton therapy facility has been evaluated by the analytic model and multi-Monte Carlo methods. We confirmed that the setting of shielding was located in well accessible area to people when the proton facility is operated.
Highlights: • Angular distribution of the proton beam in a gaseous environment. • Particle energy distribution profile and proton flux within gas-cell target with MCNPX. • Detection of the residual nuclei during the nuclear reactions. • Estimation of production yield for 22,natNe(p,x)22Na reactions. - Abstract: The 22Ne(p,n)22Na is an optimal reaction for the cyclotron production of 22Na. This work tends to monitor the proton induced production of 22Na in a gas-cell target, containing natural and enriched neon gas, using Monte Carlo method. The excitation functions of reactions are calculated by both TALYS-1.6 and ALICE/ASH codes and then the optimum energy range of projectile for the high yield production is selected. A free gaseous environment of neon at a particular pressure and temperature is prearranged and the proton beam is transported within it using Monte Carlo codes MCNPX and SRIM. The beam monitoring performed by each of these codes indicates that the gas-cell has to be designed as conical frustum to reach desired interactions. The MCNPX is also employed to calculate the energy distribution of proton in the designed target and estimation of the residual nuclei during irradiation. The production yield of 22Na in 22Ne(p,n)22Na and natNe(p,x)22Na reactions are estimated and it shows a good agreement with the experimental results. The results demonstrate that Monte Carlo makes available a beneficial manner to design and optimize the gas targets as well as calibration of detectors, which can be used for the radionuclide production purposes
MONTE: An automated Monte Carlo based approach to nuclear magnetic resonance assignment of proteins
Hitchens, T. Kevin; Lukin, Jonathan A.; Zhan Yiping; McCallum, Scott A.; Rule, Gordon S. [Carnegie Mellon University, Department of Biological Sciences (United States)], E-mail: rule@andrew.cmu.edu
2003-01-15
A general-purpose Monte Carlo assignment program has been developed to aid in the assignment of NMR resonances from proteins. By virtue of its flexible data requirements the program is capable of obtaining assignments of both heavily deuterated and fully protonated proteins. A wide variety of source data, such as inter-residue scalar connectivity, inter-residue dipolar (NOE) connectivity, and residue specific information, can be utilized in the assignment process. The program can also use known assignments from one form of a protein to facilitate the assignment of another form of the protein. This attribute is useful for assigning protein-ligand complexes when the assignments of the unliganded protein are known. The program can be also be used as an interactive research tool to assist in the choice of additional experimental data to facilitate completion of assignments. The assignment of a deuterated 45 kDa homodimeric Glutathione-S-transferase illustrates the principal features of the program.
Investigations on Monte Carlo based coupled core calculations
The present trend in advanced and next generation nuclear reactor core designs is towards increased material heterogeneity and geometry complexity. The continuous energy Monte Carlo method has the capability of modeling such core environments with high accuracy. This paper presents results from feasibility studies being performed at the Pennsylvania State University (PSU) on both accelerating Monte Carlo criticality calculations by using hybrid nodal diffusion Monte Carlo schemes and thermal-hydraulic feedback modeling in Monte Carlo core calculations. The computation process is greatly accelerated by calculating the three-dimensional (3D) distributions of fission source and thermal-hydraulics parameters with the coupled NEM/COBRA-TF code and then using coupled MCNP5/COBRA-TF code to fine tune the results to obtain an increased accuracy. The PSU NEM code employs cross-sections generated by MCNP5 for pin-cell based nodal compositions. The implementation of different code modifications facilitating coupled calculations are presented first. Then the coupled hybrid Monte Carlo based code system is applied to a 3D 2*2 pin array extracted from a Boiling Water Reactor (BWR) assembly with reflective radial boundary conditions. The obtained results are discussed and it is showed that performing Monte-Carlo based coupled core steady state calculations are feasible. (authors)
Purpose: The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the GEANT 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from GEANT 4 simulations to correlate with biological effectiveness of therapeutic protons. Methods: The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LETt and dose-averaged LET, LETd) using GEANT 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LETt and LETd of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose. Results: Simulation results showed that the step limit effect was small for LETt but significant for LETd. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in GEANT 4 can result in incorrect LETd calculation results in the dose plateau region for small step limits. The erroneous LETd results can be attributed to the algorithm to determine fluctuations in energy deposition along the
Fluctuations in the EAS radio signal derived with improved Monte Carlo simulations based on CORSIKA
Huege, T; Badea, F; Bähren, L; Bekk, K; Bercuci, A; Bertaina, M; Biermann, P L; Blumer, J; Bozdog, H; Brancus, I M; Buitink, S; Bruggemann, M; Buchholz, P; Butcher, H; Chiavassa, A; Daumiller, K; De Bruyn, A G; De Vos, C M; Di Pierro, F; Doll, P; Engel, R; Falcke, H; Gemmeke, H; Ghia, P L; Glasstetter, R; Grupen, C; Haungs, A; Heck, D; Hörandel, J R; Horneffer, A; Kampert, K H; Kant, G W; Klein, U; Kolotaev, Yu; Koopman, Y; Krömer, O; Kuijpers, J; Lafebre, S; Maier, G; Mathes, H J; Mayer, H J; Milke, J; Mitrica, B; Morello, C; Navarra, G; Nehls, S; Nigl, A; Obenland, R; Oehlschläger, J; Ostapchenko, S; Over, S; Pepping, H J; Petcu, M; Petrovic, J; Pierog, T; Plewnia, S; Rebel, H; Risse, A; Roth, M; Schieler, H; Schoonderbeek, G; Sima, O; Stumpert, M; Toma, G; Trinchero, G C; Ulrich, H; Valchierotti, S; Van Buren, J; Van Capellen, W; Walkowiak, W; Weindl, A; Wijnholds, S J; Wochele, J; Zabierowski, J; Zensus, J A; Zimmermann, D; Bowman, J D; Huege, Tim
2005-01-01
Cosmic ray air showers are known to emit pulsed radio emission which can be understood as coherent geosynchrotron radiation arising from the deflection of electron-positron pairs in the earth's magnetic field. Here, we present simulations carried out with an improved version of our Monte Carlo code for the calculation of geosynchrotron radiation. Replacing the formerly analytically parametrised longitudinal air shower development with CORSIKA-generated longitudinal profiles, we study the radio flux variations arising from inherent fluctuations between individual air showers. Additionally, we quantify the dependence of the radio emission on the nature of the primary particle by comparing the emission generated by proton- and iron-induced showers. This is only the first step in the incorporation of a more realistic air shower model into our Monte Carlo code. The inclusion of highly realistic CORSIKA-based particle energy, momentum and spatial distributions together with an analytical treatment of ionisation los...
Palmans, H. [Ghent Univ. (Belgium). Dept. of Biomedical Physics; Verhaegen, F.
1995-12-01
In the last decade, several clinical proton beam therapy facilities have been developed. To satisfy the demand for uniformity in clinical (routine) proton beam dosimetry two dosimetry protocols (ECHED and AAPM) have been published. Both protocols neglect the influence of ion chamber dependent parameters on dose determination in proton beams because of the scatter properties of these beams, although the problem has not been studied thoroughly yet. A comparison between water calorimetry and ionisation chamber dosimetry showed a discrepancy of 2.6% between the former method and ionometry following the ECHED protocol. Possibly, a small part of this difference can be attributed to chamber dependent correction factors. Indications for this possibility are found in ionometry measurements. To allow the simulation of complex geometries with different media necessary for the study of those corrections, an existing proton Monte Carlo code (PTRAN, Berger) has been modified. The original code, that applies Mollire`s multiple scattering theory and Vavilov`s energy straggling theory, calculates depth dose profiles, energy distributions and radial distributions for pencil beams in water. Comparisons with measurements and calculations reported in the literature are done to test the program`s accuracy. Preliminary results of the influence of chamber design and chamber materials on dose to water determination are presented.
In the last decade, several clinical proton beam therapy facilities have been developed. To satisfy the demand for uniformity in clinical (routine) proton beam dosimetry two dosimetry protocols (ECHED and AAPM) have been published. Both protocols neglect the influence of ion chamber dependent parameters on dose determination in proton beams because of the scatter properties of these beams, although the problem has not been studied thoroughly yet. A comparison between water calorimetry and ionisation chamber dosimetry showed a discrepancy of 2.6% between the former method and ionometry following the ECHED protocol. Possibly, a small part of this difference can be attributed to chamber dependent correction factors. Indications for this possibility are found in ionometry measurements. To allow the simulation of complex geometries with different media necessary for the study of those corrections, an existing proton Monte Carlo code (PTRAN, Berger) has been modified. The original code, that applies Mollire's multiple scattering theory and Vavilov's energy straggling theory, calculates depth dose profiles, energy distributions and radial distributions for pencil beams in water. Comparisons with measurements and calculations reported in the literature are done to test the program's accuracy. Preliminary results of the influence of chamber design and chamber materials on dose to water determination are presented
Inside an IBA proton therapy centre, secondary neutrons are produced due to nuclear interactions of the proton beam with matter mainly inside the cyclotron, the beam line, the treatment nozzle and the patient. Accurate measurements of the neutron ambient dose equivalent H*(10) in such a facility require the use of a detector that has a good sensitivity for neutrons ranging from thermal energies up to 230 MeV, such as for instance the WENDI-2 detector. WENDI-2 measurements have been performed at the Westdeutsches Protonentherapiezentrum Essen, at several positions around the cyclotron room and around a gantry treatment room operated in two different beam delivery modes: Pencil Beam Scanning and Double Scattering. These measurements are compared with Monte Carlo simulation results for the neutron H*(10) obtained with MCNPX 2.5.0 and GEANT4 9.6. In proton therapy, proton beams with energies up to typically 230 MeV are used to treat cancerous tumours very efficiently while sparing surrounding healthy tissues as much as possible. Due to nuclear interactions of the proton beams with matter, mainly inside the cyclotron, the beam line, the treatment nozzle and the patient, secondary neutrons with energies up to 230 MeV are unfortunately produced, as well as photons up to ∼10 MeV. Behind the thick concrete shielding walls which are necessary to attenuate the stray radiation fields, the total ambient dose equivalent H*(10) is very large due to the neutron component. In shielding studies for proton therapy facilities, the neutron H*(10) component is often evaluated using the Monte Carlo codes MCNPX(5), FLUKA(6) or PHITS(7). Recent benchmark simulations performed with GEANT4 have shown that this code would also be a suitable tool for the shielding studies of proton therapy centres. The experimental validation of such shielding studies requires the use of a detector with a good sensitivity for neutrons ranging from thermal energies up to 230 MeV, such as for example the
MCHIT - Monte Carlo model for proton and heavy-ion therapy
Pshenichnov, Igor; Greiner, Walter
2007-01-01
We study the propagation of nucleons and nuclei in tissue-like media within a Monte Carlo Model for Heavy-ion Therapy (MCHIT) based on the GEANT4 toolkit (version 8.2). The model takes into account fragmentation of projectile nuclei and secondary interactions of produced nuclear fragments. Model predictions are validated with available experimental data obtained for water and PMMA phantoms irradiated by monoenergetic carbon-ion beams. The MCHIT model describes well (1) the depth-dose distributions in water and PMMA, (2) the doses measured for fragments of certain charge, (3) the distributions of positron emitting nuclear fragments produced by carbon-ion beams, and (4) the energy spectra of secondary neutrons measured at different angles to the beam direction. Radial dose profiles for primary nuclei and for different projectile fragments are calculated and discussed as possible input for evaluation of biological dose distributions. It is shown that at the periphery of the transverse dose profile close to the B...
Comparison of Monte Carlo simulations with proton experiment for a thick Au absorber
Proton therapy applications deal with relatively thick targets like the human head or the trunk. Therefore, relatively small differences in the total proton stopping power given, for example, by the different models provided by GEANT4, could lead to significant disagreement in the final proton energy spectra when integrated along lengthy proton trajectories. This work presents a comparison of proton energy spectra for 49.1 MeV protons passing through a couple of Au absorbers with different thicknesses obtained by GEANT4.8.2 simulations using ICRU49, Ziegler1985 and Ziegler2000 models. The comparison was made with the experimental data of Tschalaer, with TRIM/SRIM2008 and MCNPX2.4.0 simulations, and the Payne analytical solution for the transport equation in the Fokker-Plank approximation. It is shown that the simulations reproduce the experimental spectra with some detectable contradictions. It should be noted that all the spectra lay at the proton energies significantly above 2 MeV, i.e. in the so-called 'Bethe-Bloch region'. Therefore the observed disagreements in GEANT4 results, simulated with different models, are somewhat unexpected. Further studies for a better understanding and to obtain definitive conclusions are necessary. (author)
Comparison of Monte Carlo simulations with proton experiment for a thick Au absorber
Yevseyeva, Olga; Assis, Joaquim T. de, E-mail: yevseveva@iprj.uerj.b, E-mail: joaquim@iprj.uerj.b [Universidade do Estado do Rio de Janeiro (UERJ), Nova Friburgo, RJ (Brazil). Inst. Politecnico; Evseev, Ivan; Schelin, Hugo R.; Paschuk, Sergei A.; Milhoretto, Edney; Setti, Joao A.P., E-mail: evseev@utfpr.edu.b, E-mail: schelin@utfpr.edu.b, E-mail: sergei@utfpr.edu.b, E-mail: edneymilhoretto@yahoo.co, E-mail: jsetti@gmail.co [Universidade Tecnologica Federal do Parana, Curitiba, PR (Brazil); Diaz, Katherin S., E-mail: kshtejer@infomed.sld.c [Centro de Aplicaciones Tecnologicas y Desarrollo Nuclear, Havana (Cuba); Hormaza, Joel M., E-mail: jmesa@ibb.unesp.b [UNESP, Botucatu, SP (Brazil). Inst. de Biociencias; Lopes, Ricardo T., E-mail: ricardo@lin.ufrj.b [Coordenacao dos Programas de Pos-graduacao de Engenharia (COPPE/UFRJ), Rio de Janeiro, RJ (Brazil). Lab. de Instrumentacao Nuclear
2009-07-01
Proton therapy applications deal with relatively thick targets like the human head or the trunk. Therefore, relatively small differences in the total proton stopping power given, for example, by the different models provided by GEANT4, could lead to significant disagreement in the final proton energy spectra when integrated along lengthy proton trajectories. This work presents a comparison of proton energy spectra for 49.1 MeV protons passing through a couple of Au absorbers with different thicknesses obtained by GEANT4.8.2 simulations using ICRU49, Ziegler1985 and Ziegler2000 models. The comparison was made with the experimental data of Tschalaer, with TRIM/SRIM2008 and MCNPX2.4.0 simulations, and the Payne analytical solution for the transport equation in the Fokker-Plank approximation. It is shown that the simulations reproduce the experimental spectra with some detectable contradictions. It should be noted that all the spectra lay at the proton energies significantly above 2 MeV, i.e. in the so-called 'Bethe-Bloch region'. Therefore the observed disagreements in GEANT4 results, simulated with different models, are somewhat unexpected. Further studies for a better understanding and to obtain definitive conclusions are necessary. (author)
Schmid, S.; Landry, G.; Thieke, C.; Verhaegen, F.; Ganswindt, U.; Belka, C.; Parodi, K.; Dedes, G.
2015-12-01
Proton range verification based on prompt gamma imaging is increasingly considered in proton therapy. Tissue heterogeneity normal to the beam direction or near the end of range may considerably degrade the ability of prompt gamma imaging to detect proton range shifts. The goal of this study was to systematically investigate the accuracy and precision of range detection from prompt gamma emission profiles for various fractions for intensity modulated proton therapy of prostate cancer, using a comprehensive clinical dataset of 15 different CT scans for 5 patients. Monte Carlo simulations using Geant4 were performed to generate spot-by-spot dose distributions and prompt gamma emission profiles for prostate treatment plans. The prompt gammas were scored at their point of emission. Three CT scans of the same patient were used to evaluate the impact of inter-fractional changes on proton range. The range shifts deduced from the comparison of prompt gamma emission profiles in the planning CT and subsequent CTs were then correlated to the corresponding range shifts deduced from the dose distributions for individual pencil beams. The distributions of range shift differences between prompt gamma and dose were evaluated in terms of precision (defined as half the 95% inter-percentile range IPR) and accuracy (median). In total about 1700 individual proton pencil beams were investigated. The IPR of the relative range shift differences between the dose profiles and the prompt gamma profiles varied between ±1.4 mm and ±2.9 mm when using the more robust profile shifting analysis. The median was found smaller than 1 mm. Methods to identify and reject unreliable spots for range verification due to range mixing were derived and resulted in an average 10% spot rejection, clearly improving the prompt gamma-dose correlation. This work supports that prompt gamma imaging can offer a reliable indicator of range changes due to anatomical variations and tissue heterogeneity
Impact of the material composition on proton range variation - A Monte Carlo study
Wu, S. W.; Tung, C. J.; Lee, C. C.; Fan, K. H.; Huang, H. C.; Chao, T. C.
2015-11-01
In this study, we used the Geant4 toolkit to demonstrate the impacts of the material composition of tissues on proton range variation. Bragg curves of different materials subjected to a 250 MeV mono-energy proton beam were simulated and compared. These simulated materials included adipose, heart, brain, cartilage, cortical bone and water. The results showed that there was significant proton range deviation between Bragg curves, especially for cortical bone. The R50 values for a 250 MeV proton beam were approximately 39.55 cm, 35.52 cm, 37.00 cm, 36.51 cm, 36.72 cm, 22.53 cm, and 38.52 cm in the phantoms that were composed completely of adipose, cartilage, tissue, heart, brain, cortical bone, and water, respectively. Mass density and electron density were used to scale the proton range for each material; electron density provided better range scaling. In addition, a similar comparison was performed by artificially setting all material density to 1.0 g/cm3 to evaluate the range deviation due to chemical components alone. Tissue heterogeneity effects due to density variation were more significant, and less significant for chemical composition variation unless the Z/A was very different.
Farah, J; Bonfrate, A; Donadille, L; Dubourg, N; Lacoste, V; Martinetti, F; Sayah, R; Trompier, F; Clairand, I [IRSN - Institute for Radiological Protection and Nuclear Safety, Fontenay-aux-roses (France); Caresana, M [Politecnico di Milano, Milano (Italy); Delacroix, S; Nauraye, C [Institut Curie - Centre de Protontherapie d Orsay, Orsay (France); Herault, J [Centre Antoine Lacassagne, Nice (France); Piau, S; Vabre, I [Institut de Physique Nucleaire d Orsay, Orsay (France)
2014-06-01
Purpose: Measure stray radiation inside a passive scattering proton therapy facility, compare values to Monte Carlo (MC) simulations and identify the actual needs and challenges. Methods: Measurements and MC simulations were considered to acknowledge neutron exposure associated with 75 MeV ocular or 180 MeV intracranial passively scattered proton treatments. First, using a specifically-designed high sensitivity Bonner Sphere system, neutron spectra were measured at different positions inside the treatment rooms. Next, measurement-based mapping of neutron ambient dose equivalent was fulfilled using several TEPCs and rem-meters. Finally, photon and neutron organ doses were measured using TLDs, RPLs and PADCs set inside anthropomorphic phantoms (Rando, 1 and 5-years-old CIRS). All measurements were also simulated with MCNPX to investigate the efficiency of MC models in predicting stray neutrons considering different nuclear cross sections and models. Results: Knowledge of the neutron fluence and energy distribution inside a proton therapy room is critical for stray radiation dosimetry. However, as spectrometry unfolding is initiated using a MC guess spectrum and suffers from algorithmic limits a 20% spectrometry uncertainty is expected. H*(10) mapping with TEPCs and rem-meters showed a good agreement between the detectors. Differences within measurement uncertainty (10–15%) were observed and are inherent to the energy, fluence and directional response of each detector. For a typical ocular and intracranial treatment respectively, neutron doses outside the clinical target volume of 0.4 and 11 mGy were measured inside the Rando phantom. Photon doses were 2–10 times lower depending on organs position. High uncertainties (40%) are inherent to TLDs and PADCs measurements due to the need for neutron spectra at detector position. Finally, stray neutrons prediction with MC simulations proved to be extremely dependent on proton beam energy and the used nuclear models and
Monte Carlo (MC) simulations of beam interaction and transport in matter are increasingly considered as essential tools to support several aspects of radiation therapy. Despite the vast application of MC to photon therapy and scattered proton therapy, clinical experience in scanned ion beam therapy is still scarce. This is especially the case for ions heavier than protons, which pose additional issues like nuclear fragmentation and varying biological effectiveness. In this work, we present the evaluation of a dedicated framework which has been developed at the Heidelberg Ion Beam Therapy Center to provide automated FLUKA MC simulations of clinical patient treatments with scanned proton and carbon ion beams. Investigations on the number of transported primaries and the dimension of the geometry and scoring grids have been performed for a representative class of patient cases in order to provide recommendations on the simulation settings, showing that recommendations derived from the experience in proton therapy cannot be directly translated to the case of carbon ion beams. The MC results with the optimized settings have been compared to the calculations of the analytical treatment planning system (TPS), showing that regardless of the consistency of the two systems (in terms of beam model in water and range calculation in different materials) relevant differences can be found in dosimetric quantities and range, especially in the case of heterogeneous and deep seated treatment sites depending on the ion beam species and energies, homogeneity of the traversed tissue and size of the treated volume. The analysis of typical TPS speed-up approximations highlighted effects which deserve accurate treatment, in contrast to adequate beam model simplifications for scanned ion beam therapy. In terms of biological dose calculations, the investigation of the mixed field components in realistic anatomical situations confirmed the findings of previous groups so far reported only in
Monte Carlo simulations of soft proton flares: testing the physics with XMM-Newton
Fioretti, Valentina; Malaguti, Giuseppe; Spiga, Daniele; Tiengo, Andrea
2016-01-01
Low energy protons (<100-300 keV) in the Van Allen belt and the outer regions can enter the field of view of X-ray focusing telescopes, interact with the Wolter-I optics, and reach the focal plane. The use of special filters protects the XMM-Newton focal plane below an altitude of 70000 km, but above this limit the effect of soft protons is still present in the form of sudden flares in the count rate of the EPIC instruments, causing the loss of large amounts of observing time. We try to characterize the input proton population and the physics interaction by simulating, using the BoGEMMS framework, the proton interaction with a simplified model of the X-ray mirror module and the focal plane, and comparing the result with a real observation. The analysis of ten orbits of observations of the EPIC/pn instrument show that the detection of flares in regions far outside the radiation belt is largely influenced by the different orientation of the Earth's magnetosphere respect with XMM-Newton's orbit, confirming th...
EPEWAX - a Monte Carlo generator for W production in electron proton scattering
Theuer, E. (RWTH Aachen, 1. Physikalisches Institut (Germany))
1992-04-01
EPEWAX is a Monte Carlo event generator intended to simulate the production of free single W bosons at an eP collider according to the process: e[sup -]P [yields] e[sup -]Wsub([yields]f[sub 1]f[sub 2]) X. It offers the possibility to completely simulate W production down to the particle level. Because the WW-Photon couplings are not fixed it allows to study non standard model W production. (orig.).
EPEWAX - a Monte Carlo generator for W production in electron proton scattering
EPEWAX is a Monte Carlo event generator intended to simulate the production of free single W bosons at an eP collider according to the process: e-P → e-Wsub(→f1f2) X. It offers the possibility to completely simulate W production down to the particle level. Because the WW-Photon couplings are not fixed it allows to study non standard model W production. (orig.)
Afanasiev, Alexandr; Vainio, Rami
2016-01-01
Context. Solar energetic particles observed in association with coronal mass ejections (CMEs) are produced by the CME-driven shock waves. The acceleration of particles is considered to be due to diffusive shock acceleration (DSA). Aims. We aim at a better understanding of DSA in the case of quasi-parallel shocks, in which self-generated turbulence in the shock vicinity plays a key role. Methods. We have developed and applied a new Monte Carlo simulation code for acceleration of protons in parallel coronal shocks. The code performs a self-consistent calculation of resonant interactions of particles with Alfv\\'en waves based on the quasi-linear theory. In contrast to the existing Monte Carlo codes of DSA, the new code features the full quasi-linear resonance condition of particle pitch-angle scattering. This allows us to take anisotropy of particle pitch-angle scattering into account, while the older codes implement an approximate resonance condition leading to isotropic scattering.We performed simulations with...
GPU based Monte Carlo for PET image reconstruction: detector modeling
Monte Carlo (MC) calculations and Graphical Processing Units (GPUs) are almost like the dedicated hardware designed for the specific task given the similarities between visible light transport and neutral particle trajectories. A GPU based MC gamma transport code has been developed for Positron Emission Tomography iterative image reconstruction calculating the projection from unknowns to data at each iteration step taking into account the full physics of the system. This paper describes the simplified scintillation detector modeling and its effect on convergence. (author)
Feasibility study of the neutron dose for real-time image-guided proton therapy: A Monte Carlo study
Kim, Jin Sung; Shin, Jung Suk; Kim, Daehyun; Shin, Eunhyuk; Chung, Kwangzoo; Cho, Sungkoo; Ahn, Sung Hwan; Ju, Sanggyu; Chung, Yoonsun; Jung, Sang Hoon; Han, Youngyih
2015-07-01
Two full rotating gantries with different nozzles (multipurpose nozzle with MLC, scanning dedicated nozzle) for a conventional cyclotron system are installed and being commissioned for various proton treatment options at Samsung Medical Center in Korea. The purpose of this study is to use Monte Carlo simulation to investigate the neutron dose equivalent per therapeutic dose, H/D, for X-ray imaging equipment under various treatment conditions. At first, we investigated the H/D for various modifications of the beamline devices (scattering, scanning, multi-leaf collimator, aperture, compensator) at the isocenter and at 20, 40 and 60 cm distances from the isocenter, and we compared our results with those of other research groups. Next, we investigated the neutron dose at the X-ray equipment used for real-time imaging under various treatment conditions. Our investigation showed doses of 0.07 ~ 0.19 mSv/Gy at the X-ray imaging equipment, depending on the treatment option and interestingly, the 50% neutron dose reduction was observed due to multileaf collimator during proton scanning treatment with the multipurpose nozzle. In future studies, we plan to measure the neutron dose experimentally and to validate the simulation data for X-ray imaging equipment for use as an additional neutron dose reduction method.
Monte Carlo Predictions of Proton SEE Cross-Sections from Heavy Ion Test Data
Xi, Kai; Zhang, Zhan-Gang; Hou, Ming-Dong; Sun, You-Mei; Luo, Jie; Liu, Tian-Qi; Wang, Bin; Ye, Bing; Yin, Ya-Nan; Liu, Jie
2015-01-01
The limits of previous methods promote us to design a new approach (named PRESTAGE) to predict proton single event effect (SEE) cross-sections using heavy-ion test data. To more realistically simulate the SEE mechanisms, we adopt Geant4 and the location-dependent strategy to describe the physics processes and the sensitivity of the device. Cross-sections predicted by PRESTAGE for over twenty devices are compared with the measured data. Evidences show that PRESTAGE can calculate not only single event upsets induced by proton indirect ionization, but also direct ionization effects and single event latch-ups. Most of the PRESTAGE calculated results agree with the experimental data within a factor of 2-3.
Elmekawy, Ahmed Farouk
The distal edge of therapeutic proton radiation beams was investigated by different methods. Proton beams produced at the Hampton University Proton Therapy Institute (HUPTI) were used to irradiate a Polymethylmethacrylate (PMMA) phantom for three different ranges (13.5, 17.0 and 21.0 cm) to investigate the distal slope dependence of the Bragg peak. The activation of 11 C was studied by scanning the phantom less than 10 minutes post-irradiation with a Philips Big Bore Gemini(c) PET/CT. The DICOM images were imported into the Varian Eclipse(c) Treatment Planning System (TPS) for analysis and then analyzed by ImageJ(c) . The distal slope ranged from ?0.1671 +/- 0.0036 to -0.1986 +/- 0.0052 (pixel intensity/slice number) for ranges 13.5 to 21.0 cm respectively. A realistic description of the setup was modeled using the GATE 7.0 Monte Carlo simulation tool and compared to the experiment data. The results show the distal slope ranged from -0.1158+/-0.0133 to -0.0787+/-0.002 (Gy/mm). Additionally, low activity, 11C were simulated to study the 11C reconstructed half-life dependence versus the initial activity for six ranges chosen around the previous activation study. The results of the expected/nominal half-life vs. activity ranged from -5 x 10-4 +/- 2.8104 x 10-4 to 1.6 x 10-3 +/- 9.44 x 10-4 (%diff./Bq). The comparison between two experiments with proton beams on a PMMA phantom and multi-layer ion chamber, and two GATE simulations of a proton beam incident on a water phantom and 11C PET study show that: (i) the distal fall-off variation of the steepness of the slopes are found to be similar thus validating the sensitivity of the PET technique to the range degradation and (ii) the average of the super-ratios difference between all studies observed is primarily due to the difference in the dose deposited in the media.
Validation of proton ionization cross section generators for Monte Carlo particle transport
Batic, Matej; Saracco, Paolo
2011-01-01
Three software systems, ERCS08, ISICS 2011 and \\v{S}mit's code, that implement theoretical calculations of inner shell ionization cross sections by proton impact, are validated with respect to experimental data. The accuracy of the cross sections they generate is quantitatively estimated and inter-compared through statistical methods. Updates and extensions of a cross section data library relevant to PIXE simulation with Geant4 are discussed.
Guan, Fada; Peeler, Christopher; Taleei, Reza; Randeniya, Sharmalee; Ge, Shuaiping; Mirkovic, Dragan; Mohan, Radhe; Titt, Uwe, E-mail: UTitt@mdanderson.org [Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030 (United States); Bronk, Lawrence [Department of Experimental Radiation Oncology, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030 (United States); Geng, Changran [Department of Nuclear Science and Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China and Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114 (United States); Grosshans, David [Department of Experimental Radiation Oncology, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030 and Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030 (United States)
2015-11-15
Purpose: The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the GEANT 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from GEANT 4 simulations to correlate with biological effectiveness of therapeutic protons. Methods: The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LET{sub t} and dose-averaged LET, LET{sub d}) using GEANT 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LET{sub t} and LET{sub d} of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose. Results: Simulation results showed that the step limit effect was small for LET{sub t} but significant for LET{sub d}. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in GEANT 4 can result in incorrect LET{sub d} calculation results in the dose plateau region for small step limits. The erroneous LET{sub d} results can be attributed to the algorithm to
{sup 18}F-FET-PET-based dose painting by numbers with protons
Rickhey, Mark; Moravek, Zdenek; Koelbl, Oliver; Bogner, Ludwig [Dept. of Radiotherapy, Univ. of Regensburg (Germany); Eilles, Christoph [Dept. of Nuclear Medicine, Univ. of Regensburg (Germany)
2010-06-15
Purpose: to investigate the potential of {sup 18}F-fluoroethyltyrosine-positron emission tomography-({sup 18}F-FET-PET-)based dose painting by numbers with protons. Material and methods: due to its high specificity to brain tumor cells, FET has a high potential to serve as a target for dose painting by numbers. Biological image-based dose painting might lead to an inhomogeneous dose prescription. For precise treatment planning of such a prescribed dose, an intensity-modulated radiotherapy (IMRT) algorithm including a Monte Carlo dose-calculation algorithm for spot-scanning protons was used. A linear tracer uptake to dose model was used to derive a dose prescription from the {sup 18}F-FET-PET. As a first investigation, a modified modulation transfer function (MTF) of protons was evaluated and compared to the MTF of photons. In a clinically adapted planning study, the feasibility of {sup 18}F-FET-PET-based dose painting with protons was demonstrated using three patients with glioblastome multiforme. The resulting dose distributions were evaluated by means of dose-difference and dose-volume histograms and compared to IMRT data. Results: the MTF for protons was constantly above that for photons. The standard deviations of the dose differences between the prescribed and the optimized dose were smaller in case of protons compared to photons. Furthermore, the escalation study showed that the doses within the subvolumes identified by biological imaging techniques could be escalated remarkably while the dose within the organs at risk was kept at a constant level. Conclusion: the presented investigation fortifies the feasibility of {sup 18}F-FET-PET-based dose painting with protons. (orig.)
Purpose: Predicted PET images on the basis of analytical filtering approach for proton range verification has been successful developed and validated using FLUKA Monte Carlo (MC) codes and phantom measurements. The purpose of the study is to validate the effectiveness of analytical filtering model for proton range verification on GATE/GEANT4 Monte Carlo simulation codes. Methods: In this study, we performed two experiments for validation of predicted β+-isotope by the analytical model with GATE/GEANT4 simulations. The first experiments to evaluate the accuracy of predicting β+-yields as a function of irradiated proton energies. In second experiment, we simulate homogeneous phantoms of different materials irradiated by a mono-energetic pencil-like proton beam. The results of filtered β+-yields distributions by the analytical model is compared with those of MC simulated β+-yields in proximal and distal fall-off ranges. Results: The results investigate the distribution between filtered β+-yields and MC simulated β+-yields distribution in different conditions. First, we found that the analytical filtering can be applied over the whole range of the therapeutic energies. Second, the range difference between filtered β+-yields and MC simulated β+-yields at the distal fall-off region are within 1.5mm for all materials used. The findings validated the usefulness of analytical filtering model on range verification of proton therapy on GATE Monte Carlo simulations. In addition, there is a larger discrepancy between filtered prediction and MC simulated β+-yields using GATE code, especially in proximal region. This discrepancy might Result from the absence of wellestablished theoretical models for predicting the nuclear interactions. Conclusion: Despite the fact that large discrepancies of the distributions between MC-simulated and predicted β+-yields were observed, the study prove the effectiveness of analytical filtering model for proton range verification using
Proton irradiation on textured bismuth based cuprate superconductors
Textured bulk polycrystalline samples of bismuth based cuprate superconductors have been subjected to irradiation with 15 MeV protons. In case of Bi-2212, there has been substantial increase in Tc, which may be due to proton induced knock-out of loosely bound oxygen. In case of (Bi,Pb)-2223, there has been a reduction in Tc. The difference in behaviour in these two systems towards proton irradiation has been explained. (author). 7 refs., 3 figs., 1 tab
Fuel-Cell Electrolytes Based on Organosilica Hybrid Proton Conductors
Narayan, Sri R.; Yen, Shiao-Pin S.
2008-01-01
A new membrane composite material that combines an organosilica proton conductor with perfluorinated Nafion material to achieve good proton conductivity and high-temperature performance for membranes used for fuel cells in stationary, transportation, and portable applications has been developed. To achieve high proton conductivities of the order of 10(exp -1)S/cm over a wide range of temperatures, a composite membrane based on a new class of mesoporous, proton-conducting, hydrogen-bonded organosilica, used with Nafion, will allow for water retention and high proton conductivity over a wider range of temperatures than currently offered by Nafion alone. At the time of this reporting, this innovation is at the concept level. Some of the materials and processes investigated have shown good proton conductivity, but membranes have not yet been prepared and demonstrated.
Proton Spin Based On Chiral Dynamics
Weber, H. J.
1999-01-01
Chiral spin fraction models agree with the proton spin data only when the chiral quark-Goldstone boson couplings are pure spinflip. For axial-vector coupling from soft-pion physics this is true for massless quarks but not for constituent quarks. Axial-vector quark-Goldstone boson couplings with {\\bf constituent} quarks are found to be inconsistent with the proton spin data.
Monte Carlo based weighting functions in neutron capture measurements
To determine neutron capture cross sections using C6D6 detectors, the Pulse Height Weighting Technique (PHWT) is mostly applied. The weighting function depends from the response function of the detection system in use. Therefore, the quality of the data depends on the detector response used for the calculation of the weighting function. An experimental determination of the response of C6D6 detectors is not always straightforward. We determined the detector response and, hence, the weighting function from Monte Carlo simulations, using the MCNP 4C2 code. To obtain reliable results a big effort was made in preparing geometry input file describing the experimental conditions. To validate the results of the Monte Carlo simulations we performed several experiments at GELINA. First, we measured the C6D6 detector response for standard -sources and for selected resonances in the 206Pb(n,). These responses were compared with the one based on Monte Carlo simulations. The good agreement between experimental and simulated data confirms the reliability of the Monte Carlo simulations. As a second validation exercise, we also determined the normalization factor in Ag and Au sample of different composition and thickness and the neutron width of the 1.15 keV resonance in 5 Fe using samples of different compositions. The result of this validation exercise was that the photon transport and the coupling of the photon and neutron transport must be accounted for in the determination of the weighting function. Accurate weighting functions are required for capture reactions in nuclei where the gamma cascade differs strongly from resonance to resonance, and are extremely important for neutron data related to reactor technologies where Pb-isotopes play an important role. The Monte Carlo based weighting function have been used to deduce the capture yield of 206Pb between 3 and 620 keV and of 232Th between 5 and 150 keV. This method will also be used for the analysis of other neutron capture
Quantitative Monte Carlo-based holmium-166 SPECT reconstruction
Purpose: Quantitative imaging of the radionuclide distribution is of increasing interest for microsphere radioembolization (RE) of liver malignancies, to aid treatment planning and dosimetry. For this purpose, holmium-166 (166Ho) microspheres have been developed, which can be visualized with a gamma camera. The objective of this work is to develop and evaluate a new reconstruction method for quantitative 166Ho SPECT, including Monte Carlo-based modeling of photon contributions from the full energy spectrum.Methods: A fast Monte Carlo (MC) simulator was developed for simulation of 166Ho projection images and incorporated in a statistical reconstruction algorithm (SPECT-fMC). Photon scatter and attenuation for all photons sampled from the full 166Ho energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator-detector response was modeled using precalculated convolution kernels. Phantom experiments were performed to quantitatively evaluate image contrast, image noise, count errors, and activity recovery coefficients (ARCs) of SPECT-fMC in comparison with those of an energy window-based method for correction of down-scattered high-energy photons (SPECT-DSW) and a previously presented hybrid method that combines MC simulation of photopeak scatter with energy window-based estimation of down-scattered high-energy contributions (SPECT-ppMC+DSW). Additionally, the impact of SPECT-fMC on whole-body recovered activities (Aest) and estimated radiation absorbed doses was evaluated using clinical SPECT data of six 166Ho RE patients.Results: At the same noise level, SPECT-fMC images showed substantially higher contrast than SPECT-DSW and SPECT-ppMC+DSW in spheres ≥17 mm in diameter. The count error was reduced from 29% (SPECT-DSW) and 25% (SPECT-ppMC+DSW) to 12% (SPECT-fMC). ARCs in five spherical volumes of 1.96–106.21 ml were improved from 32%–63% (SPECT-DSW) and 50%–80% (SPECT-ppMC+DSW) to 76%–103
Feasibility study of proton-based quality assurance of proton range compensator
Park, S.; Jeong, C.; Min, B. J.; Kwak, J.; Lee, J.; Cho, S.; Shin, D.; Lim, Y. K.; Park, S. Y.; Lee, S. B.
2013-06-01
All patient specific range compensators (RCs) are customized for achieving distal dose conformity of target volume in passively scattered proton therapy. Compensators are milled precisely using a computerized machine. In proton therapy, precision of the compensator is critical and quality assurance (QA) is required to protect normal tissues and organs from radiation damage. This study aims to evaluate the precision of proton-based quality assurance of range compensator. First, the geometry information of two compensators was extracted from the DICOM Radiotherapy (RT) plan. Next, RCs were irradiated on the EBT film individually by proton beam which is modulated to have a photon-like percent depth dose (PDD). Step phantoms were also irradiated on the EBT film to generate calibration curve which indicates relationship between optical density of irradiated film and perpendicular depth of compensator. Comparisons were made using the mean absolute difference (MAD) between coordinate information from DICOM RT and converted depth information from the EBT film. MAD over the whole region was 1.7, and 2.0 mm. However, MAD over the relatively flat regions on each compensator selected for comparison was within 1 mm. These results shows that proton-based quality assurance of range compensator is feasible and it is expected to achieve MAD over the whole region less than 1 mm with further correction about scattering effect of proton imaging.
Monte Carlo simulation of single spin asymmetries in pion-proton collisions
Bianconi, A; Bianconi, Andrea; Radici, Marco
2006-01-01
We present Monte Carlo simulations of both the Sivers and the Boer-Mulders effects in the polarized Drell-Yan $\\pi^\\pm p^\\uparrow \\to \\mu^+ \\mu^- X$ process at the center-of-mass energy $\\sqrt{s} \\sim 14$ GeV reachable at COMPASS with pion beams of energy 100 GeV. For the Sivers effect, we adopt two different parametrizations for the Sivers function to explore the statistical accuracy required to extract unambiguous information on this parton density. In particular, we verify the possibility of checking its predicted sign change between Semi-Inclusive Deep-Inelastic Scattering (SIDIS) and Drell-Yan processes, a crucial test of nonperturbative QCD. For the Boer-Mulders effect, because of the lack of parametrizations we can make only guesses. The goal is to explore the possibility of extracting information on the transversity distribution, the missing piece necessary to complete the knowledge of the nucleon spin structure at leading twist, and the Boer-Mulders function, which is related to the long-standing pro...
Sahoo, G. S.; Tripathy, S. P.; Molokanov, A. G.; Aleynikov, V. E.; Sharma, S. D.; Bandyopadhyay, T.
2016-05-01
In this work, we have used CR-39 detectors to estimate the LET (linear energy transfer) spectrum of secondary particles due to 171 MeV proton beam at different depths of water including the Bragg peak region. The measured LET spectra were compared with those obtained from FLUKA Monte Carlo simulation. The absorbed dose (DLET), dose equivalent (HLET) were estimated using the LET spectra. The values of DLET and HLET per incident proton fluence were found to increase with the increase in depth of water and were maximum at Bragg peak.
This study settles on a contribution to the elaboration of slow protons transport simulation. Atomic inner shell ionization is studied in the Plane Wave Born Approximations and in the Binary Encounter Approximation. BRINKMAN-KRAMER's theory and DMITRIEV's theory are used to study charge exchange phenomena. Protons slowing down is studied with the BRICE's stopping power, with the VAVILOV and SYMON's energy straggling distributions and with the MOLIERE, KEIL and MEYER's angular deflexion distributions. Transport simulation is made with Monte Carlo Method; K electrons motion is also taken into account
Park, M; Jung, H; Kim, G; Ji, Y; Kim, K; Park, S [Korea Institute of Radiological and Medical Sciences, Seoul, Nowon-gu (Korea, Republic of)
2014-06-01
Purpose: To estimate the three dimensional dose distributions in a polymer gel and a radiochromic gel by comparing with the virtual water phantom exposed to proton beams by applying Monte Carlo simulation. Methods: The polymer gel dosimeter is the compositeness material of gelatin, methacrylic acid, hydroquinone, tetrakis, and distilled water. The radiochromic gel is PRESAGE product. The densities of polymer and radiochromic gel were 1.040 and 1.0005 g/cm3, respectively. The shape of water phantom was a hexahedron with the size of 13 × 13 × 15 cm3. The proton beam energies of 72 and 116 MeV were used in the simulation. Proton beam was directed to the top of the phantom with Z-axis and the shape of beam was quadrangle with 10 × 10 cm2 dimension. The Percent depth dose and the dose distribution were evaluated for estimating the dose distribution of proton particle in two gel dosimeters, and compared with the virtual water phantom. Results: The Bragg-peak for proton particles in two gel dosimeters was similar to the virtual water phantom. Bragg-peak regions of polymer gel, radiochromic gel, and virtual water phantom were represented in the identical region (4.3 cm) for 72 MeV proton beam. For 116 MeV proton beam, the Bragg-peak regions of polymer gel, radiochromic gel, and virtual water phantom were represented in 9.9, 9.9 and 9.7 cm, respectively. The dose distribution of proton particles in polymer gel, radiochromic gel, and virtual water phantom was approximately identical in the case of 72 and 116 MeV energies. The errors for the simulation were under 10%. Conclusion: This work indicates the evaluation of three dimensional dose distributions by exposing proton particles to polymer and radiochromic gel dosimeter by comparing with the water phantom. The polymer gel and the radiochromic gel dosimeter show similar dose distributions for the proton beams.
Adaptation of GEANT4 to Monte Carlo dose calculations based on CT data
The GEANT4 Monte Carlo code provides many powerful functions for conducting particle transport simulations with great reliability and flexibility. However, as a general purpose Monte Carlo code, not all the functions were specifically designed and fully optimized for applications in radiation therapy. One of the primary issues is the computational efficiency, which is especially critical when patient CT data have to be imported into the simulation model. In this paper we summarize the relevant aspects of the GEANT4 tracking and geometry algorithms and introduce our work on using the code to conduct dose calculations based on CT data. The emphasis is focused on modifications of the GEANT4 source code to meet the requirements for fast dose calculations. The major features include a quick voxel search algorithm, fast volume optimization, and the dynamic assignment of material density. These features are ready to be used for tracking the primary types of particles employed in radiation therapy such as photons, electrons, and heavy charged particles. Re-calculation of a proton therapy treatment plan generated by a commercial treatment planning program for a paranasal sinus case is presented as an example
Clinical dosimetry in photon radiotherapy. A Monte Carlo based investigation
Practical clinical dosimetry is a fundamental step within the radiation therapy process and aims at quantifying the absorbed radiation dose within a 1-2% uncertainty. To achieve this level of accuracy, corrections are needed for calibrated and air-filled ionization chambers, which are used for dose measurement. The procedures of correction are based on cavity theory of Spencer-Attix and are defined in current dosimetry protocols. Energy dependent corrections for deviations from calibration beams account for changed ionization chamber response in the treatment beam. The corrections applied are usually based on semi-analytical models or measurements and are generally hard to determine due to their magnitude of only a few percents or even less. Furthermore the corrections are defined for fixed geometrical reference-conditions and do not apply to non-reference conditions in modern radiotherapy applications. The stochastic Monte Carlo method for the simulation of radiation transport is becoming a valuable tool in the field of Medical Physics. As a suitable tool for calculation of these corrections with high accuracy the simulations enable the investigation of ionization chambers under various conditions. The aim of this work is the consistent investigation of ionization chamber dosimetry in photon radiation therapy with the use of Monte Carlo methods. Nowadays Monte Carlo systems exist, which enable the accurate calculation of ionization chamber response in principle. Still, their bare use for studies of this type is limited due to the long calculation times needed for a meaningful result with a small statistical uncertainty, inherent to every result of a Monte Carlo simulation. Besides heavy use of computer hardware, techniques methods of variance reduction to reduce the needed calculation time can be applied. Methods for increasing the efficiency in the results of simulation were developed and incorporated in a modern and established Monte Carlo simulation environment
Purpose: To monitor the secondary dose distribution originating from a water phantom during proton therapy of prostate cancer using scintillating fibers. Methods: The Geant4 Monte Carlo toolkit version 9.6.p02 was used to simulate prostate cancer proton therapy based treatments. Two cases were studied. In the first case, 8 × 8 = 64 equally spaced fibers inside three 4 × 4 × 2.54 cmm3 DuPont™ Delrin blocks were used to monitor the emission of secondary particles in the transverse (left and right) and distal regions relative to the beam direction. In the second case, a scintillating block with a thickness of 2.54 cm and equal vertical and longitudinal dimensions as the water phantom was used. Geometrical cuts were used to extract the energy deposited in each fiber and the scintillating block. Results: The transverse dose distributions from secondary particles in both cases agree within <5% and with a very good symmetry. The energy deposited not only gradually increases as one moves from the peripheral row fibers towards the center of the block (aligned with the center of the prostate) but also decreases as one goes from the frontal to distal region of the block. The ratio of the doses from the prostate to the ones in the middle two rows of fibers showed a linear relationship with a slope (−3.55±2.26) × 10−5 MeV per treatment Gy. The distal detectors recorded a very small energy deposited due to water attenuation. Conclusion: With a good calibration and the ability to define a good correlation between the dose to the external fibers and the prostate, such fibers can be used for real time dose verification to the target
An 8-GeV Synchrotron-Based Proton Driver
In January 2002, the Fermilab Director initiated a design study for a high average power, modest energy proton facility. Such a facility is a possible candidate for a construction project in the U.S. starting in the middle of this decade. The key technical element is a new machine, dubbed the ''Proton Driver,'' as a replacement of the present Booster. The study of an 8-GeV synchrotron-based proton driver has been completed and published. This paper will give a summary report, including machine layout and performance, optics, beam dynamics issues, technical systems design, civil construction, cost estimate and schedule
Shin, Wook-Geun; Shin, Jae-Ik; Jeong, Jong Hwi; Lee, Se Byeong
2015-01-01
For the in vivo range verification in proton therapy, it has been tried to measure the spatial distribution of the prompt gammas generated by the proton-induced interactions with the close relationship with the proton dose distribution. However, the high energy of the prompt gammas and background gammas are still problematic in measuring the distribution. In this study, we suggested a new method determining the in vivo range by utilizing the time structure of the prompt gammas formed with the rotation of a range modulation wheel (RMW) in the passive scattering proton therapy. To validate the Monte Carlo code simulating the proton beam nozzle, axial percent depth doses (PDDs) were compared with the measured PDDs with the varying beam range of 4.73-24.01 cm. And the relationship between the proton dose rate and the time structure of the prompt gammas was assessed and compared in the water phantom. The results of the PDD showed accurate agreement within the relative errors of 1.1% in the distal range and 2.9% in...
Bznuni, S A; Zhamkochyan, V M; Polanski, A; Sosnin, A N; Khudaverdyan, A H
2001-01-01
Parameters of a subcritical cascade reactor driven by a proton accelerator and based on a primary lead-bismuth target, main reactor constructed analogously to the molten salt breeder (MSBR) reactor core and a booster-reactor analogous to the core of the BN-350 liquid metal cooled fast breeder reactor (LMFBR). It is shown by means of Monte-Carlo modeling that the reactor under study provides safe operation modes (k_{eff}=0.94-0.98), is apable to transmute effectively radioactive nuclear waste and reduces by an order of magnitude the requirements on the accelerator beam current. Calculations show that the maximal neutron flux in the thermal zone is 10^{14} cm^{12}\\cdot s^_{-1}, in the fast booster zone is 5.12\\cdot10^{15} cm^{12}\\cdot s{-1} at k_{eff}=0.98 and proton beam current I=2.1 mA.
Quantitative Monte Carlo-based holmium-166 SPECT reconstruction
Elschot, Mattijs; Smits, Maarten L. J.; Nijsen, Johannes F. W.; Lam, Marnix G. E. H.; Zonnenberg, Bernard A.; Bosch, Maurice A. A. J. van den; Jong, Hugo W. A. M. de [Department of Radiology and Nuclear Medicine, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht (Netherlands); Viergever, Max A. [Image Sciences Institute, University Medical Center Utrecht, Heidelberglaan 100, 3584 CX Utrecht (Netherlands)
2013-11-15
Purpose: Quantitative imaging of the radionuclide distribution is of increasing interest for microsphere radioembolization (RE) of liver malignancies, to aid treatment planning and dosimetry. For this purpose, holmium-166 ({sup 166}Ho) microspheres have been developed, which can be visualized with a gamma camera. The objective of this work is to develop and evaluate a new reconstruction method for quantitative {sup 166}Ho SPECT, including Monte Carlo-based modeling of photon contributions from the full energy spectrum.Methods: A fast Monte Carlo (MC) simulator was developed for simulation of {sup 166}Ho projection images and incorporated in a statistical reconstruction algorithm (SPECT-fMC). Photon scatter and attenuation for all photons sampled from the full {sup 166}Ho energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator-detector response was modeled using precalculated convolution kernels. Phantom experiments were performed to quantitatively evaluate image contrast, image noise, count errors, and activity recovery coefficients (ARCs) of SPECT-fMC in comparison with those of an energy window-based method for correction of down-scattered high-energy photons (SPECT-DSW) and a previously presented hybrid method that combines MC simulation of photopeak scatter with energy window-based estimation of down-scattered high-energy contributions (SPECT-ppMC+DSW). Additionally, the impact of SPECT-fMC on whole-body recovered activities (A{sup est}) and estimated radiation absorbed doses was evaluated using clinical SPECT data of six {sup 166}Ho RE patients.Results: At the same noise level, SPECT-fMC images showed substantially higher contrast than SPECT-DSW and SPECT-ppMC+DSW in spheres ≥17 mm in diameter. The count error was reduced from 29% (SPECT-DSW) and 25% (SPECT-ppMC+DSW) to 12% (SPECT-fMC). ARCs in five spherical volumes of 1.96–106.21 ml were improved from 32%–63% (SPECT-DSW) and 50%–80
Track structure Monte Carlo simulations of ionising radiation in water are often used to estimate radiation damage to DNA. For this purpose, an accurate simulation of the transport of densely ionising low-energy secondary electrons is particularly important, but is impaired by a high uncertainty of the required physical interaction cross section data of liquid water. A possible tool for the verification of the secondary electron transport in a track structure simulation has been suggested by Toburen et al. (2010), who have measured the angle-dependent energy spectra of electrons, emitted from a thin layer of amorphous solid water (ASW) upon a passage of 6 MeV protons. In this work, simulations were performed for the setup of their experiment, using the PTB Track structure code (PTra) and Geant4-DNA. To enable electron transport below the ionisation threshold, additional excitation and dissociative attachment anion states were included in PTra and activated in Geant4. Additionally, a surface potential was considered in both simulations, such that the escape probability for an electron is dependent on its energy and impact angle at the ASW/vacuum interface. For vanishing surface potential, the simulated spectra are in good agreement with the measured spectra for energies above 50 eV. Below, the simulations overestimate the yield of electrons by a factor up to 4 (PTra) or 7 (Geant4-DNA), which is still a better agreement than obtained in previous simulations of this experimental situation. The agreement of the simulations with experimental data was significantly improved by using a step-like increase of the potential energy at the ASW surface. - Highlights: ► Benchmarked electron transport in track structure simulations using liquid water. ► Simulated differential electron spectra agree with measured data. ► The agreement was improved by including a 3 eV surface potential step.
Bauer, J; Unholtz, D; Kurz, C; Parodi, K
2013-08-01
We report on the experimental campaign carried out at the Heidelberg Ion-Beam Therapy Center (HIT) to optimize the Monte Carlo (MC) modelling of proton-induced positron-emitter production. The presented experimental strategy constitutes a pragmatic inverse approach to overcome the known uncertainties in the modelling of positron-emitter production due to the lack of reliable cross-section data for the relevant therapeutic energy range. This work is motivated by the clinical implementation of offline PET/CT-based treatment verification at our facility. Here, the irradiation induced tissue activation in the patient is monitored shortly after the treatment delivery by means of a commercial PET/CT scanner and compared to a MC simulated activity expectation, derived under the assumption of a correct treatment delivery. At HIT, the MC particle transport and interaction code FLUKA is used for the simulation of the expected positron-emitter yield. For this particular application, the code is coupled to externally provided cross-section data of several proton-induced reactions. Studying experimentally the positron-emitting radionuclide yield in homogeneous phantoms provides access to the fundamental production channels. Therefore, five different materials have been irradiated by monoenergetic proton pencil beams at various energies and the induced β(+) activity subsequently acquired with a commercial full-ring PET/CT scanner. With the analysis of dynamically reconstructed PET images, we are able to determine separately the spatial distribution of different radionuclide concentrations at the starting time of the PET scan. The laterally integrated radionuclide yields in depth are used to tune the input cross-section data such that the impact of both the physical production and the imaging process on the various positron-emitter yields is reproduced. The resulting cross-section data sets allow to model the absolute level of measured β(+) activity induced in the investigated
In this work, we demonstrate inconsistencies in commonly used Monte Carlo methods of scoring linear energy transfer (LET) in proton therapy beams. In particle therapy beams, the LET is an important parameter because the relative biological effectiveness (RBE) depends on it. LET is often determined using Monte Carlo techniques. We used a realistic Monte Carlo model of a proton therapy nozzle to score proton LET in spread-out Bragg peak (SOBP) depth-dose distributions. We used three different scoring and calculation techniques to determine average LET at varying depths within a 140 MeV beam with a 4 cm SOBP and a 250 MeV beam with a 10 cm SOBP. These techniques included fluence-weighted (Φ-LET) and dose-weighted average (D-LET) LET calculations from: 1) scored energy spectra converted to LET spectra through a lookup table, 2) directly scored LET spectra and 3) accumulated LET scored ‘on-the-fly’ during simulations. All protons (primary and secondary) were included in the scoring. Φ-LET was found to be less sensitive to changes in scoring technique than D-LET. In addition, the spectral scoring methods were sensitive to low-energy (high-LET) cutoff values in the averaging. Using cutoff parameters chosen carefully for consistency between techniques, we found variations in Φ-LET values of up to 1.6% and variations in D-LET values of up to 11.2% for the same irradiation conditions, depending on the method used to score LET. Variations were largest near the end of the SOBP, where the LET and energy spectra are broader
GPU-Monte Carlo based fast IMRT plan optimization
Yongbao Li
2014-03-01
Full Text Available Purpose: Intensity-modulated radiation treatment (IMRT plan optimization needs pre-calculated beamlet dose distribution. Pencil-beam or superposition/convolution type algorithms are typically used because of high computation speed. However, inaccurate beamlet dose distributions, particularly in cases with high levels of inhomogeneity, may mislead optimization, hindering the resulting plan quality. It is desire to use Monte Carlo (MC methods for beamlet dose calculations. Yet, the long computational time from repeated dose calculations for a number of beamlets prevents this application. It is our objective to integrate a GPU-based MC dose engine in lung IMRT optimization using a novel two-steps workflow.Methods: A GPU-based MC code gDPM is used. Each particle is tagged with an index of a beamlet where the source particle is from. Deposit dose are stored separately for beamlets based on the index. Due to limited GPU memory size, a pyramid space is allocated for each beamlet, and dose outside the space is neglected. A two-steps optimization workflow is proposed for fast MC-based optimization. At first step, a rough dose calculation is conducted with only a few number of particle per beamlet. Plan optimization is followed to get an approximated fluence map. In the second step, more accurate beamlet doses are calculated, where sampled number of particles for a beamlet is proportional to the intensity determined previously. A second-round optimization is conducted, yielding the final result.Results: For a lung case with 5317 beamlets, 105 particles per beamlet in the first round, and 108 particles per beam in the second round are enough to get a good plan quality. The total simulation time is 96.4 sec.Conclusion: A fast GPU-based MC dose calculation method along with a novel two-step optimization workflow are developed. The high efficiency allows the use of MC for IMRT optimizations.--------------------------------Cite this article as: Li Y, Tian Z
A CNS calculation line based on a Monte Carlo method
Full text: The design of the moderator cell of a Cold Neutron Source (CNS) involves many different considerations regarding geometry, location, and materials. Decisions taken in this sense affect not only the neutron flux in the source neighborhood, which can be evaluated by a standard empirical method, but also the neutron flux values in experimental positions far away of the neutron source. At long distances from the neutron source, very time consuming 3D deterministic methods or Monte Carlo transport methods are necessary in order to get accurate figures. Standard and typical terminology such as average neutron flux, neutron current, angular flux, luminosity, are magnitudes very difficult to evaluate in positions located several meters away from the neutron source. The Monte Carlo method is a unique and powerful tool to transport neutrons. Its use in a bootstrap scheme appears to be an appropriate solution for this type of systems. The proper use of MCNP as the main tool leads to a fast and reliable method to perform calculations in a relatively short time with low statistical errors. The design goal is to evaluate the performance of the neutron sources, their beam tubes and neutron guides at specific experimental locations in the reactor hall as well as in the neutron or experimental hall. In this work, the calculation methodology used to design Cold, Thermal and Hot Neutron Sources and their associated Neutron Beam Transport Systems, based on the use of the MCNP code, is presented. This work also presents some changes made to the cross section libraries in order to cope with cryogenic moderators such as liquid hydrogen and liquid deuterium. (author)
GPU based Monte Carlo for PET image reconstruction: parameter optimization
This paper presents the optimization of a fully Monte Carlo (MC) based iterative image reconstruction of Positron Emission Tomography (PET) measurements. With our MC re- construction method all the physical effects in a PET system are taken into account thus superior image quality is achieved in exchange for increased computational effort. The method is feasible because we utilize the enormous processing power of Graphical Processing Units (GPUs) to solve the inherently parallel problem of photon transport. The MC approach regards the simulated positron decays as samples in mathematical sums required in the iterative reconstruction algorithm, so to complement the fast architecture, our work of optimization focuses on the number of simulated positron decays required to obtain sufficient image quality. We have achieved significant results in determining the optimal number of samples for arbitrary measurement data, this allows to achieve the best image quality with the least possible computational effort. Based on this research recommendations can be given for effective partitioning of computational effort into the iterations in limited time reconstructions. (author)
Proton radiography and proton computed tomography based on time-resolved dose measurements.
Testa, Mauro; Verburg, Joost M; Rose, Mark; Min, Chul Hee; Tang, Shikui; Bentefour, El Hassane; Paganetti, Harald; Lu, Hsiao-Ming
2013-11-21
We present a proof of principle study of proton radiography and proton computed tomography (pCT) based on time-resolved dose measurements. We used a prototype, two-dimensional, diode-array detector capable of fast dose rate measurements, to acquire proton radiographic images expressed directly in water equivalent path length (WEPL). The technique is based on the time dependence of the dose distribution delivered by a proton beam traversing a range modulator wheel in passive scattering proton therapy systems. The dose rate produced in the medium by such a system is periodic and has a unique pattern in time at each point along the beam path and thus encodes the WEPL. By measuring the time dose pattern at the point of interest, the WEPL to this point can be decoded. If one measures the time–dose patterns at points on a plane behind the patient for a beam with sufficient energy to penetrate the patient, the obtained 2D distribution of the WEPL forms an image. The technique requires only a 2D dosimeter array and it uses only the clinical beam for a fraction of second with negligible dose to patient. We first evaluated the accuracy of the technique in determining the WEPL for static phantoms aiming at beam range verification of the brain fields of medulloblastoma patients. Accurate beam ranges for these fields can significantly reduce the dose to the cranial skin of the patient and thus the risk of permanent alopecia. Second, we investigated the potential features of the technique for real-time imaging of a moving phantom. Real-time tumor tracking by proton radiography could provide more accurate validations of tumor motion models due to the more sensitive dependence of proton beam on tissue density compared to x-rays. Our radiographic technique is rapid (~100 ms) and simultaneous over the whole field, it can image mobile tumors without the problem of interplay effect inherently challenging for methods based on pencil beams. Third, we present the reconstructed p
Proton radiography and proton computed tomography based on time-resolved dose measurements
Testa, Mauro; Verburg, Joost M.; Rose, Mark; Min, Chul Hee; Tang, Shikui; Hassane Bentefour, El; Paganetti, Harald; Lu, Hsiao-Ming
2013-11-01
We present a proof of principle study of proton radiography and proton computed tomography (pCT) based on time-resolved dose measurements. We used a prototype, two-dimensional, diode-array detector capable of fast dose rate measurements, to acquire proton radiographic images expressed directly in water equivalent path length (WEPL). The technique is based on the time dependence of the dose distribution delivered by a proton beam traversing a range modulator wheel in passive scattering proton therapy systems. The dose rate produced in the medium by such a system is periodic and has a unique pattern in time at each point along the beam path and thus encodes the WEPL. By measuring the time dose pattern at the point of interest, the WEPL to this point can be decoded. If one measures the time-dose patterns at points on a plane behind the patient for a beam with sufficient energy to penetrate the patient, the obtained 2D distribution of the WEPL forms an image. The technique requires only a 2D dosimeter array and it uses only the clinical beam for a fraction of second with negligible dose to patient. We first evaluated the accuracy of the technique in determining the WEPL for static phantoms aiming at beam range verification of the brain fields of medulloblastoma patients. Accurate beam ranges for these fields can significantly reduce the dose to the cranial skin of the patient and thus the risk of permanent alopecia. Second, we investigated the potential features of the technique for real-time imaging of a moving phantom. Real-time tumor tracking by proton radiography could provide more accurate validations of tumor motion models due to the more sensitive dependence of proton beam on tissue density compared to x-rays. Our radiographic technique is rapid (˜100 ms) and simultaneous over the whole field, it can image mobile tumors without the problem of interplay effect inherently challenging for methods based on pencil beams. Third, we present the reconstructed p
Application of equivalence methods on Monte Carlo method based homogenization multi-group constants
The multi-group constants generated via continuous energy Monte Carlo method do not satisfy the equivalence between reference calculation and diffusion calculation applied in reactor core analysis. To the satisfaction of the equivalence theory, general equivalence theory (GET) and super homogenization method (SPH) were applied to the Monte Carlo method based group constants, and a simplified reactor core and C5G7 benchmark were examined with the Monte Carlo constants. The results show that the calculating precision of group constants is improved, and GET and SPH are good candidates for the equivalence treatment of Monte Carlo homogenization. (authors)
Matsumoto, Shinnosuke; Koba, Yusuke; Kohno, Ryosuke; Lee, Choonsik; Bolch, Wesley E; Kai, Michiaki
2016-04-01
Proton therapy has the physical advantage of a Bragg peak that can provide a better dose distribution than conventional x-ray therapy. However, radiation exposure of normal tissues cannot be ignored because it is likely to increase the risk of secondary cancer. Evaluating secondary neutrons generated by the interaction of the proton beam with the treatment beam-line structure is necessary; thus, performing the optimization of radiation protection in proton therapy is required. In this research, the organ dose and energy spectrum were calculated from secondary neutrons using Monte Carlo simulations. The Monte Carlo code known as the Particle and Heavy Ion Transport code System (PHITS) was used to simulate the transport proton and its interaction with the treatment beam-line structure that modeled the double scattering body of the treatment nozzle at the National Cancer Center Hospital East. The doses of the organs in a hybrid computational phantom simulating a 5-y-old boy were calculated. In general, secondary neutron doses were found to decrease with increasing distance to the treatment field. Secondary neutron energy spectra were characterized by incident neutrons with three energy peaks: 1×10, 1, and 100 MeV. A block collimator and a patient collimator contributed significantly to organ doses. In particular, the secondary neutrons from the patient collimator were 30 times higher than those from the first scatter. These results suggested that proactive protection will be required in the design of the treatment beam-line structures and that organ doses from secondary neutrons may be able to be reduced. PMID:26910030
Cell death following BNCT: A theoretical approach based on Monte Carlo simulations
In parallel to boron measurements and animal studies, investigations on radiation-induced cell death are also in progress in Pavia, with the aim of better characterisation of the effects of a BNCT treatment down to the cellular level. Such studies are being carried out not only experimentally but also theoretically, based on a mechanistic model and a Monte Carlo code. Such model assumes that: (1) only clustered DNA strand breaks can lead to chromosome aberrations; (2) only chromosome fragments within a certain threshold distance can undergo misrejoining; (3) the so-called 'lethal aberrations' (dicentrics, rings and large deletions) lead to cell death. After applying the model to normal cells exposed to monochromatic fields of different radiation types, the irradiation section of the code was purposely extended to mimic the cell exposure to a mixed radiation field produced by the 10B(n,α) 7Li reaction, which gives rise to alpha particles and Li ions of short range and high biological effectiveness, and by the 14N(n,p)14C reaction, which produces 0.58 MeV protons. Very good agreement between model predictions and literature data was found for human and animal cells exposed to X- or gamma-rays, protons and alpha particles, thus allowing to validate the model for cell death induced by monochromatic radiation fields. The model predictions showed good agreement also with experimental data obtained by our group exposing DHD cells to thermal neutrons in the TRIGA Mark II reactor of University of Pavia; this allowed to validate the model also for a BNCT exposure scenario, providing a useful predictive tool to bridge the gap between irradiation and cell death.
Cell death following BNCT: A theoretical approach based on Monte Carlo simulations
Ballarini, F., E-mail: francesca.ballarini@pv.infn.it [University of Pavia, Department of Nuclear and Theoretical Physics, via Bassi 6, Pavia (Italy)] [INFN (National Institute of Nuclear Physics)-Sezione di Pavia, via Bassi 6, Pavia (Italy); Bakeine, J. [University of Pavia, Department of Nuclear and Theoretical Physics, via Bassi 6, Pavia (Italy); Bortolussi, S. [University of Pavia, Department of Nuclear and Theoretical Physics, via Bassi 6, Pavia (Italy)] [INFN (National Institute of Nuclear Physics)-Sezione di Pavia, via Bassi 6, Pavia (Italy); Bruschi, P. [University of Pavia, Department of Nuclear and Theoretical Physics, via Bassi 6, Pavia (Italy); Cansolino, L.; Clerici, A.M.; Ferrari, C. [University of Pavia, Department of Surgery, Experimental Surgery Laboratory, Pavia (Italy); Protti, N.; Stella, S. [University of Pavia, Department of Nuclear and Theoretical Physics, via Bassi 6, Pavia (Italy)] [INFN (National Institute of Nuclear Physics)-Sezione di Pavia, via Bassi 6, Pavia (Italy); Zonta, A.; Zonta, C. [University of Pavia, Department of Surgery, Experimental Surgery Laboratory, Pavia (Italy); Altieri, S. [University of Pavia, Department of Nuclear and Theoretical Physics, via Bassi 6, Pavia (Italy)] [INFN (National Institute of Nuclear Physics)-Sezione di Pavia, via Bassi 6, Pavia (Italy)
2011-12-15
In parallel to boron measurements and animal studies, investigations on radiation-induced cell death are also in progress in Pavia, with the aim of better characterisation of the effects of a BNCT treatment down to the cellular level. Such studies are being carried out not only experimentally but also theoretically, based on a mechanistic model and a Monte Carlo code. Such model assumes that: (1) only clustered DNA strand breaks can lead to chromosome aberrations; (2) only chromosome fragments within a certain threshold distance can undergo misrejoining; (3) the so-called 'lethal aberrations' (dicentrics, rings and large deletions) lead to cell death. After applying the model to normal cells exposed to monochromatic fields of different radiation types, the irradiation section of the code was purposely extended to mimic the cell exposure to a mixed radiation field produced by the {sup 10}B(n,{alpha}) {sup 7}Li reaction, which gives rise to alpha particles and Li ions of short range and high biological effectiveness, and by the {sup 14}N(n,p){sup 14}C reaction, which produces 0.58 MeV protons. Very good agreement between model predictions and literature data was found for human and animal cells exposed to X- or gamma-rays, protons and alpha particles, thus allowing to validate the model for cell death induced by monochromatic radiation fields. The model predictions showed good agreement also with experimental data obtained by our group exposing DHD cells to thermal neutrons in the TRIGA Mark II reactor of University of Pavia; this allowed to validate the model also for a BNCT exposure scenario, providing a useful predictive tool to bridge the gap between irradiation and cell death.
Cell death following BNCT: a theoretical approach based on Monte Carlo simulations.
Ballarini, F; Bakeine, J; Bortolussi, S; Bruschi, P; Cansolino, L; Clerici, A M; Ferrari, C; Protti, N; Stella, S; Zonta, A; Zonta, C; Altieri, S
2011-12-01
In parallel to boron measurements and animal studies, investigations on radiation-induced cell death are also in progress in Pavia, with the aim of better characterisation of the effects of a BNCT treatment down to the cellular level. Such studies are being carried out not only experimentally but also theoretically, based on a mechanistic model and a Monte Carlo code. Such model assumes that: (1) only clustered DNA strand breaks can lead to chromosome aberrations; (2) only chromosome fragments within a certain threshold distance can undergo misrejoining; (3) the so-called "lethal aberrations" (dicentrics, rings and large deletions) lead to cell death. After applying the model to normal cells exposed to monochromatic fields of different radiation types, the irradiation section of the code was purposely extended to mimic the cell exposure to a mixed radiation field produced by the (10)B(n,α) (7)Li reaction, which gives rise to alpha particles and Li ions of short range and high biological effectiveness, and by the (14)N(n,p)(14)C reaction, which produces 0.58 MeV protons. Very good agreement between model predictions and literature data was found for human and animal cells exposed to X- or gamma-rays, protons and alpha particles, thus allowing to validate the model for cell death induced by monochromatic radiation fields. The model predictions showed good agreement also with experimental data obtained by our group exposing DHD cells to thermal neutrons in the TRIGA Mark II reactor of the University of Pavia; this allowed to validate the model also for a BNCT exposure scenario, providing a useful predictive tool to bridge the gap between irradiation and cell death. PMID:21481595
Purpose: Proton therapy exhibits several advantages over photon therapy due to depth-dose distributions from proton interactions within the target material. However, uncertainties associated with protons beam range in the patient limit the advantage of proton therapy applications. To quantify beam range, positron-emitting nuclei (PEN) and prompt gamma (PG) techniques have been developed. These techniques use de-excitation photons to describe the location of the beam in the patient. To develop a detector system for implementing the PG technique for range verification applications in proton therapy, we studied the yields, energy and angular distributions of the secondary particles emitted from a PMMA phantom. Methods: Proton pencil beams of various energies incident onto a PMMA phantom with dimensions of 5 x 5 x 50 cm3 were used for simulation with the Geant4 toolkit using the standard electromagnetic packages as well as the packages based on the binary-cascade nuclear model. The emitted secondary particles are analyzed . Results: For 160 MeV incident protons, the yields of secondary neutrons and photons per 100 incident protons were ~6 and ~15 respectively. Secondary photon energy spectrum showed several energy peaks in the range between 0 and 10 MeV. The energy peaks located between 4 and 6 MeV were attributed to originate from direct proton interactions with 12C (~ 4.4 MeV) and 16O (~ 6 MeV), respectively. Most of the escaping secondary neutrons were found to have energies between 10 and 100 MeV. Isotropic emissions were found for lower energy neutrons (<10 MeV) and photons for all energies, while higher energy neutrons were emitted predominantly in the forward direction. The yields of emitted photons and neutrons increased with the increase of incident proton energies. Conclusions: A detector system is currently being developed incorporating the yields, energy and angular distributions of secondary particles from proton interactions obtained from this study
Proton irradiation of liquid crystal based adaptive optical devices
To assess its radiation hardness, a liquid crystal based adaptive optical element has been irradiated using a 60 MeV proton beam. The device with the functionality of an optical beam steerer was characterised before, during and after the irradiation. A systematic set of measurements on the transmission and beam deflection angles was carried out. The measurements showed that the transmission decreased only marginally and that its optical performance degraded only after a very high proton fluence (1010p/cm2). The device showed complete annealing in the functionality as a beam steerer, which leads to the conclusion that the liquid crystal technology for optical devices is not vulnerable to proton irradiation as expected in space.
Refined Stratified Sampling for efficient Monte Carlo based uncertainty quantification
A general adaptive approach rooted in stratified sampling (SS) is proposed for sample-based uncertainty quantification (UQ). To motivate its use in this context the space-filling, orthogonality, and projective properties of SS are compared with simple random sampling and Latin hypercube sampling (LHS). SS is demonstrated to provide attractive properties for certain classes of problems. The proposed approach, Refined Stratified Sampling (RSS), capitalizes on these properties through an adaptive process that adds samples sequentially by dividing the existing subspaces of a stratified design. RSS is proven to reduce variance compared to traditional stratified sample extension methods while providing comparable or enhanced variance reduction when compared to sample size extension methods for LHS – which do not afford the same degree of flexibility to facilitate a truly adaptive UQ process. An initial investigation of optimal stratification is presented and motivates the potential for major advances in variance reduction through optimally designed RSS. Potential paths for extension of the method to high dimension are discussed. Two examples are provided. The first involves UQ for a low dimensional function where convergence is evaluated analytically. The second presents a study to asses the response variability of a floating structure to an underwater shock. - Highlights: • An adaptive process, rooted in stratified sampling, is proposed for Monte Carlo-based uncertainty quantification. • Space-filling, orthogonality, and projective properties of stratified sampling are investigated • Stratified sampling is shown to possess attractive properties for certain classes of problems. • Refined Stratified Sampling, a new sampling method is proposed that enables the adaptive UQ process. • Optimality of RSS stratum division is explored
Clinical results of proton beam therapy for skull base chordoma
Purpose: To evaluate clinical results of proton beam therapy for patients with skull base chordoma. Methods and materials: Thirteen patients with skull base chordoma who were treated with proton beams with or without X-rays at the University of Tsukuba between 1989 and 2000 were retrospectively reviewed. A median total tumor dose of 72.0 Gy (range, 63.0-95.0 Gy) was delivered. The patients were followed for a median period of 69.3 months (range, 14.6-123.4 months). Results: The 5-year local control rate was 46.0%. Cause-specific, overall, and disease-free survival rates at 5 years were 72.2%, 66.7%, and 42.2%, respectively. The local control rate was higher, without statistical significance, for those with preoperative tumors <30 mL. Partial or subtotal tumor removal did not yield better local control rates than for patients who underwent biopsy only as the latest surgery. Conclusion: Proton beam therapy is effective for patients with skull base chordoma, especially for those with small tumors. For a patient with a tumor of <30 mL with no prior treatment, biopsy without tumor removal seems to be appropriate before proton beam therapy
A global reaction route mapping-based kinetic Monte Carlo algorithm
Mitchell, Izaac; Irle, Stephan; Page, Alister J.
2016-07-01
We propose a new on-the-fly kinetic Monte Carlo (KMC) method that is based on exhaustive potential energy surface searching carried out with the global reaction route mapping (GRRM) algorithm. Starting from any given equilibrium state, this GRRM-KMC algorithm performs a one-step GRRM search to identify all surrounding transition states. Intrinsic reaction coordinate pathways are then calculated to identify potential subsequent equilibrium states. Harmonic transition state theory is used to calculate rate constants for all potential pathways, before a standard KMC accept/reject selection is performed. The selected pathway is then used to propagate the system forward in time, which is calculated on the basis of 1st order kinetics. The GRRM-KMC algorithm is validated here in two challenging contexts: intramolecular proton transfer in malonaldehyde and surface carbon diffusion on an iron nanoparticle. We demonstrate that in both cases the GRRM-KMC method is capable of reproducing the 1st order kinetics observed during independent quantum chemical molecular dynamics simulations using the density-functional tight-binding potential.
Monte Carlo-based simulation of dynamic jaws tomotherapy
Purpose: Original TomoTherapy systems may involve a trade-off between conformity and treatment speed, the user being limited to three slice widths (1.0, 2.5, and 5.0 cm). This could be overcome by allowing the jaws to define arbitrary fields, including very small slice widths (<1 cm), which are challenging for a beam model. The aim of this work was to incorporate the dynamic jaws feature into a Monte Carlo (MC) model called TomoPen, based on the MC code PENELOPE, previously validated for the original TomoTherapy system. Methods: To keep the general structure of TomoPen and its efficiency, the simulation strategy introduces several techniques: (1) weight modifiers to account for any jaw settings using only the 5 cm phase-space file; (2) a simplified MC based model called FastStatic to compute the modifiers faster than pure MC; (3) actual simulation of dynamic jaws. Weight modifiers computed with both FastStatic and pure MC were compared. Dynamic jaws simulations were compared with the convolution/superposition (C/S) of TomoTherapy in the ''cheese'' phantom for a plan with two targets longitudinally separated by a gap of 3 cm. Optimization was performed in two modes: asymmetric jaws-constant couch speed (''running start stop,'' RSS) and symmetric jaws-variable couch speed (''symmetric running start stop,'' SRSS). Measurements with EDR2 films were also performed for RSS for the formal validation of TomoPen with dynamic jaws. Results: Weight modifiers computed with FastStatic were equivalent to pure MC within statistical uncertainties (0.5% for three standard deviations). Excellent agreement was achieved between TomoPen and C/S for both asymmetric jaw opening/constant couch speed and symmetric jaw opening/variable couch speed, with deviations well within 2%/2 mm. For RSS procedure, agreement between C/S and measurements was within 2%/2 mm for 95% of the points and 3%/3 mm for 98% of the points, where dose is greater than 30% of the prescription dose (gamma analysis
MCHITS: Monte Carlo based Method for Hyperlink Induced Topic Search on Networks
Zhaoyan Jin
2013-10-01
Full Text Available Hyperlink Induced Topic Search (HITS is the most authoritative and most widely used personalized ranking algorithm on networks. The HITS algorithm ranks nodes on networks according to power iteration, and has high complexity of computation. This paper models the HITS algorithm with the Monte Carlo method, and proposes Monte Carlo based algorithms for the HITS computation. Theoretical analysis and experiments show that the Monte Carlo based approximate computing of the HITS ranking reduces computing resources a lot while keeping higher accuracy, and is significantly better than related works
Proton interaction with an exposed object material needs to be modeled with account for three basic processes: electromagnetic stopping of protons in matter, multiple coulomb scattering and nuclear interactions. Just the last type of processes is the topic of this paper. Monte Carlo codes are often used to simulate high-energy particle interaction with matter. However, nuclear interaction models implemented in these codes are rather extensive and their use in treatment planning systems requires huge computational resources. We have selected the IThMC code for its ability to reproduce experiments which measure the distribution of the projected ranges of nuclear secondary particles generated by proton beams in a multi-layer Faraday cup. The multi-layer Faraday cup detectors measure charge rather than dose and allow distinguishing between electromagnetic and nuclear interactions. The event generator used in the IThMC code is faster, but less accurate than any other used in testing. Our model of nuclear reactions demonstrates quite good agreement with experiment in the context of their effect on the Bragg peak in therapeutic applications
Unfolding an under-determined neutron spectrum using genetic algorithm based Monte Carlo
Spallation in addition to the other photon-neutron reactions in target materials and different components in accelerators may result in production of huge amount of energetic protons which further leads to the production of neutron and contributes to the main component of the total dose. For dosimetric purposes in accelerator facilities the detector measurements doesn't provide directly the actual neutron flux values but a cumulative picture. To obtain Neutron spectrum from the measured data, response functions of the measuring instrument together with the measurements are used into many unfolding techniques which are frequently used for unfolding the hidden spectral information. Here we discuss a genetic algorithm based unfolding technique which is in the process of development. Genetic Algorithm is a stochastic method based on natural selection, which mimics Darwinian theory of survival of the best. The above said method has been tested to unfold the neutron spectra obtained from a reaction carried out at an accelerator facility, with energetic carbon ions on thick silver target along with its respective neutron response of BC501A liquid scintillation detector. The problem dealt here is under-determined where the number of measurements is less than the required energy bin information. The results so obtained were compared with those obtained using the established unfolding code FERDOR, which unfolds data for completely determined problems. It is seen that the genetic algorithm based solution has a reasonable match with the results of FERDOR, when smoothening carried out by Monte Carlo is taken into consideration. This method appears to be a promising candidate for unfolding neutron spectrum in cases of under-determined and over-determined, where measurements are more. The method also has advantages of flexibility, computational simplicity and works well without need of any initial guess spectrum. (author)
This paper reports on two experiments performed at the Synchrophasotron/Nuclotron accelerator complex at JINR. Relativistic protons with energies 885 MeV and 1.5 GeV hit a massive cylindrical lead target. The spatial and energetic distributions of the neutron field produced by the spallation reactions were measured by the activation of Al, Au, Bi, Co, and Cu foils placed on the surface of the target and close to it. The yields of the radioactive nuclei produced by threshold reactions in these foils were determined by the analyses of their γ spectra. The comparison with Monte-Carlo based simulations was performed both with the LAHET+MCNP code and the MCNPX code
An empirical formula based on Monte Carlo simulation for diffuse reflectance from turbid media
Gnanatheepam, Einstein; Aruna, Prakasa Rao; Ganesan, Singaravelu
2016-03-01
Diffuse reflectance spectroscopy has been widely used in diagnostic oncology and characterization of laser irradiated tissue. However, still accurate and simple analytical equation does not exist for estimation of diffuse reflectance from turbid media. In this work, a diffuse reflectance lookup table for a range of tissue optical properties was generated using Monte Carlo simulation. Based on the generated Monte Carlo lookup table, an empirical formula for diffuse reflectance was developed using surface fitting method. The variance between the Monte Carlo lookup table surface and the surface obtained from the proposed empirical formula is less than 1%. The proposed empirical formula may be used for modeling of diffuse reflectance from tissue.
TH-A-19A-06: Site-Specific Comparison of Analytical and Monte Carlo Based Dose Calculations
Schuemann, J; Grassberger, C; Paganetti, H [Massachusetts General Hospital and Harvard Medical School, Boston, MA (United States); Dowdell, S [Illawarra Shoalhaven Local Health District, Wollongong (Australia)
2014-06-15
Purpose: To investigate the impact of complex patient geometries on the capability of analytical dose calculation algorithms to accurately predict dose distributions and to verify currently used uncertainty margins in proton therapy. Methods: Dose distributions predicted by an analytical pencilbeam algorithm were compared with Monte Carlo simulations (MCS) using TOPAS. 79 complete patient treatment plans were investigated for 7 disease sites (liver, prostate, breast, medulloblastoma spine and whole brain, lung and head and neck). A total of 508 individual passively scattered treatment fields were analyzed for field specific properties. Comparisons based on target coverage indices (EUD, D95, D90 and D50) were performed. Range differences were estimated for the distal position of the 90% dose level (R90) and the 50% dose level (R50). Two-dimensional distal dose surfaces were calculated and the root mean square differences (RMSD), average range difference (ARD) and average distal dose degradation (ADD), the distance between the distal position of the 80% and 20% dose levels (R80- R20), were analyzed. Results: We found target coverage indices calculated by TOPAS to generally be around 1–2% lower than predicted by the analytical algorithm. Differences in R90 predicted by TOPAS and the planning system can be larger than currently applied range margins in proton therapy for small regions distal to the target volume. We estimate new site-specific range margins (R90) for analytical dose calculations considering total range uncertainties and uncertainties from dose calculation alone based on the RMSD. Our results demonstrate that a reduction of currently used uncertainty margins is feasible for liver, prostate and whole brain fields even without introducing MC dose calculations. Conclusion: Analytical dose calculation algorithms predict dose distributions within clinical limits for more homogeneous patients sites (liver, prostate, whole brain). However, we recommend
Burris-Mog, Trevor J.
The interaction of intense laser light (I > 10 18 W/cm2) with a thin target foil leads to the Target Normal Sheath Acceleration mechanism (TNSA). TNSA is responsible for the generation of high current, ultra-low emittance proton beams, which may allow for the development of a compact and cost effective proton therapy system for the treatment of cancer. Before this application can be realized, control is needed over the large divergence and the 100% kinetic energy spread that are characteristic of TNSA proton beams. The work presented here demonstrates control over the divergence and energy spread using strong magnetic fields generated by a pulse power solenoid. The solenoidal field results in a parallel proton beam with a kinetic energy spread DeltaE/E = 10%. Assuming that next generation lasers will be able to operate at 10 Hz, the 10% spread in the kinetic energy along with the 23% capture efficiency of the solenoid yield enough protons per laser pulse to, for the first time, consider applications in Radiation Oncology. Current lasers can generate proton beams with kinetic energies up to 67.5 MeV, but for therapy applications, the proton kinetic energy must reach 250 MeV. Since the maximum kinetic energy Emax of the proton scales with laser light intensity as Emax ∝ I0.5, next generation lasers may very well accelerate 250 MeV protons. As the kinetic energy of the protons is increased, the magnetic field strength of the solenoid will need to increase. The scaling of the magnetic field B with the kinetic energy of the protons follows B ∝ E1/2. Therefor, the field strength of the solenoid presented in this work will need to be increased by a factor of 2.4 in order to accommodate 250 MeV protons. This scaling factor seems reasonable, even with present technology. This work not only demonstrates control over beam divergence and energy spread, it also allows for us to now perform feasibility studies to further research what a laser-based proton therapy system
Machine learning-based patient specific prompt-gamma dose monitoring in proton therapy
Online dose monitoring in proton therapy is currently being investigated with prompt-gamma (PG) devices. PG emission was shown to be correlated with dose deposition. This relationship is mostly unknown under real conditions. We propose a machine learning approach based on simulations to create optimized treatment-specific classifiers that detect discrepancies between planned and delivered dose. Simulations were performed with the Monte-Carlo platform Gate/Geant4 for a spot-scanning proton therapy treatment and a PG camera prototype currently under investigation. The method first builds a learning set of perturbed situations corresponding to a range of patient translation. This set is then used to train a combined classifier using distal falloff and registered correlation measures. Classifier performances were evaluated using receiver operating characteristic curves and maximum associated specificity and sensitivity. A leave-one-out study showed that it is possible to detect discrepancies of 5 mm with specificity and sensitivity of 85% whereas using only distal falloff decreases the sensitivity down to 77% on the same data set. The proposed method could help to evaluate performance and to optimize the design of PG monitoring devices. It is generic: other learning sets of deviations, other measures and other types of classifiers could be studied to potentially reach better performance. At the moment, the main limitation lies in the computation time needed to perform the simulations. (paper)
Eisenhart, Thomas T; Howland, William C; Dempsey, Jillian L
2016-08-18
The proton-coupled electron transfer (PCET) oxidation of p-aminophenol in acetonitrile was initiated via stopped-flow rapid-mixing and spectroscopically monitored. For oxidation by ferrocenium in the presence of 7-(dimethylamino)quinoline proton acceptors, both the electron transfer and proton transfer components could be optically monitored in the visible region; the decay of the ferrocenium absorbance is readily monitored (λmax = 620 nm), and the absorbance of the 2,4-substituted 7-(dimethylamino)quinoline derivatives (λmax = 370-392 nm) red-shifts substantially (ca. 70 nm) upon protonation. Spectral analysis revealed the reaction proceeds via a stepwise electron transfer-proton transfer process, and modeling of the kinetics traces monitoring the ferrocenium and quinolinium signals provided rate constants for elementary proton and electron transfer steps. As the pKa values of the conjugate acids of the 2,4-R-7-(dimethylamino)quinoline derivatives employed were readily tuned by varying the substituents at the 2- and 4-positions of the quinoline backbone, the driving force for proton transfer was systematically varied. Proton transfer rate constants (kPT,2 = (1.5-7.5) × 10(8) M(-1) s(-1), kPT,4 = (0.55-3.0) × 10(7) M(-1) s(-1)) were found to correlate with the pKa of the conjugate acid of the proton acceptor, in agreement with anticipated free energy relationships for proton transfer processes in PCET reactions. PMID:27500804
Incerti, S., E-mail: sebastien.incerti@tdt.edu.vn [Division of Nuclear Physics, Ton Duc Thang University, Ho Chi Minh City (Viet Nam); Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City (Viet Nam); Univ. Bordeaux, CENBG, UMR 5797, F-33170 Gradignan (France); CNRS, IN2P3, CENBG, UMR 5797, F-33170 Gradignan (France); Barberet, Ph.; Dévès, G.; Michelet, C. [Univ. Bordeaux, CENBG, UMR 5797, F-33170 Gradignan (France); CNRS, IN2P3, CENBG, UMR 5797, F-33170 Gradignan (France); Francis, Z. [Université Saint Joseph, Science Faculty, Department of Physics, Beirut (Lebanon); Ivantchenko, V. [Ecoanalytica, Moscow (Russian Federation); Geant4 Associates International Ltd, Hebden Bridge (United Kingdom); Mantero, A. [SWHARD srl, via Greto di Cornigliano 6r, 16152 Genova (Italy); El Bitar, Z. [Institut Pluridisciplinaire Hubert Curien, CNRS/IN2P3, 67037 Strasbourg Cedex (France); Bernal, M.A. [Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, SP (Brazil); Tran, H.N. [Division of Nuclear Physics, Ton Duc Thang University, Ho Chi Minh City (Viet Nam); Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City (Viet Nam); Karamitros, M.; Seznec, H. [Univ. Bordeaux, CENBG, UMR 5797, F-33170 Gradignan (France); CNRS, IN2P3, CENBG, UMR 5797, F-33170 Gradignan (France)
2015-09-01
The general purpose Geant4 Monte Carlo simulation toolkit is able to simulate radiative and non-radiative atomic de-excitation processes such as fluorescence and Auger electron emission, occurring after interaction of incident ionising radiation with target atomic electrons. In this paper, we evaluate the Geant4 modelling capability for the simulation of fluorescence spectra induced by 1.5 MeV proton irradiation of thin high-Z foils (Fe, GdF{sub 3}, Pt, Au) with potential interest for nanotechnologies and life sciences. Simulation results are compared to measurements performed at the Centre d’Etudes Nucléaires de Bordeaux-Gradignan AIFIRA nanobeam line irradiation facility in France. Simulation and experimental conditions are described and the influence of Geant4 electromagnetic physics models is discussed.
Proton radiotherapy in management of pediatric base of skull tumors
Purpose: Primary skull base tumors of the developing child are rare and present a formidable challenge to both surgeons and radiation oncologists. Gross total resection with negative margins is rarely achieved, and the risks of functional, structural, and cosmetic deficits limit the radiation dose using conventional radiation techniques. Twenty-nine children and adolescents treated with conformal proton radiotherapy (proton RT) were analyzed to assess treatment efficacy and safety. Methods and Materials: Between July 1992 and April 1999, 29 patients with mesenchymal tumors underwent fractionated proton (13 patients) or fractionated combined proton and photon (16 patients) irradiation. The age at treatment ranged from 1 to 19 years (median 12); 14 patients were male and 15 female. Tumors were grouped as malignant or benign. Twenty patients had malignant histologic findings, including chordoma (n=10), chondrosarcoma (n=3), rhabdomyosarcoma (n=4), and other sarcomas (n=3). Target doses ranged between 50.4 and 78.6 Gy/cobalt Gray equivalent (CGE), delivered at doses of 1.8-2.0 Gy/CGE per fraction. The benign histologic findings included giant cell tumors (n=6), angiofibromas (n=2), and chondroblastoma (n=1). RT doses for this group ranged from 45.0 to 71.8 Gy/CGE. Despite maximal surgical resection, 28 (97%) of 29 patients had gross disease at the time of proton RT. Follow-up after proton RT ranged from 13 to 92 months (mean 40). Results: Of the 20 patients with malignant tumors, 5 (25%) had local failure; 1 patient had failure in the surgical access route and 3 patients developed distant metastases. Seven patients had died of progressive disease at the time of analysis. Local tumor control was maintained in 6 (60%) of 10 patients with chordoma, 3 (100%) of 3 with chondrosarcoma, 4 (100%) of 4 with rhabdomyosarcoma, and 2 (66%) of 3 with other sarcomas. The actuarial 5-year local control and overall survival rate was 72% and 56%, respectively, and the overall survival
Purpose: Recent studies have demonstrated the capability of graphics processing units (GPUs) to compute dose distributions using Monte Carlo (MC) methods within clinical time constraints. However, GPUs have a rigid vectorial architecture that favors the implementation of simplified particle transport algorithms, adapted to specific tasks. Our new, fast, and multipurpose MC code, named MCsquare, runs on Intel Xeon Phi coprocessors. This technology offers 60 independent cores, and therefore more flexibility to implement fast and yet generic MC functionalities, such as prompt gamma simulations. Methods: MCsquare implements several models and hence allows users to make their own tradeoff between speed and accuracy. A 200 MeV proton beam is simulated in a heterogeneous phantom using Geant4 and two configurations of MCsquare. The first one is the most conservative and accurate. The method of fictitious interactions handles the interfaces and secondary charged particles emitted in nuclear interactions are fully simulated. The second, faster configuration simplifies interface crossings and simulates only secondary protons after nuclear interaction events. Integral depth-dose and transversal profiles are compared to those of Geant4. Moreover, the production profile of prompt gammas is compared to PENH results. Results: Integral depth dose and transversal profiles computed by MCsquare and Geant4 are within 3%. The production of secondaries from nuclear interactions is slightly inaccurate at interfaces for the fastest configuration of MCsquare but this is unlikely to have any clinical impact. The computation time varies between 90 seconds for the most conservative settings to merely 59 seconds in the fastest configuration. Finally prompt gamma profiles are also in very good agreement with PENH results. Conclusion: Our new, fast, and multi-purpose Monte Carlo code simulates prompt gammas and calculates dose distributions in less than a minute, which complies with clinical time
Ulmer, W.; Schaffner, B.
2010-01-01
We have developed a model for proton depth dose and lateral distributions based on Monte Carlo calculations (GEANT4) and an integration procedure of the Bethe-Bloch equation (BBE). The model accounts for the transport of primary and secondary protons, the creation of recoil protons and heavy recoil nuclei as well as lateral scattering of these contributions. The buildup, which is experimentally observed in higher energy depth dose curves, is modeled by inclusion of two different origins: 1. S...
Lindsay, C; Jirasek, A [University of Victoria (Australia); Blackmore, E; Hoehr, C; Schaffer, P; Trinczek, M [TRIUMF (Canada); Sossi, V [University of British Columbia (Canada)
2014-08-15
Uveal melanoma is a rare and deadly tumour of the eye with primary metastases in the liver resulting in an 8% 2-year survival rate upon detection. Large growths, or those in close proximity to the optic nerve, pose a particular challenge to the commonly employed eye-sparing technique of eye-plaque brachytherapy. In these cases external beam charged particle therapy offers improved odds in avoiding catastrophic side effects such as neuropathy or blindness. Since 1995, the British Columbia Cancer Agency in partnership with the TRIUMF national laboratory have offered proton therapy in the treatment of difficult ocular tumors. Having seen 175 patients, yielding 80% globe preservation and 82% metastasis free survival as of 2010, this modality has proven to be highly effective. Despite this success, there have been few studies into the use of the world's largest cyclotron in patient care. Here we describe first efforts of modeling the TRIUMF dose delivery system using the FLUKA Monte Carlo package. Details on geometry, estimating beam parameters, measurement of primary dose and simulation of PET isotope production are discussed. Proton depth dose in both modulated and pristine beams is successfully simulated to sub-millimeter precision in range (within limits of measurement) and 2% agreement to measurement within in a treatment volume. With the goal of using PET signals for in vivo dosimetry (alignment), a first look at PET isotope depth distribution is presented — comparing favourably to a naive method of approximating simulated PET slice activity in a Lucite phantom.
Spot-Scanning-Based Proton Therapy for Extracranial Chordoma
Staab, Adrian, E-mail: adrian.staab@psi.ch [Center for Proton Therapy, Paul Scherrer Institute, Villigen (Switzerland); Rutz, Hans Peter; Ares, Carmen; Timmermann, Beate; Schneider, Ralf; Bolsi, Alessandra; Albertini, Francesca; Lomax, Antony; Goitein, Gudrun; Hug, Eugen [Center for Proton Therapy, Paul Scherrer Institute, Villigen (Switzerland)
2011-11-15
Purpose: To evaluate effectiveness and safety of spot-scanning-based proton-radiotherapy (PT) for extracranial chordomas (ECC). Methods and Material: Between 1999-2006, 40 patients with chordoma of C-, T-, and L-spine and sacrum were treated at Paul Scherrer Institute (PSI) with PT using spot-scanning. Median patient age was 58 years (range, 10-81 years); 63% were male, and 36% were female. Nineteen patients (47%) had gross residual disease (mean 69 cc; range, 13-495 cc) before PT, and 21 patients (53%) had undergone prior titanium-based surgical stabilization (SS) and reconstruction of the axial skeleton. Proton doses were expressed as Gy(RBE). A conversion factor of 1.1 was used to account for higher relative biological effectiveness (RBE) of protons compared with photons. Mean total dose was 72.5 Gy(RBE) [range, 59.4-75.2 Gy(RBE)] delivered at 1.8-2.0 Gy(RBE) dose per fraction. Median follow-up time was 43 months. Results: In 19 patients without surgical stabilization, actuarial local control (LC) rate at 5 years was 100%. LC for patients with gross residual disease but without surgical stabilization was also 100% at 5 years. In contrast, 12 failures occurred in 21 patients with SS, yielding a significantly decreased 5-year LC rate of 30% (p = 0.0003). For the entire cohort, 5-year LC rates were 62%, disease-free survival rates were 57%, and overall survival rates were 80%. Rates were 100% for patients without SS. No other factor, including dosimetric parameters (V95, V80) were predictive for tumor control on univariate analysis. Conclusion: Spot-scanning-based PT at PSI delivered subsequently to function-preserving surgery for tumor debulking, decompression of spinal cord, or biopsy only is safe and highly effective in patients with ECC without major surgical instrumentation even in view of large, unresectable disease.
Pentanol-based target material with polarized protons
1-pentanol is a promising material for a target with polarized protons owing to its high resistance to radiation damage. To develop the target, the solutions of 1-pentanol or 2-pentanol with complexes of pentavalent chromium ware investigated. The material based EHBA-Cr(V) solution in a glass-like matrix, consisting of 1-pentanol, 3-pentanol and 1,2-propanediol, was proposed as a target material. It was investigated by the electron paramagnetic resonance and differential scanning calorimetry methods. 24 refs.; 3 figs.; 1 tab
Excited States of Proton-bound DNA/RNA Base Homo-dimers: Pyrimidines
Féraud, Géraldine; Dedonder, Claude; Jouvet, Christophe; Pino, Gustavo A
2015-01-01
We are presenting the electronic photo fragment spectra of the protonated pyrimidine DNA bases homo-dimers. Only the thymine dimer exhibits a well structured vibrational progression, while protonated monomer shows broad vibrational bands. This shows that proton bonding can block some non radiative processes present in the monomer.
Monte-Carlo based uncertainty analysis: Sampling efficiency and sampling convergence
Monte Carlo analysis has become nearly ubiquitous since its introduction, now over 65 years ago. It is an important tool in many assessments of the reliability and robustness of systems, structures or solutions. As the deterministic core simulation can be lengthy, the computational costs of Monte Carlo can be a limiting factor. To reduce that computational expense as much as possible, sampling efficiency and convergence for Monte Carlo are investigated in this paper. The first section shows that non-collapsing space-filling sampling strategies, illustrated here with the maximin and uniform Latin hypercube designs, highly enhance the sampling efficiency, and render a desired level of accuracy of the outcomes attainable with far lesser runs. In the second section it is demonstrated that standard sampling statistics are inapplicable for Latin hypercube strategies. A sample-splitting approach is put forward, which in combination with a replicated Latin hypercube sampling allows assessing the accuracy of Monte Carlo outcomes. The assessment in turn permits halting the Monte Carlo simulation when the desired levels of accuracy are reached. Both measures form fairly noncomplex upgrades of the current state-of-the-art in Monte-Carlo based uncertainty analysis but give a substantial further progress with respect to its applicability.
Ulmer, W
2010-01-01
We have developed a model for proton depth dose and lateral distributions based on Monte Carlo calculations (GEANT4) and an integration procedure of the Bethe-Bloch equation (BBE). The model accounts for the transport of primary and secondary protons, the creation of recoil protons and heavy recoil nuclei as well as lateral scattering of these contributions. The buildup, which is experimentally observed in higher energy depth dose curves, is modeled by inclusion of two different origins: 1. Secondary reaction protons with a contribution of ca. 65 % of the buildup (for monoenergetic protons). 2. Landau tails as well as Gaussian type of fluctuations for range straggling effects. All parameters of the model for initially monoenergetic proton beams have been obtained from Monte Carlo calculations or checked by them. Furthermore, there are a few parameters, which can be obtained by fitting the model to measured depth dose curves in order to describe individual characteristics of the beamline - the most important b...
Obot, I. B.; Kaya, Savaş; Kaya, Cemal; Tüzün, Burak
2016-06-01
DFT and Monte Carlo simulation were performed on three Schiff bases namely, 4-(4-bromophenyl)-N‧-(4-methoxybenzylidene)thiazole-2-carbohydrazide (BMTC), 4-(4-bromophenyl)-N‧-(2,4-dimethoxybenzylidene)thiazole-2-carbohydrazide (BDTC), 4-(4-bromophenyl)-N‧-(4-hydroxybenzylidene)thiazole-2-carbohydrazide (BHTC) recently studied as corrosion inhibitor for steel in acid medium. Electronic parameters relevant to their inhibition activity such as EHOMO, ELUMO, Energy gap (ΔE), hardness (η), softness (σ), the absolute electronegativity (χ), proton affinity (PA) and nucleophilicity (ω) etc., were computed and discussed. Monte Carlo simulations were applied to search for the most stable configuration and adsorption energies for the interaction of the inhibitors with Fe (110) surface. The theoretical data obtained are in most cases in agreement with experimental results.
Implementation of a Monte Carlo based inverse planning model for clinical IMRT with MCNP code
He, Tongming Tony
In IMRT inverse planning, inaccurate dose calculations and limitations in optimization algorithms introduce both systematic and convergence errors to treatment plans. The goal of this work is to practically implement a Monte Carlo based inverse planning model for clinical IMRT. The intention is to minimize both types of error in inverse planning and obtain treatment plans with better clinical accuracy than non-Monte Carlo based systems. The strategy is to calculate the dose matrices of small beamlets by using a Monte Carlo based method. Optimization of beamlet intensities is followed based on the calculated dose data using an optimization algorithm that is capable of escape from local minima and prevents possible pre-mature convergence. The MCNP 4B Monte Carlo code is improved to perform fast particle transport and dose tallying in lattice cells by adopting a selective transport and tallying algorithm. Efficient dose matrix calculation for small beamlets is made possible by adopting a scheme that allows concurrent calculation of multiple beamlets of single port. A finite-sized point source (FSPS) beam model is introduced for easy and accurate beam modeling. A DVH based objective function and a parallel platform based algorithm are developed for the optimization of intensities. The calculation accuracy of improved MCNP code and FSPS beam model is validated by dose measurements in phantoms. Agreements better than 1.5% or 0.2 cm have been achieved. Applications of the implemented model to clinical cases of brain, head/neck, lung, spine, pancreas and prostate have demonstrated the feasibility and capability of Monte Carlo based inverse planning for clinical IMRT. Dose distributions of selected treatment plans from a commercial non-Monte Carlo based system are evaluated in comparison with Monte Carlo based calculations. Systematic errors of up to 12% in tumor doses and up to 17% in critical structure doses have been observed. The clinical importance of Monte Carlo based
Reactor physics analysis method based on Monte Carlo homogenization
Background: Many new concepts of nuclear energy systems with complicated geometric structures and diverse energy spectra have been put forward to meet the future demand of nuclear energy market. The traditional deterministic neutronics analysis method has been challenged in two aspects: one is the ability of generic geometry processing; the other is the multi-spectrum applicability of the multi-group cross section libraries. The Monte Carlo (MC) method predominates the suitability of geometry and spectrum, but faces the problems of long computation time and slow convergence. Purpose: This work aims to find a novel scheme to take the advantages of both methods drawn from the deterministic core analysis method and MC method. Methods: A new two-step core analysis scheme is proposed to combine the geometry modeling capability and continuous energy cross section libraries of MC method, as well as the higher computational efficiency of deterministic method. First of all, the MC simulations are performed for assembly, and the assembly homogenized multi-group cross sections are tallied at the same time. Then, the core diffusion calculations can be done with these multi-group cross sections. Results: The new scheme can achieve high efficiency while maintain acceptable precision. Conclusion: The new scheme can be used as an effective tool for the design and analysis of innovative nuclear energy systems, which has been verified by numeric tests. (authors)
Application of Photon Transport Monte Carlo Module with GPU-based Parallel System
Park, Chang Je [Sejong University, Seoul (Korea, Republic of); Shon, Heejeong [Golden Eng. Co. LTD, Seoul (Korea, Republic of); Lee, Donghak [CoCo Link Inc., Seoul (Korea, Republic of)
2015-05-15
In general, it takes lots of computing time to get reliable results in Monte Carlo simulations especially in deep penetration problems with a thick shielding medium. To mitigate such a weakness of Monte Carlo methods, lots of variance reduction algorithms are proposed including geometry splitting and Russian roulette, weight windows, exponential transform, and forced collision, etc. Simultaneously, advanced computing hardware systems such as GPU(Graphics Processing Units)-based parallel machines are used to get a better performance of the Monte Carlo simulation. The GPU is much easier to access and to manage when comparing a CPU cluster system. It also becomes less expensive these days due to enhanced computer technology. There, lots of engineering areas adapt GPU-bases massive parallel computation technique. based photon transport Monte Carlo method. It provides almost 30 times speedup without any optimization and it is expected almost 200 times with fully supported GPU system. It is expected that GPU system with advanced parallelization algorithm will contribute successfully for development of the Monte Carlo module which requires quick and accurate simulations.
Application of Photon Transport Monte Carlo Module with GPU-based Parallel System
In general, it takes lots of computing time to get reliable results in Monte Carlo simulations especially in deep penetration problems with a thick shielding medium. To mitigate such a weakness of Monte Carlo methods, lots of variance reduction algorithms are proposed including geometry splitting and Russian roulette, weight windows, exponential transform, and forced collision, etc. Simultaneously, advanced computing hardware systems such as GPU(Graphics Processing Units)-based parallel machines are used to get a better performance of the Monte Carlo simulation. The GPU is much easier to access and to manage when comparing a CPU cluster system. It also becomes less expensive these days due to enhanced computer technology. There, lots of engineering areas adapt GPU-bases massive parallel computation technique. based photon transport Monte Carlo method. It provides almost 30 times speedup without any optimization and it is expected almost 200 times with fully supported GPU system. It is expected that GPU system with advanced parallelization algorithm will contribute successfully for development of the Monte Carlo module which requires quick and accurate simulations
Acceptance and implementation of a system of planning computerized based on Monte Carlo
It has been done the acceptance for use clinical Monaco computerized planning system, based on an on a virtual model of the energy yield of the head of the linear electron Accelerator and that performs the calculation of the dose with an algorithm of x-rays (XVMC) based on Monte Carlo algorithm. (Author)
A Monte-Carlo-Based Network Method for Source Positioning in Bioluminescence Tomography
Zhun Xu; Xiaolei Song; Xiaomeng Zhang; Jing Bai
2007-01-01
We present an approach based on the improved Levenberg Marquardt (LM) algorithm of backpropagation (BP) neural network to estimate the light source position in bioluminescent imaging. For solving the forward problem, the table-based random sampling algorithm (TBRS), a fast Monte Carlo simulation method ...
Human anatomical models have been indispensable to radiation protection dosimetry using Monte Carlo calculations. Existing MIRD-based mathematical models are easy to compute and standardize, but they are simplified and crude compared to human anatomy. This article describes the development of an image-based whole-body model, called VIP-Man, using transversal color photographic images obtained from the National Library of Medicine's Visible Human Project for Monte Carlo organ dose calculations involving photons, electron, neutrons, and protons. As the first of a series of papers on dose calculations based on VIP-Man, this article provides detailed information about how to construct an image-based model, as well as how to adopt it into well-tested Monte Carlo codes, EGS4, MCNP4B, and MCNPX
Proton therapy for tumors of the skull base
Munzenrider, J.E.; Liebsch, N.J. [Dept. of Radiation Oncology, Harvard Univ. Medical School, Boston, MA (United States)
1999-06-01
Charged particle beams are ideal for treating skull base and cervical spine tumors: dose can be focused in the target, while achieving significant sparing of the brain, brain stem, cervical cord, and optic nerves and chiasm. For skull base tumors, 10-year local control rates with combined proton-photon therapy are highest for chondrosarcomas, intermediate for male chordomas, and lowest for female chordomas (94%, 65%, and 42%, respectively). For cervical spine tumors, 10-year local control rates are not significantly different for chordomas and chondrosarcomas (54% and 48%, respectively), nor is there any difference in local control between males and females. Observed treatment-related morbidity has been judged acceptable, in view of the major morbidity and mortality which accompany uncontrolled tumor growth. (orig.)
Parodi, K; Kraemer, M; Sommerer, F; Naumann, J; Mairani, A; Brons, S
2010-01-01
Scanned ion beam delivery promises superior flexibility and accuracy for highly conformal tumour therapy in comparison to the usage of passive beam shaping systems. The attainable precision demands correct overlapping of the pencil-like beams which build up the entire dose distribution in the treatment field. In particular, improper dose application due to deviations of the lateral beam profiles from the nominal planning conditions must be prevented via appropriate beam monitoring in the beamline, prior to the entrance in the patient. To assess the necessary tolerance thresholds of the beam monitoring system at the Heidelberg Ion Beam Therapy Center, Germany, this study has investigated several worst-case scenarios for a sensitive treatment plan, namely scanned proton and carbon ion delivery to a small target volume at a shallow depth. Deviations from the nominal lateral beam profiles were simulated, which may occur because of misaligned elements or changes of the beam optic in the beamline. Data have been an...
Lopez-Tarjuelo, J.; Garcia-Molla, R.; Suan-Senabre, X. J.; Quiros-Higueras, J. Q.; Santos-Serra, A.; Marco-Blancas, N.; Calzada-Feliu, S.
2013-07-01
It has been done the acceptance for use clinical Monaco computerized planning system, based on an on a virtual model of the energy yield of the head of the linear electron Accelerator and that performs the calculation of the dose with an algorithm of x-rays (XVMC) based on Monte Carlo algorithm. (Author)
Makovicka, L.; Vasseur, A.; Sauget, M.; Martin, E.; Gschwind, R.; Henriet, J. [Universite de Franche-Comte, Equipe IRMA/ENISYS/FEMTO-ST, UMR6174 CNRS, 25 - Montbeliard (France); Vasseur, A.; Sauget, M.; Martin, E.; Gschwind, R.; Henriet, J.; Salomon, M. [Universite de Franche-Comte, Equipe AND/LIFC, 90 - Belfort (France)
2009-01-15
Monte Carlo codes, precise but slow, are very important tools in the vast majority of specialities connected to Radiation Physics, Radiation Protection and Dosimetry. A discussion about some other computing solutions is carried out; solutions not only based on the enhancement of computer power, or on the 'biasing'used for relative acceleration of these codes (in the case of photons), but on more efficient methods (A.N.N. - artificial neural network, C.B.R. - case-based reasoning - or other computer science techniques) already and successfully used for a long time in other scientific or industrial applications and not only Radiation Protection or Medical Dosimetry. (authors)
Response matrix Monte Carlo based on a general geometry local calculation for electron transport
A Response Matrix Monte Carlo (RMMC) method has been developed for solving electron transport problems. This method was born of the need to have a reliable, computationally efficient transport method for low energy electrons (below a few hundred keV) in all materials. Today, condensed history methods are used which reduce the computation time by modeling the combined effect of many collisions but fail at low energy because of the assumptions required to characterize the electron scattering. Analog Monte Carlo simulations are prohibitively expensive since electrons undergo coulombic scattering with little state change after a collision. The RMMC method attempts to combine the accuracy of an analog Monte Carlo simulation with the speed of the condensed history methods. Like condensed history, the RMMC method uses probability distributions functions (PDFs) to describe the energy and direction of the electron after several collisions. However, unlike the condensed history method the PDFs are based on an analog Monte Carlo simulation over a small region. Condensed history theories require assumptions about the electron scattering to derive the PDFs for direction and energy. Thus the RMMC method samples from PDFs which more accurately represent the electron random walk. Results show good agreement between the RMMC method and analog Monte Carlo. 13 refs., 8 figs
Convex-based void filling method for CAD-based Monte Carlo geometry modeling
Highlights: • We present a new void filling method named CVF for CAD based MC geometry modeling. • We describe convex based void description based and quality-based space subdivision. • The results showed improvements provided by CVF for both modeling and MC calculation efficiency. - Abstract: CAD based automatic geometry modeling tools have been widely applied to generate Monte Carlo (MC) calculation geometry for complex systems according to CAD models. Automatic void filling is one of the main functions in the CAD based MC geometry modeling tools, because the void space between parts in CAD models is traditionally not modeled while MC codes such as MCNP need all the problem space to be described. A dedicated void filling method, named Convex-based Void Filling (CVF), is proposed in this study for efficient void filling and concise void descriptions. The method subdivides all the problem space into disjointed regions using Quality based Subdivision (QS) and describes the void space in each region with complementary descriptions of the convex volumes intersecting with that region. It has been implemented in SuperMC/MCAM, the Multiple-Physics Coupling Analysis Modeling Program, and tested on International Thermonuclear Experimental Reactor (ITER) Alite model. The results showed that the new method reduced both automatic modeling time and MC calculation time
Pietrzak, Robert; Konefał, Adam; Sokół, Maria; Orlef, Andrzej
2016-08-01
The success of proton therapy depends strongly on the precision of treatment planning. Dose distribution in biological tissue may be obtained from Monte Carlo simulations using various scientific codes making it possible to perform very accurate calculations. However, there are many factors affecting the accuracy of modeling. One of them is a structure of objects called bins registering a dose. In this work the influence of bin structure on the dose distributions was examined. The MCNPX code calculations of Bragg curve for the 60 MeV proton beam were done in two ways: using simple logical detectors being the volumes determined in water, and using a precise model of ionization chamber used in clinical dosimetry. The results of the simulations were verified experimentally in the water phantom with Marcus ionization chamber. The average local dose difference between the measured relative doses in the water phantom and those calculated by means of the logical detectors was 1.4% at first 25 mm, whereas in the full depth range this difference was 1.6% for the maximum uncertainty in the calculations less than 2.4% and for the maximum measuring error of 1%. In case of the relative doses calculated with the use of the ionization chamber model this average difference was somewhat greater, being 2.3% at depths up to 25 mm and 2.4% in the full range of depths for the maximum uncertainty in the calculations of 3%. In the dose calculations the ionization chamber model does not offer any additional advantages over the logical detectors. The results provided by both models are similar and in good agreement with the measurements, however, the logical detector approach is a more time-effective method.
Monte-Carlo based prediction of radiochromic film response for hadrontherapy dosimetry
A model has been developed to calculate MD-55-V2 radiochromic film response to ion irradiation. This model is based on photon film response and film saturation by high local energy deposition computed by Monte-Carlo simulation. We have studied the response of the film to photon irradiation and we proposed a calculation method for hadron beams.
New memory devices based on the proton transfer process
Wierzbowska, Małgorzata
2016-01-01
Memory devices operating due to the fast proton transfer (PT) process are proposed by the means of first-principles calculations. Writing information is performed using the electrostatic potential of scanning tunneling microscopy (STM). Reading information is based on the effect of the local magnetization induced at the zigzag graphene nanoribbon (Z-GNR) edge—saturated with oxygen or the hydroxy group—and can be realized with the use of giant magnetoresistance (GMR), a magnetic tunnel junction or spin-transfer torque devices. The energetic barriers for the hop forward and backward processes can be tuned by the distance and potential of the STM tip; this thus enables us to tailor the non-volatile logic states. The proposed system enables very dense packing of the logic cells and could be used in random access and flash memory devices.
New memory devices based on the proton transfer process.
Wierzbowska, Małgorzata
2016-01-01
Memory devices operating due to the fast proton transfer (PT) process are proposed by the means of first-principles calculations. Writing information is performed using the electrostatic potential of scanning tunneling microscopy (STM). Reading information is based on the effect of the local magnetization induced at the zigzag graphene nanoribbon (Z-GNR) edge-saturated with oxygen or the hydroxy group-and can be realized with the use of giant magnetoresistance (GMR), a magnetic tunnel junction or spin-transfer torque devices. The energetic barriers for the hop forward and backward processes can be tuned by the distance and potential of the STM tip; this thus enables us to tailor the non-volatile logic states. The proposed system enables very dense packing of the logic cells and could be used in random access and flash memory devices. PMID:26596910
Proton beam micromachining on strippable aqueous base developable negative resist
Rajta, I. [Institute of Nuclear Research of the Hungarian Academy of Sciences, H-4001 Debrecen, P.O. Box 51 (Hungary)]. E-mail: rajta@atomki.hu; Baradacs, E. [University of Debrecen, Department of Environmental Physics, H-4026 Debrecen, Poroszlay u. 6 (Hungary); Chatzichristidi, M. [Institute of Microelectronics, NCSR-' Demokritos' , POB 62230, 153 10 Ag. Paraskevi (Greece); Valamontes, E.S. [Department of Electronics Technological Educational Institute of Athens, 12210 Aegaleo (Greece); Uzonyi, I. [Institute of Nuclear Research of the Hungarian Academy of Sciences, H-4001 Debrecen, P.O. Box 51 (Hungary); Raptis, I. [Institute of Microelectronics, NCSR-' Demokritos' , POB 62230, 153 10 Ag. Paraskevi (Greece)
2005-04-01
Nowadays a significant amount of research effort is devoted to the development of technologies for the fabrication of microcomponents and microsystems worldwide. In certain applications of micromachining high aspect ratio (HAR) structures are required. However, the resist materials used in HAR technologies are usually not compatible with the IC fabrication, either because they cannot be stripped away or because they are developed in organic solvents. In the present work the application of a novel chemically amplified resist for proton beam micromachining is presented. The resist based on epoxy and polyhydroxystyrene polymers is developed in the IC standard aqueous developers. The exposed areas can be stripped away using conventional organic stripping solutions. In order to test the exposure dose sensitivity and the lateral resolution, various test structures were irradiated. Using this formulation 5-8 {mu}m wide lines with aspect ratio 4-6 were resolved.
Proton beam micromachining on strippable aqueous base developable negative resist
Nowadays a significant amount of research effort is devoted to the development of technologies for the fabrication of microcomponents and microsystems worldwide. In certain applications of micromachining high aspect ratio (HAR) structures are required. However, the resist materials used in HAR technologies are usually not compatible with the IC fabrication, either because they cannot be stripped away or because they are developed in organic solvents. In the present work the application of a novel chemically amplified resist for proton beam micromachining is presented. The resist based on epoxy and polyhydroxystyrene polymers is developed in the IC standard aqueous developers. The exposed areas can be stripped away using conventional organic stripping solutions. In order to test the exposure dose sensitivity and the lateral resolution, various test structures were irradiated. Using this formulation 5-8 μm wide lines with aspect ratio 4-6 were resolved
A Muon Source Proton Driver at JPARC-based Parameters
Neuffer, David [Fermilab
2016-06-01
An "ultimate" high intensity proton source for neutrino factories and/or muon colliders was projected to be a ~4 MW multi-GeV proton source providing short, intense proton pulses at ~15 Hz. The JPARC ~1 MW accelerators provide beam at parameters that in many respects overlap these goals. Proton pulses from the JPARC Main Ring can readily meet the pulsed intensity goals. We explore these parameters, describing the overlap and consider extensions that may take a JPARC-like facility toward this "ultimate" source. JPARC itself could serve as a stage 1 source for such a facility.
Magro, G; Mairani, A; Mirandola, A; Panizza, D; Russo, S; Ferrari, A; Valvo, F; Fossati, P; Ciocca, M
2015-01-01
This study was performed to evaluate the accuracy of a commercial treatment planning system (TPS), in optimising proton pencil beam dose distributions for small targets of different sizes (5–30 mm side) located at increasing depths in water. The TPS analytical algorithm was benchmarked against experimental data and the FLUKA Monte Carlo (MC) code, previously validated for the selected beam-line. We tested the Siemens syngo® TPS plan optimisation module for water cubes fixing the configurable parameters at clinical standards, with homogeneous target coverage to a 2 Gy (RBE) dose prescription as unique goal. Plans were delivered and the dose at each volume centre was measured in water with a calibrated PTW Advanced Markus® chamber. An EBT3® film was also positioned at the phantom entrance window for the acquisition of 2D dose maps. Discrepancies between TPS calculated and MC simulated values were mainly due to the different lateral spread modeling and resulted in being related to the field-to-spot size r...
Development of Monte Carlo-based pebble bed reactor fuel management code
Highlights: • A new Monte Carlo-based fuel management code for OTTO cycle pebble bed reactor was developed. • The double-heterogeneity was modeled using statistical method in MVP-BURN code. • The code can perform analysis of equilibrium and non-equilibrium phase. • Code-to-code comparisons for Once-Through-Then-Out case were investigated. • Ability of the code to accommodate the void cavity was confirmed. - Abstract: A fuel management code for pebble bed reactors (PBRs) based on the Monte Carlo method has been developed in this study. The code, named Monte Carlo burnup analysis code for PBR (MCPBR), enables a simulation of the Once-Through-Then-Out (OTTO) cycle of a PBR from the running-in phase to the equilibrium condition. In MCPBR, a burnup calculation based on a continuous-energy Monte Carlo code, MVP-BURN, is coupled with an additional utility code to be able to simulate the OTTO cycle of PBR. MCPBR has several advantages in modeling PBRs, namely its Monte Carlo neutron transport modeling, its capability of explicitly modeling the double heterogeneity of the PBR core, and its ability to model different axial fuel speeds in the PBR core. Analysis at the equilibrium condition of the simplified PBR was used as the validation test of MCPBR. The calculation results of the code were compared with the results of diffusion-based fuel management PBR codes, namely the VSOP and PEBBED codes. Using JENDL-4.0 nuclide library, MCPBR gave a 4.15% and 3.32% lower keff value compared to VSOP and PEBBED, respectively. While using JENDL-3.3, MCPBR gave a 2.22% and 3.11% higher keff value compared to VSOP and PEBBED, respectively. The ability of MCPBR to analyze neutron transport in the top void of the PBR core and its effects was also confirmed
Relativistic protons with energies 0.7-1.5 GeV interacting with a thick, cylindrical, lead target, surrounded by a uranium blanket and a polyethylene moderator, produced spallation neutrons. The spatial and energetic distributions of the produced neutron field were measured by the Activation Analysis Method using Al, Au, Bi, and Co radio-chemical sensors. The experimental yields of isotopes induced in the sensors were compared with Monte-Carlo calculations performed with the MCNPX 2.4.0 code
Jeraj, Robert; Keall, Paul
2000-12-01
The effect of the statistical uncertainty, or noise, in inverse treatment planning for intensity modulated radiotherapy (IMRT) based on Monte Carlo dose calculation was studied. Sets of Monte Carlo beamlets were calculated to give uncertainties at Dmax ranging from 0.2% to 4% for a lung tumour plan. The weights of these beamlets were optimized using a previously described procedure based on a simulated annealing optimization algorithm. Several different objective functions were used. It was determined that the use of Monte Carlo dose calculation in inverse treatment planning introduces two errors in the calculated plan. In addition to the statistical error due to the statistical uncertainty of the Monte Carlo calculation, a noise convergence error also appears. For the statistical error it was determined that apparently successfully optimized plans with a noisy dose calculation (3% 1σ at Dmax ), which satisfied the required uniformity of the dose within the tumour, showed as much as 7% underdose when recalculated with a noise-free dose calculation. The statistical error is larger towards the tumour and is only weakly dependent on the choice of objective function. The noise convergence error appears because the optimum weights are determined using a noisy calculation, which is different from the optimum weights determined for a noise-free calculation. Unlike the statistical error, the noise convergence error is generally larger outside the tumour, is case dependent and strongly depends on the required objectives.
Proton beam micromachining on strippable aqueous base developable negative resist
Complete text of publication follows. Proton Beam Micromachining (PBM, also known as P-beam writing), a novel direct- write process for the production of 3D microstructures, can be used to make multilevel structures in a single layer of resist by varying the ion energy. The interaction between the bombarding ions and the target material is mainly ionization, and very few ions suffer high angle nuclear collisions, therefore structures made with PBM have smooth near vertical side walls. The most commony applied resists in PBM are the positive, conventional, polymethyl methacrylate (PMMA); and the negative, chemically amplified, SU-8 (Micro Chem Corp). SU-8 is an epoxy based resist suitable also for LIGA and UV-LIGA processes, it offers good sensitivity, good process latitude, very high aspect ratio and therefore it dominates in the high aspect ratio micromachining applications. SU-8 requires 30 nC/mm2 fluence for PBM irradiations at 2 MeV protons. Its crosslinking chemistry is based on the eight epoxy rings in the polymer chain, which provide a very dense three dimensional network in the presence of suitably activated photo acid generators (PAGs) which is very difficult to be stripped away after development. Thus, stripping has to be assisted with plasma processes or with special liquid removers. Moreover, the SU-8 developer is organic, propylene glycol methyl ether acetate (PGMEA), and thus environmentally non-friendly. To overcome the SU-8 stripping limitations, design of a negative resist system where solubility change is not based solely on cross- linking but also on the differentiation of hydrophilicity between exposed and non-exposed areas is desirable. A new resist formulation, fulfilling the above specifications has been developed recently [1]. This formulation is based on a specific grade epoxy novolac (EP) polymer, a partially hydrogenated poly-4-hydroxy styrene (PHS) polymer, and an onium salt as photoacid generator (PAG), and has been successfully applied
Feasibility Study of Accumulator and Compressor for the 6-bunches SPL based Proton Driver
Aiba, M
2008-01-01
Feasibility of the accumulator and the compressor ring for the SPL based proton driver have been studied for a future neutrino factory. The scenario retained for the SPL proton driver uses six bunches, with 10^14 protons in total at 50 Hz. Possible lattices for the accumulator and the compressor are presented. The beam injection/accumulation and the bunch compression are delicate issues and discussed in detail in this note. Throughout the presented study, these difficulties are disclosed not to be critical issues, and together with a discussion on the focusing towards production target, the feasibility of the 6-bunches SPL based proton driver has been confirmed.
Feasibility study of accumulator and compressor for the 6-bunches SPL based proton driver
Aiba, M
2008-01-01
Feasibility studies of the accumulator and the compressor rings for the SPL based proton driver have been studied for a future neutrino factory. The scenario retained for the SPL proton driver uses six bunches, with 10^14 protons in total at 50 Hz. Possible lattices for the accumulator and the compressor are presented. The beam injection/accumulation and the bunch compression are delicate issues and discussed in detail in this note. Throughout the presented study, these difficulties are disclosed not to be critical issues, and together with a discussion on the focusing towards production target, the feasibility of the 6-bunches SPL based proton driver has been confirmed
Monte Carlo vs. Pencil Beam based optimization of stereotactic lung IMRT
Weinmann Martin; Söhn Matthias; Muzik Jan; Sikora Marcin; Alber Markus
2009-01-01
Abstract Background The purpose of the present study is to compare finite size pencil beam (fsPB) and Monte Carlo (MC) based optimization of lung intensity-modulated stereotactic radiotherapy (lung IMSRT). Materials and methods A fsPB and a MC algorithm as implemented in a biological IMRT planning system were validated by film measurements in a static lung phantom. Then, they were applied for static lung IMSRT planning based on three different geometrical patient models (one phase static CT, ...
Monte Carlo vs. Pencil Beam based optimization of stereotactic lung IMRT
Sikora, Marcin; Muzik, Jan; Söhn, Matthias; Weinmann, Martin; Alber, Markus
2009-01-01
Background The purpose of the present study is to compare finite size pencil beam (fsPB) and Monte Carlo (MC) based optimization of lung intensity-modulated stereotactic radiotherapy (lung IMSRT). Materials and methods A fsPB and a MC algorithm as implemented in a biological IMRT planning system were validated by film measurements in a static lung phantom. Then, they were applied for static lung IMSRT planning based on three different geometrical patient models (one phase static CT, density o...
Monte Carlo vs. Pencil Beam based optimization of stereotactic lung IMRT
Sikora, Marcin Pawel; Muzik, Jan; Söhn, Matthias; Weinmann, Martin; Alber, Markus
2009-01-01
Background: The purpose of the present study is to compare finite size pencil beam (fsPB) and Monte Carlo (MC) based optimization of lung intensity-modulated stereotactic radiotherapy (lung IMSRT). Materials and methods: A fsPB and a MC algorithm as implemented in a biological IMRT planning system were validated by film measurements in a static lung phantom. Then, they were applied for static lung IMSRT planning based on three different geometrical patient models (one phase ...
Grządziel Małgorzata
2014-07-01
Full Text Available Verification of calculations of the depth-dose distributions in water, using GEANT4 (version of 4.9.3 and MCNPX (version of 2.7.0 Monte Carlo codes, was performed for the scatterer-phantom system used in the dosimetry measurements in the proton therapy of eye tumours. The simulated primary proton beam had the energy spectra distributed according to the Gauss distribution with the cut at energy greater than that related to the maximum of the spectrum. The energy spectra of the primary protons were chosen to get the possibly best agreement between the measured relative depth-dose distributions along the central-axis of the proton beam in a water phantom and that derived from the Monte Carlo calculations separately for the both tested codes. The local depth-dose differences between results from the calculations and the measurements were mostly less than 5% (the mean value of 2.1% and 3.6% for the MCNPX and GEANT4 calculations. In the case of the MCNPX calculations, the best fit to the experimental data was obtained for the spectrum with maximum at 60.8 MeV (more probable energy, FWHM of the spectrum of 0.4 MeV and the energy cut at 60.85 MeV whereas in the GEANT4 calculations more probable energy was 60.5 MeV, FWHM of 0.5 MeV, the energy cut at 60.7 MeV. Thus, one can say that the results obtained by means of the both considered Monte Carlo codes are similar but they are not the same. Therefore the agreement between the calculations and the measurements has to be verified before each application of the MCNPX and GEANT4 codes for the determination of the depth-dose curves for the therapeutic protons.
A Markov Chain Monte Carlo Based Method for System Identification
Glaser, R E; Lee, C L; Nitao, J J; Hanley, W G
2002-10-22
This paper describes a novel methodology for the identification of mechanical systems and structures from vibration response measurements. It combines prior information, observational data and predictive finite element models to produce configurations and system parameter values that are most consistent with the available data and model. Bayesian inference and a Metropolis simulation algorithm form the basis for this approach. The resulting process enables the estimation of distributions of both individual parameters and system-wide states. Attractive features of this approach include its ability to: (1) provide quantitative measures of the uncertainty of a generated estimate; (2) function effectively when exposed to degraded conditions including: noisy data, incomplete data sets and model misspecification; (3) allow alternative estimates to be produced and compared, and (4) incrementally update initial estimates and analysis as more data becomes available. A series of test cases based on a simple fixed-free cantilever beam is presented. These results demonstrate that the algorithm is able to identify the system, based on the stiffness matrix, given applied force and resultant nodal displacements. Moreover, it effectively identifies locations on the beam where damage (represented by a change in elastic modulus) was specified.
Effects of CT based Voxel Phantoms on Dose Distribution Calculated with Monte Carlo Method
Chen Chaobin; Huang Qunying; Wu Yican
2005-01-01
A few CT-based voxel phantoms were produced to investigate the sensitivity of Monte Carlo simulations of X-ray beam and electron beam to the proportions of elements and the mass densities of the materials used to express the patient's anatomical structure. The human body can be well outlined by air, lung, adipose, muscle, soft bone and hard bone to calculate the dose distribution with Monte Carlo method. The effects of the calibration curves established by using various CT scanners are not clinically significant based on our investigation. The deviation from the values of cumulative dose volume histogram derived from CT-based voxel phantoms is less than 1% for the given target.
Effects of CT based Voxel Phantoms on Dose Distribution Calculated with Monte Carlo Method
Chen, Chaobin; Huang, Qunying; Wu, Yican
2005-04-01
A few CT-based voxel phantoms were produced to investigate the sensitivity of Monte Carlo simulations of x-ray beam and electron beam to the proportions of elements and the mass densities of the materials used to express the patient's anatomical structure. The human body can be well outlined by air, lung, adipose, muscle, soft bone and hard bone to calculate the dose distribution with Monte Carlo method. The effects of the calibration curves established by using various CT scanners are not clinically significant based on our investigation. The deviation from the values of cumulative dose volume histogram derived from CT-based voxel phantoms is less than 1% for the given target.
Caldwell, Allen
2015-01-01
Based on current CERN infrastructure, an electron-proton collider is proposed at a centre-of-mass energy of about 9 TeV. A 7 TeV LHC bunch is used as the proton driver to create a plasma wakefield which then accelerates electrons to 3 TeV, these then colliding with the other 7 TeV LHC proton beam. The basic parameters of the collider are presented, which although of very high energy, has integrated luminosities of the order of 1 pb$^{-1}$/year. For such a collider, with a centre-of-mass energy 30 times greater than HERA, parton momentum fractions, $x$, down to about $10^{-8}$ are accessible for $Q^2$ of 1 GeV$^2$ and could lead to effects of saturation or some other breakdown of DGLAP being observed. The total photon-proton cross section can be measured up to very high energies and also at different energies as the possibility of varying the electron beam energy is assumed; this could have synergy with cosmic-ray physics. Other physics which can be pursued at such a collider are contact interaction searches, ...
Relative motion between a tumor and a scanning proton beam results in a degradation of the dose distribution (interplay effect). This study investigates the relationship between beam scanning parameters and the interplay effect, with the goal of finding parameters that minimize interplay. 4D Monte Carlo simulations of pencil beam scanning proton therapy treatments were performed using the 4DCT geometry of five lung cancer patients of varying tumor size (50.4–167.1 cc) and motion amplitude (2.9–30.1 mm). Treatments were planned assuming delivery in 35 × 2.5 Gy(RBE) fractions. The spot size, time to change the beam energy (τes), time required for magnet settling (τss), initial breathing phase, spot spacing, scanning direction, scanning speed, beam current and patient breathing period were varied for each of the five patients. Simulations were performed for a single fraction and an approximation of conventional fractionation. For the patients considered, the interplay effect could not be predicted using the superior–inferior motion amplitude alone. Larger spot sizes (σ ∼ 9–16 mm) were less susceptible to interplay, giving an equivalent uniform dose (EUD) of 99.0 ± 4.4% (1 standard deviation) in a single fraction compared to 86.1 ± 13.1% for smaller spots (σ ∼ 2–4 mm). The smaller spot sizes gave EUD values as low as 65.3% of the prescription dose in a single fraction. Reducing the spot spacing improved the target dose homogeneity. The initial breathing phase can have a significant effect on the interplay, particularly for shorter delivery times. No clear benefit was evident when scanning either parallel or perpendicular to the predominant axis of motion. Longer breathing periods decreased the EUD. In general, longer delivery times led to lower interplay effects. Conventional fractionation showed significant improvement in terms of interplay, giving a EUD of at least 84.7% and 100.0% of the prescription dose for the small and larger spot sizes
Biniam Tesfamicael
2016-03-01
Full Text Available Purpose: The main purpose of this study was to monitor the secondary dose distribution originating from a water phantom during proton therapy of prostate cancer using scintillating fibers.Methods: The Geant4 Monte Carlo toolkit version 9.6.p02 was used to simulate a proton therapy of prostate cancer. Two cases were studied. In the first case, 8 × 8 = 64 equally spaced fibers inside three 4 × 4 × 2.54 cm3 Delrin® blocks were used to monitor the emission of secondary particles in the transverse (left and right and distal regions relative to the beam direction. In the second case, a scintillating block with a thickness of 2.54 cm and equal vertical and longitudinal dimensions as the water phantom was used. Geometrical cuts were implemented to extract the energy deposited in each fiber and inside the scintillating block.Results: The transverse dose distributions from the detected secondary particles in both cases are symmetric and agree to within <3.6%. The energy deposited gradually increases as one moves from the peripheral row of fibers towards the center of the block (aligned with the center of the prostate by a factor of approximately 5. The energy deposited was also observed to decrease as one goes from the frontal to distal region of the block. The ratio of the energy deposited in the prostate to the energy deposited in the middle two rows of fibers showed a linear relationship with a slope of (-3.55±2.26 × 10-5 MeV per treatment Gy delivered. The distal detectors recorded a negligible amount of energy deposited due to higher attenuation of the secondary particles by the water in that direction.Conclusion: With a good calibration and with the ability to define a good correlation between the radiation flux recorded by the external fibers and the dose delivered to the prostate, such fibers can be used for real time dose verification to the target. The system was also observed to respond to the series of Bragg Peaks used to generate the
Dowdell, S.; Grassberger, C.; Sharp, G. C.; Paganetti, H.
2013-06-01
Relative motion between a tumor and a scanning proton beam results in a degradation of the dose distribution (interplay effect). This study investigates the relationship between beam scanning parameters and the interplay effect, with the goal of finding parameters that minimize interplay. 4D Monte Carlo simulations of pencil beam scanning proton therapy treatments were performed using the 4DCT geometry of five lung cancer patients of varying tumor size (50.4-167.1 cc) and motion amplitude (2.9-30.1 mm). Treatments were planned assuming delivery in 35 × 2.5 Gy(RBE) fractions. The spot size, time to change the beam energy (τes), time required for magnet settling (τss), initial breathing phase, spot spacing, scanning direction, scanning speed, beam current and patient breathing period were varied for each of the five patients. Simulations were performed for a single fraction and an approximation of conventional fractionation. For the patients considered, the interplay effect could not be predicted using the superior-inferior motion amplitude alone. Larger spot sizes (σ ˜ 9-16 mm) were less susceptible to interplay, giving an equivalent uniform dose (EUD) of 99.0 ± 4.4% (1 standard deviation) in a single fraction compared to 86.1 ± 13.1% for smaller spots (σ ˜ 2-4 mm). The smaller spot sizes gave EUD values as low as 65.3% of the prescription dose in a single fraction. Reducing the spot spacing improved the target dose homogeneity. The initial breathing phase can have a significant effect on the interplay, particularly for shorter delivery times. No clear benefit was evident when scanning either parallel or perpendicular to the predominant axis of motion. Longer breathing periods decreased the EUD. In general, longer delivery times led to lower interplay effects. Conventional fractionation showed significant improvement in terms of interplay, giving a EUD of at least 84.7% and 100.0% of the prescription dose for the small and larger spot sizes respectively
Development of 3d reactor burnup code based on Monte Carlo method and exponential Euler method
Burnup analysis plays a key role in fuel breeding, transmutation and post-processing in nuclear reactor. Burnup codes based on one-dimensional and two-dimensional transport method have difficulties in meeting the accuracy requirements. A three-dimensional burnup analysis code based on Monte Carlo method and Exponential Euler method has been developed. The coupling code combines advantage of Monte Carlo method in complex geometry neutron transport calculation and FISPACT in fast and precise inventory calculation, meanwhile resonance Self-shielding effect in inventory calculation can also be considered. The IAEA benchmark text problem has been adopted for code validation. Good agreements were shown in the comparison with other participants' results. (authors)
Proton pump inhibitors in cirrhosis: Tradition or evidence based practice?
Francesca Lodato; Francesco Azzaroli; Maria Di Girolamo; Valentina Feletti; Paolo Cecinato; Andrea Lisotti; Davide Festi; Enrico Roda; Giuseppe Mazzella
2008-01-01
Proton Pump Inhibitors (PPI) are very effective in inhibiting acid secretion and are extensively used in many acid related diseases. They are also often used in patients with cirrhosis sometimes in the absence of a specific acid related disease, with the aim of preventing peptic complications in patients with variceal or hypertensive gastropathic bleeding receiving multidrug treatment. Contradicting reports support their use in cirrhosis and evidence of their efficacy in this condition is poor. Moreover there are convincing papers suggesting that acid secretion is reduced in patients with liver cirrhosis. With regard to H pylori infection, its prevalence in patients with cirrhosis is largely variable among different studies, and it seems that H pylori eradication does not prevent gastro-duodenal ulcer formation and bleeding. With regard to the prevention and treatment of oesophageal complications after banding or sclerotherapy of oesophageal varices, there is little evidence for a protective role of PPI. Moreover, due to liver metabolism of PPI, the dose of most available PPIs should be reduced in cirrhotics. In conclusion, the use of this class of drugs seems more habit related than evidence-based eventually leading to an increase in health costs.
New approach to polarized proton scattering based on Dirac dynamics
The Dirac impulse approximation has to date provided dramatic improvement in our ability to predict, with no free parameters, spin observables in proton-nucleus elastic scattering at intermediate energies. The key ingredients of this approach are Dirac propagation and the nucleon-nucleon invariant amplitudes. So far, local approximations to the NN amplitudes have been used. The standard NN representation in terms of Dirac scalar, vector, and so on, parts which is free of kinematical singularities seems to naturally predict the correct coupling to negative energy states for energies above 300 MeV. At low energy, this coupling is subject to an ambiguity between pseudoscalar and pseudovector πN coupling mechanisms and it is evident that the pseudoscalar coupling treated in a local approximation causes too much scalar-vector difference and thus too large pair contributions. Once this problem is remedied, the Dirac optical potential is expected to be calculable from a nucleon-nucleon quasi-potential over the range 0 to 1000 MeV. For the energy region above about 300 MeV, the large scalar and vector potentials of Dirac phenomenology are seen to be accurately predicted by the impulse approximation. Work by Shakin and collaborators provides complementary results at low energy based on a nuclear matter g-matrix. A basic conclusion is that relativistic spin effects cannot be neglected in nuclear physics. 36 references
Dowdell, Stephen; Paganetti, Harald [Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114 (United States); Grassberger, Clemens [Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts 02114 and Centre for Proton Therapy, Paul Scherrer Institut, 5232 Villigen-PSI (Switzerland)
2013-12-15
Purpose: To compare motion effects in intensity modulated proton therapy (IMPT) lung treatments with different levels of intensity modulation.Methods: Spot scanning IMPT treatment plans were generated for ten lung cancer patients for 2.5Gy(RBE) and 12Gy(RBE) fractions and two distinct energy-dependent spot sizes (σ∼8–17 mm and ∼2–4 mm). IMPT plans were generated with the target homogeneity of each individual field restricted to <20% (IMPT{sub 20%}). These plans were compared to full IMPT (IMPT{sub full}), which had no restriction on the single field homogeneity. 4D Monte Carlo simulations were performed upon the patient 4DCT geometry, including deformable image registration and incorporating the detailed timing structure of the proton delivery system. Motion effects were quantified via comparison of the results of the 4D simulations (4D-IMPT{sub 20%}, 4D-IMPT{sub full}) with those of a 3D Monte Carlo simulation (3D-IMPT{sub 20%}, 3D-IMPT{sub full}) upon the planning CT using the equivalent uniform dose (EUD), V{sub 95} and D{sub 1}-D{sub 99}. The effects in normal lung were quantified using mean lung dose (MLD) and V{sub 90%}.Results: For 2.5Gy(RBE), the mean EUD for the large spot size is 99.9%± 2.8% for 4D-IMPT{sub 20%} compared to 100.1%± 2.9% for 4D-IMPT{sub full}. The corresponding values are 88.6%± 8.7% (4D-IMPT{sub 20%}) and 91.0%± 9.3% (4D-IMPT{sub full}) for the smaller spot size. The EUD value is higher in 69.7% of the considered deliveries for 4D-IMPT{sub full}. The V{sub 95} is also higher in 74.7% of the plans for 4D-IMPT{sub full}, implying that IMPT{sub full} plans experience less underdose compared to IMPT{sub 20%}. However, the target dose homogeneity is improved in the majority (67.8%) of plans for 4D-IMPT{sub 20%}. The higher EUD and V{sub 95} suggests that the degraded homogeneity in IMPT{sub full} is actually due to the introduction of hot spots in the target volume, perhaps resulting from the sharper in-target dose gradients. The
We briefly present our atomistic kinetic Monte Carlo approach to model the diffusion of point-defects in Fe-based alloys, and therefore to simulate diffusion induced mass transport and subsequent nano-structural and microchemical changes. This methodology has been hitherto successfully applied to the simulation of thermal annealing experiments. We here present our achievements in the generalization of this method to the simulation of neutron irradiation damage. (authors)
Validation of a Monte Carlo Based Depletion Methodology Using HFIR Post-Irradiation Measurements
Chandler, David [ORNL; Maldonado, G Ivan [ORNL; Primm, Trent [ORNL
2009-11-01
Post-irradiation uranium isotopic atomic densities within the core of the High Flux Isotope Reactor (HFIR) were calculated and compared to uranium mass spectrographic data measured in the late 1960s and early 70s [1]. This study was performed in order to validate a Monte Carlo based depletion methodology for calculating the burn-up dependent nuclide inventory, specifically the post-irradiation uranium
Laser-based detection and tracking moving objects using data-driven Markov chain Monte Carlo
Vu, Trung-Dung; Aycard, Olivier
2009-01-01
We present a method of simultaneous detection and tracking moving objects from a moving vehicle equipped with a single layer laser scanner. A model-based approach is introduced to interpret the laser measurement sequence by hypotheses of moving object trajectories over a sliding window of time. Knowledge of various aspects including object model, measurement model, motion model are integrated in one theoretically sound Bayesian framework. The data-driven Markov chain Monte Carlo (DDMCMC) tech...
Development of a space radiation Monte Carlo computer simulation based on the FLUKA and ROOT codes
Pinsky, L; Ferrari, A; Sala, P; Carminati, F; Brun, R
2001-01-01
This NASA funded project is proceeding to develop a Monte Carlo-based computer simulation of the radiation environment in space. With actual funding only initially in place at the end of May 2000, the study is still in the early stage of development. The general tasks have been identified and personnel have been selected. The code to be assembled will be based upon two major existing software packages. The radiation transport simulation will be accomplished by updating the FLUKA Monte Carlo program, and the user interface will employ the ROOT software being developed at CERN. The end-product will be a Monte Carlo-based code which will complement the existing analytic codes such as BRYNTRN/HZETRN presently used by NASA to evaluate the effects of radiation shielding in space. The planned code will possess the ability to evaluate the radiation environment for spacecraft and habitats in Earth orbit, in interplanetary space, on the lunar surface, or on a planetary surface such as Mars. Furthermore, it will be usef...
The influence of air cavities within the PTV on Monte Carlo-based IMRT optimization
Smedt, Bart de [Department of Medical Physics, Ghent University, Gent (Belgium); Vanderstraeten, Barbara [Department of Medical Physics, Ghent University, Gent (Belgium); Reynaert, Nick [Department of Medical Physics, Ghent University, Gent (Belgium); Gersem, Werner de [Department of Radiotherapy, Ghent University Hospital, Gent (Belgium); Neve, Wilfried de [Department of Radiotherapy, Ghent University Hospital, Gent (Belgium); Thierens, Hubert [Department of Medical Physics, Ghent University, Gent (Belgium)
2007-06-15
Integrating Monte Carlo calculated dose distributions into an iterative aperture-based IMRT optimization process can improve the final treatment plan. However, the influence of large air cavities in the planning target volume (PTV) on the outcome of the optimization process should not be underestimated. To study this influence, the treatment plan of an ethmoid sinus cancer patient, which has large air cavities included in the PTV, is iteratively optimized in two different situations, namely when the large air cavities are included in the PTV and when these air cavities are excluded from the PTV. Two optimization methods were applied to integrate the Monte Carlo calculated dose distributions into the optimization process, namely the 'Correction-method' and the 'Per Segment-method'. The 'Correction-method' takes the Monte Carlo calculated global dose distribution into account in the optimization process by means of a correction matrix, which is in fact a dose distribution that is equal to the difference between the Monte Carlo calculated global dose distribution and the global dose distribution calculated by a conventional dose calculation algorithm. The 'Per Segment-method' uses directly the Monte Carlo calculated dose distributions of the individual segments in the optimization process. Both methods tend to converge whether or not large air cavities are excluded from the PTV during the optimization process. However, the 'Per Segment-method' performs better than the 'Correction-method' in both situations and the 'Per Segment-method' in the case where the large air cavities are excluded from the PTV leads to a better treatment plan then when these air cavities are included. Therefore we advise to exclude large air cavities and to apply the 'Per Segment-method' to integrate the Monte Carlo dose calculations into an iterative aperture-based optimization process. Nevertheless, the &apos
Supanitsky, A. D.
2016-02-01
Gamma rays and neutrinos are produced as a result of proton-proton interactions that occur in different astrophysical contexts. The detection of these two types of messengers is of great importance for the study of different physical phenomena, related to nonthermal processes, taking place in different astrophysical scenarios. Therefore, the knowledge of the energy spectrum of these two types of particles, as a function of the incident proton energy, is essential for the interpretation of the observational data. In this paper, parametrizations of the energy spectra of gamma rays and neutrinos, originated in proton-proton collisions, are presented. The energy range of the incident protons considered extends from 102 to 108 GeV . The parametrizations are based on Monte Carlo simulations of proton-proton interactions performed with the hadronic interaction models QGSJET-II-04 and EPOS-LHC, which have recently been updated with the data taken by the Large Hadron Collider.
Protonation Equilibrium of Linear Homopolyacids
Požar J.
2015-07-01
Full Text Available The paper presents a short summary of investigations dealing with protonation equilibrium of linear homopolyacids, in particularly those of high charge density. Apart from the review of experimental results which can be found in the literature, a brief description of theoretical models used in processing the dependence of protonation constants on monomer dissociation degree and ionic strength is given (cylindrical model based on Poisson-Boltzmann equation, cylindrical Stern model, the models according to Ising, Högfeldt, Mandel and Katchalsky. The applicability of these models regarding the polyion charge density, electrolyte concentration and counterion type is discussed. The results of Monte Carlo simulations of protonation equilibrium are also briefly mentioned. In addition, frequently encountered errors connected with calibration of of glass electrode and the related unreliability of determined protonation constants are pointed out.
ERSN-OpenMC, a Java-based GUI for OpenMC Monte Carlo code
Jaafar EL Bakkali
2016-07-01
Full Text Available OpenMC is a new Monte Carlo transport particle simulation code focused on solving two types of neutronic problems mainly the k-eigenvalue criticality fission source problems and external fixed fission source problems. OpenMC does not have any Graphical User Interface and the creation of one is provided by our java-based application named ERSN-OpenMC. The main feature of this application is to provide to the users an easy-to-use and flexible graphical interface to build better and faster simulations, with less effort and great reliability. Additionally, this graphical tool was developed with several features, as the ability to automate the building process of OpenMC code and related libraries as well as the users are given the freedom to customize their installation of this Monte Carlo code. A full description of the ERSN-OpenMC application is presented in this paper.
Polarization imaging of multiply-scattered radiation based on integral-vector Monte Carlo method
A new integral-vector Monte Carlo method (IVMCM) is developed to analyze the transfer of polarized radiation in 3D multiple scattering particle-laden media. The method is based on a 'successive order of scattering series' expression of the integral formulation of the vector radiative transfer equation (VRTE) for application of efficient statistical tools to improve convergence of Monte Carlo calculations of integrals. After validation against reference results in plane-parallel layer backscattering configurations, the model is applied to a cubic container filled with uniformly distributed monodispersed particles and irradiated by a monochromatic narrow collimated beam. 2D lateral images of effective Mueller matrix elements are calculated in the case of spherical and fractal aggregate particles. Detailed analysis of multiple scattering regimes, which are very similar for unpolarized radiation transfer, allows identifying the sensitivity of polarization imaging to size and morphology.
Monte Carlo Methods in Materials Science Based on FLUKA and ROOT
Pinsky, Lawrence; Wilson, Thomas; Empl, Anton; Andersen, Victor
2003-01-01
A comprehensive understanding of mitigation measures for space radiation protection necessarily involves the relevant fields of nuclear physics and particle transport modeling. One method of modeling the interaction of radiation traversing matter is Monte Carlo analysis, a subject that has been evolving since the very advent of nuclear reactors and particle accelerators in experimental physics. Countermeasures for radiation protection from neutrons near nuclear reactors, for example, were an early application and Monte Carlo methods were quickly adapted to this general field of investigation. The project discussed here is concerned with taking the latest tools and technology in Monte Carlo analysis and adapting them to space applications such as radiation shielding design for spacecraft, as well as investigating how next-generation Monte Carlos can complement the existing analytical methods currently used by NASA. We have chosen to employ the Monte Carlo program known as FLUKA (A legacy acronym based on the German for FLUctuating KAscade) used to simulate all of the particle transport, and the CERN developed graphical-interface object-oriented analysis software called ROOT. One aspect of space radiation analysis for which the Monte Carlo s are particularly suited is the study of secondary radiation produced as albedoes in the vicinity of the structural geometry involved. This broad goal of simulating space radiation transport through the relevant materials employing the FLUKA code necessarily requires the addition of the capability to simulate all heavy-ion interactions from 10 MeV/A up to the highest conceivable energies. For all energies above 3 GeV/A the Dual Parton Model (DPM) is currently used, although the possible improvement of the DPMJET event generator for energies 3-30 GeV/A is being considered. One of the major tasks still facing us is the provision for heavy ion interactions below 3 GeV/A. The ROOT interface is being developed in conjunction with the
Proton exchange in acid–base complexes induced by reaction coordinates with heavy atom motions
Highlights: ► Proton exchange in acid–base complexes is studied. ► The structures, binding energies, and normal mode vibrations are calculated. ► Transition state structures of proton exchange mechanism are determined. ► In the complexes studied, the reaction coordinate involves heavy atom rocking. ► The reaction coordinate is not simply localized in the proton movements. - Abstract: We extend previous work on nitric acid–ammonia and nitric acid–alkylamine complexes to illustrate that proton exchange reaction coordinates involve the rocking motion of the base moiety in many double hydrogen-bonded gas phase strong acid–strong base complexes. The complexes studied involve the biologically and atmospherically relevant glycine, formic, acetic, propionic, and sulfuric acids with ammonia/alkylamine bases. In these complexes, the magnitude of the imaginary frequencies associated with the proton exchange transition states are −1. This contrasts with widely studied proton exchange reactions between symmetric carboxylic acid dimers or asymmetric DNA base pair and their analogs where the reaction coordinate is localized in proton motions and the magnitude of the imaginary frequencies for the transition states are >1100 cm−1. Calculations on complexes of these acids with water are performed for comparison. Variations of normal vibration modes along the reaction coordinate in the complexes are described.
Espel, Federico Puente
The main objective of this PhD research is to develop a high accuracy modeling tool using a Monte Carlo based coupled system. The presented research comprises the development of models to include the thermal-hydraulic feedback to the Monte Carlo method and speed-up mechanisms to accelerate the Monte Carlo criticality calculation. Presently, deterministic codes based on the diffusion approximation of the Boltzmann transport equation, coupled with channel-based (or sub-channel based) thermal-hydraulic codes, carry out the three-dimensional (3-D) reactor core calculations of the Light Water Reactors (LWRs). These deterministic codes utilize nuclear homogenized data (normally over large spatial zones, consisting of fuel assembly or parts of fuel assembly, and in the best case, over small spatial zones, consisting of pin cell), which is functionalized in terms of thermal-hydraulic feedback parameters (in the form of off-line pre-generated cross-section libraries). High accuracy modeling is required for advanced nuclear reactor core designs that present increased geometry complexity and material heterogeneity. Such high-fidelity methods take advantage of the recent progress in computation technology and coupled neutron transport solutions with thermal-hydraulic feedback models on pin or even on sub-pin level (in terms of spatial scale). The continuous energy Monte Carlo method is well suited for solving such core environments with the detailed representation of the complicated 3-D problem. The major advantages of the Monte Carlo method over the deterministic methods are the continuous energy treatment and the exact 3-D geometry modeling. However, the Monte Carlo method involves vast computational time. The interest in Monte Carlo methods has increased thanks to the improvements of the capabilities of high performance computers. Coupled Monte-Carlo calculations can serve as reference solutions for verifying high-fidelity coupled deterministic neutron transport methods
Oxide-based protonic conductors: Point defects and transport properties
Bonanos, N.
, hydrogen pumps, fuel cells, etc. The extent to which protonic defects form depends mainly on the partial pressure of water vapour, temperature and basicity of the constituent oxides, while their mobility depends, among other factors, on the metal-oxygen bond length and bond energy. The defect equilibria...... that determine the protonic concentrations are considered, with emphasis on the regime of low oxygen partial pressure. The measurement of the thermoelectric power (TEP) and of the H+/D+ isotope effect in conductivity are discussed as a means of characterising the conduction process. (C) 2001 Elsevier...
Oxide-based protonic conductors: Point defects and transport properties
Bonanos, N.
2001-01-01
, hydrogen pumps, fuel cells, etc. The extent to which protonic defects form depends mainly on the partial pressure of water vapour, temperature and basicity of the constituent oxides, while their mobility depends, among other factors, on the metal-oxygen bond length and bond energy. The defect equilibria...... that determine the protonic concentrations are considered, with emphasis on the regime of low oxygen partial pressure. The measurement of the thermoelectric power (TEP) and of the H+/D+ isotope effect in conductivity are discussed as a means of characterising the conduction process. (C) 2001 Elsevier...
Comparing analytical and Monte Carlo optical diffusion models in phosphor-based X-ray detectors
Kalyvas, N.; Liaparinos, P.
2014-03-01
Luminescent materials are employed as radiation to light converters in detectors of medical imaging systems, often referred to as phosphor screens. Several processes affect the light transfer properties of phosphors. Amongst the most important is the interaction of light. Light attenuation (absorption and scattering) can be described either through "diffusion" theory in theoretical models or "quantum" theory in Monte Carlo methods. Although analytical methods, based on photon diffusion equations, have been preferentially employed to investigate optical diffusion in the past, Monte Carlo simulation models can overcome several of the analytical modelling assumptions. The present study aimed to compare both methodologies and investigate the dependence of the analytical model optical parameters as a function of particle size. It was found that the optical photon attenuation coefficients calculated by analytical modeling are decreased with respect to the particle size (in the region 1- 12 μm). In addition, for particles sizes smaller than 6μm there is no simultaneous agreement between the theoretical modulation transfer function and light escape values with respect to the Monte Carlo data.
Jones, Kevin C.; Solberg, Timothy D.; Avery, Stephen, E-mail: Stephen.Avery@uphs.upenn.edu [Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania 19104 (United States); Vander Stappen, François; Janssens, Guillaume; Prieels, Damien [Ion Beam Applications SA, Louvain-la-Neuve 1348 (Belgium); Bawiec, Christopher R.; Lewin, Peter A. [School of Biomedical Engineering, Drexel University, Philadelphia, Pennsylvania 19104 (United States); Sehgal, Chandra M. [Department of Radiology, University of Pennsylvania, Philadelphia, Pennsylvania 19104 (United States)
2015-12-15
Purpose: To measure the acoustic signal generated by a pulsed proton spill from a hospital-based clinical cyclotron. Methods: An electronic function generator modulated the IBA C230 isochronous cyclotron to create a pulsed proton beam. The acoustic emissions generated by the proton beam were measured in water using a hydrophone. The acoustic measurements were repeated with increasing proton current and increasing distance between detector and beam. Results: The cyclotron generated proton spills with rise times of 18 μs and a maximum measured instantaneous proton current of 790 nA. Acoustic emissions generated by the proton energy deposition were measured to be on the order of mPa. The origin of the acoustic wave was identified as the proton beam based on the correlation between acoustic emission arrival time and distance between the hydrophone and proton beam. The acoustic frequency spectrum peaked at 10 kHz, and the acoustic pressure amplitude increased monotonically with increasing proton current. Conclusions: The authors report the first observation of acoustic emissions generated by a proton beam from a hospital-based clinical cyclotron. When modulated by an electronic function generator, the cyclotron is capable of creating proton spills with fast rise times (18 μs) and high instantaneous currents (790 nA). Measurements of the proton-generated acoustic emissions in a clinical setting may provide a method for in vivo proton range verification and patient monitoring.
Purpose: To measure the acoustic signal generated by a pulsed proton spill from a hospital-based clinical cyclotron. Methods: An electronic function generator modulated the IBA C230 isochronous cyclotron to create a pulsed proton beam. The acoustic emissions generated by the proton beam were measured in water using a hydrophone. The acoustic measurements were repeated with increasing proton current and increasing distance between detector and beam. Results: The cyclotron generated proton spills with rise times of 18 μs and a maximum measured instantaneous proton current of 790 nA. Acoustic emissions generated by the proton energy deposition were measured to be on the order of mPa. The origin of the acoustic wave was identified as the proton beam based on the correlation between acoustic emission arrival time and distance between the hydrophone and proton beam. The acoustic frequency spectrum peaked at 10 kHz, and the acoustic pressure amplitude increased monotonically with increasing proton current. Conclusions: The authors report the first observation of acoustic emissions generated by a proton beam from a hospital-based clinical cyclotron. When modulated by an electronic function generator, the cyclotron is capable of creating proton spills with fast rise times (18 μs) and high instantaneous currents (790 nA). Measurements of the proton-generated acoustic emissions in a clinical setting may provide a method for in vivo proton range verification and patient monitoring
Fast GPU-based Monte Carlo simulations for LDR prostate brachytherapy
Bonenfant, Éric; Magnoux, Vincent; Hissoiny, Sami; Ozell, Benoît; Beaulieu, Luc; Després, Philippe
2015-07-01
The aim of this study was to evaluate the potential of bGPUMCD, a Monte Carlo algorithm executed on Graphics Processing Units (GPUs), for fast dose calculations in permanent prostate implant dosimetry. It also aimed to validate a low dose rate brachytherapy source in terms of TG-43 metrics and to use this source to compute dose distributions for permanent prostate implant in very short times. The physics of bGPUMCD was reviewed and extended to include Rayleigh scattering and fluorescence from photoelectric interactions for all materials involved. The radial and anisotropy functions were obtained for the Nucletron SelectSeed in TG-43 conditions. These functions were compared to those found in the MD Anderson Imaging and Radiation Oncology Core brachytherapy source registry which are considered the TG-43 reference values. After appropriate calibration of the source, permanent prostate implant dose distributions were calculated for four patients and compared to an already validated Geant4 algorithm. The radial function calculated from bGPUMCD showed excellent agreement (differences within 1.3%) with TG-43 accepted values. The anisotropy functions at r = 1 cm and r = 4 cm were within 2% of TG-43 values for angles over 17.5°. For permanent prostate implants, Monte Carlo-based dose distributions with a statistical uncertainty of 1% or less for the target volume were obtained in 30 s or less for 1 × 1 × 1 mm3 calculation grids. Dosimetric indices were very similar (within 2.7%) to those obtained with a validated, independent Monte Carlo code (Geant4) performing the calculations for the same cases in a much longer time (tens of minutes to more than a hour). bGPUMCD is a promising code that lets envision the use of Monte Carlo techniques in a clinical environment, with sub-minute execution times on a standard workstation. Future work will explore the use of this code with an inverse planning method to provide a complete Monte Carlo-based planning solution.
State-of-the-art in Comprehensive Cascade Control Approach through Monte-Carlo Based Representation
A.H. Mazinan
2015-10-01
Full Text Available The research relies on the comprehensive cascade control approach to be developed in the area of spacecraft, as long as Monte-Carlo based representation is taken into real consideration with respect to state-of-the-art. It is obvious that the conventional methods do not have sufficient merit to be able to deal with such a process under control, constantly, provided that a number of system parameters variations are to be used in providing real situations. It is to note that the new insights in the area of the research’s topic are valuable to outperform a class of spacecrafts performance as the realizations of the acquired results are to be addressed in both real and academic environments. In a word, there are a combination of double closed loop based upon quaternion based control approach in connection with Euler based control approach to handle the three-axis rotational angles and its rates, synchronously, in association with pulse modulation analysis and control allocation, where the dynamics and kinematics of the present system under control are analyzed. A series of experiments are carried out to consider the approach performance in which the aforementioned Monte-Carlo based representation is to be realized in verifying the investigated outcomes.
Design study of Be-target for proton accelerator based neutron source with 13MeV cyclotron
There is a cyclotron named KIRAMS-13 in Pusan National University, Busan, Korea, which has the proton energy of 13MeV and the beam current of 0.05mA. Originally, it was developed for producing medical radioisotopes and nuclear physics research. To improve the utilization of the facility, we are considering the possibilities of installing a neutron generation target in it. The Beryllium target has been considered and neutrons can be generated by 9Be(p,n)9B reaction above the threshold proton energy of 2.057MeV. In this presentation, we suggest candidate materials and structures, thicknesses, metal layers and cooling systems of target, which is optimal for the KIRAMS-13. We chose the Beryllium material of 1.14mm thick, which is calculated by stopping power of Beryllium, based on PSTAR, NIST. As for the cooling system, we chose to use water as a coolant, which will also act as a moderator. As protons pass through the target, hydrogen ions continue to pile up in the material and this makes the material brittle. To solve this problem, we chose Vanadium material because it has high hydrogen diffusion coefficient and short half-life isotope after being activated by neutrons. We simulated the neutron characteristics by the Monte Carlo simulation code, Geant4, CERN and performed thermal analysis on the target. The design of target system is very important to produce neutrons for the desired purposes. There are several other existing facilities in Korea, in addition to the cyclotron facility considered in this study, where new neutron target system can be installed and neutrons can be generated. Two prominent facilities are KOMAC, Gyeongju and RFT-30, Jeongeup and we are planning to do study on the possibilities of utilizing the accelerators for neutron generation.
Proton conductive membranes based on doped sulfonated polytriazole
Boaventura, M.; Brandao, L.; Mendes, A. [Laboratorio de Engenharia de Processos, Ambiente e Energia (LEPAE), Faculdade de Engenharia da Universidade do Porto, Rua Roberto Frias, 4200-465 Porto (Portugal); Ponce, M.L.; Nunes, S.P. [GKSS Research Centre Geesthacht GmbH, Max Planck Str. 1, D-21502, Geesthacht (Germany)
2010-11-15
This work reports the preparation and characterization of proton conducting sulfonated polytriazole membranes doped with three different agents: 1H-benzimidazole-2-sulfonic acid, benzimidazole and phosphoric acid. The modified membranes were characterized by scanning electron microscopy (SEM), infrared spectra, thermogravimetric analysis (TGA), dynamical mechanical thermal analysis (DMTA) and electrochemical impedance spectroscopy (EIS). The addition of doping agents resulted in a decrease of the glass transition temperature. For membranes doped with 85 wt.% phosphoric acid solution proton conductivity increased up to 2.10{sup -3} S cm{sup -1} at 120 C and at 5% relative humidity. The performance of the phosphoric acid doped membranes was evaluated in a fuel cell set-up at 120 C and 2.5% relative humidity. (author)
Channel capacity of ocean water is limited by propagation distance and optical properties. Previous studies on this problem are based on water-tank experiments with different amounts of Maalox antacid. However, propagation distance is limited by the experimental set-up and the optical properties are different from ocean water. Therefore, the experiment result is not accurate for the physical design of underwater wireless communications links. This letter developed a Monte Carlo model to study channel capacity of underwater optical communications. Moreover, this model can flexibly configure various parameters of transmitter, receiver and channel, and is suitable for physical underwater optical communications links design. (paper)
A new Monte-Carlo based simulation for the CryoEDM experiment
Raso-Barnett, Matthew
2015-01-01
This thesis presents a new Monte-Carlo based simulation of the physics of ultra-cold neutrons (UCN) in complex geometries and its application to the CryoEDM experiment. It includes a detailed description of the design and performance of this simulation along with its use in a project to study the magnetic depolarisation time of UCN within the apparatus due to magnetic impurities in the measurement cell, which is a crucial parameter in the sensitivity of a neutron electricdipole-moment (nEDM) ...
CARMEN: a system Monte Carlo based on linear programming from direct openings
The use of Monte Carlo (MC) has shown an improvement in the accuracy of the calculation of the dose compared to other analytics algorithms installed on the systems of business planning, especially in the case of non-standard situations typical of complex techniques such as IMRT and VMAT. Our treatment planning system called CARMEN, is based on the complete simulation, both the beam transport in the head of the accelerator and the patient, and simulation designed for efficient operation in terms of the accuracy of the estimate and the required computation times. (Author)
Krosigk, B. von; Zuber, K. [Technische Universitaet Dresden, Institut fuer Kern- und Teilchenphysik, Dresden (Germany); Neumann, L. [Technische Universitaet Dresden, Institut fuer Kern- und Teilchenphysik, Dresden (Germany); Karlsruher Institut fuer Technologie, Institut fuer Experimentelle Kernphysik, Karlsruhe (Germany); Nolte, R.; Roettger, S. [Physikalisch-Technische Bundesanstalt (PTB), Braunschweig (Germany)
2013-04-15
The proton light output function in electron-equivalent energy of various scintillators based on linear alkylbenzene (LAB) has been measured in the energy range from 1 MeV to 17.15 MeV for the first time. The measurement was performed at the Physikalisch-Technische Bundesanstalt (PTB) using a neutron beam with continuous energy distribution. The proton light output data is extracted from proton recoil spectra originating from neutron-proton scattering in the scintillator. The functional behavior of the proton light output is described successfully by Birks' law with a Birks constant kB between (0.0094{+-}0.0002) cm MeV{sup -1} and (0.0098{+-}0.0003) cm MeV{sup -1} for the different LAB solutions. The constant C, parameterizing the quadratic term in the generalized Birks law, is consistent with zero for all investigated scintillators with an upper limit (95 % CL) of about 10{sup -7} cm{sup 2} MeV{sup -2}. The resulting quenching factors are especially important for future planned supernova neutrino detection based on the elastic scattering of neutrinos on protons. The impact of proton quenching on the supernova event yield from neutrino-proton scattering is discussed. (orig.)