Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Commissioning a CT-compatible LDR tandem and ovoid applicator using Monte Carlo calculation and 3D dosimetry

Purpose: To determine the geometric and dose attenuation characteristics of a new commercially available CT-compatible LDR tandem and ovoid (T and O) applicator using Monte Carlo calculation and 3D dosimetry. Methods: For geometric characterization, we quantified physical dimensions and investigated a systematic difference found to exist between nominal ovoid angle and the angle at which the afterloading buckets fall within the ovoid. For dosimetric characterization, we determined source attenuation through asymmetric gold shielding in the buckets using Monte Carlo simulations and 3D dosimetry. Monte Carlo code MCNP5 was used to simulate 1.5 Multiplication-Sign 10{sup 9} photon histories from a {sup 137}Cs source placed in the bucket to achieve statistical uncertainty of 1% at a 6 cm distance. For 3D dosimetry, the distribution about an unshielded source was first measured to evaluate the system for {sup 137}Cs, after which the distribution was measured about sources placed in each bucket. Cylindrical PRESAGE{sup Registered-Sign} dosimeters (9.5 cm diameter, 9.2 cm height) with a central channel bored for source placement were supplied by Heuris Inc. The dosimeters were scanned with the Duke Large field of view Optical CT-Scanner before and after delivering a nominal dose at 1 cm of 5-8 Gy. During irradiation the dosimeter was placed in a water phantom to provide backscatter. Optical CT scan time lasted 15 min during which 720 projections were acquired at 0.5 Degree-Sign increments, and a 3D distribution was reconstructed with a (0.05 cm){sup 3} isotropic voxel size. The distributions about the buckets were used to calculate a 3D distribution of transmission rate through the bucket, which was applied to a clinical CT-based T and O implant plan. Results: The systematic difference in bucket angle relative to the nominal ovoid angle (105 Degree-Sign ) was 3.1 Degree-Sign -4.7 Degree-Sign . A systematic difference in bucket angle of 1 Degree-Sign , 5 Degree-Sign , and 10 Degree-Sign caused a 1%{+-} 0.1%, 1.7%{+-} 0.4%, and 2.6%{+-} 0.7% increase in rectal dose, respectively, with smaller effect to dose to Point A, bladder, sigmoid, and bowel. For 3D dosimetry, 90.6% of voxels had a 3D {gamma}-index (criteria = 0.1 cm, 3% local signal) below 1.0 when comparing measured and expected dose about the unshielded source. Dose transmission through the gold shielding at a radial distance of 1 cm was 85.9%{+-} 0.2%, 83.4%{+-} 0.7%, and 82.5%{+-} 2.2% for Monte Carlo, and measurement for left and right buckets, respectively. Dose transmission was lowest at oblique angles from the bucket with a minimum of 56.7%{+-} 0.8%, 65.6%{+-} 1.7%, and 57.5%{+-} 1.6%, respectively. For a clinical T and O plan, attenuation from the buckets leads to a decrease in average Point A dose of {approx}3.2% and decrease in D{sub 2cc} to bladder, rectum, bowel, and sigmoid of 5%, 18%, 6%, and 12%, respectively. Conclusions: Differences between dummy and afterloading bucket position in the ovoids is minor compared to effects from asymmetric ovoid shielding, for which rectal dose is most affected. 3D dosimetry can fulfill a novel role in verifying Monte Carlo calculations of complex dose distributions as are common about brachytherapy sources and applicators.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-T-303: Spot Scanning Dose Delivery with Rapid Cycling Proton Beams From RCMS

Purpose: A rapid cycling proton beam has several distinct characteristics superior to a slow extraction synchrotron: The beam energy and energy spread, beam intensity and spot size can be varied spot by spot. The feasibility of using a spot scanning beam from a rapidc-ycling-medical-synchrotron (RCMS) at 10 Hz repetition frequency is investigated in this study for its application in proton therapy. Methods: The versatility of the beam is illustrated by two examples in water phantoms: (1) a cylindrical PTV irradiated by a single field and (2) a spherical PTV irradiated by two parallel opposed fields. A uniform dose distribution is to be delivered to the volumes. Geant4 Monte Carlo code is used to validate the dose distributions in each example. Results: Transverse algorithms are developed to produce uniform distributions in each transverseplane in the two examples with a cylindrical and a spherical PTV respectively. Longitudinally, different proton energies are used in successive transverse planes toproduce the SOBP required to cover the PTVs. In general, uniformity of dosedistribution within 3% is obtained for the cylinder and 3.5% for the sphere. The transversealgorithms requires only few hundred beam spots for each plane The algorithms may beapplied to larger volumes by increasing the intensity spot by spot for the same deliverytime of the same dose. The treatment time can be shorter than 1 minute for any fieldconfiguration and tumor shape. Conclusion: The unique beam characteristics of a spot scanning beam from a RCMS at 10 Hz repetitionfrequency are used to design transverse and longitudinal algorithms to produce uniformdistribution for any arbitrary shape and size of targets. The proposed spot scanning beam ismore versatile than existing spot scanning beams in proton therapy with better beamcontrol and lower neutron dose. This work is supported in part by grants from the US Department of Energy under contract; DE-FG02-12ER41800 and the National Science Foundation NSF PHY-1205431.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

SU-F-T-360: Dosimetric Impacts On the Mucosa and Bone in Radiotherapy with Unflattened Photon Beams

Purpose: This study investigated the dosimetric impacts on the mucosa and bone when using the unflattened photon beams in radiotherapy. Dose calculations were carried out by Monte Carlo simulation. Methods: Heterogeneous phantoms containing water (soft tissue and mucosa), air and bone, with mucosa thicknesses varying from 0.5 – 3 mm were irradiated by the 6 MV unflattened and flattened photon beams (field size = 10 × 10 cm{sup 2}), produced by a Varian TrueBEAM linear accelerator. The photon energy spectra of the beams, mean bone and mucosal doses with different mucosa thicknesses were calculated using the EGSnrc Monte Carlo code. Results: It is found that the flattened photon beams had higher mean bone doses (1.3% and 2% for upper and lower bone regarding the phantom geometry, respectively) than the unflattened beams, and the mean bone doses of both beams did not vary significantly with the mucosa thickness. Similarly, flattened photon beams had higher mucosal dose (0.9% and 1.6% for upper and lower mucosa, respectively) than the unflattened beams. This is due to the larger slope of the depth dose for the unflattened photon beams compared to the flattened. The mucosal doses of both beams were found increased with the mucosa thickness. Moreover, the mucosal dose differences between the unflattened and flattened beams increased with the mucosa thickness. For photon energy spectra on the mucosal layers, it is found that the unflattened photon beams contained a larger portion of lowenergy photons than the flattened beams. The photon energy spectra did not change significantly with the mucosa thickness. Conclusion: It is concluded that the mucosal and bone dose for the unflattened photon beams were not more than 2% lower than the flattened beams, though the flattening filter free beams contained larger portion of low-energy photons than the flattened beams.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Evaluation of S-values and dose distributions for {sup 90}Y, {sup 131}I, {sup 166}Ho, and {sup 188}Re in seven lobes of the rat liver

Purpose: Rats have been widely used in radionuclide therapy research for the treatment of hepatocellular carcinoma (HCC). This has created the need to assess rat liver absorbed radiation dose. In most dose estimation studies, the rat liver is considered as a homogeneous integrated target organ with a tissue composition assumed to be similar to that of human liver tissue. However, the rat liver is composed of several lobes having different anatomical and chemical characteristics. To assess the overall impact on rat liver dose calculation, the authors use a new voxel-based rat model with identified suborgan regions of the liver. Methods: The liver in the original cryosectional color images was manually segmented into seven individual lobes and subsequently integrated into a voxel-based computational rat model. Photon and electron particle transport was simulated using the MCNPX Monte Carlo code to calculate absorbed fractions and S-values for {sup 90}Y, {sup 131}I, {sup 166}Ho, and {sup 188}Re for the seven liver lobes. The effect of chemical composition on organ-specific absorbed dose was investigated by changing the chemical composition of the voxel filling liver material. Radionuclide-specific absorbed doses at the voxel level were further assessed for a small spherical hepatic tumor. Results: The self-absorbed dose for different liver lobes varied depending on their respective masses. A maximum difference of 3.5% was observed for the liver self-absorbed fraction between rat and human tissues for photon energies below 100 keV. {sup 166}Ho and {sup 188}Re produce a uniformly distributed high dose in the tumor and relatively low absorbed dose for surrounding tissues. Conclusions: The authors evaluated rat liver radiation doses from various radionuclides used in HCC treatments using a realistic computational rat model. This work contributes to a better understanding of all aspects influencing radiation transport in organ-specific radiation dose evaluation for preclinical therapy studies, from tissue composition to organ morphology and activity distribution.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-T-605: A New Design for a Rotating Gamma Knife. Monte Carlo Simulation

Purpose: to determine the characteristics of the 60Co beam emerging from a new design of a rotating Gamma Knife system and to calculate 3D dose distributions at the isocenter for different source configurations and collimator openings. Methods: We employed the BEAM-Monte Carlo code to realistically model the geometry design, including 30 60Co source capsules, two circular primary collimators (diameter of 6.6mm and 6.1mm) and four different changeable collimators. The shielding of the head design was also simulated. The sources (2.8mm diameter) are distributed in six groups in the spherical geometry. Each source is individually collimated to obtained four different circular fields at the isocenter (3mm, 3.5mm, 6mm and 8mm).The phase-space particles reaching the scoring plane below the primary collimation assembly were recorded and the BEAMDP code was used to determine the fluence and energy spectra of the particles emerging from each source-collimator configuration. The dose distributions at the isocenter plane (397.6mm from the source) were calculated in a spherical component module for the circular field sizes studied. Results: The energy spectra below the head assembly and primary collimators have been obtained, which exhibited the typical 60Co peaks and a small low-energy tail due to scattered photons (from about 200keV to 1MeV). The scattered component of the spectra represents about 8 % of the total number of photons reaching the scoring plane. The radial photon fluence does not vary significantly inside the collimator openings. The spectra of particles from different source groups are compared. Conclusion: The 60Co beam emerging from each source configuration was characterized, which can be used to establish a generic source model for all the sources for fast MC dose calculation. Further investigations are needed to determine the dose variations as a result of partial switching on/off different groups of sources for advanced Gamma Knife SRS/SBRT planning.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

SU-E-T-409: Evaluation of Tissue Composition Effect On Dose Distribution in Radiotherapy with 6 MV Photon Beam of a Medical Linac

Purpose: The aim of this study is to evaluate soft tissue composition effect on dose distribution for various soft tissues and various depths in radiotherapy with 6 MV photon beam of a medical linac. Methods: A phantom and Siemens Primus linear accelerator were simulated using MCNPX Monte Carlo code. In a homogeneous cubic phantom, six types of soft tissue and three types of tissue-equivalent materials were defined separately. The soft tissues were muscle (skeletal), adipose tissue, blood (whole), breast tissue, soft tissue (9-component) and soft tissue (4-component). The tissue-equivalent materials included: water, A-150 tissue-equivalent plastic and perspex. Photon dose relative to dose in 9-component soft tissue at various depths on the beam’s central axis was determined for the 6 MV photon beam. The relative dose was also calculated and compared for various MCNPX tallies including,F8, F6 and,F4. Results: The results of the relative photon dose in various materials relative to dose in 9-component soft tissue and using different tallies are reported in the form of tabulated data. Minor differences between dose distributions in various soft tissues and tissue-equivalent materials were observed. The results from F6 and F4 were practically the same but different with,F8 tally. Conclusion: Based on the calculations performed, the differences in dose distributions in various soft tissues and tissue-equivalent materials are minor but they could be corrected in radiotherapy calculations to upgrade the accuracy of the dosimetric calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MO-F-CAMPUS-J-04: Radiation Heat Load On the MR System of the Elekta Atlantic System

Purpose: The Elekta Atlantic system combines a digital linear accelerator system with a 1.5T Philips MRI machine.This study aimed to assess the energy deposited within the cryostat system when the radiation beam passes through the cryostat. The cryocooler on the magnet has a cooling capacity which is about 1 Watt in excess of the cryogenic heat leak into the magnet’s cold mass. A pressure-controlled heater inside the magnet balances the excess refrigeration power such that the helium pressure in the tank is kept slightly above ambient air pressure. If radiation power is deposited in the cold mass then this heater will need less power to maintain pressure equilibrium and if the radiation heat load exceeds the excess cryocooler capacity the pressure will rise. Methods: An in-house CAD based Monte Carlo code based on Penelope was used to model the entire MR-Linac system to quantify the heat load on the magnet’s cold mass. These results were then compared to experimental results obtained from an Elekta Atlantic system installed in UMC-Utrecht. Results: For a field size of 25 cm x 22 cm and a dose rate of 107 mu.min-1, the energy deposited by the radiation beam led to a reduction in heater power from 1.16 to 0.73 W. Simulations predicted a reduction to 0.69 W which is in good agreement. For the worst case field size (largest) and maximum dose rate the cryostat cooler capacity was exceeded. This resulted in a pressure rise within the system but was such that continuous irradiation for over 12 hours would be required before the magnet would start blowing off helium. Conclusion: The study concluded that the Atlantic system does not have to be duty cycle restricted, even for the worst case non-clinical scenario and that there are no adverse effects on the MR system. Stephen Towe and David Roberts Both work for Elekta; Ezra Van Lanen works for Philips Healthcare; Johan Overweg works for Philips Innovative Technologies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-EF-207-11: Energy and Depth Response of Thermoluminescent Dosimeters and Ionization Chambers in Water for Kilovoltage X-Ray Beams

Purpose: To assess the effects of changes in beam quality on detector response in the kilovoltage energy range by modulating the x-ray tube voltage and the measurement depth in water. Methods: Measurements were performed with TLD-100 and TLD-100H thermoluminescent dosimeters and an A12 farmer-type ionization chamber. To assess the energy response of the detectors, irradiations were performed at a depth of 3 cm in a custom-built thin-window water phantom using the moderately filtered x-ray beams at the UWADCL (20 kVp-250 kVp) and a Co-60 beam.The x-ray beams and detectors were modeled using the EGSnrc Monte Carlo code. The model was validated by simulating dose to the collecting volume of an A12 farmer chamber and comparing it with measured A12 signal as a function of depth. Dose was tallied to each detector and to water for comparison with measurements. Simulations were used to calculate the predicted energy response, which was compared to the measured response of each detector. Dose to each detector and dose to water as a function of depth were also simulated. Results: Detector output per dose to water was found to deviate by up to 15%, 20% and 30% as a function of energy relative to Co-60 for the A12, TLD-100H and TLD-100, respectively. The EGSnrc simulations produced results similar to the measurements for ionization chambers, but discrepancies of up to 30% were observed for TLD-100H. Simulated detector response as a function of depth was found to vary by up to 3%. Conclusion: These results suggest that changes in beam quality in kilovoltage x-ray beams can have a significant impact on detector response. In-water detector response was found to differ from the previously investigated in-air response. Deviations in detector response as a function of depth were less significant, but could potentially cause dosimetric errors if ignored.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Reoptimization of Intensity Modulated Proton Therapy Plans Based on Linear Energy Transfer

Purpose: We describe a treatment plan optimization method for intensity modulated proton therapy (IMPT) that avoids high values of linear energy transfer (LET) in critical structures located within or near the target volume while limiting degradation of the best possible physical dose distribution. Methods and Materials: To allow fast optimization based on dose and LET, a GPU-based Monte Carlo code was extended to provide dose-averaged LET in addition to dose for all pencil beams. After optimizing an initial IMPT plan based on physical dose, a prioritized optimization scheme is used to modify the LET distribution while constraining the physical dose objectives to values close to the initial plan. The LET optimization step is performed based on objective functions evaluated for the product of LET and physical dose (LET×D). To first approximation, LET×D represents a measure of the additional biological dose that is caused by high LET. Results: The method is effective for treatments where serial critical structures with maximum dose constraints are located within or near the target. We report on 5 patients with intracranial tumors (high-grade meningiomas, base-of-skull chordomas, ependymomas) in whom the target volume overlaps with the brainstem and optic structures. In all cases, high LET×D in critical structures could be avoided while minimally compromising physical dose planning objectives. Conclusion: LET-based reoptimization of IMPT plans represents a pragmatic approach to bridge the gap between purely physical dose-based and relative biological effectiveness (RBE)-based planning. The method makes IMPT treatments safer by mitigating a potentially increased risk of side effects resulting from elevated RBE of proton beams near the end of range.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

SU-G-201-08: Energy Response of Thermoluminescent Microcube Dosimeters in Water for Kilovoltage X-Ray Beams

Purpose: To characterize the energy dependence for TLD-100 microcubes in water at kilovoltage energies. Methods: TLD-100 microcubes with dimensions of (1 × 1 × 1) mm{sup 3} were irradiated with kilovoltage x-rays in a custom-built thin-window liquid water phantom. The TLD-100 microcubes were held in Virtual Water™ probes and aligned at a 2 cm depth in water. Irradiations were performed using the M-series x-ray beams of energies ranging from 50-250 kVp and normalized to a {sup 60}Co beam located at the UWADCL. Simulations using the EGSnrc Monte Carlo Code System were performed to model the x-ray beams, the {sup 60}Co beam, the water phantom and the dosimeters in the phantom. The egs-chamber user code was used to tally the dose to the TLDs and the dose to water. The measurements and calculations were used to determine the intrinsic energy dependence, absorbed-dose energy dependence, and absorbed-dose sensitivity. These values were compared to TLD-100 chips with dimensions of (3.2 × 0.9 × 0.9) mm{sup 3}. Results: The measured TLD-100 microcube response per dose to water among all investigated x-ray energies had a maximum percent difference of 61% relative to {sup 60}Co. The simulated ratio of dose to water to the dose to TLD had a maximum percent difference of 29% relative to {sup 60}Co. The ratio of dose to TLD to the TLD output had a maximum percent difference of 13% relative to {sup 60}Co. The maximum percent difference for the absorbed-dose sensitivity was 15% more than the used value of 1.41. Conclusion: These results confirm that differences in beam quality have a significant effect on TLD response when irradiated in water. These results also indicated a difference in TLD-100 response between microcube and chip geometries. The intrinsic energy dependence and the absorbed-dose energy dependence deviated up to 10% between TLD-100 microcubes and chips.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

SU-G-TeP3-13: The Role of Nanoscale Energy Deposition in the Development of Gold Nanoparticle-Enhanced Radiotherapy

Purpose: Gold nanoparticles (GNPs) can enhance radiotherapy effects. The high photoelectric cross section of gold relative to tissue, particularly at lower energies, leads to localized dose enhancement. However in a clinical context, photon energies must also be sufficient to reach a target volume at a given depth. These properties must be balanced to optimize such a therapy. Given that nanoscale energy deposition patterns around GNPs play a role in determining biological outcomes, in this work we seek to establish their role in this optimization process. Methods: The PENELOPE Monte Carlo code was used to generate spherical dose deposition kernels in 1000 nm diameter spheres around 50 nm diameter GNPs in response to monoenergetic photons incident on the GNP. Induced “lesions” were estimated by either a local effect model (LEM) or a mean dose model (MDM). The ratio of these estimates was examined for a range of photon energies (10 keV to 2 MeV), for three sets of linear-quadratic parameters. Results: The models produce distinct differences in expected lesion values, the lower the alpha-beta ratio, the greater the difference. The ratio of expected lesion values remained constant within 5% for energies of 40 keV and above across all parameter sets and rose to a difference of 35% for lower energies only for the lowest alpha-beta ratio. Conclusion: Consistent with other work, these calculations suggest nanoscale energy deposition patterns matter in predicting biological response to GNP-enhanced radiotherapy. However the ratio of expected lesions between the different models is largely independent of energy, indicating that GNP-enhanced radiotherapy scenarios can be optimized in photon energy without consideration of the nanoscale patterns. Special attention may be warranted for energies of 20 keV or below and low alpha-beta ratios.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Monte Carlo simulation for the electron cascade due to gamma rays in semiconductor radiation detectors

A Monte Carlo code was developed for simulating the electron cascade in radiation detector materials. The electron differential scattering cross sections were derived from measured electron energy-loss and optical spectra, making the method applicable for a wide range of materials. The detector resolution in a simplified model system shows dependence on the bandgap, the plasmon strength and energy, and the valence band width. In principle, these parameters could be optimized to improve detector performance. The intrinsic energy resolution was calculated for three semiconductors: silicon (Si), gallium arsenide (GaAs), and zinc telluride (ZnTe). Setting the ionization thresholds for electrons and holes is identified as a critical issue, as this strongly affects both the average electron-hole pair energy w and the Fano factor F. Using an ionization threshold from impact ionization calculations as an effective bandgap yields pair energies that are well matched to measured values. Fano factors of 0.091 (Si), 0.100 (GaAs), and 0.075 (ZnTe) were calculated. The Fano factor calculated for silicon using this model was lower than some results from past simulations and experiments. This difference could be attributed to problems in simulating inter-band transitions and the scattering of low-energy electrons.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

SU-E-T-748: Theoretical Investigation On Using High Energy Proton Beam for Total-Body-Irradiation

Purpose: The broad-slow-rising entrance dose region proximal to the Bragg peak made by a mono-energetic proton beam could potentially be used for total body irradiation (TBI). Due to the quasi-uniform dose deposition, customized thickness compensation may not be required to deliver a uniform dose to patients with varied thickness. We investigated the possibility, efficacy, and hardware requirement to use such proton beam for TBI. Methods: A wedge shaped water phantom with thickness varying from 2 cm to 40 cm was designed to mimic a patient. Geant4 based Monte Carlo code was used to simulate broad mono-energetic proton beams with energy ranging from 250 MeV to 300 MeV radiating the phantom. A 6 MV photon with 1 cm water equivalent build-up used for conventional TBI was also calculated. A paired-opposing beam arrangement with no thickness compensation was used to generate TBI plans for all beam energies. Dose from all particles were scored on a grid size of 2 mm{sup 3}. Dose uniformity across the phantom was calculated to evaluate the plan. The field size limit and the dose uniformity of Mevion S250 proton system was examined by using radiochromic films placed at extended treatment distance with the open large applicator and 90° gantry angle. Results: To achieve a maximum ± 7.5% dose variation, the largest patient thickness variation allowed for 250 MeV, 275 MeV, and 300 MeV proton beams were 27.0 cm, 34.9 cm and 36.7 cm. The value for 6 MV photon beam was only 8.0 cm to achieve the same dose variation. With open gantry, Mevion S250 system allows 5 m source-to-surface distance producing an expected 70 cm{sup 2} field size. Conclusion: Energetic proton beam can potentially be used to deliver TBI. Treatment planning and delivery would be much simple since no thickness compensation is required to achieve a uniform dose distribution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Detailed high-accuracy megavoltage transmission measurements: A sensitive experimental benchmark of EGSnrc

Purpose: There are three goals for this study: (a) to perform detailed megavoltage transmission measurements in order to identify the factors that affect the measurement accuracy, (b) to use the measured data as a benchmark for the EGSnrc system in order to identify the computational limiting factors, and (c) to provide data for others to benchmark Monte Carlo codes. Methods: Transmission measurements are performed at the National Research Council Canada on a research linac whose incident electron parameters are independently known. Automated transmission measurements are made on-axis, down to a transmission value of {approx}1.7%, for eight beams between 10 MV (the lowest stable MV beam on the linac) and 30 MV, using fully stopping Be, Al, and Pb bremsstrahlung targets and no fattening filters. To diversify energy differentiation, data are acquired for each beam using low-Z and high-Z attenuators (C and Pb) and Farmer chambers with low-Z and high-Z buildup caps. Experimental corrections are applied for beam drifts (2%), polarity (2.5% typical maximum, 6% extreme), ion recombination (0.2%), leakage (0.3%), and room scatter (0.8%)-the values in parentheses are the largest corrections applied. The experimental setup and the detectors are modeled using EGSnrc, with the newly added photonuclear attenuation included (up to a 5.6% effect). A detailed sensitivity analysis is carried out for the measured and calculated transmission data. Results: The developed experimental protocol allows for transmission measurements with 0.4% uncertainty on the smallest signals. Suggestions for accurate transmission measurements are provided. Measurements and EGSnrc calculations agree typically within 0.2% for the sensitivity of the transmission values to the detector details, to the bremsstrahlung target material, and to the incident electron energy. Direct comparison of the measured and calculated transmission data shows agreement better than 2% for C (3.4% for the 10 MV beam) and typically better than 1% for Pb. The differences can be explained by acceptable photon cross section changes of Less-Than-Or-Slanted-Equal-To 0.4%. Conclusions: Accurate transmission measurements require accounting for a number of influence quantities which, if ignored, can collectively introduce errors larger than 10%. Accurate transmission calculations require the use of the most accurate data and physics options available in EGSnrc, particularly the more accurate bremsstrahlung angular sampling option and the newly added modeling of photonuclear attenuation. Comparison between measurements and calculations implies that EGSnrc is accurate within 0.2% for relative ion chamber response calculations. Photon cross section uncertainties are the ultimate limiting factor for the accuracy of the calculated transmission data (Monte Carlo or analytical).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Improved proton computed tomography by dual modality image reconstruction

Purpose: Proton computed tomography (CT) is a promising image modality for improving the stopping power estimates and dose calculations for particle therapy. However, the finite range of about 33 cm of water of most commercial proton therapy systems limits the sites that can be scanned from a full 360° rotation. In this paper the authors propose a method to overcome the problem using a dual modality reconstruction (DMR) combining the proton data with a cone-beam x-ray prior. Methods: A Catphan 600 phantom was scanned using a cone beam x-ray CT scanner. A digital replica of the phantom was created in the Monte Carlo code Geant4 and a 360° proton CT scan was simulated, storing the entrance and exit position and momentum vector of every proton. Proton CT images were reconstructed using a varying number of angles from the scan. The proton CT images were reconstructed using a constrained nonlinear conjugate gradient algorithm, minimizing total variation and the x-ray CT prior while remaining consistent with the proton projection data. The proton histories were reconstructed along curved cubic-spline paths. Results: The spatial resolution of the cone beam CT prior was retained for the fully sampled case and the 90° interval case, with the MTF = 0.5 (modulation transfer function) ranging from 5.22 to 5.65 linepairs/cm. In the 45° interval case, the MTF = 0.5 dropped to 3.91 linepairs/cm For the fully sampled DMR, the maximal root mean square (RMS) error was 0.006 in units of relative stopping power. For the limited angle cases the maximal RMS error was 0.18, an almost five-fold improvement over the cone beam CT estimate. Conclusions: Dual modality reconstruction yields the high spatial resolution of cone beam x-ray CT while maintaining the improved stopping power estimation of proton CT. In the case of limited angles, the use of prior image proton CT greatly improves the resolution and stopping power estimate, but does not fully achieve the quality of a 360° proton CT scan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-F-303-15: Ion Chamber Dose Response in Magnetic Fields as a Function of Incident Photon Energy

Purpose: In considering the continued development of synergetic MRI-radiation therapy machines, we seek to quantify the variability of ion chamber response per unit dose in the presence of magnetic fields of varying strength as a function of incident photon beam quality and geometric configuration. Methods: To account for the effect of magnetic fields on the trajectory of charged particles a new algorithm was introduced into the EGSnrc Monte Carlo code. In the egs-chamber user code the dose to the cavity of an NE2571 ion chamber is calculated in two configurations, in 0 to 2 T magnetic fields, with an incoming parallel 10×10 cm{sup 2} photon beam with energies ranging between 0.5 MeV and 8 MeV. In the first, the photon beam is incident on the long-axis of the ion chamber (config-1), and in the second the beam is parallel to the long-axis and incident from the conical end of the chamber (config-2). For both, the magnetic field is perpendicular to the direction of the beam and the long axis of the chamber. Results: The ion chamber response per unit dose to water at the same point is determined as a function of magnetic field and is normalized to the 0T case for each of incoming photon energies. For both configurations, accurate modeling of the ion chamber yielded closer agreement with the experimental results obtained by Meijsing et. al (2009). Config-1 yields a gradual increase in response with increasing field strength to a maximum of 13.4% and 1.4% for 1 MeV and 8 MeV photon beams, respectively. Config-2 produced a decrease in response of up to 6% and 13% for 0.5 MeV and 8 MeV beams, respectively. Conclusion: These results provide further support for ion chamber calibration in MRI-radiotherapy coupled systems and demonstrates noticeable energy dependence for clinically relevant fields.

60 APPLIED LIFE SCIENCES↗