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

SU-D-213-01: Transparent Photon Detector For The Online Monitoring Of IMRT Beams

Purpose: An innovative Transparent Detector for Radiotherapy (TraDeRa) has been developed. The detector aims at real-time monitoring of modulated beam ahead of the patient during delivery sessions, with a field cover up to 40×40 cm {sup 2}. Methods: TraDeRa consists in a pixelated matrix of ionization chambers with a patented electrodes design. An in-house designed specific integrated circuit allows to extract the signal and provides a real-time map of beam intensity and shape, at the linac pulse-scale. The measurements under irradiation are made with a 6 MV clinical X-Ray beam. Dose calculations are performed with the Monte Carlo code PENELOPE, modeling the full accelerator head and the TraDeRa detector. Results: A 2 % attenuation of the beam was measured in the presence of TraDeRa and the PENELOPE dosimetric study showed no significant modification of the photon beam properties. TraDeRa detects error leaf position as small as 1 mm compared to a reference field, for both static and modulated fields. In addition, measurements are accurate over a large dynamic range from low intensity signals, as inter-leaves leaks, to very high intensities as obtained on the medical line of the European Synchrotron Radiation Facility. The detector is fully operational for conventional and high dose rate beams as FFF modes (up to 2400 MU/min). Conclusion: The current version of TraDeRa shows promising results for IMRT quality assurance (QA), allowing pulse-scale monitoring of the beam and high sensitivity for errors detection. The attenuation is small enough not to hinder the irradiation while keeping the beam upstream of the patient under constant control. A final prototype under development will include 1600 independent electrodes, half of them with a high resolution centered on the beam axis. This compact detector provides an independent set of measurements for a better QA. Funding support : This work was supported by the LABEX PRIMES (ANR-11-LABX-0063) of Universite de Lyon, within the program “Investissements d’Avenir” (ANR-11-IDEX-0007) operated by the French National Research Agency (ANR) and within the project “INSPIRA” operated by the OSEO institution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-129: Impact of Radial Fluctuations in RBE for Therapeutic Proton Beams

Purpose: To evaluate the off axis relative biological effectiveness (RBE) for actively scanned proton beams and determine if a constant radial RBE can be assumed. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary protons were simulated for a 0.6cm diameter pencil beam. Beam energies corresponding to Bragg Peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially for annuli of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm outer radius. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model, for human submandibular gland (HSG) cells, to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. Results: RBE was generally seen to increase as distance from the central axis (CAX) increased. However, this increase was only seen in low dose regions and its overall effects on the transverse biological dose remains low. In the entrance region of the phantom (10mm depth), minimum and maximum calculated RBEs varied between 15.22 and 18.88% for different energies. At the Bragg peak, this difference ranged from 3.15 to 26.77%. Despite these rather large variations the dose-weighted RBE and the CAX RBE varied by less than 0.14% at 10mm depth and less than 0.16% at the Bragg peak. Similarly small variations were found at all depths proximal of the Bragg peak. Conclusion: Although proton RBE does vary radially, its overall effect on biological dose is minimal and the use of a radially constant RBE in treatment planning for scanned proton beams would not produce large errors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-D-209-07: Monte Carlo Assessment of Dose to the Lens of the Eye of Radiologist Using Realistic Phantoms and Eyeglass Models

Purpose: To study how eyeglass design features and postures of the interventional radiologist affect the radiation dose to the lens of the eye. Methods: A mesh-based deformable phantom, consisting of an ultra-fine eye model, was used to simulate postures of a radiologist in fluoroscopically guided interventional procedure (facing the patient, 45 degree to the left, and 45 degree to the right). Various eyewear design features were studied, including the shape, lead-equivalent thickness, and separation from the face. The MCNPX Monte Carlo code was used to simulate the X-ray source used for the transcatheter arterial chemoembolization procedure (The X-ray tube is located 35 cm from the ground, emitting X-rays toward to the ceiling; Field size is 40cm X 40cm; X-ray tube voltage is 90 kVp). Experiments were also performed using dosimeter placed on a physical phantom behind eyeglasses. Results: Without protective eyewear, the radiologist’s eye lens can receive an annual dose equivalent of about 80 mSv. When wearing a pair of lead eyeglasses with lead-equivalent of 0.5-mm Pb, the annual dose equivalent of the eye lens is reduced to 31.47 mSv, but both exceed the new ICRP limit of 20 mSv. A face shield with a lead-equivalent of 0.125-mm Pb in the shape of a semi-cylinder (13cm in radius and 20-cm in height) would further reduce the exposure to the lens of the eye. Examination of postures and eyeglass features reveal surprising information, including that the glass-to-eye separation also plays an important role in the dose to the eye lens from scattered X-ray from underneath and the side. Results are in general agreement with measurements. Conclusion: There is an urgent need to further understand the relationship between the radiation environment and the radiologist’s eyewear and posture in order to provide necessary protection to the interventional radiologists under newly reduced dose limits.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Radiological Characterisation of VLLW Produced at High-Energy Particle Accelerators - 19118

In this paper we describe the characterisation method of very-low-level radioactive waste (VLLW) produced at the European Organization for Nuclear Research (CERN). The method is based on the calculation of a preliminary radionuclide inventory using either analytical or Monte Carlo codes, on the measurement of γ-emitters and on the quantification of pure-α, pure-β and low-energy X-ray emitters, called difficult-to-measure (DTM) radionuclides, using the so-called scaling factor (SF) technique. The first stage of the characterisation process is the calculation of the radionuclide inventory, which is a list of radionuclides that can potentially be produced via the interaction of primary and secondary particles with the surrounding media. Once the preliminary radionuclide inventory is obtained, the γ-emitting radionuclides are measured via γ-ray spectrometry on each package of the waste population. The major γemitter, called key nuclide (KN), is also identified. The scaling factor method estimates the activity of DTM radionuclides by checking for a consistent and repeated relationship between the key nuclide and the activity of the difficult to measure radionuclides from samples. Between 2016 and 2018, several hundred samples were collected at CERN and over 1000 radiochemical analyses were performed to estimate robust statistics to quantify scaling factors for the most common DTM radionuclides, such as H-3, Fe-55 and Ni-63. This paper illustrates the sample collection strategy, the tests performed to choose the most appropriate statistical models describing the behaviour of scaling factors and a complete quantification of the uncertainties associated with the characterisation process. Finally, we present the results obtained during the characterisation of over 2000 m{sup 3} radioactive metals and cables disposed of at the French National repository for VLLW. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Using Rose’s metal alloy as a pinhole collimator material in preclinical small-animal imaging: A Monte Carlo evaluation

Purpose: Pinhole collimation is the most common method of high-resolution preclinical single photon emission computed tomography imaging. The collimators are usually constructed from dense materials with high atomic numbers, such as gold and platinum, which are expensive and not always flexible in the fabrication step. In this work, the authors have investigated the properties of a fusible alloy called Rose’s metal and its potential in pinhole preclinical imaging. When compared to current standard pinhole materials such as gold and platinum, Rose’s metal has a lower density and a relatively low effective atomic number. However, it is inexpensive, has a low melting point, and does not contract when solidifying. Once cast, the piece can be machined with high precision. The aim of this study was to evaluate the imaging properties for Rose’s metal and compare them with those of standard materials. Methods: After validating their Monte Carlo code by comparing its results with published data and the results from analytical calculations, they investigated different pinhole geometries by varying the collimator material, acceptance angle, aperture diameter, and photon incident angle. The penetration-to-scatter and penetration-to-total component ratios, sensitivity, and the spatial resolution were determined for gold, tungsten, and Rose’s metal for two radionuclides, {sup 99}Tc{sup m} and {sup 125}I. Results: The Rose’s metal pinhole-imaging simulations show higher penetration/total and scatter/total ratios. For example, the penetration/total is 50% for gold and 75% for Rose’s metal when simulating {sup 99}Tc{sup m} with a 0.3 mm aperture diameter and a 60° acceptance angle. However, the degradation in spatial resolution remained below 10% relative to the spatial resolution for gold for acceptance angles below 40° and aperture diameters larger than 0.5 mm. Conclusions: Extra penetration and scatter associated with Rose’s metal contribute to degradation in the spatial resolution, but this degradation is not always substantial. The most important factor besides the collimator material was the acceptance angle. This should be kept to a minimum to prevent unnecessary scatter and penetration. For {sup 125}I, the difference in spatial resolution between gold and Rose’s metal is very small, 2.2% in the worst-case scenario. Based on these results, the authors conclude that Rose’s metal is an alternative to standard materials not only for low-energy photon imaging but also for medium-energy applications that require low-cost, flexible pinhole configurations and designs, and that can tolerate a degraded spatial resolution.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2012

FORMATION OF MASSIVE BLACK HOLES IN DENSE STAR CLUSTERS. II. INITIAL MASS FUNCTION AND PRIMORDIAL MASS SEGREGATION

A promising mechanism to form intermediate-mass black holes is the runaway merger in dense star clusters, where main-sequence stars collide and form a very massive star (VMS), which then collapses to a black hole (BH). In this paper, we study the effects of primordial mass segregation and the importance of the stellar initial mass function (IMF) on the runaway growth of VMSs using a dynamical Monte Carlo code for N-body systems with N as high as 10{sup 6} stars. Our code now includes an explicit treatment of all stellar collisions. We place special emphasis on the possibility of top-heavy IMFs, as observed in some very young massive clusters. We find that both primordial mass segregation and the shape of the IMF affect the rate of core collapse of star clusters and thus the time of the runaway. When we include primordial mass segregation, we generally see a decrease in core-collapse time (t{sub cc}). Although for smaller degrees of primordial mass segregation this decrease in t{sub cc} is mostly due to the change in the density profile of the cluster, for highly mass-segregated (primordial) clusters, it is the increase in the average mass in the core which reduces the central relaxation time decreasing t{sub cc}. The final mass of the VMS formed is always close to {approx}10{sup -3} of the total cluster mass, in agreement with previous studies and is reminiscent of the observed correlation between the central BH mass and the bulge mass of the galaxies. As the degree of primordial mass segregation is increased, the mass of the VMS increases at most by a factor of three. Flatter IMFs generally increase the average mass in the whole cluster, which increases t{sub cc}. For the range of IMFs investigated in this paper, this increase in t{sub cc} is to some degree balanced by stellar collisions, which accelerate core collapse. Thus, there is no significant change in t{sub cc} for the somewhat flatter global IMFs observed in very young massive clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-131: No Increase in Biological Effectiveness Through Collimator Scattered Low Energy Protons

Purpose: To reduce the lateral penumbra of low-energy proton beams, brass collimators are often used in spot-scanning proton therapy (SSPT). This study investigates the increase in biological effectiveness through collimator scattered protons in SSPT. Methods: The SSPT system of the Hokkaido University Hospital Proton Beam Therapy Center, which consists of a scanning nozzle, a 2-cm thick brass collimator, and a 4-cm thick energy absorber, was simulated with our validated Geant4 Monte Carlo code (ver. 9.3). A water phantom was irradiated with proton pencil beams of 76, 110, and 143 MeV. The tested collimator opening areas (COA) were 5×5, 10×10, and 15×15 cm{sup 2}. Comparisons were made among the dose-averaged LET values of protons that hit the collimators (LETDColl), protons that did not hit the collimators (LETDNoColl), and all protons (LETDTotal). X-ray equivalent doses (Deq) were calculated using the linear-quadratic model with LETDNoColl and LETDTotal, and their maximum difference was determined over regions where the physical dose was greater than 10% of the peak dose of 2 Gy. Results: The ratio of the dose contribution of collimator scattered protons to that of all protons, defined as λ, was large at high proton energies and large COAs. The maximum λ value ranged from 3% (76 MeV, 5×5 cm{sup 2}) to 29% (143 MeV, 15×15 cm{sup 2}). Moreover, a large difference between LETDColl and LETDNoColl was only found in regions where λ was below 20% (ΔLETD > 2 keV/µm) and 8% (ΔLETD > 5 keV/µm). Consequently, the maximum difference between LETDNoColl and LETDTotal was as small as 0.8 keV/µm in all simulated voxels, and the difference of Deq reached a maximum of 1.5% that of the peak dose obtained at the water surface with a 76 MeV beam. Conclusion: Although collimator scattered protons have high LET, they only increase the physical dose, not the biological effectiveness.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2012

THE MERGER RATE OF BINARY WHITE DWARFS IN THE GALACTIC DISK

We use multi-epoch spectroscopy of {approx}4000 white dwarfs in the Sloan Digital Sky Survey to constrain the properties of the Galactic population of binary white dwarf systems and calculate their merger rate. With a Monte Carlo code, we model the distribution of {Delta}RV{sub max}, the maximum radial velocity shift between exposures of the same star, as a function of the binary fraction within 0.05 AU, f{sub bin}, and the power-law index in the separation distribution at the end of the common-envelope phase, {alpha}. Although there is some degeneracy between f{sub bin} and {alpha}, the 15 high-{Delta}RV{sub max} systems that we find constrain the combination of these parameters, which determines a white dwarf merger rate per unit stellar mass of 1.4{sup +3.4}{sub -1.0} Multiplication-Sign 10{sup -13} yr{sup -1} M{sup -1}{sub Sun} (1{sigma} limits). This is remarkably similar to the measured rate of Type Ia supernovae (SNe Ia) per unit stellar mass in Milky-Way-like Sbc galaxies. The rate of super-Chandrasekhar mergers is only 1.0{sup +1.6}{sub -0.6} Multiplication-Sign 10{sup -14} yr{sup -1} M{sup -1}{sub Sun }. We conclude that there are not enough close binary white dwarf systems to reproduce the observed SN Ia rate in the 'classic' double degenerate super-Chandrasekhar scenario. On the other hand, if sub-Chandrasekhar mergers can lead to SNe Ia, as has been recently suggested by some studies, they could make a major contribution to the overall SN Ia rate. Although unlikely, we cannot rule out contamination of our sample by M-dwarf binaries or non-Gaussian errors. These issues will be clarified in the near future by completing the follow-up of all 15 high-{Delta}RV{sub max} systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

WE-H-BRA-09: Application of a Modified Microdosimetric-Kinetic Model to Analyze Relative Biological Effectiveness of Ions Relevant to Light Ion Therapy Using the Particle Heavy Ion Transport System

Purpose: To evaluate the dose and biological effectiveness of various ions that could potentially be used for actively scanned particle therapy. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary particles for 1H beams and ten million particles for 4He, 7Li, 10B, 12C, 14N, 16O, and 20Ne were simulated for 0.6cm diameter pencil beams. Beam energies corresponding to Bragg peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially in annuli with outer radius of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model for five different cell types to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. The product of the calculated RBEs and the simulated physical dose was taken to create biological dose and comparisons were then made between the various ions. Results: Transversely, the 10B beam was seen to minimize relative biological dose in both the constant and accelerated dose change regions, proximal to the Bragg Peak, for all beams traveling greater than 50mm. For the 50mm beam, 7Li was seen to provide the most optimal biological dose profile. Radially small fluctuations (<4.2%) were seen in RBE while physical dose was greater than 1% for all beams. Conclusion: Even with the growing usage of 12C, it may not be the most optimal ion in all clinical situations. Boron was calculated to have slightly enhanced RBE characteristics, leading to lower relative biological doses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Absolute Dosimetry Model of the University of Washington Clinical Neutron Therapy System (CNTS)

The University of Washington (UW) Clinical Neutron Therapy System (CNTS) generates a fast neutron beam using 50.5 MeV protons incident on a beryllium target. The CNTS is mainly used for fast neutron therapy in the treatment of cancer but is also used for radiation effects testing and the production of rare isotopes, or isotopes with a short half-life, for medical imaging and therapy. The accelerator and treatment head is modeled using the MCNP6 Monte Carlo code to enable beam characterization (neutron and photon fluence, linear energy transfer, secondary particle production), study the behavior of neutron dosimeters and ionization chamber (IC) responses, and to conduct first principle studies of the relative biological effectiveness (RBE) of fast neutrons compared to {sup 60}Co γ-rays. Measurements of neutron dose in selected CNTS fields were made using a model-17 ion chamber (IC-17) flowing methane-based, tissue-equivalent (TE) gas. As a first approximation, we compared IC-17 measurements to the MCNP calculated dose in a thin layer of water at selected locations along the central axis of the beam. To examine the impact of the composition of the chamber fill-gas and the chamber design characteristics, we also developed a high-resolution (∼ 50 μm) MCNP model of the IC-17 from a computerized tomography (CT) scan. Initial results indicate that the dose response of the IC-17 ion chamber is up to 30% larger (smaller?) compared to the dose to water for neutrons with energies >20 MeV. Differences in chamber and water response as large as 800% occur for neutrons with kinetic energies below 1 keV. The simulations of IC-17 response (and dose to water) for CNTS fast neutron spectra match measured doses within 4.6% for all field sizes and irradiation geometries considered. This reported studies provide additional evidence establishing use of the IC-17 ion chamber model or a thin layer of water as a reasonable surrogate for CNTS fast neutron dose calculations.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Dosimetric characterization of a microDiamond detector in clinical scanned carbon ion beams

Purpose: To investigate for the first time the dosimetric properties of a new commercial synthetic diamond detector (PTW microDiamond) in high-energy scanned clinical carbon ion beams generated by a synchrotron at the CNAO facility. Methods: The detector response was evaluated in a water phantom with actively scanned carbon ion beams ranging from 115 to 380 MeV/u (30–250 mm Bragg peak depth in water). Homogeneous square fields of 3 × 3 and 6 × 6 cm{sup 2} were used. Short- and medium-term (2 months) detector response stability, dependence on beam energy as well as ion type (carbon ions and protons), linearity with dose, and directional and dose-rate dependence were investigated. The depth dose curve of a 280 MeV/u carbon ion beam, scanned over a 3 × 3 cm{sup 2} area, was measured with the microDiamond detector and compared to that measured using a PTW Advanced Markus ionization chamber, and also simulated using FLUKA Monte Carlo code. The detector response in two spread-out-Bragg-peaks (SOBPs), respectively, centered at 9 and 21 cm depths in water and calculated using the treatment planning system (TPS) used at CNAO, was measured. Results: A negligible drift of detector sensitivity within the experimental session was seen, indicating that no detector preirradiation was needed. Short-term response reproducibility around 1% (1 standard deviation) was found. Only 2% maximum variation of microDiamond sensitivity was observed among all the evaluated proton and carbon ion beam energies. The detector response showed a good linear behavior. Detector sensitivity was found to be dose-rate independent, with a variation below 1.3% in the evaluated dose-rate range. A very good agreement between measured and simulated Bragg curves with both microDiamond and Advanced Markus chamber was found, showing a negligible LET dependence of the tested detector. A depth dose curve was also measured by positioning the microDiamond with its main axis oriented orthogonally to the beam direction. A strong distortion in Bragg peak measurement was observed, confirming manufacturer recommendation on avoiding such configuration. Very good results were obtained for SOBP measurements, with a difference below 1% between measured and TPS-calculated doses. The stability of detector sensitivity in the observation period was within the experimental uncertainty. Conclusions: Dosimetric characterization of a PTW microDiamond detector in high-energy scanned carbon ion beams was performed. The results of the present study showed that this detector is suitable for dosimetry of clinical carbon ion beams, with a negligible LET and dose-rate dependence.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2012

Kilovoltage Imaging Doses in the Radiotherapy of Pediatric Cancer Patients

Purpose: To investigate doses induced by kilovoltage cone-beam computed tomography (kVCBCT) to pediatric cancer patients undergoing radiotherapy, as well as strategies for dose reduction. Methods and Materials: An EGS4 Monte Carlo code was used to calculate three-dimensional dose deposition due to kVCBCT on 4 pediatric cancer patients. Absorbed doses to various organs were analyzed for both half-fan and full-fan modes. Clinical conditions, such as distance from organ at risk (OAR) to CBCT field border, kV peak energy, and testicular shielding, were studied. Results: The mean doses induced by one CBCT scan operated at 125 kV in half-fan mode to testes, liver, kidneys, femoral heads, spinal cord, brain, eyes, lens, and optical nerves were 2.9, 4.7, 7.7, 10.5, 8.8, 7.6, 7.7, 7.8, and 7.2 cGy, respectively. Increasing the distances from OARs to CBCT field border greatly reduced the doses to OARs, ranging from 33% reduction for spinal cord to 2300% reduction for testes. As photon beam energy increased from 60 to 125 kV, the dose increase due to kVCBCT ranged from 170% for lens to 460% for brain and spinal cord. A testicular shielding made of 1-cm cerrobend could reduce CBCT doses down to 31%, 51%, 68%, and 82%, respectively, for 60, 80, 100, and 125 kV when the testes lay within the CBCT field. Conclusions: Generally speaking, kVCBCT deposits much larger doses to critical structures in children than in adults, usually by a factor of 2 to 3. Increasing the distances from OARs to CBCT field border greatly reduces doses to OARs. Depending on OARs, kVCBCT-induced doses increase linearly or exponentially with photon beam energy. Testicular shielding works more efficiently at lower kV energies. On the basis of our study, it is essential to choose an appropriate scanning protocol when kVCBCT is applied to pediatric cancer patients routinely.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

TH-CD-BRA-07: MRI-Linac Dosimetry: Parameters That Change in a Magnetic Field

Purpose: In MRI-linac integrated systems, the presence of the magnetic (B-)field has a large impact of the dose-distribution and the dose-responses of detectors; yet established protocols and previous experience may lead to assumptions about the commissioning process that are no longer valid. This study quantifies parameters that change when performing dosimetry with an MRI-linac including beam quality specifiers and the effective-point-of-measurement (EPOM) of ionization chambers. Methods: We used the Geant4 Monte Carlo code for this work with physics parameters that pass the Fano cavity test to within 0.1% for the simulated conditions with and without a 1.5 T B-field. A point source model with the energy distribution of an MRI-linac beam was used with and without the B-field to calculate the beam quality specifiers %dd(10)× and TPR{sup 20}{sub 10}, the variation of chamber response with orientation and the how the B-field affects the EPOM of ionization chambers by comparing depth-dose curves calculated in water to those generated by a model PTW30013 Farmer chamber. Results: The %dd(10)× changes by over 2% in the presence of the B-field while the TPR{sup 20}{sub 10} is unaffected. Ionization chamber dose-response is known to depend on the orientation w.r.t. the B-field, but two alternative perpendicular orientations (anti-parallel to each other) also differ in dose-response by over 1%. The B-field shifts the EPOM downstream (closer to the chamber center) but it is also shifted laterally by 0.27 times the chamber’s cavity radius. Conclusion: The EPOM is affected by the B-field and it even shifts laterally. The relationship between %dd(10)× and the Spencer-Attix stopping powers is also changed. Care must be taken when using chambers perpendicular to the field as the dose-response changes depending on which perpendicular orientation is used. All of these effects must be considered when performing dosimetry in B-fields and should be accounted for in future dosimetry protocols. This project was partially funded by Elekta Ltd.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Dose response of selected ion chambers in applied homogeneous transverse and longitudinal magnetic fields

Purpose: The magnetic fields of an integrated MR-Linac system will alter the paths of electrons that produce ions in the ionization chambers. The dose response of selected ion chambers is evaluated in the presence of varying transverse and longitudinal magnetic fields. The investigation is useful in calibration of therapeutic x-ray beams associated with MR-Linac systems. Methods: The Monte Carlo code PENELOPE was used to model the irradiation of NE2571, and PR06C ionization chambers in the presence of a transverse and longitudinal (with respect to the photon beam) magnetic fields of varying magnitude. The long axis of each chamber was simulated both parallel and perpendicular to the incident photon beam for each magnetic field case. The dose deposited in each chamber for each case was compared to the case with zero magnetic field by means of a ratio. The PR06C chamber's response was measured in the presence of a transverse magnetic field with field strengths ranging from 0.0 to 0.2 T to compare to simulated results. Results: The simulations and measured data show that in the presence of a transverse magnetic field there is a considerable dose response (maximum of 11% near 1.0 T in the ion chambers investigated, which depends on the magnitude of magnetic field, and relative orientation of the magnetic field, radiation beam, and ion chamber. Measurements made with the PR06C chamber verify these results in the region of measurement. In contrast, a longitudinal magnetic field produces only a slight increase in dose response (2% at 1.5 T) that rises slowly with increasing magnetic field and is seemingly independent of chamber orientation. Response trends were similar for the two ion chambers and relative orientations considered, but slight variations are present from chamber to chamber. Conclusions: Care must be taken when making ion chamber measurements in a transverse magnetic field. Ion chamber responses vary not only with transverse field strength, but with chamber orientation and type, and can be considerable. Longitudinal magnetic fields influence ion chamber responses relatively little (2% at 1.5 T), and only at field strengths in excess of 1.0 T.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-AB-201-09 [Medical Physics, Jun 2016, v. 43(6)]

Purpose: To determine the k-(Q-msr,Q)^(f-msr,f-ref ) factor, introduced in the small field formalism for five common type chambers used in the calibration of Leksell Gamma-Knife Perfexion model over a range of different phantom electron densities. Methods: Five chamber types including Exradin-A16, A14SL, A14, A1SL and IBA-CC04 are modeled in EGSnrc and PENELOPE Monte Carlo codes using the blueprints provided by the manufacturers. The chambers are placed in a previously proposed water-filled phantom and four 16-cm diameter spherical phantoms made of liquid water, Solid Water, ABS and polystyrene. Dose to the cavity of the chambers and a small water volume are calculated using EGSnrc/PENELOPE codes. The calculations are performed over a range of phantom electron densities for two chamber orientations. Using the calculated dose-ratio in reference and machine specific reference field, the k-(Q-msr,Q)^(f-msr,f-ref ) factor can be determined. Results: When chambers are placed along the symmetry axis of the collimator block (z-axis), the CC04 requires the smallest correction followed by A1SL and A16. However, when detectors are placed perpendicular to z-axis, A14SL needs the smallest and A16 the largest correction. Moreover, an increase in the phantom electron density results in a linear increase in the k-(Q-msr,Q)^(f-msr,f-ref ). Depending on the chambers, the agreement between this study and a previous study performed varies between 0.05–0.70% for liquid water, 0.07–0.85% for Solid Water and 0.00–0.60% for ABS phantoms. After applying the EGSnrc-calculated k-(Q-msr,Q)^(f-msr,f-ref ) factors for A16 to the previously measured dose-rates in liquid water, Solid Water and ABS normalized to the dose-rate measured with TG-21 protocol and ABS phantom, the dose-rate ratios are found to be 1.004±0.002, 0.996±0.002 and 0.998±0.002 (3σ) respectively. Conclusion: Knowing the electron density of the phantoms, the calculated k-(Q-msr,Q)^(f-msr,f-ref ) values in this work will enable users to apply the appropriate correction for their own specific phantom material. LM acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MO-F-CAMPUS-I-02: Accuracy in Converting the Average Breast Dose Into the Mean Glandular Dose (MGD) Using the F-Factor in Cone Beam Breast CT- a Monte Carlo Study Using Homogeneous and Quasi-Homogeneous Phantoms

Purpose: To investigate the accuracy in estimating the mean glandular dose (MGD) for homogeneous breast phantoms by converting from the average breast dose using the F-factor in cone beam breast CT. Methods: EGSnrc-based Monte Carlo codes were used to estimate the MGDs. 13-cm in diameter, 10-cm high hemi-ellipsoids were used to simulate pendant-geometry breasts. Two different types of hemi-ellipsoidal models were employed: voxels in quasi-homogeneous phantoms were designed as either adipose or glandular tissue while voxels in homogeneous phantoms were designed as the mixture of adipose and glandular tissues. Breast compositions of 25% and 50% volume glandular fractions (VGFs), defined as the ratio of glandular tissue voxels to entire breast voxels in the quasi-homogeneous phantoms, were studied. These VGFs were converted into glandular fractions by weight and used to construct the corresponding homogeneous phantoms. 80 kVp x-rays with a mean energy of 47 keV was used in the simulation. A total of 109 photons were used to image the phantoms and the energies deposited in the phantom voxels were tallied. Breast doses in homogeneous phantoms were averaged over all voxels and then used to calculate the MGDs using the F-factors evaluated at the mean energy of the x-rays. The MGDs for quasi-homogeneous phantoms were computed directly by averaging the doses over all glandular tissue voxels. The MGDs estimated for the two types of phantoms were normalized to the free-in-air dose at the iso-center and compared. Results: The normalized MGDs were 0.756 and 0.732 mGy/mGy for the 25% and 50% VGF homogeneous breasts and 0.761 and 0.733 mGy/mGy for the corresponding quasi-homogeneous breasts, respectively. The MGDs estimated for the two types of phantoms were similar within 1% in this study. Conclusion: MGDs for homogeneous breast models may be adequately estimated by converting from the average breast dose using the F-factor.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Proton radiography in three dimensions: A proof of principle of a new technique

Purpose: Monte Carlo simulations were used to investigate a range of phantom configurations to establish enabling three-dimensional proton radiographic techniques.Methods: A large parameter space of stacked phantom geometries composed of tissue inhomogeneity materials such as lung, bone, and cartilage inserted within water background were simulated using a purposefully modified version of TOPAS, an application running on top of the GEANT4 Monte Carlo code. The phantoms were grouped in two classes, one with the inhomogeneity inserted only half-way in the lateral direction and another with complete inhomogeneity insertion. The former class was used to calculate the track count and the energy fluence of the protons as they exit the phantoms either having traversed the inhomogeneity or not. The latter class was used to calculate one yield value accounting for loss of protons due to physical processes only and another yield value accounting for deliberately discarded protons due to large scattering angles. A graphical fingerprinting method was developed to determine the inhomogeneity thickness and location within the phantom based on track count and energy fluence information. Two additional yield values extended this method to the general case which also determines the inhomogeneity material and the phantom thickness.Results: The graphical fingerprinting method was manually validated for two, and automatically tested for all, tissue materials using an exhaustive set of inhomogeneity geometries for 16 cm thick phantoms. Unique recognition of test phantom configurations was achieved in the large majority of cases. The method in the general case was further tested using an exhaustive set of inhomogeneity and phantom tissues and geometries where the phantom thicknesses ranged between 8 and 24 cm. Unique recognition of the test phantom configurations was achieved only for part of the phantom parameter space. The correlations between the remaining false positive recognitions were analyzed.Conclusions: The concept of 3D proton radiography for tissue inhomogeneities of simple geometries was established with the current work. In contrast to conventional 2D proton radiography, the main objective of the demonstrated 3D technique is not proton range. Rather, it is to measure the depth and thickness of an inhomogeneity located in an imaged geometry. Further work is needed to extend and apply the method to more complex geometries.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

3D inpatient dose reconstruction from the PET-CT imaging of {sup 90}Y microspheres for metastatic cancer to the liver: Feasibility study

Purpose: The introduction of radioembolization with microspheres represents a significant step forward in the treatment of patients with metastatic disease to the liver. This technique uses semiempirical formulae based on body surface area or liver and target volumes to calculate the required total activity for a given patient. However, this treatment modality lacks extremely important information, which is the three-dimensional (3D) dose delivered by microspheres to different organs after their administration. The absence of this information dramatically limits the clinical efficacy of this modality, specifically the predictive power of the treatment. Therefore, the aim of this study is to develop a 3D dose calculation technique that is based on the PET imaging of the infused microspheres.Methods: The Fluka Monte Carlo code was used to calculate the voxel dose kernel for {sup 90}Y source with voxel size equal to that of the PET scan. The measured PET activity distribution was converted to total activity distribution for the subsequent convolution with the voxel dose kernel to obtain the 3D dose distribution. In addition, dose-volume histograms were generated to analyze the dose to the tumor and critical structures.Results: The 3D inpatient dose distribution can be reconstructed from the PET data of a patient scanned after the infusion of microspheres. A total of seven patients have been analyzed so far using the proposed reconstruction method. Four patients underwent treatment with SIR-Spheres for liver metastases from colorectal cancer and three patients were treated with Therasphere for hepatocellular cancer. A total of 14 target tumors were contoured on post-treatment PET-CT scans for dosimetric evaluation. Mean prescription activity was 1.7 GBq (range: 0.58–3.8 GBq). The resulting mean maximum measured dose to targets was 167 Gy (range: 71–311 Gy). Mean minimum dose to 70% of target (D70) was 68 Gy (range: 25–155 Gy). Mean minimum dose to 90% of target (D90) was 53 Gy (range: 13–125 Gy).Conclusions: A three-dimensional inpatient dose reconstruction method has been developed that is based on the PET/CT data of a patient treated with {sup 90}Y microspheres. It allows for a complete description of the absorbed dose by the tumor and critical structures. It represents the first step in building predictive models for treatment outcomes for patients receiving this therapeutic modality as well as it allows for better analysis of patients' dose response and will ultimately improve future treatment administration.

61 RADIATION PROTECTION AND DOSIMETRY↗