Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Concomitant Imaging Dose and Cancer Risk in Image Guided Thoracic Radiation Therapy

Purpose: Kilovoltage cone beam computed tomography (CT) (kVCBCT) imaging guidance improves the accuracy of radiation therapy but imposes an extra radiation dose to cancer patients. This study aimed to investigate concomitant imaging dose and associated cancer risk in image guided thoracic radiation therapy. Methods and Materials: The planning CT images and structure sets of 72 patients were converted to CT phantoms whose chest circumferences (C{sub chest}) were calculated retrospectively. A low-dose thorax protocol on a Varian kVCBCT scanner was simulated by a validated Monte Carlo code. Computed doses to organs and cardiac substructures (for 5 selected patients of various dimensions) were regressed as empirical functions of C{sub chest}, and associated cancer risk was calculated using the published models. The exposures to nonthoracic organs in children were also investigated. Results: The structural mean doses decreased monotonically with increasing C{sub chest}. For all 72 patients, the median doses to the heart, spinal cord, breasts, lungs, and involved chest were 1.68, 1.33, 1.64, 1.62, and 1.58 cGy/scan, respectively. Nonthoracic organs in children received 0.6 to 2.8 cGy/scan if they were directly irradiated. The mean doses to the descending aorta (1.43 ± 0.68 cGy), left atrium (1.55 ± 0.75 cGy), left ventricle (1.68 ± 0.81 cGy), and right ventricle (1.85 ± 0.84 cGy) were significantly different (P<.05) from the heart mean dose (1.73 ± 0.82 cGy). The blade shielding alleviated the exposure to nonthoracic organs in children by an order of magnitude. Conclusions: As functions of patient size, a series of models for personalized estimation of kVCBCT doses to thoracic organs and cardiac substructures have been proposed. Pediatric patients received much higher doses than did the adults, and some nonthoracic organs could be irradiated unexpectedly by the default scanning protocol. Increased cancer risks and disease adverse events in the thorax were strongly related to higher imaging doses and smaller chest dimensions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-T-238: Monte Carlo Estimation of Cerenkov Dose for Photo-Dynamic Radiotherapy

Purpose: Estimation of Cerenkov dose from high-energy megavoltage photon and electron beams in tissue and its impact on the radiosensitization using Protoporphyrine IX (PpIX) for tumor targeting enhancement in radiotherapy. Methods: The GEPTS Monte Carlo code is used to generate dose distributions from 18MV Varian photon beam and generic high-energy (45-MV) photon and (45-MeV) electron beams in a voxel-based tissueequivalent phantom. In addition to calculating the ionization dose, the code scores Cerenkov energy released in the wavelength range 375–425 nm corresponding to the pick of the PpIX absorption spectrum (Fig. 1) using the Frank-Tamm formula. Results: The simulations shows that the produced Cerenkov dose suitable for activating PpIX is 4000 to 5500 times lower than the overall radiation dose for all considered beams (18MV, 45 MV and 45 MeV). These results were contradictory to the recent experimental studies by Axelsson et al. (Med. Phys. 38 (2011) p 4127), where Cerenkov dose was reported to be only two orders of magnitude lower than the radiation dose. Note that our simulation results can be corroborated by a simple model where the Frank and Tamm formula is applied for electrons with 2 MeV/cm stopping power generating Cerenkov photons in the 375–425 nm range and assuming these photons have less than 1mm penetration in tissue. Conclusion: The Cerenkov dose generated by high-energy photon and electron beams may produce minimal clinical effect in comparison with the photon fluence (or dose) commonly used for photo-dynamic therapy. At the present time, it is unclear whether Cerenkov radiation is a significant contributor to the recently observed tumor regression for patients receiving radiotherapy and PpIX versus patients receiving radiotherapy only. The ongoing study will include animal experimentation and investigation of dose rate effects on PpIX response.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

STELLAR COLLISIONS AND BLUE STRAGGLER STARS IN DENSE GLOBULAR CLUSTERS

Blue straggler stars (BSSs) are abundantly observed in all Galactic globular clusters (GGCs) where data exist. However, observations alone cannot reveal the relative importance of various formation channels or the typical formation times for this well-studied population of anomalous stars. Using a state-of-the-art Hénon-type Monte Carlo code that includes all relevant physical processes, we create 128 models with properties typical of the observed GGCs. These models include realistic numbers of single and binary stars, use observationally motivated initial conditions, and span large ranges in central density, concentration, binary fraction, and mass. Their properties can be directly compared with those of observed GGCs. We can easily identify the BSSs in our models and determine their formation channels and birth times. We find that for central densities above ∼10{sup 3} M{sub ☉} pc{sup –3}, the dominant formation channel is stellar collisions, while for lower density clusters, mass transfer in binaries provides a significant contribution (up to 60% in our models). The majority of these collisions are binary-mediated, occurring during three-body and four-body interactions. As a result, a strong correlation between the specific frequency of BSSs and the binary fraction in a cluster can be seen in our models. We find that the number of BSSs in the core shows only a weak correlation with the collision rate estimator Γ traditionally used by observers, in agreement with the latest Hubble Space Telescope Advanced Camera for Surveys data. Using an idealized 'full mixing' prescription for collision products, our models indicate that the BSSs observed today may have formed several Gyr ago. However, denser clusters tend to have younger (∼1 Gyr) BSSs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A New On-the-Fly Sampling Method for Incoherent Inelastic Thermal Neutron Scattering Data in MCNP6

At thermal energies, the scattering of neutrons in a system is complicated by the comparable velocities of the neutron and target, resulting in competing upscattering and downscattering events. The neutron wavelength is also similar in size to the target's interatomic spacing making the scattering process a quantum mechanical problem. Because of the complicated nature of scattering at low energies, the thermal data files in ACE format used in continuous-energy Monte Carlo codes are quite large { on the order of megabytes for a single temperature and material. In this paper, a new storage and sampling method is introduced that is orders of magnitude less in size and is used to sample scattering parameters at any temperature on-the-fly. In addition to the reduction in storage, the need to pre-generate thermal scattering data tables at fine temperatures has been eliminated. This is advantageous for multiphysics simulations which may involve temperatures not known in advance. A new module was written for MCNP6 that bypasses the current S(α,β) table lookup in favor of the new format. The new on-the-fly sampling method was tested for graphite for two benchmark problems at ten temperatures: 1) an eigenvalue test with a fuel compact of uranium oxycarbide fuel homogenized into a graphite matrix, 2) a surface current test with a \broomstick" problem with a monoenergetic point source. The largest eigenvalue difference was 152pcm for T= 1200K. For the temperatures and incident energies chosen for the broomstick problem, the secondary neutron spectrum showed good agreement with the traditional S(α,β) sampling method. These preliminary results show that sampling thermal scattering data on-the-fly is a viable option to eliminate both the storage burden of keeping thermal data at discrete temperatures and the need to know temperatures before simulation runtime.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

SU-F-BRD-14: Dose Weighted Linear Energy Transfer Analysis of Critical Structures in Proton Therapy of Pediatric Brain Tumor Patients

Purpose: Dose weighted Linear Energy Transfer (LETd) analysis of critical structures may be useful in understanding the side effects of the proton therapy. The objective is to analyze the differences between LETd and dose distributions in brain tumor patients receiving double scattering proton therapy, to quantify LETd variation in critical organs, and to identify beam arrangements contributing to high LETd in critical organs. Methods: Monte Carlo simulations of 9 pediatric brain tumor patients were performed. The treatment plans were reconstructed with the TOPAS Monte Carlo code to calculate LETd and dose. The beam data were reconstructed proximal to the aperture of the double scattering nozzle. The dose and LETd to target and critical organs including brain stem, optic chiasm, lens, optic nerve, pituitary gland, and hypothalamus were computed for each beam. Results: Greater variability in LETd compared to dose was observed in the brainstem for patients with a variety of tumor types including 5 patients with tumors located in the posterior fossa. Approximately 20%–44% brainstem volume received LETd of 5kev/µm or greater from beams within gantry angles 180°±30° for 5 patients treated with a 3 beam arrangement. Critical organs received higher LETd when located in the vicinity of the beam distal edge. Conclusion: This study presents a novel strategy in the evaluation of the proton treatment impact on critical organs. While the dose to critical organs is confined below the required limits, the LETd may have significant variation. Critical organs in the vicinity of beam distal edge receive higher LETd and depended on beam arrangement, e.g. in posterior fossa tumor treatment, brainstem receive higher LETd from posterior-anterior beams. This study shows importance of the LETd analysis of the radiation impact on the critical organs in proton therapy and may be used to explain clinical imaging observations after therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-399: Evaluation of Selection Criteria for Computational Human Phantoms for Use in Out-Of-Field Organ Dosimetry for Radiotherapy Patients

Purpose: To quantify the dosimetric uncertainty due to organ position errors when using height and weight as phantom selection criteria in the UF/NCI Hybrid Phantom Library for the purpose of out-of-field organ dose reconstruction. Methods: Four diagnostic patient CT images were used to create 7-field IMRT plans. For each patient, dose to the liver, right lung, and left lung were calculated using the XVMC Monte Carlo code. These doses were taken to be the ground truth. For each patient, the phantom with the most closely matching height and weight was selected from the body size dependent phantom library. The patient plans were then transferred to the computational phantoms and organ doses were recalculated. Each plan was also run on 4 additional phantoms with reference heights and or weights. Maximum and mean doses for the three organs were computed, and the DVHs were extracted and compared. One sample t-tests were performed to compare the accuracy of the height and weight matched phantoms against the additional phantoms in regards to both maximum and mean dose. Results: For one of the patients, the height and weight matched phantom yielded the most accurate results across all three organs for both maximum and mean doses. For two additional patients, the matched phantom yielded the best match for one organ only. In 13 of the 24 cases, the matched phantom yielded better results than the average of the other four phantoms, though the results were only statistically significant at the .05 level for three cases. Conclusion: Using height and weight matched phantoms does yield better results in regards to out-of-field dosimetry than using average phantoms. Height and weight appear to be moderately good selection criteria, though this selection criteria failed to yield any better results for one patient.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-240: Monte Carlo Modelling of SMC Proton Nozzles Using TOPAS

Purpose: To simulate the 60Co beam from a novel Gamma—Tomo SBRT system and to calculate the output factors and dose rates for different source configurations and collimator sizes. Methods: The BEAM-Monte Carlo code is used to realistically model the system geometry, including 32 source capsules, source housing, primary collimator and 4 different changeable collimators. The sources (6.22 mm of diameter) are located in 2 rows with different angles and distances to the longitudinal axis and the beams can be collimated to obtain four different circular field sizes at the isocenter (35 mm, 16 mm, 7mm, and 3.5mm). Since the dose distributions will be determined in a spherical polystyrene phantom of 8cm radius centered at the isocenter, a new geometry module is designed for the dose calculation. Results: : The caracteristics of the particle spectra emerging from each source are determined and the influence of the source position on the spectra reaching the isocenter is found to be minimal (<2%). We found that the scattered photons represent about 15% of the total number of photons reaching the scoring plane immediately (1 mm) below the source and primary collimator. The radial photon fluence does not vary significantly inside the collimator opening and decreases quickly beyond the 3cm radius (projected at the isocenter plane). Conclusion: Our preliminary results indicate that the position of the source and primary collimator assembly does not affect the energy spectra and fluence distributions reaching the isocenter significantly. This allows the use of a well-represented source model for all the sources to perform the dose calculations for all the changeable collimators for this novel Gamma-Tomo SBRT system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-T-448: On the Perturbation Factor P-cav of the Markus Parallel Plate Ion Chambers in Clinical Electron Beams, Monte Carlo Based Reintegration of An Historical Experiment

Purpose: All present dosimetry protocols recommend well-guarded parallel-plate ion chambers for electron dosimetry. For the guard-less Markus chamber an energy dependent fluence perturbation correction pcav is given. This perturbation correction was experimentally determined by van der Plaetsen by comparison of the read-out of a Markus and a NACP chamber, which was assumed to be “perturbation-free”. Aim of the present study is a Monte Carlo based reiteration of this experiment. Methods: Detailed models of four parallel-plate chambers (Roos, Markus, NACP and Advanced Markus) were designed using the Monte Carlo code EGSnrc and placed in a water phantom. For all chambers the dose to the active volume filled with low density water was calculated for 13 clinical electron spectra (E{sub 0}=6-21 MeV) at the depth of maximum and at the reference depth under reference conditions. In all cases the chamber's reference point was positioned at the depth of measurement. Moreover, the dose to water DW was calculated in a small water voxel positioned at the same depth. Results: The calculated dose ratio D{sub NACP}/D{sub Markus}, which according to van der Plaetsen reflects the fluence perturbation correction of the Markus chamber, deviates less from unity than the values given by van der Plaetsen's but exhibits a similar energy dependence. The same holds for the dose ratios of the other well guarded chambers. But, in comparison to water, the Markus chamber reveals the smallest overall perturbation correction which is nearly energy independent at both investigated depths. Conclusion: The simulations principally confirm the energy dependence of the dose ratio D{sub NACP}/D{sub Markus} as published by van der Plaetsen. But, as shown by our simulations of the ratio D{sub W}/D{sub Markus}, the conclusion drawn in all dosimetry protocols is questionable: in contrast to all well-guarded chambers the guard-less Markus chamber reveals the smallest overall perturbation correction and also the smallest energy dependence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-F-J-149: Beam and Cryostat Scatter Characteristics of the Elekta MR-Linac

Purpose: The Elekta MR-Linac combines a digital linear accelerator system with a 1.5T Philips MRI machine. This study aimed to determine key characteristic information regarding the MR-Linac beam and in particular it evaluated the effect of the MR cryostat on the out of field scatter dose. Methods: Tissue phantom ratios, profiles and depth doses were acquired in plastic water with an IC-profiler or with an MR compatible water tank using multiple system configurations (Full (B0= 1.5T), Full (B0=0T) and No cryostat). Additionally, an in-house CAD based Monte Carlo code based on Penelope was used to provide comparative data. Results: With the cryostat in place and B0=0T, the measured TPR for the MR Linac system was 0.702, indicating an energy of around 7MV. Without the cryostat, the measured TPR was 0.669. For the Full (B0=0T) case, out of field dose at a depth of 10 cm in the isocentric plane, 5 cm from the field edge was 0.8%, 3.1% and 5.4% for 3×3 cm{sup 2}, 10×10 cm{sup 2} and 20×20 cm{sup 2} fields respectively.The out of field dose (averaged between 5 cm and 10 cm beyond the field edges) in the “with cryostat” case is 0.78% (absolute difference) higher than without the cryostat for clinically relevant field sizes (i.e. 10×10 cm{sup 2}) and comparable to measured conventional 6MV treatment beams at a depth of 10 cm (within 0.1% between 5 cm and 6 cm from field edge). At dose maximum and at 5 cm from the field edge, the “with cryostat” out of field scatter for a 10×10 cm{sup 2} field is 1.5% higher than “without cryostat', with a modest increase (0.9%) compared to Agility 6MV in the same conditions. Conclusion: The study has presented typical characteristics of the MR-Linac beam and determined that out of field dose is comparable to conventional treatment beams. All authors are employed by Elekta Ltd., who are developing an MR-Linac.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Microdosimetric study on influence of low energy photons on relative biological effectiveness under therapeutic conditions using 6 MV linac

Purpose: Microdosimetry has been developed for the evaluation of radiation quality, and single-event dose-mean lineal energy y{sub D} is well-used to represent the radiation quality. In this study, the changes of the relative biological effectiveness (RBE) values under the therapeutic conditions using a 6 MV linac were investigated with a microdosimetric method. Methods: The y{sub D} values under the various irradiation conditions for x-rays from a 6 MV linac were measured with a tissue-equivalent proportional counter (TEPC) at an extremely low dose rate of a few tens of {mu}Gy/min by decreasing the gun grid voltage of the linac. According to the microdosimetric kinetic model (MK model), the RBE{sub MK} values for cell killing of the human salivary gland (HSG) tumor cells can be derived if the y{sub D} values are obtained from TEPC measurements. The Monte Carlo code GEANT4 was also used to calculate the photon energy distributions and to investigate the changes of the y{sub D} values under the various conditions. Results: The changes of the y{sub D} values were less than approximately 10% when the field size and the depth in a phantom varied. However, in the measurements perpendicular to a central beam axis, large changes were observed between the y{sub D} values inside the field and those outside the field. The maximum increase of approximately 50% in the y{sub D} value outside the field was obtained compared with those inside the field. The GEANT4 calculations showed that there existed a large relative number of low energy photons outside of the field as compared with inside of the field. The percentages of the photon fluences below 200 keV outside the field were approximately 40% against approximately 8% inside the field. By using the MK model, the field size and the depth dependence of the RBE{sub MK} values were less than approximately 2% inside the field. However, the RBE{sub MK} values outside the field were 6.6% higher than those inside the field. Conclusions: The increase of the RBE{sub MK} values by 6.6% outside the field was observed. This increase is caused by the change of the photon energy distributions, especially the increase of the relative number of low energy photons outside the field.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

SU-E-J-51: Dose Response of Common Solid State Detectors in Homogeneous Transverse and Longitudinal Magnetic Fields

Purpose: Solid state radiation detectors are often used for dose profiles and percent depth dose measurements. The dose response of selected solid state detectors is evaluated in varying transverse and longitudinal magnetic fields for eventual use in MR-Linac devices. Methods: A PTW 60003 and IBA PFD detector were modeled in the Monte Carlo code PENELOPE, incorporating a magnetic field which was varied in strength and oriented both transversely and longitudinally with respect to the incident photon beam. The detectors' long axis was in turn oriented either parallel or perpendicular to the photon beam. Dose to the active volume of each detector was scored, and its ratio to dose with zero magnetic field strength (dose response) was determined. Accuracy of the simulations was evaluated by measurements using both chambers taken at low field with a small electromagnet. Simulations were also performed in a water phantom to compare to the in air results. Results: Significant dose response was found in transverse field geometries, nearing 20% at 1.5T. The response is highly dependent on relative orientations to the magnetic field and photon beam, and on detector composition. Low field measurements confirm these results. In the presence of longitudinal magnetic fields, the detectors exhibit little dose response, reaching 0.5–1% at 1.5T regardless of detector orientation. Water tank simulations compared well to the in air simulations when not at the beam periphery, where in transverse magnetic fields only, the water tank simulations differed from the in air results. Conclusion: Transverse magnetic fields can cause large deviations in dose response, and are highly position orientation dependent. Comparatively, longitudinal magnetic fields exhibit little to no dose response in each detector as a function of magnetic field strength. Water tank simulations show longitudinal fields are generally easier to work with, but each detector must be evaluated separately.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2013

A neutronic feasibility study of the AP1000 design loaded with fully ceramic micro-encapsulated fuel

A neutronic feasibility study is performed to evaluate the utilization of fully ceramic microencapsulated (FCM) fuel in the AP1000 reactor design. The widely used Monte Carlo code MCNP is employed to perform the full core analysis at the beginning of cycle (BOC). Both the original AP1000 design and the modified design with the replacement of uranium dioxide fuel pellets with FCM fuel compacts are modeled and simulated for comparison. To retain the original excess reactivity, ranges of fuel particle packing fraction and fuel enrichment in the FCM fuel design are first determined. Within the determined ranges, the reactor control mechanism employed by the original design is directly used in the modified design and the utilization feasibility is evaluated. The worth of control of each type of fuel burnable absorber (discrete/integral fuel burnable absorbers and soluble boron in primary coolant) is calculated for each design and significant differences between the two designs are observed. Those differences are interpreted by the fundamental difference of the fuel form used in each design. Due to the usage of silicon carbide as the matrix material and the fuel particles fuel form in FCM fuel design, neutron slowing down capability is increased in the new design, leading to a much higher thermal spectrum than the original design. This results in different reactivity and fission power density distributions in each design. We conclude that a direct replacement of fuel pellets by the FCM fuel in the AP1000 cannot retain the original optimum reactor core performance. Necessary modifications of the core design should be done and the original control mechanism needs to be re-designed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Validation for {sup 233}U-Fueled Systems in KENO V.a in SCALE 6.2

Validation of the KENO V.a Monte Carlo code has been performed in recent years using benchmark experiments contained in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The models used in these recent validation efforts have been drawn from the Verified, Archived Library of Inputs and Data (VALID). Models must be generated and reviewed by qualified individuals to be added to the VALID library, which leads to high confidence in the models and reduces the probability that apparent validation issues are a result of modeling errors. The contents of the VALID library are somewhat limited because of the effort involved in generating and reviewing models, so no {sup 233}U-fueled experiments are included in VALID. A large number of models of {sup 233}U experiments have been modeled in the past in support of the criticality safety of {sup 233}U in storage at Oak Ridge National Laboratory (ORNL). In June, 2016, two midshipmen from the United States Naval Academy were assigned to ORNL, and they reviewed and updated these models for execution in SCALE 6.2. The validation of {sup 233}U is of interest because some discrepancies in SCALE results were identified by a beta tester at the Nuclear Energy Agency modeling these systems. These discrepancies were related to a data thinning procedure in the continuous-energy nuclear data process, and had been resolved prior to the release of SCALE 6.2. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-E-T-350: Effective Point of Measurement and Total Perturbation Correction P for Parallel-Plate Ion Chambers in High-Energy Photon Beams

Purpose: This paper aims to determine the effective point of measurement and the total perturbation correction p of parallel-plate chambers for clinical photon dosimetry. Methods: The effective point of measurement (EPOM) was calculated using the EGSnrc Monte Carlo code system with the EGSnrc user code egs- chamber. Depth dose curves of the ionization chambers were calculated in a water phantom for several high energy photon spectra (4, 6, 10, 15, 18 MV-X). Different normalization criterions (normalization to the maximum of the depth dose curve and normalization to the value in 10 cm depth) have been applied. The EPOM was determined by shifting the normalized depth dose curve of a small water voxel against the depth ionization curve until the disagreement (calculated by the root mean square deviation) reaches a minimum. In addition, the total perturbation correction p was calculated by the ratio of the dose to water and the product of the dose determined in the chamber and the water to air stopping power ratio. Results: The EPOM varied slightly depending on the chosen normalization criterion. For all chambers the necessary shift of the EPOM decreased linearly with increasing beam quality specifier TPR{sub 20/10}. For the Roos and NACP chamber, the results were positive suggesting that the chambers need to be shifted towards the focus. For the Markus chamber, the required shift was negative and for the Advanced Markus chamber partly negative and partly positive. The total perturbation correction p was almost independent of the depth. Only for regions below 1 cm the perturbation correction deviated significantly from unity. Conclusion: In the present study, the effective point of measurement and the total perturbation correction p was determined for four parallel-plate ionization chambers and five clinical relevant photon spectra. Applying the calculated EPOM, the residual perturbation correction p was mostly depth independent.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Modeling the molecular gas in NGC 6240

We present the first observations of H{sup 13}CN (1−0), H{sup 13}CO{sup +}(1−0), and SiO (2−1)in NGC 6240, obtained with the IRAM Plateau de Bure Interferometer. Combining a Markov Chain Monte Carlo code with Large Velocity Gradient (LVG) modeling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass−luminosity conversion factors. We find ∼10{sup 10}M{sub ⊙} of dense molecular gas in cold, dense clouds (T{sub k}∼10 K, n{sub H{sub 2}}∼10{sup 6} cm{sup −3}) with a volume filling factor <0.002, embedded in a shock heated molecular medium (T{sub k}∼2000 K, n{sub H{sub 2}}∼10{sup 3.6} cm{sup −3}), both surrounded by an extended diffuse phase (T{sub k}∼200 K, n{sub H{sub 2}}∼10{sup 2.5} cm{sup −3}). We derive a global α{sub CO}=1.5{sub 1.1}{sup 7.1} with gas masses log{sub 10}(M/[M{sub ⊙}])=10.1{sub 10.0}{sup 10.8}, dominated by the dense gas. We also find α{sub HCN}=32{sub 13}{sup 89}, which traces the cold, dense gas. The [{sup 12}C]/[{sup 13}C] ratio is only slightly elevated (98{sub 65}{sup 230}), contrary to the very high [CO]/[{sup 13}CO] ratio (300–500) reported in the literature. However, we find very high [HCN]/[H{sup 13}CN] and [HCO{sup +}]/[H{sup 13}CO{sup +}] abundance ratios (300{sub 200}{sup 500}) which we attribute to isotope fractionation in the cold, dense clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-J-56: The Connection Between Cherenkov Light Emission and Radiation Absorbed Dose in Proton Irradiated Phantoms

Purpose: Range verification in proton therapy is of great importance. Cherenkov light follows the photon and electron energy deposition in water phantom. The purpose of this study is to investigate the connection between Cherenkov light generation and radiation absorbed dose in a water phantom irradiated with proton beams. Methods: Monte Carlo simulation was performed by employing FLUKA Monte Carlo code to stochastically simulate radiation transport, ionizing radiation dose deposition, and Cherenkov radiation in water phantoms. The simulations were performed for proton beams with energies in the range 50–600 MeV to cover a wide range of proton energies. Results: The mechanism of Cherenkov light production depends on the initial energy of protons. For proton energy with 50–400 MeV energy that is below the threshold (∼483 MeV in water) for Cherenkov light production directly from incident protons, Cherenkov light is produced mainly from the secondary electrons liberated as a result of columbic interactions with the incident protons. For proton beams with energy above 500 MeV, in the initial depth that incident protons have higher energy than the Cherenkov light production threshold, the light has higher intensity. As the slowing down process results in lower energy protons in larger depths in the water phantom, there is a knee point in the Cherenkov light curve vs. depth due to switching the Cherenkov light production mechanism from primary protons to secondary electrons. At the end of the depth dose curve the Cherenkov light intensity does not follow the dose peak because of the lack of high energy protons to produce Cherenkov light either directly or through secondary electrons. Conclusion: In contrast to photon and electron beams, Cherenkov light generation induced by proton beams does not follow the proton energy deposition specially close to the end of the proton range near the Bragg peak.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-462: Impact of the Radiochromic Film Energy Response On Dose Measurements of Low Energy Electronic Brachytherapy Sources

Purpose: We investigated the effect of the EBT3 GafChromicTM film model absorbed dose energy response when used for percent depth dose (PDD) measurements in low-energy photon beams. Methods: We measured PDDs in water from a Xoft 50 kVp source using EBT3 film, and compared them to PDD measurements acquired with a PTW-TN34013 parallel-plate ionization chamber. For the x-ray source, we simulated spectra using the EGSnrc (BEAMnrc) Monte Carlo code, and calculated Half Value Layer (HVL) at different distances from the source in water. Pieces of EBT3 film were irradiated in air and calibration curves were created in terms of air-kerma in air ((Kair)air) for different beam qualities. Pieces of EBT3 film were positioned at distances of 2–6 cm from the Xoft source in a water phantom using a custom-made holder, and irradiated at the same time. As scatter is incorporated in the measured film signal in water, measured (Kair)wat was subsequently converted into absorbed dose to water by the ratio of mass energy absorption coefficients following the AAPM TG-61 dosimetry protocol. Results: Our results show that film calibration curves obtained at beam qualities near the effective energy of the Xoft 50 kVp source in water lead to variation in absorbed dose energy dependence of the response of around 3%. However, if the calibration curve was established at MV beam quality, the error in absorbed dose could be as large as 15%. We observed agreement within 1% between PDD measurements using EBT3 film model (using a calibration curve obtained at 80 kVp, HVL=2.18 mm Al, Eeff=29.5 keV) and the parallel-plate ionization chamber. Conclusion: Accurate dose measurements using radiochromic films at low photon energies require that the radiochromic film dosimetry system be calibrated at corresponding low energies, as large absorbed dose errors are expected for calibrations performed at MV beam qualities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-F-J-147: Magnetic Field Dose Response Considerations for a Linac Monitor Chamber

Purpose: The impact of magnetic fields on the readings of a linac monitor chamber have not yet been investigated. Herein we examine the total dose response as well as any deviations in the beam parameters of flatness and symmetry when a Varian monitor chamber is irradiated within an applied magnetic field. This work has direct application to the development of Linac-MR systems worldwide. Methods: A Varian monitor chamber was modeled in the Monte Carlo code PENELOPE and irradiated in the presence of a magnetic field with a phase space generated from a model of a Linac-MR prototype system. The magnetic field strength was stepped from 0 to 3.0T in both parallel and perpendicular directions with respect to the normal surface of the phase space. Dose to each of the four regions in the monitor chamber were scored separately for every magnetic field adaptation to evaluate the effect of the magnetic field on flatness and symmetry. Results: When the magnetic field is perpendicular to the phase space normal we see a change in dose response with a maximal deviation (10–25% depending on the chamber region) near 0.75T. In the direction of electron deflection we expectedly see opposite responses in chamber regions leading to a measured asymmetry. With a magnetic field parallel to the phase space normal we see no measured asymmetries, however there is a monotonic rise in dose response leveling off at about +12% near 2.5T. Conclusion: Attention must be given to correct for the strength and direction of the magnetic field at the location of the linac monitor chamber in hybrid Linac-MR devices. Elsewise the dose sampled by these chambers may not represent the actual dose expected at isocentre; additionally there may be a need to correct for the symmetry of the beam recorded by the monitor chamber. Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗