Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dose 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 433 records · Page 24

PNNL Richland Campus Radionuclide Air Emissions Report for Calendar Year 2021

This report documents radionuclide air emissions that result in the 2021 highest effective dose equivalent (EDE) to an offsite member of the public, referred to as the maximally exposed individual (MEI). The report has been prepared in compliance with the Code of Federal Regulations, Title 40, Protection of the Environment, Part 61, National Emission Standards for Hazardous Air Pollutants, Subpart H, “National Emission Standards for Emissions of Radionuclides Other than Radon from Department of Energy Facilities” and Washington Administrative Code Chapter 246-247, “Radiation Protection–Air Emissions.” The dose to the PNNL Richland Campus MEI from routine major and minor point source emissions in 2021 from PNNL Richland Campus sources is 1.7E-5 mrem (1.7E-7 mSv) EDE. The dose from all fugitive sources is 1.2E-6 mrem (1.2E-8 mSv) EDE. The dose from radon emissions is 2.1E-9 mrem (2.1E-11 mSv) EDE. No nonroutine emissions occurred in 2021. The total radiological dose to the MEI from all PNNL Richland Campus radionuclide emissions, including fugitive emissions and radon, is 1.8E-5 mrem (1.8E-7 mSv) EDE, or more than 100,000 times less than the federal and state standard of 10 mrem/yr, with which the PNNL Richland Campus is in compliance.

40 CFR 61 Subpart H↗

PNNL-Richland Campus Radionuclide Air Emissions Report for Calendar Year 2022

This report documents radionuclide air emissions that result in the 2022 highest effective dose equivalent (EDE) to an offsite member of the public, referred to as the maximally exposed individual (MEI). The report has been prepared in compliance with the Code of Federal Regulations, Title 40, Protection of the Environment, Part 61, National Emission Standards for Hazardous Air Pollutants, Subpart H, “National Emission Standards for Emissions of Radionuclides Other than Radon from Department of Energy Facilities” and Washington Administrative Code Chapter 246-247, “Radiation Protection–Air Emissions.” The dose to the PNNL Richland Campus MEI from routine major and minor point source emissions in 2022 from PNNL Richland Campus sources is 2.1E-5 mrem (2.1E-7 mSv) EDE. The dose from all fugitive sources is 1.2E-6 mrem (1.2E-8 mSv) EDE. The dose from radon emissions is 9.2E-9 mrem (9.2E-11 mSv) EDE. No nonroutine emissions occurred in 2022. The total radiological dose to the MEI from all PNNL Richland Campus radionuclide emissions, including fugitive emissions and radon, is 2.3E-5 mrem (2.3E-7 mSv) EDE, or more than 100,000 times less than the federal and state standard of 10 mrem/yr, with which the PNNL Richland Campus is in compliance.

40 CFR 61 Subpart↗

PNNL-Richland Campus Radionuclide Air Emissions Report for Calendar Year 2023

This report documents radionuclide air emissions that result in the 2023 highest effective dose equivalent (EDE) to an offsite member of the public, referred to as the maximally exposed individual (MEI). The report has been prepared in compliance with the Code of Federal Regulations, Title 40, Protection of the Environment, Part 61, National Emission Standards for Hazardous Air Pollutants, Subpart H, “National Emission Standards for Emissions of Radionuclides Other than Radon from Department of Energy Facilities” and Washington Administrative Code Chapter 246-247, “Radiation Protection–Air Emissions.” The dose to the PNNL-Richland campus MEI from routine emissions sources, excluding radon, in 2023 from campus sources is 2.0E-5 mrem (2.0E-7 mSv) EDE. The dose from radon emissions is 4.0E-7 mrem (4.0E-09 mSv) EDE. No nonroutine emissions occurred in 2023. The total radiological dose to the MEI from all PNNL-Richland campus radionuclide emissions, including fugitive emissions and radon, is 2.1E-5 mrem (2.1E-7 mSv) EDE, or more than 100,000 times less than the federal and state standard of 10 mrem/yr, with which the PNNL-Richland campus is in compliance.

40 CFR 61 Subpart H↗

SU-F-T-302: Implementation of the RadCalc Image Analysis Tool for IMRT QA

Purpose: To assess RadCalc as a means for analyzing fluence maps of IMRT QA patient plans. Eclipse generated fluence maps were compared to RadCalc fluence maps generated from both Eclipse parameters and from dynalog trajectory files. Methods: Six IMRT plans consisting of fifty individual fields were compared both field-by-field and as composite plans. These plans were exported from Eclipse to RadCalc. Each plan was then delivered in QA mode on a Varian Clinac iX, while recording and saving the dynalog trajectory files. These files were then imported into RadCalc, providing three sources of generating fluence maps: Eclipse, RadCalc, and dynalog files. Using the image analysis tool in RadCalc, a gamma analysis was then performed as a point-by-point comparison for each IMRT field. The Eclipse fluence map was used as the baseline for each comparison. Results: Based on the manufacturer’s recommendations, all fields were normalized to the maximum pixel value within each Eclipse field. All data was analyzed using a gamma index of 3mm/3% with passing criteria of 90%. For both dynalog analysis and RadCalc fluence analysis, 48/50 created fluence maps passed at greater than 90% when compared with the Eclipse baseline. In analyzing the composite of each of the six patients’ plans, all plans passed over 90%. Conclusion: Using the dynalog trajectory files in combination with RadCalc Version 6.3 image analysis is a promising metric for verifying IMRT QA passing rates. Notably, if a field failed, it did so for both dynalog and RadCalc compared to Eclipse. This suggests RadCalc can accurately simulate fluence maps, with similar results to dynalog comparisons. Limitations within the RadCalc software were discovered in the analytical process such as pixel resolution and the inability to set minimum threshold values. Comparisons could be extended to include dose map and distance to agreement analysis by expanding software capabilities.

60 APPLIED LIFE SCIENCES↗

Position-probability-sampled Monte Carlo calculation of VMAT, 3DCRT, step-shoot IMRT, and helical tomotherapy dose distributions using BEAMnrc/DOSXYZnrc

Purpose: The commercial release of volumetric modulated arc therapy techniques using a conventional linear accelerator and the growing number of helical tomotherapy users have triggered renewed interest in dose verification methods, and also in tools for exploring the impact of machine tolerance and patient motion on dose distributions without the need to approximate time-varying parameters such as gantry position, MLC leaf motion, or patient motion. To this end we have developed a Monte Carlo-based calculation method capable of simulating a wide variety of treatment techniques without the need to resort to discretization approximations. Methods: The ability to perform complete position-probability-sampled Monte Carlo dose calculations was implemented in the BEAMnrc/DOSXZYnrc user codes of EGSnrc. The method includes full accelerator head simulations of our tomotherapy and Elekta linacs, and a realistic representation of continous motion via the sampling of a time variable. The functionality of this algorithm was tested via comparisons with both measurements and treatment planning dose distributions for four types of treatment techniques: 3D conformal, step-shoot intensity modulated radiation therapy, helical tomotherapy, and volumetric modulated arc therapy. Results: For static fields, the absolute dose agreement between the EGSnrc Monte Carlo calculations and measurements is within 2%/1 mm. Absolute dose agreement between Monte Carlo calculations and treatment planning system for the four different treatment techniques is within 3%/3 mm. Discrepancies with the tomotherapy TPS on the order of 10%/5 mm were observed for the extreme example of a small target located 15 cm off-axis and planned with a low modulation factor. The increase in simulation time associated with using position-probability sampling, as opposed to the discretization approach, was less than 2% in most cases. Conclusions: A single Monte Carlo simulation method can be used to calculate patient dose distribution for various types of treatment techniques delivered with either tomotherapy or a conventional linac. The method simplifies the simulation process, improves dose calculation accuracy, and involves an acceptably small change in computation time.

60 APPLIED LIFE SCIENCES↗

A generic high-dose rate {sup 192}Ir brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism

Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) {sup 192}Ir source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR {sup 192}Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic {sup 192}Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra{sup ®} Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS{sup TM}]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201){sup 3} voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR {sup 192}Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR {sup 192}Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

07 ISOTOPE AND RADIATION SOURCES↗

Calculations vs. measurements of remnant dose rates for SNS spent structures

Residual dose rate measurements were conducted on target vessel #13 and proton beam window #5 after extraction from their service locations. These measurements were used to verify calculation methods of radionuclide inventory assessment that are typically performed for nuclear waste characterization and transportation of these structures. Neutronics analyses for predicting residual dose rates were carried out using the transport code MCNPX and the transmutation code CINDER90. For transport analyses complex and rigorous geometry model of the structures and their surrounding are applied. The neutronics analyses were carried out using Bertini and CEM high energy physics models for simulating particles interaction. Obtained preliminary calculational results were analysed and compared to the measured dose rates and overall are showing good agreement with in 40% in average.

Popova, Irina I.↗

A modern Monte Carlo investigation of the TG-43 dosimetry parameters for an {sup 125}I seed already having AAPM consensus data

Purpose: To investigate potential causes for differences in TG-43 brachytherapy dosimetry parameters in the existent literature for the model IAI-125A{sup 125}I seed and to propose new standard dosimetry parameters. Methods: The MCNP5 code was used for Monte Carlo (MC) simulations. Sensitivity of dose distributions, and subsequently TG-43 dosimetry parameters, was explored to reproduce historical methods upon which American Association of Physicists in Medicine (AAPM) consensus data are based. Twelve simulation conditions varying{sup 125}I coating thickness, coating mass density, photon interaction cross-section library, and photon emission spectrum were examined. Results: Varying{sup 125}I coating thickness, coating mass density, photon cross-section library, and photon emission spectrum for the model IAI-125A seed changed the dose-rate constant by up to 0.9%, about 1%, about 3%, and 3%, respectively, in comparison to the proposed standard value of 0.922 cGy h{sup −1} U{sup −1}. The dose-rate constant values by Solberg et al. [“Dosimetric parameters of three new solid core {sup 125}I brachytherapy sources,” J. Appl. Clin. Med. Phys. 3, 119–134 (2002)], Meigooni et al. [“Experimental and theoretical determination of dosimetric characteristics of IsoAid ADVANTAGE™ {sup 125}I brachytherapy source,” Med. Phys. 29, 2152–2158 (2002)], and Taylor and Rogers [“An EGSnrc Monte Carlo-calculated database of TG-43 parameters,” Med. Phys. 35, 4228–4241 (2008)] for the model IAI-125A seed and Kennedy et al. [“Experimental and Monte Carlo determination of the TG-43 dosimetric parameters for the model 9011 THINSeed™ brachytherapy source,” Med. Phys. 37, 1681–1688 (2010)] for the model 6711 seed were +4.3% (0.962 cGy h{sup −1} U{sup −1}), +6.2% (0.98 cGy h{sup −1} U{sup −1}), +0.3% (0.925 cGy h{sup −1} U{sup −1}), and −0.2% (0.921 cGy h{sup −1} U{sup −1}), respectively, in comparison to the proposed standard value. Differences in the radial dose functions between the current study and both Solberg et al. and Meigooni et al. were <10% for r ≤ 5 cm, and increased for r > 5 cm with a maximum difference of 29% at r = 9 cm. In comparison to Taylor and Rogers, these differences were lower (maximum of 2% at r = 9 cm). For the similarly designed model 6711 {sup 125}I seed, differences did not exceed 0.5% for 0.5 ≤ r ≤ 10 cm. Radial dose function values varied by 1% as coating thickness and coating density were changed. Varying the cross-section library and source spectrum altered the radial dose function by 25% and 12%, respectively, but these differences occurred at r = 10 cm where the dose rates were very low. The 2D anisotropy function results were most similar to those of Solberg et al. and most different to those of Meigooni et al. The observed order of simulation condition variables from most to least important for influencing the 2D anisotropy function was spectrum, coating thickness, coating density, and cross-section library. Conclusions: Several MC radiation transport codes are available for calculation of the TG-43 dosimetry parameters for brachytherapy seeds. The physics models in these codes and their related cross-section libraries have been updated and improved since publication of the 2007 AAPM TG-43U1S1 report. Results using modern data indicated statistically significant differences in these dosimetry parameters in comparison to data recommended in the TG-43U1S1 report. Therefore, it seems that professional societies such as the AAPM should consider reevaluating the consensus data for this and others seeds and establishing a process of regular evaluations in which consensus data are based upon methods that remain state-of-the-art.

61 RADIATION PROTECTION AND DOSIMETRY↗

Computation of Cosmic Ray Ionization and Dose at Mars: a Comparison of HZETRN and Planetocosmics for Proton and Alpha Particles

The ability to evaluate the cosmic ray environment at Mars is of interest for future manned exploration. To support exploration, tools must be developed to accurately access the radiation environment in both free space and on planetary surfaces. The primary tool NASA uses to quantify radiation exposure behind shielding materials is the space radiation transport code, HZETRN. In order to build confidence in HZETRN, code benchmarking against Monte Carlo radiation transport codes is often used. This work compares the dose calculations at Mars by HZETRN and the Geant4 application Planetocosmics. The dose at ground and the energy deposited in the atmosphere by galactic cosmic ray protons and alpha particles has been calculated for the Curiosity landing conditions. In addition, this work has considered Solar Energetic Particle events, allowing for the comparison of varying input radiation environments. The results for protons and alpha particles show very good agreement between HZETRN and Planetocosmics.

Gronoff, Guillaume↗

SU-F-T-618: Evaluation of a Mono-Isocentric Treatment Planning Software for Stereotactic Radiosurgery of Multiple Brain Metastases

Purpose: To evaluate planning performance of an automated treatment planning software (BrainLAB; Elements) for stereotactic radiosurgery (SRS) of multiple brain metastases. Methods: Brainlab’s Multiple Metastases Elements (MME) uses single isocentric technique to treat up to 10 cranial planning target volumes (PTVs). The planning algorithm of the MME accounts for multiple PTVs overlapping with one another on the beam eyes view (BEV) and automatically selects a subset of all overlapping PTVs on each arc for sparing normal tissues in the brain. The algorithm also optimizes collimator angles, margins between multi-leaf collimators (MLCs) and PTVs, as well as monitor units (MUs) using minimization of conformity index (CI) for all targets. Planning performance was evaluated by comparing the MME-calculated treatment plan parameters with the same parameters calculated with the Volumetric Modulated Arc Therapy (VMAT) optimization on Varian’s Eclipse platform. Results: Figures 1 to 3 compare several treatment plan outcomes calculated between the MME and VMAT for 5 clinical multi-targets SRS patient plans. Prescribed target dose was volume-dependent and defined based on the RTOG recommendation. For a total number of 18 PTV’s, mean values for the CI, PITV, and GI were comparable between the MME and VMAT within one standard deviation (σ). However, MME-calculated MDPD was larger than the same VMAT-calculated parameter. While both techniques delivered similar maximum point doses to the critical cranial structures and total MU’s for the 5 patient plans, the MME required less treatment planning time by an order of magnitude compared to VMAT. Conclusion: The MME and VMAT produce similar plan qualities in terms of MUs, target dose conformation, and OAR dose sparing. While the selective use of PTVs for arc-optimization with the MME reduces significantly the total planning time in comparison to VMAT, the target dose homogeneity was also compromised due to its simplified inverse planning algorithm used.

60 APPLIED LIFE SCIENCES↗

Implementation and verification of PyNE R2S with DAG-OpenMC

The mesh-based Rigorous-Two-Step (R2S) method has been widely used in the accurate estimation of the shutdown dose rate (SDR) of fusion systems. Several mesh-based R2S code has been implemented based on Monte Carlo particle transport code MCNP5 or DAG-MCNP5 and inventory calculation code such as FISPACT-II, ALARA or ACAB. Recently, the OpenMC, a community-developed Monte Carlo particle transport code, implemented the photon particle transport capability, which enables OpenMC to be used in R2S. In this paper, PyNE R2S is extended with the support of DAG-OpenMC, i.e., OpenMC with DAGMC. Modifications that allow PyNE to read the flux result in the state point file of OpenMC has been implemented to PyNE R2S workflow. A custom source routine of OpenMC that reads the photon source file from PyNE R2S and sampling photon particle in OpenMC photon transport has been implemented in OpenMC. With these modifications, PyNE R2S is now capable of performing R2S calculations with both DAG-MCNP5 and DAG-OpenMC. The ITER FNG dose rate benchmark problem has been used to validate the code. The FNG neutron source has been modified for neutron transport with DAG-OpenMC. The shutdown dose rate of 19 cooling times has been calculated and compared with the experimental data and the computational results of other R2S codes. The results of PyNE R2S with DAG-OpenMC show satisfactory agreement with the experimental data and other computational results. Therefore, we consider PyNE R2S with DAG-OpenMC is a reliable code to calculate the SDR of fusion systems.

fusion↗

Minibeam radiation therapy for the management of osteosarcomas: A Monte Carlo study

Purpose: Minibeam radiation therapy (MBRT) exploits the well-established tissue-sparing effect provided by the combination of submillimetric field sizes and a spatial fractionation of the dose. The aim of this work is to evaluate the feasibility and potential therapeutic gain of MBRT, in comparison with conventional radiotherapy, for osteosarcoma treatments. Methods: Monte Carlo simulations (PENELOPE/PENEASY code) were used as a method to study the dose distributions resulting from MBRT irradiations of a rat femur and a realistic human femur phantoms. As a figure of merit, peak and valley doses and peak-to-valley dose ratios (PVDR) were assessed. Conversion of absorbed dose to normalized total dose (NTD) was performed in the human case. Several field sizes and irradiation geometries were evaluated. Results: It is feasible to deliver a uniform dose distribution in the target while the healthy tissue benefits from a spatial fractionation of the dose. Very high PVDR values (⩾20) were achieved in the entrance beam path in the rat case. PVDR values ranged from 2 to 9 in the human phantom. NTD{sub 2.0} of 87 Gy might be reached in the tumor in the human femur while the healthy tissues might receive valley NTD{sub 2.0} lower than 20 Gy. The doses in the tumor and healthy tissues might be significantly higher and lower than the ones commonly delivered used in conventional radiotherapy. Conclusions: The obtained dose distributions indicate that a gain in normal tissue sparing might be expected. This would allow the use of higher (and potentially curative) doses in the tumor. Biological experiments are warranted.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-E-T-786: Utility of Gold Wires to Optimize Intensity Modulation Capacity of a Novel Directional Modulated Brachytherapy Tandem Applicator for Image Guided Cervical Cancer Brachytherapy

Purpose: To evaluate the impact of using gold wires to differentially fill various channels on plan quality compared with conventional T&R applicator, inside a novel directional modulated brachytherapy (DMBT) tandem applicator for cervical cancer brachytherapy. Materials and Methods: The novel DMBT tandem applicator has a 5.4-mm diameter MR-compatible tungsten alloy enclosed in a 0.3-mm thick plastic tubing that wraps around the tandem. To modulate the radiation intensity, 6 symmetric peripheral holes of 1.3-mm diameter are grooved along the tungsten alloy rod. These grooved holes are differentially filled with gold wires to generate various degrees of directional beams. For example, three different fill patterns of 1) all void, 2) all filled except the hole containing the 192-Ir source, and 3) two adjacent holes to the 192-Ir source filled were Monte Carlo simulated. The resulting 3D dose distributions were imported into an in-house-coded inverse optimization planning system to generate HDR brachytherapy clinical plans for 19 patient cases. All plans generated were normalized to the same D90 as the clinical plans and D2cc doses of OARs were evaluated. Prescription ranged between 15 and 17.5Gy. Results: In general, the plans in case 1) resulted in the highest D2cc doses for the OARs with 11.65±2.30Gy, 7.47±3.05Gy, and 9.84±2.48Gy for bladder, rectum, and sigmoid, respectively, although the differences were small. For the case 2), D2cc doses were 11.61±2.29Gy, 7.41±3.07Gy, and 9.75±2.45Gy, respectively. And, for the case 3), D2cc doses were 11.60±2.28Gy, 7.41±3.05Gy, and 9.74±2.45Gy, respectively. Difference between 1) and 2) cases were small with the average D2cc difference of <0.64%. Difference between 1) and 3) cases were even smaller with the average D2cc difference of <0.1%. Conclusions: There is a minimal clinical benefit by differentially filling grooved holes in the novel DMBT tandem applicator for image guided cervical cancer brachytherapy.

60 APPLIED LIFE SCIENCES↗

Commissioning of a Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Monte Carlo modelling of a linear accelerator is the first and most important step in Monte Carlo dose calculations in radiotherapy. Monte Carlo is considered today to be the most accurate and detailed calculation method in different fields of medical physics. In this research, we developed a photon beam model for Varian Clinac iX 6 MV equipped with MilleniumMLC120 for dose calculation purposes using BEAMnrc/DOSXYZnrc Monte Carlo system based on the underlying EGSnrc particle transport code. Monte Carlo simulation for this commissioning head LINAC divided in two stages are design head Linac model using BEAMnrc, characterize this model using BEAMDP and analyze the difference between simulation and measurement data using DOSXYZnrc. In the first step, to reduce simulation time, a virtual treatment head LINAC was built in two parts (patient-dependent component and patient-independent component). The incident electron energy varied 6.1 MeV, 6.2 MeV and 6.3 MeV, 6.4 MeV, and 6.6 MeV and the FWHM (full width at half maximum) of source is 1 mm. Phase-space file from the virtual model characterized using BEAMDP. The results of MC calculations using DOSXYZnrc in water phantom are percent depth doses (PDDs) and beam profiles at depths 10 cm were compared with measurements. This process has been completed if the dose difference of measured and calculated relative depth-dose data along the central-axis and dose profile at depths 10 cm is ≤ 5%. The effect of beam width on percentage depth doses and beam profiles was studied. Results of the virtual model were in close agreement with measurements in incident energy electron 6.4 MeV. Our results showed that photon beam width could be tuned using large field beam profile at the depth of maximum dose. The Monte Carlo model developed in this study accurately represents the Varian Clinac iX with millennium MLC 120 leaf and can be used for reliable patient dose calculations. In this commissioning process, the good criteria of dose difference in PDD and dose profiles were achieve using incident electron energy 6.4 MeV.

61 RADIATION PROTECTION AND DOSIMETRY↗

Radiation doses in volume-of-interest breast computed tomography—A Monte Carlo simulation study

Purpose: Cone beam breast computed tomography (breast CT) with true three-dimensional, nearly isotropic spatial resolution has been developed and investigated over the past decade to overcome the problem of lesions overlapping with breast anatomical structures on two-dimensional mammographic images. However, the ability of breast CT to detect small objects, such as tissue structure edges and small calcifications, is limited. To resolve this problem, the authors proposed and developed a volume-of-interest (VOI) breast CT technique to image a small VOI using a higher radiation dose to improve that region’s visibility. In this study, the authors performed Monte Carlo simulations to estimate average breast dose and average glandular dose (AGD) for the VOI breast CT technique. Methods: Electron–Gamma-Shower system code-based Monte Carlo codes were used to simulate breast CT. The Monte Carlo codes estimated were validated using physical measurements of air kerma ratios and point doses in phantoms with an ion chamber and optically stimulated luminescence dosimeters. The validated full cone x-ray source was then collimated to simulate half cone beam x-rays to image digital pendant-geometry, hemi-ellipsoidal, homogeneous breast phantoms and to estimate breast doses with full field scans. 13-cm in diameter, 10-cm long hemi-ellipsoidal homogeneous phantoms were used to simulate median breasts. Breast compositions of 25% and 50% volumetric glandular fractions (VGFs) were used to investigate the influence on breast dose. The simulated half cone beam x-rays were then collimated to a narrow x-ray beam with an area of 2.5 × 2.5 cm{sup 2} field of view at the isocenter plane and to perform VOI field scans. The Monte Carlo results for the full field scans and the VOI field scans were then used to estimate the AGD for the VOI breast CT technique. Results: The ratios of air kerma ratios and dose measurement results from the Monte Carlo simulation to those from the physical measurements were 0.97 ± 0.03 and 1.10 ± 0.13, respectively, indicating that the accuracy of the Monte Carlo simulation was adequate. The normalized AGD with VOI field scans was substantially reduced by a factor of about 2 over the VOI region and by a factor of 18 over the entire breast for both 25% and 50% VGF simulated breasts compared with the normalized AGD with full field scans. The normalized AGD for the VOI breast CT technique can be kept the same as or lower than that for a full field scan with the exposure level for the VOI field scan increased by a factor of as much as 12. Conclusions: The authors’ Monte Carlo estimates of normalized AGDs for the VOI breast CT technique show that this technique can be used to markedly increase the dose to the breast and thus the visibility of the VOI region without increasing the dose to the breast. The results of this investigation should be helpful for those interested in using VOI breast CT technique to image small calcifications with dose concern.

60 APPLIED LIFE SCIENCES↗

SU-F-BRA-05: Utility of the Combined Use of Two Types of HDR Sources with the Direction Modulation Brachytherapy (DMBT) Tandem Applicator for Cervical Cancer Treatment

Purpose: To maximize the dose to HRCTV while minimizing dose to the OARs, the combination of two HDR brachytherapy sources, 192-Ir and 169-Yb, used in combination with the recently-proposed novel direction modulated brachytherapy (DMBT) tandem applicator were examined. Methods: The DMBT tandem, made from nonmagnetic tungsten-alloy rod, with diameter of 5.4mm, has 6 symmetric peripheral holes of 1.3mm diameter. The 0.3mm thick bio-compatible plastic tubing wraps the tandem. MCNPX v.2.6 was used to simulate the mHDR 192-Ir V2 and 4140 HDR 169-Yb sources inside the DMBT applicator. Thought was by combining the higher energy 192-Ir (380keV) and lower energy 169-Yb (92.7keV) sources could create unprecedented level of dose conformality when combined with the high-degree intensity modulation capable DMBT tandem applicator. 3D dose matrices, with 1 mm3 resolution, were imported into an in-house-coded inverse optimization planning system to evaluate plan quality of 19 clinical patient cases. Prescription dose was 15Gy. All plans were normalized to receive the same HRCTV D90. Results: Generally, the use of dual sources produced better plans than using either of the sources alone, with significantly better performance in some patients. The mean D2cc for bladder, rectum, and sigmoid were 11.65±2.30Gy, 7.47±3.05Gy, and 9.84±2.48Gy for 192-Ir-only, respectively. For 169 -Yb-only, they were 11.67±2.26Gy, 7.44±3.02Gy, and 9.83±2.38Gy, respectively. The corresponding data for the dual sources were 11.51±2.24Gy, 7.30±3.00Gy, and 9.68 ±2.39Gy, respectively. The HRCTV D98 and V100 were 16.37±1.86Gy and 97.37±1.92Gy for Ir-192-only, respectively. For 169-Yb-only, they were 16.43±1.86Gy, and 97.51±1.91Gy, respectively. For the dual source, they were 16.42±1.87Gy and 97.47±1.93Gy, respectively. Conclusion: The plan quality improves, in some cases quite significantly, for when dual 192-Ir and 169-Yb sources are used in combination with highly intensity modulation capable DMBT tandem applicator for image guided cervical cancer brachytherapy.

60 APPLIED LIFE SCIENCES↗

Technical Note: Study of the electron transport parameters used in PENELOPE for the Monte Carlo simulation of Linac targets

Purpose: The Monte Carlo simulation of electron transport in Linac targets using the condensed history technique is known to be problematic owing to a potential dependence of absorbed dose distributions on the electron step length. In the PENELOPE code, the step length is partially determined by the transport parameters C1 and C2. The authors have investigated the effect on the absorbed dose distribution of the values given to these parameters in the target. Methods: A monoenergetic 6.26 MeV electron pencil beam from a point source was simulated impinging normally on a cylindrical tungsten target. Electrons leaving the tungsten were discarded. Radial absorbed dose profiles were obtained at 1.5 cm of depth in a water phantom located at 100 cm for values of C1 and C2 in the target both equal to 0.1, 0.01, or 0.001. A detailed simulation case was also considered and taken as the reference. Additionally, lateral dose profiles were estimated and compared with experimental measurements for a 6 MV photon beam of a Varian Clinac 2100 for the cases of C1 and C2 both set to 0.1 or 0.001 in the target. Results: On the central axis, the dose obtained for the case C1 = C2 = 0.1 shows a deviation of (17.2% ± 1.2%) with respect to the detailed simulation. This difference decreases to (3.7% ± 1.2%) for the case C1 = C2 = 0.01. The case C1 = C2 = 0.001 produces a radial dose profile that is equivalent to that of the detailed simulation within the reached statistical uncertainty of 1%. The effect is also appreciable in the crossline dose profiles estimated for the realistic geometry of the Linac. In another simulation, it was shown that the error made by choosing inappropriate transport parameters can be masked by tuning the energy and focal spot size of the initial beam. Conclusions: The use of large path lengths for the condensed simulation of electrons in a Linac target with PENELOPE conducts to deviations of the dose in the patient or phantom. Based on the results obtained in this work, values of C1 and C2 larger than 0.001 should not be used in Linac targets without further investigation.

07 ISOTOPE AND RADIATION SOURCES↗

Monte Carlo calculations of k{sub Q}, the beam quality conversion factor

Purpose: To use EGSnrc Monte Carlo simulations to directly calculate beam quality conversion factors, k{sub Q}, for 32 cylindrical ionization chambers over a range of beam qualities and to quantify the effect of systematic uncertainties on Monte Carlo calculations of k{sub Q}. These factors are required to use the TG-51 or TRS-398 clinical dosimetry protocols for calibrating external radiotherapy beams. Methods: Ionization chambers are modeled either from blueprints or manufacturers' user's manuals. The dose-to-air in the chamber is calculated using the EGSnrc user-code egs{sub c}hamber using 11 different tabulated clinical photon spectra for the incident beams. The dose to a small volume of water is also calculated in the absence of the chamber at the midpoint of the chamber on its central axis. Using a simple equation, k{sub Q} is calculated from these quantities under the assumption that W/e is constant with energy and compared to TG-51 protocol and measured values. Results: Polynomial fits to the Monte Carlo calculated k{sub Q} factors as a function of beam quality expressed as %dd(10){sub x} and TPR{sub 10}{sup 20} are given for each ionization chamber. Differences are explained between Monte Carlo calculated values and values from the TG-51 protocol or calculated using the computer program used for TG-51 calculations. Systematic uncertainties in calculated k{sub Q} values are analyzed and amount to a maximum of one standard deviation uncertainty of 0.99% if one assumes that photon cross-section uncertainties are uncorrelated and 0.63% if they are assumed correlated. The largest components of the uncertainty are the constancy of W/e and the uncertainty in the cross-section for photons in water. Conclusions: It is now possible to calculate k{sub Q} directly using Monte Carlo simulations. Monte Carlo calculations for most ionization chambers give results which are comparable to TG-51 values. Discrepancies can be explained using individual Monte Carlo calculations of various correction factors which are more accurate than previously used values. For small ionization chambers with central electrodes composed of high-Z materials, the effect of the central electrode is much larger than that for the aluminum electrodes in Farmer chambers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗