Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toolkit”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

A software-based x-ray scatter correction method for breast tomosynthesis

Purpose: To develop a software-based scatter correction method for digital breast tomosynthesis (DBT) imaging and investigate its impact on the image quality of tomosynthesis reconstructions of both phantoms and patients. Methods: A Monte Carlo (MC) simulation of x-ray scatter, with geometry matching that of the cranio-caudal (CC) view of a DBT clinical prototype, was developed using the Geant4 toolkit and used to generate maps of the scatter-to-primary ratio (SPR) of a number of homogeneous standard-shaped breasts of varying sizes. Dimension-matched SPR maps were then deformed and registered to DBT acquisition projections, allowing for the estimation of the primary x-ray signal acquired by the imaging system. Noise filtering of the estimated projections was then performed to reduce the impact of the quantum noise of the x-ray scatter. Three dimensional (3D) reconstruction was then performed using the maximum likelihood-expectation maximization (MLEM) method. This process was tested on acquisitions of a heterogeneous 50/50 adipose/glandular tomosynthesis phantom with embedded masses, fibers, and microcalcifications and on acquisitions of patients. The image quality of the reconstructions of the scatter-corrected and uncorrected projections was analyzed by studying the signal-difference-to-noise ratio (SDNR), the integral of the signal in each mass lesion (integrated mass signal, IMS), and the modulation transfer function (MTF). Results: The reconstructions of the scatter-corrected projections demonstrated superior image quality. The SDNR of masses embedded in a 5 cm thick tomosynthesis phantom improved 60%-66%, while the SDNR of the smallest mass in an 8 cm thick phantom improved by 59% (p < 0.01). The IMS of the masses in the 5 cm thick phantom also improved by 15%-29%, while the IMS of the masses in the 8 cm thick phantom improved by 26%-62% (p < 0.01). Some embedded microcalcifications in the tomosynthesis phantoms were visible only in the scatter-corrected reconstructions. The visibility of the findings in two patient images was also improved by the application of the scatter correction algorithm. The MTF of the images did not change after application of the scatter correction algorithm, indicating that spatial resolution was not adversely affected. Conclusions: Our software-based scatter correction algorithm exhibits great potential in improving the image quality of DBT acquisitions of both phantoms and patients. The proposed algorithm does not require a time-consuming MC simulation for each specific case to be corrected, making it applicable in the clinical realm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Dedicated breast CT: radiation dose for circle-plus-line trajectory

Purpose: Dedicated breast CT prototypes used in clinical investigations utilize single circular source trajectory and cone-beam geometry with flat-panel detectors that do not satisfy data-sufficiency conditions and could lead to cone beam artifacts. Hence, this work investigated the glandular dose characteristics of a circle-plus-line trajectory that fulfills data-sufficiency conditions for image reconstruction in dedicated breast CT. Methods: Monte Carlo-based computer simulations were performed using the GEANT4 toolkit and was validated with previously reported normalized glandular dose coefficients for one prototype breast CT system. Upon validation, Monte Carlo simulations were performed to determine the normalized glandular dose coefficients as a function of x-ray source position along the line scan. The source-to-axis of rotation distance and the source-to-detector distance were maintained constant at 65 and 100 cm, respectively, in all simulations. The ratio of the normalized glandular dose coefficient at each source position along the line scan to that for the circular scan, defined as relative normalized glandular dose coefficient (RD{sub g}N), was studied by varying the diameter of the breast at the chest wall, chest-wall to nipple distance, skin thickness, x-ray beam energy, and glandular fraction of the breast. Results: The RD{sub g}N metric when stated as a function of source position along the line scan, relative to the maximum length of line scan needed for data sufficiency, was found to be minimally dependent on breast diameter, chest-wall to nipple distance, skin thickness, glandular fraction, and x-ray photon energy. This observation facilitates easy estimation of the average glandular dose of the line scan. Polynomial fit equations for computing the RD{sub g}N and hence the average glandular dose are provided. Conclusions: For a breast CT system that acquires 300-500 projections over 2{pi} for the circular scan, the addition of a line trajectory with equal source spacing and constant x-ray beam quality (kVp and HVL) and mAs matched to the circular scan, will result in less than 0.18% increase in average glandular dose to the breast per projection along the line scan.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Estimation of organ and effective dose due to Compton backscatter security scans

Purpose: To estimate organ and effective radiation doses due to backscatter security scanners using Monte Carlo simulations and a voxelized phantom set. Methods: Voxelized phantoms of male and female adults and children were used with the GEANT4 toolkit to simulate a backscatter security scan. The backscatter system was modeled based on specifications available in the literature. The simulations modeled a 50 kVp spectrum with 1.0 mm-aluminum-equivalent filtration and a previously measured exposure of approximately 4.6 {mu}R at 30 cm from the source. Photons and secondary interactions were tracked from the source until they reached zero kinetic energy or exited from the simulation's boundaries. The energy deposited in the phantoms' respective organs was tallied and used to calculate total organ dose and total effective dose for frontal, rear, and full scans with subjects located 30 and 75 cm from the source. Results: For a full screen, all phantoms' total effective doses were below the established 0.25 {mu}Sv standard, with an estimated maximum total effective dose of 0.07 {mu}Sv for full screen of a male child. The estimated maximum organ dose due to a full screen was 1.03 {mu}Gy, deposited in the adipose tissue of the male child phantom when located 30 cm from the source. All organ dose estimates had a coefficient of variation of less than 3% for a frontal scan and less than 11% for a rear scan. Conclusions: Backscatter security scanners deposit dose in organs beyond the skin. The effective dose is below recommended standards set by the Health Physics Society (HPS) and the American National Standards Institute (ANSI) assuming the system provides a maximum exposure of approximately 4.6 {mu}R at 30 cm.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Development of a randomized 3D cell model for Monte Carlo microdosimetry simulations

Purpose: The objective of the current work was to develop an algorithm for growing a macroscopic tumor volume from individual randomized quasi-realistic cells. The major physical and chemical components of the cell need to be modeled. It is intended to import the tumor volume into GEANT4 (and potentially other Monte Carlo packages) to simulate ionization events within the cell regions. Methods: A MATLAB Copyright-Sign code was developed to produce a tumor coordinate system consisting of individual ellipsoidal cells randomized in their spatial coordinates, sizes, and rotations. An eigenvalue method using a mathematical equation to represent individual cells was used to detect overlapping cells. GEANT4 code was then developed to import the coordinate system into GEANT4 and populate it with individual cells of varying sizes and composed of the membrane, cytoplasm, reticulum, nucleus, and nucleolus. Each region is composed of chemically realistic materials. Results: The in-house developed MATLAB Copyright-Sign code was able to grow semi-realistic cell distributions ({approx}2 Multiplication-Sign 10{sup 8} cells in 1 cm{sup 3}) in under 36 h. The cell distribution can be used in any number of Monte Carlo particle tracking toolkits including GEANT4, which has been demonstrated in this work. Conclusions: Using the cell distribution and GEANT4, the authors were able to simulate ionization events in the individual cell components resulting from 80 keV gamma radiation (the code is applicable to other particles and a wide range of energies). This virtual microdosimetry tool will allow for a more complete picture of cell damage to be developed.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

Water and tissue equivalence of a new PRESAGE{sup Registered-Sign} formulation for 3D proton beam dosimetry: A Monte Carlo study

Purpose: To evaluate the water and tissue equivalence of a new PRESAGE{sup Registered-Sign} 3D dosimeter for proton therapy. Methods: The GEANT4 software toolkit was used to calculate and compare total dose delivered by a proton beam with mean energy 62 MeV in a PRESAGE{sup Registered-Sign} dosimeter, water, and soft tissue. The dose delivered by primary protons and secondary particles was calculated. Depth-dose profiles and isodose contours of deposited energy were compared for the materials of interest. Results: The proton beam range was found to be Almost-Equal-To 27 mm for PRESAGE{sup Registered-Sign }, 29.9 mm for soft tissue, and 30.5 mm for water. This can be attributed to the lower collisional stopping power of water compared to soft tissue and PRESAGE{sup Registered-Sign }. The difference between total dose delivered in PRESAGE{sup Registered-Sign} and total dose delivered in water or tissue is less than 2% across the entire water/tissue equivalent range of the proton beam. The largest difference between total dose in PRESAGE{sup Registered-Sign} and total dose in water is 1.4%, while for soft tissue it is 1.8%. In both cases, this occurs at the distal end of the beam. Nevertheless, the authors find that PRESAGE{sup Registered-Sign} dosimeter is overall more tissue-equivalent than water-equivalent before the Bragg peak. After the Bragg peak, the differences in the depth doses are found to be due to differences in primary proton energy deposition; PRESAGE{sup Registered-Sign} and soft tissue stop protons more rapidly than water. The dose delivered by secondary electrons in the PRESAGE{sup Registered-Sign} differs by less than 1% from that in soft tissue and water. The contribution of secondary particles to the total dose is less than 4% for electrons and Almost-Equal-To 1% for protons in all the materials of interest. Conclusions: These results demonstrate that the new PRESAGE{sup Registered-Sign} formula may be considered both a tissue- and water-equivalent 3D dosimeter for a 62 MeV proton beam. The results further suggest that tissue-equivalent thickness may provide better dosimetric and geometric accuracy than water-equivalent thickness for 3D dosimetry of this proton beam.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Accuracy of volume measurement using 3D ultrasound and development of CT-3D US image fusion algorithm for prostate cancer radiotherapy

Purpose: To evaluate the accuracy of measuring volumes using three-dimensional ultrasound (3D US), and to verify the feasibility of the replacement of CT-MR fusion images with CT-3D US in radiotherapy treatment planning. Methods: Phantoms, consisting of water, contrast agent, and agarose, were manufactured. The volume was measured using 3D US, CT, and MR devices. A CT-3D US and MR-3D US image fusion software was developed using the Insight Toolkit library in order to acquire three-dimensional fusion images. The quality of the image fusion was evaluated using metric value and fusion images. Results: Volume measurement, using 3D US, shows a 2.8 {+-} 1.5% error, 4.4 {+-} 3.0% error for CT, and 3.1 {+-} 2.0% error for MR. The results imply that volume measurement using the 3D US devices has a similar accuracy level to that of CT and MR. Three-dimensional image fusion of CT-3D US and MR-3D US was successfully performed using phantom images. Moreover, MR-3D US image fusion was performed using human bladder images. Conclusions: 3D US could be used in the volume measurement of human bladders and prostates. CT-3D US image fusion could be used in monitoring the target position in each fraction of external beam radiation therapy. Moreover, the feasibility of replacing the CT-MR image fusion to the CT-3D US in radiotherapy treatment planning was verified.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Digital tomosynthesis (DTS) for verification of target position in early stage lung cancer patients

Purpose: The ability to verify intrafraction tumor position is clinically useful for hypofractionated treatments. Short arc kV digital tomosynthesis (DTS) could facilitate more frequent target verification. The authors used DTS combined with triangulation to determine the mean temporal position of small-volume lung tumor targets treated with stereotactic radiotherapy. DTS registration results were benchmarked against online clinical localization using registration between free-breathing cone-beam computed tomography (CBCT) and the average intensity projection (AvIP) of the planning 4DCT.Methods: In this retrospective study, 76 sets of kV-projection images from online CBCT scans of 13 patients were used to generate DTS image slices (CB-DTS) with nonclinical research software (DTS Toolkit, Varian Medical Systems). Three-dimensional tumor motion was 1.3–4 mm in six patients and 6.1–25.4 mm in seven patients on 4DCT (significant difference in the mean of the groups, P < 0.01). The 4DCT AvIP was used to digitally reconstruct the Reference-DTS. DTS registration and DTS registration combined with triangulation were investigated. Progressive shortening of total DTS arc lengths from 95° to 35° around 0° gantry position was evaluated for different scenarios: DTS registration using the entire arc; DTS registration plus triangulation using two nonoverlapping arcs; and for 55° and 45° total gantry rotation, DTS registration plus triangulation using two overlapping arcs. Finally, DTS registration plus triangulation performed at eight gantry angles, each separated by 45° was evaluated using full fan kV projection data for one patient with an immobile tumor and five patients with mobile tumors.Results: For DTS registration alone, shortening arc length did not influence accuracy in X- and Y-directions, but in Z-direction, mean deviations from online CBCT localization systematically increased for shorter arc length (P < 0.05). For example, using a 95° arc mean DTS-CBCT difference was 0.8 mm (1 SD = 0.6 mm) and for a 35° arc the mean was 2.4 mm (1 SD = 1.7 mm). DTS plus triangulation using nonoverlapping-arcs increased accuracy in Z-direction for tested arc lengths ≤55° (P < 0.01). Overlapping arcs increased accuracy in Y-direction for tumors with motion >4 mm (P < 0.02) but increased Z-direction accuracy was only observed with 55° total gantry rotation. The 95th percentile deviations with this overlapping technique in X-, Y-, and Z-directions were 1.3, 2.0, and 2.5 mm, respectively. For the five patients with mobile tumors where DTS + triangulation was performed with 45° intervals, the pooled deviation from online CBCT correction showed, for X-, Y-, and Z-directions, mean of 1.1 mm, standard deviations (SD) of 0.9, 1.0, and 0.9 mm, respectively. The mean + 2 SD was <3 mm for each direction.Conclusions: Short-arc DTS verification of time averaged lung tumor position is feasible using free-breathing kV projection data and the AvIP of the 4DCT as a reference. Observed differences between DTS and online CBCT registration with AvIP were ≤3 mm (mean + 2 SD), however, the increased temporal resolution of DTS + triangulation also identified short period deviations from the average target position on the CBCT. Short-arc DTS appears promising for intrafraction tumor position monitoring during stereotactic lung radiotherapy delivered with a rotational technique.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Investigation of realistic PET simulations incorporating tumor patient's specificity using anthropomorphic models: Creation of an oncology database

Purpose: The GATE Monte Carlo simulation toolkit is used for the implementation of realistic PET simulations incorporating tumor heterogeneous activity distributions. The reconstructed patient images include noise from the acquisition process, imaging system's performance restrictions and have limited spatial resolution. For those reasons, the measured intensity cannot be simply introduced in GATE simulations, to reproduce clinical data. Investigation of the heterogeneity distribution within tumors applying partial volume correction (PVC) algorithms was assessed. The purpose of the present study was to create a simulated oncology database based on clinical data with realistic intratumor uptake heterogeneity properties.Methods: PET/CT data of seven oncology patients were used in order to create a realistic tumor database investigating the heterogeneity activity distribution of the simulated tumors. The anthropomorphic models (NURBS based cardiac torso and Zubal phantoms) were adapted to the CT data of each patient, and the activity distribution was extracted from the respective PET data. The patient-specific models were simulated with the Monte Carlo Geant4 application for tomography emission (GATE) in three different levels for each case: (a) using homogeneous activity within the tumor, (b) using heterogeneous activity distribution in every voxel within the tumor as it was extracted from the PET image, and (c) using heterogeneous activity distribution corresponding to the clinical image following PVC. The three different types of simulated data in each case were reconstructed with two iterations and filtered with a 3D Gaussian postfilter, in order to simulate the intratumor heterogeneous uptake. Heterogeneity in all generated images was quantified using textural feature derived parameters in 3D according to the ground truth of the simulation, and compared to clinical measurements. Finally, profiles were plotted in central slices of the tumors, across lines with heterogeneous activity distribution for visual assessment.Results: The accuracy of the simulated database was assessed against the original clinical images. The PVC simulated images matched the clinical ones best. Local, regional, and global features extracted from the PVC simulated images were closest to the clinical measurements, with the exception of the size zone variability and the mean intensity values, where heterogeneous tumors showed better reproducibility. The profiles on PVC simulated tumors after postfiltering seemed to represent the more realistic heterogeneous regions with respect to the clinical reference.Conclusions: In this study, the authors investigated the input activity map heterogeneity in the GATE simulations of tumors with heterogeneous activity distribution. The most realistic heterogeneous tumors were obtained by inserting PVC activity distributions from the clinical image into the activity map of the simulation. Partial volume effect (PVE) can play a crucial role in the quantification of heterogeneity within tumors and have an important impact on applications such as patient follow-up during treatment and assessment of tumor response to therapy. The development of such a database incorporating patient anatomical and functional variability can be used to evaluate new image processing or analysis algorithms, while providing control of the ground truth, which is not available when dealing with clinical datasets. The database includes all images used and generated in this study, as well as the sinograms and the attenuation phantoms for further investigation. It is freely available to the interested reader of the journal at http://www.med.upatras.gr/oncobase/.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Toward an organ based dose prescription method for the improved accuracy of murine dose in orthovoltage x-ray irradiators

Purpose: Accurate dosimetry is essential when irradiating mice to ensure that functional and molecular endpoints are well understood for the radiation dose delivered. Conventional methods of prescribing dose in mice involve the use of a single dose rate measurement and assume a uniform average dose throughout all organs of the entire mouse. Here, the authors report the individual average organ dose values for the irradiation of a 12, 23, and 33 g mouse on a 320 kVp x-ray irradiator and calculate the resulting error from using conventional dose prescription methods. Methods: Organ doses were simulated in the Geant4 application for tomographic emission toolkit using the MOBY mouse whole-body phantom. Dosimetry was performed for three beams utilizing filters A (1.65 mm Al), B (2.0 mm Al), and C (0.1 mm Cu + 2.5 mm Al), respectively. In addition, simulated x-ray spectra were validated with physical half-value layer measurements. Results: Average doses in soft-tissue organs were found to vary by as much as 23%–32% depending on the filter. Compared to filters A and B, filter C provided the hardest beam and had the lowest variation in soft-tissue average organ doses across all mouse sizes, with a difference of 23% for the median mouse size of 23 g. Conclusions: This work suggests a new dose prescription method in small animal dosimetry: it presents a departure from the conventional approach of assigninga single dose value for irradiation of mice to a more comprehensive approach of characterizing individual organ doses to minimize the error and uncertainty. In human radiation therapy, clinical treatment planning establishes the target dose as well as the dose distribution, however, this has generally not been done in small animal research. These results suggest that organ dose errors will be minimized by calibrating the dose rates for all filters, and using different dose rates for different organs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Influence of polarization and a source model for dose calculation in MRT

Purpose: Microbeam Radiation Therapy (MRT), an alternative preclinical treatment strategy using spatially modulated synchrotron radiation on a micrometer scale, has the great potential to cure malignant tumors (e.g., brain tumors) while having low side effects on normal tissue. Dose measurement and calculation in MRT is challenging because of the spatial accuracy required and the arising high dose differences. Dose calculation with Monte Carlo simulations is time consuming and their accuracy is still a matter of debate. In particular, the influence of photon polarization has been discussed in the literature. Moreover, it is controversial whether a complete knowledge of phase space trajectories, i.e., the simulation of the machine from the wiggler to the collimator, is necessary in order to accurately calculate the dose. Methods: With Monte Carlo simulations in the Geant4 toolkit, the authors investigate the influence of polarization on the dose distribution and the therapeutically important peak to valley dose ratios (PVDRs). Furthermore, the authors analyze in detail phase space information provided byMartínez-Rovira et al. [“Development and commissioning of a Monte Carlo photon model for the forthcoming clinical trials in microbeam radiation therapy,” Med. Phys. 39(1), 119–131 (2012)] and examine its influence on peak and valley doses. A simple source model is developed using parallel beams and its applicability is shown in a semiadjoint Monte Carlo simulation. Results are compared to measurements and previously published data. Results: Polarization has a significant influence on the scattered dose outside the microbeam field. In the radiation field, however, dose and PVDRs deduced from calculations without polarization and with polarization differ by less than 3%. The authors show that the key consequences from the phase space information for dose calculations are inhomogeneous primary photon flux, partial absorption due to inclined beam incidence outside the field center, increased beam width and center to center distance due to the beam propagation from the collimator to the phantom surface and imperfect absorption in the absorber material of the Multislit Collimator. These corrections have an effect of approximately 10% on the valley dose and suffice to describe doses in MRT within the measurement uncertainties of currently available dosimetry techniques. Conclusions: The source for the first clinical pet trials in MRT is characterized with respect to its phase space and the photon polarization. The results suggest the use of a presented simplified phase space model in dose calculations and hence pave the way for alternative and fast dose calculation algorithms. They also show that the polarization is of minor importance for the clinical important peak and valley doses inside the microbeam field.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

A review of the use and potential of the GATE Monte Carlo simulation code for radiation therapy and dosimetry applications

In this paper, the authors' review the applicability of the open-source GATE Monte Carlo simulation platform based on the GEANT4 toolkit for radiation therapy and dosimetry applications. The many applications of GATE for state-of-the-art radiotherapy simulations are described including external beam radiotherapy, brachytherapy, intraoperative radiotherapy, hadrontherapy, molecular radiotherapy, and in vivo dose monitoring. Investigations that have been performed using GEANT4 only are also mentioned to illustrate the potential of GATE. The very practical feature of GATE making it easy to model both a treatment and an imaging acquisition within the same frameworkis emphasized. The computational times associated with several applications are provided to illustrate the practical feasibility of the simulations using current computing facilities.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Concept and simulation study of a novel localization method for robotic endoscopic capsules using multiple positron emission markers

Purpose: Over the last decade, wireless capsule endoscope has been the tool of choice for noninvasive inspection of the gastrointestinal tract, especially in the small intestine. However, the latest clinical products have not been equipped with a sufficiently accurate localization system which makes it difficult to determine the location of intestinal abnormalities, and to apply follow-up interventions such as biopsy or drug delivery. In this paper, the authors present a novel localization method based on tracking three positron emission markers embedded inside an endoscopic capsule. Methods: Three spherical {sup 22}Na markers with diameters of less than 1 mm are embedded in the cover of the capsule. Gamma ray detectors are arranged around a patient body to detect coincidence gamma rays emitted from the three markers. The position of each marker can then be estimated using the collected data by the authors’ tracking algorithm which consists of four consecutive steps: a method to remove corrupted data, an initialization method, a clustering method based on the Fuzzy C-means clustering algorithm, and a failure prediction method. Results: The tracking algorithm has been implemented inMATLAB utilizing simulation data generated from the Geant4 Application for Emission Tomography toolkit. The results show that this localization method can achieve real-time tracking with an average position error of less than 0.4 mm and an average orientation error of less than 2°. Conclusions: The authors conclude that this study has proven the feasibility and potential of the proposed technique in effectively determining the position and orientation of a robotic endoscopic capsule.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-J-175: Proton Dose Calculation On Scatter-Corrected CBCT Image: Feasibility Study for Adaptive Proton Therapy

Purpose: To demonstrate feasibility of proton dose calculation on scattercorrected CBCT images for the purpose of adaptive proton therapy. Methods: Two CBCT image sets were acquired from a prostate cancer patient and a thorax phantom using an on-board imaging system of an Elekta infinity linear accelerator. 2-D scatter maps were estimated using a previously introduced CT-based technique, and were subtracted from each raw projection image. A CBCT image set was then reconstructed with an open source reconstruction toolkit (RTK). Conversion from the CBCT number to HU was performed by soft tissue-based shifting with reference to the plan CT. Passively scattered proton plans were simulated on the plan CT and corrected/uncorrected CBCT images using the XiO treatment planning system. For quantitative evaluation, water equivalent path length (WEPL) was compared in those treatment plans. Results: The scatter correction method significantly improved image quality and HU accuracy in the prostate case where large scatter artifacts were obvious. However, the correction technique showed limited effects on the thorax case that was associated with fewer scatter artifacts. Mean absolute WEPL errors from the plans with the uncorrected and corrected images were 1.3 mm and 5.1 mm in the thorax case and 13.5 mm and 3.1 mm in the prostate case. The prostate plan dose distribution of the corrected image demonstrated better agreement with the reference one than that of the uncorrected image. Conclusion: A priori CT-based CBCT scatter correction can reduce the proton dose calculation error when large scatter artifacts are involved. If scatter artifacts are low, an uncorrected CBCT image is also promising for proton dose calculation when it is calibrated with the soft-tissue based shifting.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-247: A Simulation of X-Ray Emission with Gold Nanoparticle Irradiated by Energetic Proton Beam

Purpose: To investigate the proton induced X-ray emissions in gold-water mixture materials. Methods: In this study a Monte Carlo simulation was created using the GEANT4 toolkit (version 4.9.6). The geometry in this setup includes a 2 cm × 2 cm × 2 cm target, a scoring sphere (radius = 10 cm) and a 65 MeV planar proton source (2 cm × 2 cm). Four concentrations of a gold-water solution were irradiated with 5×10{sup 5} incident protons at a distance of 0.5 cm perpendicular to the surface of the target. The solutions of gold-water mixture had 10%, 5%, 1% and 0.5% of gold by mass, respectively. The number of photon emitting for the target was counted in the scoring sphere for the energy range of 0-86.0 keV in 0.1 keV bins. For this study the reference physics list PhysListEmStandard was used together with the x-ray fluorescence, Auger electron and PIXE (particle induced xray emission) options enabled. The range cuts for photons and electrons were set at 0.5 mm and 1.0 mm, respectively. Results: In the energy spectra of emitting X-ray fluorescence, peaks from gold K shell characteristic x-rays (68.8 and 66.9 keV) were observed. The number of counts under the peaks of Ka1 and Ka2 was found to increase with the increasing of the gold concentrations in the mixture materials. The X-ray yields (for both Ka1 and Ka2) when fitted with least-square method as a function of gold concentration demonstrate a linear dependency with R{sup 2} > 0.96. The Ka1yield per incident proton was found to be 0.0016 for 10% gold-water mixture solutions. Conclusion: This preliminary study with PIXE technique with gold nanoparticle has demonstrated potentials for its utilization in the development of range and dose verification methodology that is currently of great interest in the field of proton radiation therapy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-41: Analysis of GI Dose Variability Due to Intrafraction Setup Variance

Purpose: Proton SBRT (stereotactic body radiation therapy) can be an effective modality for treatment of gastrointestinal tumors, but limited in practice due to sensitivity with respect to variation in the RPL (radiological path length). Small, intrafractional shifts in patient anatomy can lead to significant changes in the dose distribution. This study describes a tool designed to visualize uncertainties in radiological depth in patient CT's and aid in treatment plan design. Methods: This project utilizes the Shadie toolkit, a GPU-based framework that allows for real-time interactive calculations for volume visualization. Current SBRT simulation practice consists of a serial CT acquisition for the assessment of inter- and intra-fractional motion utilizing patient specific immobilization systems. Shadie was used to visualize potential uncertainties, including RPL variance and changes in gastric content. Input for this procedure consisted of two patient CT sets, contours of the desired organ, and a pre-calculated dose. In this study, we performed rigid registrations between sets of 4DCT's obtained from a patient with varying setup conditions. Custom visualizations are written by the user in Shadie, permitting one to create color-coded displays derived from a calculation along each ray. Results: Serial CT data acquired on subsequent days was analyzed for variation in RPB and gastric content. Specific shaders were created to visualize clinically relevant features, including RPL (radiological path length) integrated up to organs of interest. Using pre-calculated dose distributions and utilizing segmentation masks as additional input allowed us to further refine the display output from Shadie and create tools suitable for clinical usage. Conclusion: We have demonstrated a method to visualize potential uncertainty for intrafractional proton radiotherapy. We believe this software could prove a useful tool to guide those looking to design treatment plans least insensitive to motion for patients undergoing proton SBRT in the GI tract.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-T-78: Comparison of Dose-Averaged Linear Energy Transfer Calculation Methods Used in Monte Carlo Simulations of Clinical Proton Beams

Purpose: To evaluate the differences in dose-averaged linear energy transfer (LETd) maps calculated in water by means of different strategies found in the literature in proton therapy Monte Carlo simulations and to compare their values with dose-mean lineal energy microdosimetry calculations. Methods: The Geant4 toolkit (version 9.6.2) was used. Dose and LETd maps in water were scored for primary protons with cylindrical voxels defined around the beam axis. Three LETd calculation methods were implemented. First, the LETd values were computed by calculating the unrestricted linear energy transfer (LET) associated to each single step weighted by the energy deposition (including delta-rays) along the step. Second, the LETd was obtained for each voxel by computing the LET along all the steps simulated for each proton track within the voxel, weighted by the energy deposition of those steps. Third, the LETd was scored as the quotient between the second momentum of the LET distribution, calculated per proton track, over the first momentum. These calculations were made with various voxel thicknesses (0.2 – 2.0 mm) for a 160 MeV proton beamlet and spread-out Bragg Peaks (SOBP). The dose-mean lineal energy was calculated in a uniformly-irradiated water sphere, 0.005 mm radius. Results: The value of the LETd changed systematically with the voxel thickness due to delta-ray emission and the enlargement of the LET distribution spread, especially at shallow depths. Differences of up to a factor 1.8 were found at the depth of maximum dose, leading to similar differences at the central and distal depths of the SOBPs. The third LETd calculation method gave better agreement with microdosimetry calculations around the Bragg Peak. Conclusion: Significant differences were found between LETd map Monte Carlo calculations due to both the calculation strategy and the voxel thickness used. This could have a significant impact in radiobiologically-optimized proton therapy treatments.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

SU-E-T-186: Feasibility Study of Glass Cherenkov Detector for Prompt Gamma Detection in Proton Therapy

Purpose: To simulate a Cherenkov glass detector system utilizing prompt gamma (PG) technique to quantify range uncertainties in proton radiation therapy. Methods: A simulation of high energy photons typically produced in proton interactions with materials incident onto a block of Cherenkov glass was performed with the Geant4 toolkit. The standard electromagnetic package was used along with several decay modules (G4Decay, G4DecayPhysics, and G4RadioactiveDecayPhysics) and the optical photon components (G4OpticalPhysics). Our setup included a pencil beam consisting of a hundred thousand 6 MeV photons (approximately the deexcitation energy released from 16O) incident onto a 2.5 ⊗ 2.5 ⊗ 1.5 cm3 of a Cherenkov glass (7.2 g of In2O3 + 90 g cladding, density of 2.82 g/cm3, Zeff = 33.7, index of refraction 1.56). The energy deposited from incident 6 MeV photons as well as secondary electrons and resulting optical photons were recorded. Results: The energy deposited by 6 MeV photons in glass material showed several peaks that included the photoelectric, the single and double escape peaks. About 11% of incident photons interacted with glass material to deposit energy. Most of the photons collected were in the region of double escape peak (approximately 4.98 MeV). The secondary electron spectrum produced from incident photons showed a high energy peak located near 6 MeV and a sharp peak located ∼120 keV with a continuous distribution between these two points. The resulting Cherenkov photons produced showed a continuous energy distribution between 2 and 5 eV with a slight increase in yield beginning about 3 eV. The amount of Cherenkov photons produced per interacting incident 6 MeV photon was ∼240.7. Conclusion: This study suggests the viability of utilizing the Cherenkov glass material as a possible prompt gamma photon detection device. Future work will include optimization of the detector system to maximize photon detection efficiency.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-T-241: Monte Carlo Simulation Study About the Prediction of Proton-Induced DNA Strand Breakage On the Double Helix Structure

Purpose: In particle therapy and radiobiology, the investigation of mechanisms leading to the death of target cancer cells induced by ionising radiation is an active field of research. Recently, several studies based on Monte Carlo simulation codes have been initiated in order to simulate physical interactions of ionising particles at cellular scale and in DNA. Geant4-DNA is the one of them; it is an extension of the general purpose Geant4 Monte Carlo simulation toolkit for the simulation of physical interactions at sub-micrometre scale. In this study, we present Geant4-DNA Monte Carlo simulations for the prediction of DNA strand breakage using a geometrical modelling of DNA structure. Methods: For the simulation of DNA strand breakage, we developed a specific DNA geometrical structure. This structure consists of DNA components, such as the deoxynucleotide pairs, the DNA double helix, the nucleosomes and the chromatin fibre. Each component is made of water because the cross sections models currently available in Geant4-DNA for protons apply to liquid water only. Also, at the macroscopic-scale, protons were generated with various energies available for proton therapy at the National Cancer Center, obtained using validated proton beam simulations developed in previous studies. These multi-scale simulations were combined for the validation of Geant4-DNA in radiobiology. Results: In the double helix structure, the deposited energy in a strand allowed to determine direct DNA damage from physical interaction. In other words, the amount of dose and frequency of damage in microscopic geometries was related to direct radiobiological effect. Conclusion: In this report, we calculated the frequency of DNA strand breakage using Geant4- DNA physics processes for liquid water. This study is now on-going in order to develop geometries which use realistic DNA material, instead of liquid water. This will be tested as soon as cross sections for DNA material become available in Geant4-DNA.

60 APPLIED LIFE SCIENCES↗