Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Radiological Data”

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 361 records · Page 20

Does administering iodine in radiological procedures increase patient doses?

Purpose: The authors investigated the changes in the pattern of energy deposition in tissue equivalent phantoms following the introduction of iodinated contrast media. Methods: The phantom consisted of a small “contrast sphere,” filled with water or iodinated contrast, located at the center of a 28 cm diameter water sphere. Monte Carlo simulations were performed using MCNP5 codes, validated by simulating irradiations with analytical solutions. Monoenergetic x-rays ranging from 35 to 150 keV were used to simulate exposures to spheres containing contrast agent with iodine concentrations ranging from 1 to 100 mg/ml. Relative values of energy imparted to the contrast sphere, as well as to the whole phantom, were calculated. Changes in patterns of energy deposition around the contrast sphere were also investigated. Results: Small contrast spheres can increase local absorbed dose by a factor of 13, but the corresponding increase in total energy absorbed was negligible (<1%). The highest localized dose increases were found to occur at incident photon energies of about 60 keV. For a concentration of about 10 mg/ml, typical of clinical practice, localized absorbed doses were generally increased by about a factor of two. At this concentration of 10 mg/ml, the maximum increase in total energy deposition in the phantom was only 6%. These simulations demonstrated that increases in contrast sphere doses were offset by corresponding dose reductions at distal and posterior locations. Conclusions: Adding iodine can result in values of localized absorbed dose increasing by more than an order of magnitude, but the total energy deposition is generally very modest (i.e., <10%). Their data show that adding iodine primarily changes the pattern of energy deposition in the irradiated region, rather than increasing patient doses per se.

60 APPLIED LIFE SCIENCES↗

Determining the effects of microsphere and surrounding material composition on {sup 90}Y dose kernels using egsnrc and mcnp5

Purpose: Recent advances in the imaging of {sup 90}Y using positron emission tomography (PET) and improved uncertainty in the branching ratio for the internal pair production component of {sup 90}Y decay allow for a more accurate determination of the activity distribution of {sup 90}Y microspheres within a patient. This improved activity distribution can be convolved with the dose kernel of {sup 90}Y to calculate the dose distribution within a patient. This work investigates the effects of microsphere and surrounding material composition on {sup 90}Y dose kernels using egsnrc and mcnp5 and compares the results of these two transport codes. Methods: Monte Carlo simulations were performed with egsnrc and mcnp5 to calculate the dose rate at multiple radial distances around various {sup 90}Y sources. Point source simulations were completed with mcnp5 to determine the optimal electron transport settings for this work. After determining the optimal settings, point source simulations were completed using egsnrc (user code edknrc) and mcnp5 in water and liver [as defined by the International Commission on Radiation Units and Measurements (ICRU) Report 44]. The results were compared to ICRU Report 72 reference data. Point source simulations were also completed in water with a density of 1.06 g{center_dot}cm{sup -3} to evaluate the effect of the density of the surrounding material. Glass and resin microsphere simulations were performed with average and maximum diameter and density values (based on values given in the literature) in water and in liver. The results were compared to point source simulation results using the same transport code and in the same surrounding material. All simulations had statistical uncertainties less than 1%. Results: The optimal transport settings in mcnp5 for this work included using the energy-and step-specific algorithm (DBCN 17J 2) and ESTEP set to 10. These settings were used for all subsequent simulations with mcnp5. The point source simulations in water for both egsnrc and mcnp5 were found to agree within 2% of the ICRU 72 reference data over the investigated range. Point source simulations in liver had large differences relative to ICRU 72, approaching -60% near the maximum range of {sup 90}Y. These differences are mostly attributed to the difference in density between water (1.0 g{center_dot}cm{sup -3}) and liver (1.06 g{center_dot}cm{sup -3}). Glass and resin microsphere simulations showed a slight decrease in the dose rate near the maximum range of {sup 90}Y relative to the point source simulations. The largest relative differences were approximately -4.2% and -2.8% for the glass and resin microspheres, respectively. Agreement between the egsnrc and mcnp5 simulations results was generally good. Conclusions: The presence of the microsphere material causes slight differences in the {sup 90}Y dose kernel compared to those calculated with point sources. Large differences were seen between simulations in water and those in liver. For the most accurate calculation of the dose distribution, the density of the patient's liver should be accounted for in the calculation of the dose kernel. Lastly, due to the need to determine the optimal transport settings with mcnp5, electron transport with this code should be used with caution.

61 RADIATION PROTECTION AND DOSIMETRY↗

Predicting diagnostic error in radiology via eye-tracking and image analytics: Preliminary investigation in mammography

Purpose: The primary aim of the present study was to test the feasibility of predicting diagnostic errors in mammography by merging radiologists’ gaze behavior and image characteristics. A secondary aim was to investigate group-based and personalized predictive models for radiologists of variable experience levels.Methods: The study was performed for the clinical task of assessing the likelihood of malignancy of mammographic masses. Eye-tracking data and diagnostic decisions for 40 cases were acquired from four Radiology residents and two breast imaging experts as part of an IRB-approved pilot study. Gaze behavior features were extracted from the eye-tracking data. Computer-generated and BIRADS images features were extracted from the images. Finally, machine learning algorithms were used to merge gaze and image features for predicting human error. Feature selection was thoroughly explored to determine the relative contribution of the various features. Group-based and personalized user modeling was also investigated.Results: Machine learning can be used to predict diagnostic error by merging gaze behavior characteristics from the radiologist and textural characteristics from the image under review. Leveraging data collected from multiple readers produced a reasonable group model [area under the ROC curve (AUC) = 0.792 ± 0.030]. Personalized user modeling was far more accurate for the more experienced readers (AUC = 0.837 ± 0.029) than for the less experienced ones (AUC = 0.667 ± 0.099). The best performing group-based and personalized predictive models involved combinations of both gaze and image features.Conclusions: Diagnostic errors in mammography can be predicted to a good extent by leveraging the radiologists’ gaze behavior and image content.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Development of a prototype gantry system for preclinical x-ray phase-contrast computed tomography

Purpose: To explore the potential of grating-based x-ray phase-contrast imaging for clinical applications, a first compact gantry system was developed. It is designed such that it can be implemented into an in-vivo small-animal phase-contrast computed tomography (PC-CT) scanner. The purpose of the present study is to assess the accuracy and quantitativeness of the described gantry in both absorption and phase-contrast. Methods: A phantom, containing six chemically well-defined liquids, was constructed. A tomography scan with cone-beam reconstruction of this phantom was performed yielding the spatial distribution of the linear attenuation coefficient {mu} and decrement {delta} of the complex refractive index. Theoretical values of {mu} and {delta} were calculated for each liquid from tabulated data and compared with the experimentally measured values. Additionally, a color-fused image representation is proposed to display the complementary absorption and phase-contrast information in a single image. Results: Experimental and calculated data of the phantom agree well confirming the quantitativeness and accuracy of the reconstructed spatial distributions of {mu} and {delta}. The proposed color-fused image representation, which combines the complementary absorption and phase information, considerably helps in distinguishing the individual substances. Conclusions: The concept of grating-based phase-contrast computed tomography (CT) can be implemented into a compact, cone-beam geometry gantry setup. The authors believe that this work represents an important milestone in translating phase-contrast x-ray imaging from previous proof-of-principle experiments to first preclinical biomedical imaging applications on small-animal models.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Assessment of interpatient heterogeneity in tumor radiosensitivity for nonsmall cell lung cancer using tumor-volume variation data

Purpose: In our previous work, the authors showed that a distribution of cell surviving fractionsS{sub 2} in a heterogeneous group of patients could be derived from tumor-volume variation curves during radiotherapy for head and neck cancer. In this research study, the authors show that this algorithm can be applied to other tumors, specifically in nonsmall cell lung cancer. This new application includes larger patient volumes and includes comparison of data sets obtained at independent institutions. Methods: Our analysis was based on two data sets of tumor-volume variation curves for heterogeneous groups of 17 patients treated for nonsmall cell lung cancer with conventional dose fractionation. The data sets were obtained previously at two independent institutions by using megavoltage computed tomography. Statistical distributions of cell surviving fractionsS{sub 2} and clearance half-lives of lethally damaged cells T{sub 1/2} have been reconstructed in each patient group by using a version of the two-level cell population model of tumor response and a simulated annealing algorithm. The reconstructed statistical distributions of the cell surviving fractions have been compared to the distributions measured using predictive assays in vitro. Results: Nonsmall cell lung cancer presents certain difficulties for modeling surviving fractions using tumor-volume variation curves because of relatively large fractional hypoxic volume, low gradient of tumor-volume response, and possible uncertainties due to breathing motion. Despite these difficulties, cell surviving fractionsS{sub 2} for nonsmall cell lung cancer derived from tumor-volume variation measured at different institutions have similar probability density functions (PDFs) with mean values of 0.30 and 0.43 and standard deviations of 0.13 and 0.18, respectively. The PDFs for cell surviving fractions S{sub 2} reconstructed from tumor volume variation agree with the PDF measured in vitro. Conclusions: The data obtained in this work, when taken together with the data obtained previously for head and neck cancer, suggests that the cell surviving fractionsS{sub 2} can be reconstructed from the tumor volume variation curves measured during radiotherapy with conventional fractionation. The proposed method can be used for treatment evaluation and adaptation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

3D modeling of effects of increased oxygenation and activity concentration in tumors treated with radionuclides and antiangiogenic drugs

Purpose: Formation of new blood vessels (angiogenesis) in response to hypoxia is a fundamental event in the process of tumor growth and metastatic dissemination. However, abnormalities in tumor neovasculature often induce increased interstitial pressure (IP) and further reduce oxygenation (pO{sub 2}) of tumor cells. In radiotherapy, well-oxygenated tumors favor treatment. Antiangiogenic drugs may lower IP in the tumor, improving perfusion, pO{sub 2} and drug uptake, by reducing the number of malfunctioning vessels in the tissue. This study aims to create a model for quantifying the effects of altered pO{sub 2}-distribution due to antiangiogenic treatment in combination with radionuclide therapy. Methods: Based on experimental data, describing the effects of antiangiogenic agents on oxygenation of GlioblastomaMultiforme (GBM), a single cell based 3D model, including 10{sup 10} tumor cells, was developed, showing how radionuclide therapy response improves as tumor oxygenation approaches normal tissue levels. The nuclides studied were {sup 90}Y, {sup 131}I, {sup 177}Lu, and {sup 211}At. The absorbed dose levels required for a tumor control probability (TCP) of 0.990 are compared for three different log-normal pO{sub 2}-distributions: {mu}{sub 1} = 2.483, {sigma}{sub 1} = 0.711; {mu}{sub 2} = 2.946, {sigma}{sub 2} = 0.689; {mu}{sub 3} = 3.689, and {sigma}{sub 3} = 0.330. The normal tissue absorbed doses will, in turn, depend on this. These distributions were chosen to represent the expected oxygen levels in an untreated hypoxic tumor, a hypoxic tumor treated with an anti-VEGF agent, and in normal, fully-oxygenated tissue, respectively. The former two are fitted to experimental data. The geometric oxygen distributions are simulated using two different patterns: one Monte Carlo based and one radially increasing, while keeping the log-normal volumetric distributions intact. Oxygen and activity are distributed, according to the same pattern. Results: As tumor pO{sub 2} approaches normal tissue levels, the therapeutic effect is improved so that the normal tissue absorbed doses can be decreased by more than 95%, while retaining TCP, in the most favorable scenario and by up to about 80% with oxygen levels previously achieved in vivo, when the least favourable oxygenation case is used as starting point. The major difference occurs in poorly oxygenated cells. This is also where the pO{sub 2}-dependence of the oxygen enhancement ratio is maximal. Conclusions: Improved tumor oxygenation together with increased radionuclide uptake show great potential for optimising treatment strategies, leaving room for successive treatments, or lowering absorbed dose to normal tissues, due to increased tumor response. Further studies of the concomitant use of antiangiogenic drugs and radionuclide therapy therefore appear merited.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Molecular Breast Imaging Using Emission Tomosynthesis

Purpose: Tour objective is to design a novel SPECT system for molecular breast imaging (MBI) and evaluate its performance. The limited angle SPECT system, or emission tomosynthesis, is designed to achieve 3D images of the breast with high spatial resolution/sensitivity. The system uses a simplified detector motion and is conducive to on-board biopsy and mult-modal imaging with mammography. Methods: The novel feature of the proposed gamma camera is a variable-angle, slant-hole (VASH) collimator, which is well suited for limited angle SPECT of a mildly compressed breast. The collimator holes change slant angle while the camera surface remains flush against the compression paddle. This allows the camera to vary the angular view ({+-}30{degrees}, {+-}45{degrees}) for tomographic imaging while keeping the camera close to the object for high spatial resolution and/or sensitivity. Theoretical analysis and Monte Carlo simulations were performed assuming a point source and isolated breast phantom. Spatial resolution, sensitivity, contrast and SNR were measured. Results were compared to single-view, planar images and conventional SPECT. For both conventional SPECT and VASH, data were reconstructed using iterative algorithms. Finally, a proof-of-concept VASH collimator was constructed for experimental evaluation. Results: Measured spatial resolution/sensitivity with VASH showed good agreement with theory including depth-of-interaction (DOI) effects. The DOI effect diminished the depth resolution by approximately 2 mm. Increasing the slant angle range from {+-}30{degrees} to {+-}45{degrees} resulted in an approximately 1 mm improvement in the depth resolution. In the breast phantom images, VASH showed improved contrast and SNR over conventional SPECT and improved contrast over planar scintimmammography. Reconstructed images from the proof-of-concept VASH collimator demonstrated reasonable depth resolution capabilities using limited angle projection data. Conclusion: We conclude that this limited angle SPECT approach using a VASH collimator can achieve superior 3D image quality for MBI in a design that has attractive characteristics for clinical imaging.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Effects of sparse sampling schemes on image quality in low-dose CT

Purpose: Various scanning methods and image reconstruction algorithms are actively investigated for low-dose computed tomography (CT) that can potentially reduce a health-risk related to radiation dose. Particularly, compressive-sensing (CS) based algorithms have been successfully developed for reconstructing images from sparsely sampled data. Although these algorithms have shown promises in low-dose CT, it has not been studied how sparse sampling schemes affect image quality in CS-based image reconstruction. In this work, the authors present several sparse-sampling schemes for low-dose CT, quantitatively analyze their data property, and compare effects of the sampling schemes on the image quality.Methods: Data properties of several sampling schemes are analyzed with respect to the CS-based image reconstruction using two measures: sampling density and data incoherence. The authors present five different sparse sampling schemes, and simulated those schemes to achieve a targeted dose reduction. Dose reduction factors of about 75% and 87.5%, compared to a conventional scan, were tested. A fully sampled circular cone-beam CT data set was used as a reference, and sparse sampling has been realized numerically based on the CBCT data.Results: It is found that both sampling density and data incoherence affect the image quality in the CS-based reconstruction. Among the sampling schemes the authors investigated, the sparse-view, many-view undersampling (MVUS)-fine, and MVUS-moving cases have shown promising results. These sampling schemes produced images with similar image quality compared to the reference image and their structure similarity index values were higher than 0.92 in the mouse head scan with 75% dose reduction.Conclusions: The authors found that in CS-based image reconstructions both sampling density and data incoherence affect the image quality, and suggest that a sampling scheme should be devised and optimized by use of these indicators. With this strategic approach, one can acquire optimally sampled sparse data so that the CS-based algorithms can best perform in terms of image quality.

61 RADIATION PROTECTION AND DOSIMETRY↗

Specific recommendations for accurate and direct use of PET-CT in PET guided radiotherapy for head and neck sites

Purpose: To provide specific experience-based guidance and recommendations for centers wishing to develop, validate, and implement an accurate and efficient process for directly using positron emission tomography-computed tomography (PET-CT) for the radiotherapy planning of head and neck cancer patients. Methods: A PET-CT system was modified with hard-top couch, external lasers and radiotherapy immobilization and indexing devices and was subject to a commissioning and quality assurance program. PET-CT imaging protocols were developed specifically for radiotherapy planning and the image quality and pathway tested using phantoms and five patients recruited into an in-house study. Security and accuracy of data transfer was tested throughout the whole data pathway. The patient pathway was fully established and tested ready for implementation in a PET-guided dose-escalation trial for head and neck cancer patients. Results: Couch deflection was greater than for departmental CT simulator machines. An area of high attenuation in the couch generated image artifacts and adjustments were made accordingly. Using newly developed protocols CT image quality was suitable to maintain delineation and treatment accuracy. Upon transfer of data to the treatment planning system a half pixel offset between PET and CT was observed and corrected. By taking this into account, PET to CT alignment accuracy was maintained below 1 mm in all systems in the data pathway. Transfer of structures delineated in the PET fusion software to the radiotherapy treatment planning system was validated. Conclusions: A method to perform direct PET-guided radiotherapy planning was successfully validated and specific recommendations were developed to assist other centers. Of major concern is ensuring that the quality of PET and CT data is appropriate for radiotherapy treatment planning and on-treatment verification. Couch movements can be compromised, bore-size can be a limitation for certain immobilization techniques, laser positioning may affect setup accuracy and couch deflection may be greater than scanners dedicated to radiotherapy. The full set of departmental commissioning and routine quality assurance tests applied to radiotherapy CT simulators must be carried out on the PET-CT scanner. CT image quality must be optimized for radiotherapy planning whilst understanding that the appearance will differ between scanners and may affect delineation. PET-CT quality assurance schedules will need to be added to and modified to incorporate radiotherapy quality assurance. Methods of working for radiotherapy and PET staff will change to take into account considerations of both parties. PET to CT alignment must be subject to quality control on a loaded and unloaded couch preferably using a suitable emission phantom, and tested throughout the whole data pathway. Data integrity must be tested throughout the whole pathway and a system included to verify that delineated structures are transferred correctly. Excellent multidisciplinary team communication and working is vital, and key staff members on both sides should be specifically dedicated to the project. Patient pathway should be clearly devised to optimize patient care and the resources of all departments. Recruitment of a cohort of patients into a methodology study is valuable to test the quality assurance methods and pathway.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-J-02: 4D Digital Tomosynthesis Based On Algebraic Image Reconstruction and Total-Variation Minimization for the Improvement of Image Quality

Purpose: In this paper, we implemented the four-dimensional (4D) digital tomosynthesis (DTS) imaging based on algebraic image reconstruction technique and total-variation minimization method in order to compensate the undersampled projection data and improve the image quality. Methods: The projection data were acquired as supposed the cone-beam computed tomography system in linear accelerator by the Monte Carlo simulation and the in-house 4D digital phantom generation program. We performed 4D DTS based upon simultaneous algebraic reconstruction technique (SART) among the iterative image reconstruction technique and total-variation minimization method (TVMM). To verify the effectiveness of this reconstruction algorithm, we performed systematic simulation studies to investigate the imaging performance. Results: The 4D DTS algorithm based upon the SART and TVMM seems to give better results than that based upon the existing method, or filtered-backprojection. Conclusion: The advanced image reconstruction algorithm for the 4D DTS would be useful to validate each intra-fraction motion during radiation therapy. In addition, it will be possible to give advantage to real-time imaging for the adaptive radiation therapy. This research was supported by Leading Foreign Research Institute Recruitment Program (Grant No.2009-00420) and Basic Atomic Energy Research Institute (BAERI); (Grant No. 2009-0078390) through the National Research Foundation of Korea(NRF) funded by the Ministry of Science, ICT and Future Planning (MSIP)

62 RADIOLOGY AND NUCLEAR MEDICINE↗

The Potential of Computational Fluid Dynamics Simulation on Serial Monitoring of Hemodynamic Change in Type B Aortic Dissection

PurposeWe aimed to assess the potential of computational fluid dynamics simulation (CFD) in detecting changes in pressure and flow velocity in response to morphological changes in type B aortic dissection.Materials and MethodsPressure and velocity in four morphological models of type B aortic dissection before and after closure of the entry tear were calculated with CFD and analyzed for changes among the different scenarios. The control model (Model 1) was patient specific and built from the DICOM data of CTA, which bore one entry tear and three re-entry tears. Models 2–4 were modifications of Model 1, with two re-entry tears less in Model 2, one re-entry tear more in Model 3, and a larger entry tear in Model 4.ResultsThe pressure and velocity pertaining to each of the morphological models were unique. Changes in pressure and velocity findings were accountable by the changes in morphological features of the different models. There was no blood flow in the false lumen across the entry tear after its closure, the blood flow direction across the re-entry tears was reversed after closure of the entry tear.ConclusionCFD simulation is probably useful to detect hemodynamic changes in the true and false lumens of type B aortic dissection in response to morphological changes, it may potentially be developed into a non-invasive and patient-specific tool for serial monitoring of hemodynamic changes of type B aortic dissection before and after treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Relationship between electron density and effective densities of body tissues for stopping, scattering, and nuclear interactions of proton and ion beams

Purpose: In treatment planning of charged-particle radiotherapy, patient heterogeneity is conventionally modeled as variable-density water converted from CT images to best reproduce the stopping power, which may lead to inaccuracies in the handling of multiple scattering and nuclear interactions. Although similar conversions can be defined for these individual interactions, they would be valid only for specific CT systems and would require additional tasks for clinical application. This study aims to improve the practicality of the interaction-specific heterogeneity correction. Methods: The authors calculated the electron densities and effective densities for stopping power, multiple scattering, and nuclear interactions of protons and ions, using the standard elemental-composition data for body tissues to construct the invariant conversion functions. The authors also simulated a proton beam in a lung-like geometry and a carbon-ion beam in a prostate-like geometry to demonstrate the procedure and the effects of the interaction-specific heterogeneity correction. Results: Strong correlations were observed between the electron density and the respective effective densities, with which the authors formulated polyline conversion functions. Their effects amounted to 10% differences in multiple-scattering angle and nuclear interaction mean free path for bones compared to those in the conventional heterogeneity correction. Although their realistic effect on patient dose distributions would be generally small, it could be at the level of a few percent when a carbon-ion beam traverses a large bone. Conclusions: The present conversion functions are invariant and may be incorporated in treatment planning systems with a common function relating CT number to electron density. This will enable improved beam dose calculation while minimizing initial setup and quality management of the user's specific system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Validation of a Monte Carlo model for the MLC-Elekta Precise with MCNP6 - Paper 76

Cancer treatment using radiotherapy techniques needs an accurate prediction of the dose delivered to the patient at each position. Treatment planning systems (TPS) allow the optimization of the dose delivered to the tumor volume while minimizing the dose absorbed by the surrounding healthy tissue. Commercial treatment planning systems use empirical deterministic algorithms that do not correctly estimate the dose in the surroundings of heterogeneities. This lack of accuracy is taken into account as bigger security margins in the treatment, reducing the potential effect of the treatment. However, Monte Carlo (MC) simulation techniques have proved to predict the dose distribution more accurately and also representing experimental data better. Nevertheless, the clinical application of these techniques has been limited by the long calculation times they require. As high performance computers decrease in cost MC simulations become more feasible to be used in TPS. The present work improves a simplified model of the MLC-Elekta Precise linear accelerator (linac or unit from now on) previously developed and correctly validated by the group. A multileaf collimator (MLC) has been modelled according to the manufacturer specifications and added to the previous model to allow higher conformity in the studies. The full geometry has been entered to the MCNP6 Monte Carlo code to simulate the transport of particles through the linac. The geometry also includes a cubic water phantom for dosimetry studies and validation with experimental data. The experimental data employed for the dose distribution comparison has been measured at the linac from the linac Hospital Universitari Clinic de Valencia (Spain). Additionally, particles are stored in a phase space surface before their entrance in the water to be analyzed and for performance optimization. Given the results presented, we conclude that this first approach to validate a model of the Elekta Precise linac from the Hospital Cl nic Universitari de Valancia represents well the depth-dose distribution, but fails at representing the dose profiles. It has been shown that they predict too much penumbra, probably caused by the filtration between the leaves. In order to solve this problems, further work is needed to improve the simulated model. This work will mainly focus on the multileaf collimator. Also, the codes developed by the group to analyze phase space files produced by MCNP6 have proven highly useful. As work progress and the codes are debugged, they will be published so that the community can take advantage of them.

61 RADIATION PROTECTION AND DOSIMETRY↗

Application of Geant4 Monte Carlo simulation in dose calculations for small radiosurgical fields

The Geant4 toolkit was used to develop a Monte Carlo (MC)-based engine for accurate dose calculations in small radiation field sizes. The Geant4 toolkit (version 10.1.p02) was used to simulate 6-MV photon beam of a Varian2100C linear accelerator that is being used for stereotactic radiosurgery (SRS) treatment with small radiation fields. Geometric models of 3 in-house designed radiosurgical divergent cones, with the diameters of their projections at the isocenter being 10, 20, and 30 mm, were simulated. The accuracy of the MC simulation technique was examined by reproducing several different simulated dosimetric parameters of the primary beams with the experimental data. The dose distributions are first checked for single beams for each cone, then standard multiple field (SMF) techniques are applied. A sample set of DICOM files from computed tomography (CT) scan imaging of a patient's head was converted to the Geant4 geometry format to implement MC-based engine for a clinical test. To validate the accuracy of the MC-based calculations for SMF arrangements, the isodose lines from MC simulation in water phantom were compared with the measured isodose lines using EBT3 Gafchromic film in Solid Water phantoms. Agreements between measured and simulated depth dose values and beam profiles for SRS cones were generally within 2%/2 mm. For output factors, the largest discrepancy was observed for 10 mm SRS cone, which was 1.7%. For SMF techniques, in SRS cones, the MC simulation and EBT3 Gafchromic film dosimetry were in acceptable agreement (5%/5 mm). Excellent agreement between the results of the MC-based and measured dose values for both single and SMF techniques in SRS cones indicates the ability of the Geant4 toolkit to be applied as the platform for treatment planning of advanced radiotherapy techniques.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

MO-E-17A-08: Attenuation-Based Size Adjusted, Scanner-Independent Organ Dose Estimates for Head CT Exams: TG 204 for Head CT

Purpose: AAPM Task Group 204 described size specific dose estimates (SSDE) for body scans. The purpose of this work is to use a similar approach to develop patient-specific, scanner-independent organ dose estimates for head CT exams using an attenuation-based size metric. Methods: For eight patient models from the GSF family of voxelized phantoms, dose to brain and lens of the eye was estimated using Monte Carlo simulations of contiguous axial scans for 64-slice MDCT scanners from four major manufacturers. Organ doses were normalized by scannerspecific 16 cm CTDIvol values and averaged across all scanners to obtain scanner-independent CTDIvol-to-organ-dose conversion coefficients for each patient model. Head size was measured at the first slice superior to the eyes; patient perimeter and effective diameter (ED) were measured directly from the GSF data. Because the GSF models use organ identification codes instead of Hounsfield units, water equivalent diameter (WED) was estimated indirectly. Using the image data from 42 patients ranging from 2 weeks old to adult, the perimeter, ED and WED size metrics were obtained and correlations between each metric were established. Applying these correlations to the GSF perimeter and ED measurements, WED was calculated for each model. The relationship between the various patient size metrics and CTDIvol-to-organ-dose conversion coefficients was then described. Results: The analysis of patient images demonstrated the correlation between WED and ED across a wide range of patient sizes. When applied to the GSF patient models, an exponential relationship between CTDIvol-to-organ-dose conversion coefficients and the WED size metric was observed with correlation coefficients of 0.93 and 0.77 for the brain and lens of the eye, respectively. Conclusion: Strong correlation exists between CTDIvol normalized brain dose and WED. For the lens of the eye, a lower correlation is observed, primarily due to surface dose variations. Funding Support: Siemens-UCLA Radiology Master Research Agreement; Disclosures - Michael McNitt-Gray: Institutional Research Agreement, Siemens AG; Research Support, Siemens AG; Consultant, Flaherty Sensabaugh Bonasso PLLC; Consultant, Fulbright and Jaworski.

60 APPLIED LIFE SCIENCES↗

Imaging of moving fiducial markers during radiotherapy using a fast, efficient active pixel sensor based EPID

Purpose: The purpose of this work was to investigate the use of an experimental complementary metal-oxide-semiconductor (CMOS) active pixel sensor (APS) for tracking of moving fiducial markers during radiotherapy. Methods: The APS has an active area of 5.4 x 5.4 cm and maximum full frame read-out rate of 20 frame s{sup -1}, with the option to read out a region-of-interest (ROI) at an increased rate. It was coupled to a 4 mm thick ZnWO4 scintillator which provided a quantum efficiency (QE) of 8% for a 6 MV x-ray treatment beam. The APS was compared with a standard iViewGT flat panel amorphous Silicon (a-Si) electronic portal imaging device (EPID), with a QE of 0.34% and a frame-rate of 2.5 frame s{sup -1}. To investigate the ability of the two systems to image markers, four gold cylinders of length 8 mm and diameter 0.8, 1.2, 1.6, and 2 mm were placed on a motion-platform. Images of the stationary markers were acquired using the APS at a frame-rate of 20 frame s{sup -1}, and a dose-rate of 143 MU min{sup -1} to avoid saturation. EPID images were acquired at the maximum frame-rate of 2.5 frame s{sup -1}, and a reduced dose-rate of 19 MU min{sup -1} to provide a similar dose per frame to the APS. Signal-to-noise ratio (SNR) of the background signal and contrast-to-noise ratio (CNR) of the marker signal relative to the background were evaluated for both imagers at doses of 0.125 to 2 MU. Results: Image quality and marker visibility was found to be greater in the APS with SNR {approx}5 times greater than in the EPID and CNR up to an order of magnitude greater for all four markers. To investigate the ability to image and track moving markers the motion-platform was moved to simulate a breathing cycle with period 6 s, amplitude 20 mm and maximum speed 13.2 mm s{sup -1}. At the minimum integration time of 50 ms a tracking algorithm applied to the APS data found all four markers with a success rate of {>=}92% and positional error {<=}90 {mu}m. At an integration time of 400 ms the smallest marker became difficult to detect when moving. The detection of moving markers using the a-Si EPID was difficult even at the maximum dose-rate of 592 MU min{sup -1} due to the lower QE and longer integration time of 400 ms. Conclusions: This work demonstrates that a fast read-out, high QE APS may be useful in the tracking of moving fiducial markers during radiotherapy. Further study is required to investigate the tracking of markers moving in 3D in a treatment beam attenuated by moving patient anatomy. This will require a larger sensor with ROI read-out to maintain speed and a manageable data-rate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

62 RADIOLOGY AND NUCLEAR MEDICINE↗