Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ray trace”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Two-Layer Anti-Reflection Coatings with Optimized Sub-Bandgap Reflection for Solar Modules

The efficiency of a crystalline silicon solar module decreases as its operating temperature rises. Module cooling is possible via selective reflection of sub-bandgap photons so that they are not parasitically absorbed. Selecting from a library of common dielectrics, we numerically optimize the design of two-layer mirrors at the outer glass surface of a crystalline Si solar cell module. The mirrors are designed to maximize the annual energy yield of a module by both reflecting light below the bandgap and enhancing the transmission of light above the bandgap. Combined ray-tracing and finite element simulations determine the power output and temperature of the module over time. Since any two-layer mirror would replace a conventional single-layer glass anti-reflection coating on the module glass, we study the ability of a two-layer structure to improve on a single-layer coating. The best two-layer designs improve the annual energy yield over a module with a glass anti-reflection coating and reduce the module operating temperature.

14 SOLAR ENERGY↗
DOE OSTI2019

Upgrade to the SHARP EUV mask microscope

© 2019 SPIE. The Sharp High-NA Actinic Reticle review Project (SHARP) is a synchrotron-based, extreme ultraviolet (EUV) microscope dedicated to photomask research. A potential upgrade to the SHARP microscope is presented. The upgrade includes changing the light path in the instrument from its current off-Axis configuration to an on-Axis configuration. This change allows for an increased working distance of 2.5 mm or more. A central obscuration, added to the zoneplate aperture, blocks stray light from reaching the central part of the image, thus improving the image contrast. The imaging performance of the two configurations is evaluated by means of ray tracing.

Benk, M↗
DOE OSTI2020

Analytical expressions of the surface shape of ‘diaboloid’ mirrors

Modern deterministic polishing processes allow fabrication of x-ray optics with almost any arbitrary aspherical surface shape. Among these optics, the so called "diaboloid"mirror is of special interest. The diaboloid mirror that converts a cylindrical wave to a spherical wave would improve focusing in x-ray beamlines implementing a diffraction element between a parabolic cylinder and a toroidal mirror. The replacement of the toroidal mirror in existing beamlines by the diaboloid mirror would mitigate aberrations. The shape of the diaboloid mirror is usually calculated numerically based on a truncated polynomial solution of the optical path problem. Here, we present an exact analytical solution for the shape of a diaboloid mirror as a function of the conjugate parameters of the mirror placed in a beamline. The derived analytical expressions for the diaboloid mirror in both the canonical and mirror-based coordinate systems are implemented in ray-tracing simulations to verify the beamline performances.

Yashchuk, Valeriy V.↗
DOE OSTI2012

Wide-Field-of-View Millimeter-Wave Telescope Design with Ultra-Low Cross-Polarization

As millimeter-wave arrays become available, off-axis imaging performance of the fore optics increases in importance due to the relatively large physical extent of the arrays. Typically, simple optical telescope designs are adapted to millimeter-wave imaging but single-mirror spherical or classic conic designs cannot deliver adequate image quality except near the optical axis. Since most millimeter-wave designs are quasi-optical, optical ray tracing and commercial design software can be used to optimize designs to improve off-axis imaging as well as minimize cross-polarization. Methods that obey the Dragone-Mizuguchi condition for the design of reflective millimeter-wave telescopes with low cross-polarization also provide additional degrees of freedom that offer larger fields of view than possible with single-reflector designs. Dragone’s graphical design method does not lend itself readily to computer-based optical design approaches, but subsequent authors expanded on Dragone’s geometric design approach with analytic expressions that describe the location, shape, off-axis height and tilt of the telescope elements that satisfy Dragone’s design rules and can be used as a first-order design for subsequent computer-based design and optimization. We investigate two design variants that obey the Dragone-Mizuguchi conditions that exhibit ultra-low polarization crosstalk and a large diffraction-limited field of view well suited to millimeter-wave imaging arrays.

Millimeter-Wave Optical Design↗
DOE OSTI2013

Experimental evaluations of the accuracy of 3D and 4D planning in robotic tracking stereotactic body radiotherapy for lung cancers

Purpose: Due to the complexity of 4D target tracking radiotherapy, the accuracy of this treatment strategy should be experimentally validated against established standard 3D technique. This work compared the accuracy of 3D and 4D dose calculations in respiration tracking stereotactic body radiotherapy (SBRT). Methods: Using the 4D planning module of the CyberKnife treatment planning system, treatment plans for a moving target and a static off-target cord structure were created on different four-dimensional computed tomography (4D-CT) datasets of a thorax phantom moving in different ranges. The 4D planning system used B-splines deformable image registrations (DIR) to accumulate dose distributions calculated on different breathing geometries, each corresponding to a static 3D-CT image of the 4D-CT dataset, onto a reference image to compose a 4D dose distribution. For each motion, 4D optimization was performed to generate a 4D treatment plan of the moving target. For comparison with standard 3D planning, each 4D plan was copied to the reference end-exhale images and a standard 3D dose calculation was followed. Treatment plans of the off-target structure were first obtained by standard 3D optimization on the end-exhale images. Subsequently, they were applied to recalculate the 4D dose distributions using DIRs. All dose distributions that were initially obtained using the ray-tracing algorithm with equivalent path-length heterogeneity correction (3D{sub EPL} and 4D{sub EPL}) were recalculated by a Monte Carlo algorithm (3D{sub MC} and 4D{sub MC}) to further investigate the effects of dose calculation algorithms. The calculated 3D{sub EPL}, 3D{sub MC}, 4D{sub EPL}, and 4D{sub MC} dose distributions were compared to measurements by Gafchromic EBT2 films in the axial and coronal planes of the moving target object, and the coronal plane for the static off-target object based on the {gamma} metric at 5%/3mm criteria ({gamma}{sub 5%/3mm}). Treatment plans were considered acceptable if the percentage of pixels passing {gamma}{sub 5%/3mm} (P{sub {gamma}<1}) {>=} 90%. Results: The averaged P{sub {gamma}<1} values of the 3D{sub EPL}, 3D{sub MC}, 4D{sub EPL}, and 4D{sub MC} dose calculation methods for the moving target plans are 95%, 95%, 94%, and 95% for reproducible motion, and 95%, 96%, 94%, and 93% for nonreproducible motion during actual treatment delivery. The overall measured target dose distributions are in better agreement with the 3D{sub MC} dose distributions than the 4D{sub MC} dose distributions. Conversely, measured dose distributions agree much better with the 4D{sub EPL/MC} than the 3D{sub EPL/MC} dose distributions in the static off-target structure, resulting in higher P{sub {gamma}<1} values with 4D{sub EPL/MC} (91%) vs 3D{sub EPL} (24%) and 3D{sub MC} (25%). Systematic changes of target motion reduced the averaged P{sub {gamma}<1} to 47% and 53% for 4D{sub EPL} and 4D{sub MC} dose calculations, and 22% for 3D{sub EPL/MC} dose calculations in the off-target films. Conclusions: In robotic tracking SBRT, 4D treatment planning was found to yield better prediction of the dose distributions in the off-target structure, but not necessarily in the moving target, compared to standard 3D treatment planning, for reproducible and nonreproducible target motion. It is important to ensure on a patient-by-patient basis that the cumulative uncertainty associated with the 4D-CT artifacts, deformable image registration, and motion variability is significantly smaller than the cumulative uncertainty occurred in standard 3D planning in order to make 4D planning a justified option.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

HDRMC, an accelerated Monte Carlo dose calculator for high dose rate brachytherapy with CT-compatible applicators

Purpose: To present a new accelerated Monte Carlo code for CT-based dose calculations in high dose rate (HDR) brachytherapy. The new code (HDRMC) accounts for both tissue and nontissue heterogeneities (applicator and contrast medium). Methods: HDRMC uses a fast ray-tracing technique and detailed physics algorithms to transport photons through a 3D mesh of voxels representing the patient anatomy with applicator and contrast medium included. A precalculated phase space file for the{sup 192}Ir source is used as source term. HDRM is calibrated to calculated absolute dose for real plans. A postprocessing technique is used to include the exact density and composition of nontissue heterogeneities in the 3D phantom. Dwell positions and angular orientations of the source are reconstructed using data from the treatment planning system (TPS). Structure contours are also imported from the TPS to recalculate dose-volume histograms. Results: HDRMC was first benchmarked against the MCNP5 code for a single source in homogenous water and for a loaded gynecologic applicator in water. The accuracy of the voxel-based applicator model used in HDRMC was also verified by comparing 3D dose distributions and dose-volume parameters obtained using 1-mm{sup 3} versus 2-mm{sup 3} phantom resolutions. HDRMC can calculate the 3D dose distribution for a typical HDR cervix case with 2-mm resolution in 5 min on a single CPU. Examples of heterogeneity effects for two clinical cases (cervix and esophagus) were demonstrated using HDRMC. The neglect of tissue heterogeneity for the esophageal case leads to the overestimate of CTV D90, CTV D100, and spinal cord maximum dose by 3.2%, 3.9%, and 3.6%, respectively. Conclusions: A fast Monte Carlo code for CT-based dose calculations which does not require a prebuilt applicator model is developed for those HDR brachytherapy treatments that use CT-compatible applicators. Tissue and nontissue heterogeneities should be taken into account in modern HDR brachytherapy planning.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-J-60: Efficient Monte Carlo Dose Calculation On CPU-GPU Heterogeneous Systems

Purpose: It is well-known that the performance of GPU-based Monte Carlo dose calculation implementations is bounded by memory bandwidth. One major cause of this bottleneck is the random memory writing patterns in dose deposition, which leads to several memory efficiency issues on GPU such as un-coalesced writing and atomic operations. We propose a new method to alleviate such issues on CPU-GPU heterogeneous systems, which achieves overall performance improvement for Monte Carlo dose calculation. Methods: Dose deposition is to accumulate dose into the voxels of a dose volume along the trajectories of radiation rays. Our idea is to partition this procedure into the following three steps, which are fine-tuned for CPU or GPU: (1) each GPU thread writes dose results with location information to a buffer on GPU memory, which achieves fully-coalesced and atomic-free memory transactions; (2) the dose results in the buffer are transferred to CPU memory; (3) the dose volume is constructed from the dose buffer on CPU. We organize the processing of all radiation rays into streams. Since the steps within a stream use different hardware resources (i.e., GPU, DMA, CPU), we can overlap the execution of these steps for different streams by pipelining. Results: We evaluated our method using a Monte Carlo Convolution Superposition (MCCS) program and tested our implementation for various clinical cases on a heterogeneous system containing an Intel i7 quad-core CPU and an NVIDIA TITAN GPU. Comparing with a straightforward MCCS implementation on the same system (using both CPU and GPU for radiation ray tracing), our method gained 2-5X speedup without losing dose calculation accuracy. Conclusion: The results show that our new method improves the effective memory bandwidth and overall performance for MCCS on the CPU-GPU systems. Our proposed method can also be applied to accelerate other Monte Carlo dose calculation approaches. This research was supported in part by NSF under Grants CCF-1217906, and also in part by a research contract from the Sandia National Laboratories.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-215: Onboard Split PET Including the Effects of Attenuation, Scatter, and Random Events

Purpose: Positron emission tomography (PET) using split PET geometry was investigated as an on-board system for functional imaging and PET marker tracking, specifically with tomotherapy. The open dual ring PET would allow measurement of both inter and intra-fractional variation, improving the delineation of tumor volume at any stage in the radiation treatment delivery process. We present results from data obtained using Monte Carlo simulations including the effects of attenuation, random events, and scatter. Methods: PET design was accomplished via Monte Carlo simulations with GATE, the Geant4 Application for Emission Tomography. Images were reconstructed using Software for Image Reconstruction (STIR) with a fully 3D ordered subset expectation maximization (OSEM) image reconstruction technique. Results: Monte Carlo simulations of the split PET geometry indicate that including physical factors that degrade image quality such as attenuation, random events, and scatter still prove the feasibility of near real-time PET for both inter and intra-fractional radiation delivery. The image quality of scan times under 1 minute reveals that it is possible to utilize PET scanning and reconstruction during the treatment session intrafractionally. GATE also simulates the depth information but does not correct obliqueness of the path of line of response so the data is much more realistic than the data obtained with ray tracing. Conclusion: Onboard split PET with TomoTherapy can generate quality images under 1 minute scan times without the need to correct for attenuation, scatter, random events, or depth information of the interaction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-T-294: Simulations to Investigate the Feasibility of ‘dry’ Optical-CT Imaging for 3D Dosimetry

Purpose: To perform simulations investigating the feasibility of “dry” optical-CT, and determine optimal design and scanning parameters for a novel dry tank telecentric optical-CT 3D dosimetry system. Such a system would have important advantages in terms of practical convenience and reduced cost. Methods: A Matlab based ray-tracing simulation platform, ScanSim, was used to model a telecentric system with a polyurethane dry tank, cylindrical dosimeter, and surrounding fluid. This program's capabilities were expanded for the geometry and physics of dry scanning. To categorize the effects of refractive index (RI) mismatches, simulations were run for several dosimeter (RI = 1.5−1.48) and fluid (RI = 1.55−1.33) combinations. Additional simulations examined the effect of increasing gap size (1–5mm) between the dosimeter and tank wall, and of changing the telecentric lens tolerance (0.5°−5°). The evaluation metric is the usable radius; the distance from the dosimeter center where the measured and true doses differ by less than 2%. Results: As the tank/dosimeter RI mismatch increases from 0–0.02, the usable radius decreases from 97.6% to 50.2%. The fluid RI for matching is lower than either the tank or dosimeter RI. Changing gap sizes has drastic effects on the usable radius, requiring more closely matched fluid at large gap sizes. Increasing the telecentric tolerance through a range from 0.5°–5.0° improved the usable radius for every combination of media. Conclusion: Dry optical-CT with telecentric lenses is feasible when the dosimeter and tank RIs are closely matched (<0.01 difference), or when data in the periphery is not required. The ScanSim tool proved very useful in situations when the tank and dosimeter have slight differences in RI by enabling estimation of the optimal choice of RI of the small amount of fluid still required. Some spoiling of the telecentric beam and increasing the tolerance helps recover the usable radius.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Prospective estimation of organ dose in CT under tube current modulation

Purpose: Computed tomography (CT) has been widely used worldwide as a tool for medical diagnosis and imaging. However, despite its significant clinical benefits, CT radiation dose at the population level has become a subject of public attention and concern. In this light, optimizing radiation dose has become a core responsibility for the CT community. As a fundamental step to manage and optimize dose, it may be beneficial to have accurate and prospective knowledge about the radiation dose for an individual patient. In this study, the authors developed a framework to prospectively estimate organ dose for chest and abdominopelvic CT exams under tube current modulation (TCM). Methods: The organ dose is mainly dependent on two key factors: patient anatomy and irradiation field. A prediction process was developed to accurately model both factors. To model the anatomical diversity and complexity in the patient population, the authors used a previously developed library of computational phantoms with broad distributions of sizes, ages, and genders. A selected clinical patient, represented by a computational phantom in the study, was optimally matched with another computational phantom in the library to obtain a representation of the patient’s anatomy. To model the irradiation field, a previously validated Monte Carlo program was used to model CT scanner systems. The tube current profiles were modeled using a ray-tracing program as previously reported that theoretically emulated the variability of modulation profiles from major CT machine manufacturers Li et al., [Phys. Med. Biol. 59, 4525–4548 (2014)]. The prediction of organ dose was achieved using the following process: (1) CTDI{sub vol}-normalized-organ dose coefficients (h{sub organ}) for fixed tube current were first estimated as the prediction basis for the computational phantoms; (2) each computation phantom, regarded as a clinical patient, was optimally matched with one computational phantom in the library; (3) to account for the effect of the TCM scheme, a weighted organ-specific CTDI{sub vol} [denoted as (CTDI{sub vol}){sub organ,weighted}] was computed for each organ based on the TCM profile and the anatomy of the “matched” phantom; (4) the organ dose was predicted by multiplying the weighted organ-specific CTDI{sub vol} with the organ dose coefficients (h{sub organ}). To quantify the prediction accuracy, each predicted organ dose was compared with the corresponding organ dose simulated from the Monte Carlo program with the TCM profile explicitly modeled. Results: The predicted organ dose showed good agreements with the simulated organ dose across all organs and modulation profiles. The average percentage error in organ dose estimation was generally within 20% across all organs and modulation profiles, except for organs located in the pelvic and shoulder regions. For an average CTDI{sub vol} of a CT exam of 10 mGy, the average error at full modulation strength (α = 1) across all organs was 0.91 mGy for chest exams, and 0.82 mGy for abdominopelvic exams. Conclusions: This study developed a quantitative model to predict organ dose for clinical chest and abdominopelvic scans. Such information may aid in the design of optimized CT protocols in relation to a targeted level of image quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Investigations into the feasibility of optical-CT 3D dosimetry with minimal use of refractively matched fluids

Purpose: In optical-CT, the use of a refractively matched polyurethane solid-tank in place of a fluid bath has the potential to greatly increase practical convenience, reduce cost, and possibly improve the efficacy of flood corrections. This work investigates the feasibility of solid-tank optical-CT imaging for 3D dosimetry through computer simulation. Methods: A MATLAB ray-tracing simulation platform, ScanSim, was used to model a parallel-source telecentric optical-CT imaging system through a polyurethane solid-tank containing a central cylindrical hollow into which PRESAGE radiochromic dosimeters can be placed. A small amount of fluid fills the 1–5 mm gap between the dosimeter and the walls of the tank. The use of the solid-tank reduces the required amount of fluid by approximately 97%. To characterize the efficacy of solid-tank, optical-CT scanning simulations investigated sensitivity to refractive index (RI) mismatches between dosimeter, solid-tank, and fluid, for a variety of dosimeter (RI = 1.5–1.47) and fluid (RI = 1.55–1.0) combinations. Efficacy was evaluated through the usable radius (r{sub u}) metric, defined as the fraction of the radius of the dosimeter where measured dose is predicted to be within 2% of the ground truth entered into the simulation. Additional simulations examined the effect of increasing gap size (1–5 mm) between the dosimeter and solid-tank well. The effects of changing the lens tolerance (0.5°–5.0°) were also investigated. Results: As the RI mismatch between the dosimeter and solid-tank increased from 0 to 0.02, the usable radius decreased from 97.6% to 50.2%. The optimal fluid RI decreased nonlinearly from 1.5 to 1.34 as the mismatch increased and was up to 9% lower than the tank. Media mismatches between the dosimeter and solid-tank also exacerbate the effects of changing the gap size, with no easily quantifiable relationship with usable radius. Generally, the optimal fluid RI value increases as gap size increases and is closely matched to the dosimeter at large gap sizes (>3 mm). Increasing the telecentric lens tolerance increases the usable radius for all refractive media combinations and improves the maximum usable radius of mismatched media to that of perfectly matched media for tolerances >5.0°. The maximum usable radius can be improved up to a factor of 2 when lens tolerances are small (<1.0°). Conclusions: Dry solid-tank optical-CT imaging in a telecentric system is feasible if the dosimeter RI is a close match with the solid-tank (<0.01 difference), providing accurate dose measurements within ±2% of true dose to over 80% of the dosimeter volume. In order to achieve accurate measurements over 96% of the dosimeter volume (representing out to 2 mm from the dosimeter edge), the dosimeter-tank RI mismatch must be less than 0.005. Optimal results occur when the RI of the dosimeter and tank is the same, in which case the fluid will have the same RI. If mismatches between the tank and dosimeter RI occur, the RI of the matching fluid needs to be fine tuned to achieve the highest usable radius.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

SU-E-T-226: Correction of a Standard Model-Based Dose Calculator Using Measurement Data

Purpose: To propose a hybrid method that combines advantages of the model-based and measurement-based method for independent dose calculation. Modeled-based dose calculation, such as collapsed-cone-convolution/superposition (CCCS) or the Monte-Carlo method, models dose deposition in the patient body accurately; however, due to lack of detail knowledge about the linear accelerator (LINAC) head, commissioning for an arbitrary machine is tedious and challenging in case of hardware changes. On the contrary, the measurement-based method characterizes the beam property accurately but lacks the capability of dose disposition modeling in heterogeneous media. Methods: We used a standard CCCS calculator, which is commissioned by published data, as the standard model calculator. For a given machine, water phantom measurements were acquired. A set of dose distributions were also calculated using the CCCS for the same setup. The difference between the measurements and the CCCS results were tabulated and used as the commissioning data for a measurement based calculator. Here we used a direct-ray-tracing calculator (ΔDRT). The proposed independent dose calculation consists of the following steps: 1. calculate D-model using CCCS. 2. calculate D-ΔDRT using ΔDRT. 3. combine Results: D=D-model+D-ΔDRT. Results: The hybrid dose calculation was tested on digital phantoms and patient CT data for standard fields and IMRT plan. The results were compared to dose calculated by the treatment planning system (TPS). The agreement of the hybrid and the TPS was within 3%, 3 mm for over 98% of the volume for phantom studies and lung patients. Conclusion: The proposed hybrid method uses the same commissioning data as those for the measurement-based method and can be easily extended to any non-standard LINAC. The results met the accuracy, independence, and simple commissioning criteria for an independent dose calculator.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-252: Developing a Pencil Beam Dose Calculation Algorithm for CyberKnife System

Purpose: Currently there are two dose calculation algorithms available in the Cyberknife planning system: ray-tracing and Monte Carlo, which is either not accurate or time-consuming for irregular field shaped by the MLC that was recently introduced. The purpose of this study is to develop a fast and accurate pencil beam dose calculation algorithm which can handle irregular field. Methods: A pencil beam dose calculation algorithm widely used in Linac system is modified. The algorithm models both primary (short range) and scatter (long range) components with a single input parameter: TPR{sub 20}/{sub 10}. The TPR{sub 20}/{sub 20}/{sub 10} value was first estimated to derive an initial set of pencil beam model parameters (PBMP). The agreement between predicted and measured TPRs for all cones were evaluated using the root mean square of the difference (RMSTPR), which was then minimized by adjusting PBMPs. PBMPs are further tuned to minimize OCR RMS (RMSocr) by focusing at the outfield region. Finally, an arbitrary intensity profile is optimized by minimizing RMSocr difference at infield region. To test model validity, the PBMPs were obtained by fitting to only a subset of cones (4) and applied to all cones (12) for evaluation. Results: With RMS values normalized to the dmax and all cones combined, the average RMSTPR at build-up and descending region is 2.3% and 0.4%, respectively. The RMSocr at infield, penumbra and outfield region is 1.5%, 7.8% and 0.6%, respectively. Average DTA in penumbra region is 0.5mm. There is no trend found in TPR or OCR agreement among cones or depths. Conclusion: We have developed a pencil beam algorithm for Cyberknife system. The prediction agrees well with commissioning data. Only a subset of measurements is needed to derive the model. Further improvements are needed for TPR buildup region and OCR penumbra. Experimental validations on MLC shaped irregular field needs to be performed. This work was partially supported by the National Natural Science Foundation of China (61171005) and the China Scholarship Council (CSC)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-404: Evaluation of the Effect of Spine Hardware for CyberKnife Spinal Stereotactic Radiosurgery

Purpose: Spine hardware made of high-Z materials such as titanium has the potential to affect the dose distribution around the metal rods in CyberKnife spinal stereotactic radiosurgery (SRS) treatments. The purpose of this work was to evaluate the magnitude of such effect retrospectively for clinical CyberKnife plans. Methods: The dose calculation was performed within the MultiPlan treatment planning system using the ray tracing (RT) and Monte Carlo (MC) method. A custom density model was created by extending the CT-to-Density table to titanium density of 4.5 g/cm3 with the CT number of 4095. To understand the dose perturbation caused by the titanium rod, a simple beam setup (7.5 mm IRIS collimator) was used to irradiate a mimic rod (5 mm) with overridden high density. Five patient spinal SRS cases were found chronologically from 2010 to 2015 in our institution. For each case, the hardware was contoured manually. The original plan was re-calculated using both RT and MC methods with and without rod density override without changing clinical beam parameters. Results: The simple beam irradiation shows that there is 10% dose increase at the interface because of electron backscattering and 7% decrease behind the rod because of photon attenuation. For actual clinical plans, the iso-dose lines and DVHs are almost identical (<2%) for calculations with and without density override for both RT and MC methods. However, there is a difference of more than 10% for D90 between RT and MC method. Conclusion: Although the dose perturbation around the metal rods can be as large as 10% for a single beam irradiation, for clinical treatments with complex beam composition the effect of spinal hardware to the PTV and spinal dose is minimal. As such, the MC dose algorithm without rod density override for CyberKnife spinal SRS is acceptable.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Population of 224 realistic human subject-based computational breast phantoms

Purpose: To create a database of highly realistic and anatomically variable 3D virtual breast phantoms based on dedicated breast computed tomography (bCT) data. Methods: A tissue classification and segmentation algorithm was used to create realistic and detailed 3D computational breast phantoms based on 230 + dedicated bCT datasets from normal human subjects. The breast volume was identified using a coarse three-class fuzzy C-means segmentation algorithm which accounted for and removed motion blur at the breast periphery. Noise in the bCT data was reduced through application of a postreconstruction 3D bilateral filter. A 3D adipose nonuniformity (bias field) correction was then applied followed by glandular segmentation using a 3D bias-corrected fuzzy C-means algorithm. Multiple tissue classes were defined including skin, adipose, and several fractional glandular densities. Following segmentation, a skin mask was produced which preserved the interdigitated skin, adipose, and glandular boundaries of the skin interior. Finally, surface modeling was used to produce digital phantoms with methods complementary to the XCAT suite of digital human phantoms. Results: After rejecting some datasets due to artifacts, 224 virtual breast phantoms were created which emulate the complex breast parenchyma of actual human subjects. The volume breast density (with skin) ranged from 5.5% to 66.3% with a mean value of 25.3% ± 13.2%. Breast volumes ranged from 25.0 to 2099.6 ml with a mean value of 716.3 ± 386.5 ml. Three breast phantoms were selected for imaging with digital compression (using finite element modeling) and simple ray-tracing, and the results show promise in their potential to produce realistic simulated mammograms. Conclusions: This work provides a new population of 224 breast phantoms based on in vivo bCT data for imaging research. Compared to previous studies based on only a few prototype cases, this dataset provides a rich source of new cases spanning a wide range of breast types, volumes, densities, and parenchymal patterns.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-F-J-201: Validation Study of Proton Radiography Against CT Data for Quantitative Imaging of Anatomical Changes in Head and Neck Patients

Purpose: In clinical pencil-beam-scanned (PBS) proton therapy, the advantage of the characteristic sharp dose fall-off after the Bragg Peak (BP) becomes a disadvantage if the BP positions of a plan’s constituent pencil beams are shifted, eg.due to anatomical changes. Thus, for fractionated PBS proton therapy, accurate knowledge of the water equivalent path length (WEPL) of the traversed anatomy is critical. In this work we investigate the feasibility of using 2D proton range maps (proton radiography, PR) with the active-scanning gantry at PSI. Methods: We simulated our approach using Monte Carlo methods (MC) to simulate proton beam interactions in patients using clinical imaging data. We selected six head and neck cases having significant anatomical changes detected in per-treatment CTs.PRs (two at 0°/90°) were generated from MC simulations of low-dose pencil beams at 230MeV. Each beam’s residual depth-dose was propagated through the patient geometry (from CT) and detected on exiting the patient anatomy in an ideal depth-resolved detector (eg. range telescope). Firstly, to validate the technique, proton radiographs were compared to the ground truth, which was the WEPL from ray-tracing in the patient CT at the pencil beam location. Secondly, WEPL difference maps (per-treatment – planning imaging timepoints) were then generated to locate the anatomical changes, both in the CT (ground truth) and in the PRs. Binomial classification was performed to evaluate the efficacy of the technique relative to CT. Results: Over the projections simulated over all six patients, 70%, 79% and 95% of the grid points agreed with the ground truth proton range to within ±0.5%, ±1%, and ±3% respectively. The sensitivity, specificity, precision and accuracy were high (mean±1σ, 83±8%, 87±13%, 95±10%, 83±7% respectively). Conclusion: We show that proton-based radiographic images can accurately monitor patient positioning and in vivo range verification, while providing equivalent WEPL information to a CT scan, with the advantage of a much lower imaging dose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-204: Carbon Digitally Reconstructed Radiography (CDRR): A GPU Based Tool for Fast and Versatile Carbonimaging Simulation

Purpose: Research in carbon imaging has been growing over the past years, as a way to increase treatment accuracy and patient positioning in carbon therapy. The purpose of this tool is to allow a fast and flexible way to generate CDRR data without the need to use Monte Carlo (MC) simulations. It can also be used to predict future clinically measured data. Methods: A python interface has been developed, which uses information from CT or 4DCT and thetreatment calibration curve to compute the Water Equivalent Path Length (WEPL) of carbon ions. A GPU based ray tracing algorithm computes the WEPL of each individual carbon traveling through the CT voxels. A multiple peak detection method to estimate high contrast margin positioning has been implemented (described elsewhere). MC simulations have been used to simulate carbons depth dose curves in order to simulate the response of a range detector. Results: The tool allows the upload of CT or 4DCT images. The user has the possibility to selectphase/slice of interested as well as position, angle…). The WEPL is represented as a range detector which can be used to assess range dilution and multiple peak detection effects. The tool also provides knowledge of the minimum energy that should be considered for imaging purposes. The multiple peak detection method has been used in a lung tumor case, showing an accuracy of 1mm in determine the exact interface position. Conclusion: The tool offers an easy and fast way to simulate carbon imaging data. It can be used for educational and for clinical purposes, allowing the user to test beam energies and angles before real acquisition. An analysis add-on is being developed, where the used will have the opportunity to select different reconstruction methods and detector types (range or energy). Fundacao para a Ciencia e a Tecnologia (FCT), PhD Grant number SFRH/BD/85749/2012.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-569: Implementation of a Patient Specific QA Method Using EBT-XD for CyberKnife SRS/SBRT Plans

Purpose: To implement the new EBT-XD Gafchromic films for accurate dosimetric and geometric validation of stereotactic radiosurgery (SRS) and stereotactic body radiation therapy (SBRT) CyberKnife (CK) patient specific QA. Methods: Film calibration was performed using a triplechannel film analysis on an Epson 10000XL scanner. Calibration films were irradiated using a Varian Clinac 21EX flattened beam (0 to 20 Gy), to ensure sufficient dose homogeneity. Films were scanned to a resolution of 0.3 mm, 24 hours post irradiation following a well-defined protocol. A set of 12 QA was performed for several types of CK plans: trigeminal neuralgia, brain metastasis, prostate and lung tumors. A custom made insert for the CK head phantom has been manufactured to yield an accurate measured to calculated dose registration. When the high dose region was large enough, absolute dose was also measured with an ionization chamber. Dose calculation is performed using MultiPlan Ray-tracing algorithm for all cases since the phantom is mostly made from near water-equivalent plastic. Results: Good agreement (<2%) was found between the dose to the chamber and the film, when a chamber measurement was possible The average dose difference and standard deviations between film measurements and TPS calculations were respectively 1.75% and 3%. The geometric accuracy has been estimated to be <1 mm, combining robot positioning uncertainty and film registration to calculated dose. Conclusion: Patient specific QA measurements using EBT-XD films yielded a full 2D dose plane with high spatial resolution and acceptable dose accuracy. This method is particularly promising for trigeminal neuralgia plan QA, where the positioning of the spatial dose distribution is equally or more important than the absolute delivered dose to achieve clinical goals.

60 APPLIED LIFE SCIENCES↗