Risk of Second Malignant Neoplasms Following Proton Arc Therapy and Volumetric Modulated Arc Therapy for Prostate Cancer

The risk of second malignant neoplasms (SMNs) following prostate radiotherapy is a concern due to the large population of survivors and decreasing age at diagnosis. It is known that parallel-opposed beam proton therapy carries a lower risk than photon IMRT. However, a comparison of SMN risk following proton and photon arc therapies has not previously been reported. The purpose of this study was to predict the ratio of excess relative risk (RRR) of SMN incidence following proton arc therapy to that after volumetric modulated arc therapy (VMAT). Additionally, we investigated the impact of margin size and the effect of risk-minimized proton beam weighting on predicted RRR. Physician-approved treatment plans were created for both modalities for three patients. Therapeutic dose was obtained with differential dose-volume histograms from the treatment planning system, and stray dose was estimated from the literature or calculated with Monte Carlo simulations. Then, various risk models were applied to the total dose. Additional treatment plans were also investigated with varying margin size and risk-minimized proton beam weighting. The mean RRR ranged from 0.74 to 0.99, depending on risk model. The additional treatment plans revealed that the RRR remained approximately constant with varying margin size, and that the predicted RRR was reduced by 12% using a risk-minimized proton arc therapy planning technique. In conclusion, proton arc therapy was found to provide an advantage over VMAT in regard to predicted risk of SMN following prostate radiotherapy. This advantage was independent of margin size and was amplified with risk-optimized proton beam weighting.

Clinical Impact of Couch Top and Rails on IMRT and Arc Therapy

Purpose: To evaluate the clinical impact of the Varian Exact Couch on dose and volume coverage to targets and critical structures and tumor control probability (TCP) for 6-MV IMRT and Arc Therapy. Methods: Five clinical prostate patients were planned with both, 6-MV 8-field IMRT and 6-MV 2-field RapidArc using the Eclipse treatment planning system (TPS). These plans neglected treatment couch attenuation, as is standard clinical practice. Dose distributions were then recalculated in Eclipse with the inclusion of the Varian Exact Couch (imaging couch top) and the rails in varying configurations. The changes in dose and coverage were evaluated using the DVHs from each plan iteration. We used a tumor control probability (TCP) model to calculate losses in tumor control resulting from not accounting for the couch top and rails. We also verified dose measurements in a phantom. Results: Failure to account for the treatment couch and rails resulted in clinically unacceptable dose and volume coverage losses to the target for both IMRT and RapidArc. The couch caused average dose losses (relative to plans that ignored the couch) to the prostate of 4.2% and 2.0% for IMRT with the rails out and in, respectively, and 3.2% and 2.9% for RapidArc with the rails out and in, respectively. On average, the percentage of the target covered by the prescribed dose dropped to 35% and 84% for IMRT (rails out and in, respectively) and to 18% and 17% for RapidArc (rails out and in, respectively). The TCP was also reduced by as much as 10.5% (6.3% on average). Dose and volume coverage losses for IMRT plans were primarily due to the rails, while the imaging couch top contributed most to losses for RapidArc. Both the couch top and rails contribute to dose and coverage losses that can render plans clinically unacceptable. A follow-up study we performed found that the less attenuating unipanel mesh couch top available with the Varian Exact couch does not cause a clinically impactful loss of dose or coverage for IMRT but still causes an unacceptable loss for RapidArc. Conclusions: Both the imaging couch top and rails contribute to dose and coverage loss to a degree that, if included, would prevent the plan from meeting clinical planning criteria. Therefore, the imaging and mesh couch tops and rails should be accounted for in Arc Therapy and the imaging couch and rails only in IMRT treatment planning.

Development and implementation of a dynamic heterogeneous proton equivalent anthropomorphic thorax phantom for the assessment of scanned proton beam therapy

DEVELOPMENT AND IMPLEMENTATION OF A DYNAMIC HETEROGENEOUS PROTON EQUIVALENT ANTHROPOMORPHIC THORAX PHANTOM FOR THE ASSESSMENT OF SCANNED PROTON BEAM THERAPY by James Leroy Neihart, B.S. APPROVED: ______________________________David Followill, Ph.D. ______________________________Peter Balter, Ph.D. ______________________________Narayan Sahoo, Ph.D. ______________________________Kenneth Hess, Ph.D. ______________________________Paige Summers, M.S. APPROVED: ____________________________ Dean, The University of Texas Graduate School of Biomedical Sciences at Houston DEVELOPMENT AND IMPLEMENTATION OF A DYNAMIC HETEROGENEOUS PROTON EQUIVALENT ANTHROPOMORPHIC THORAX PHANTOM FOR THE ASSESSMENT OF SCANNED PROTON BEAM THERAPY A THESIS Presented to the Faculty of The University of Texas Health Science Center at Houston andThe University of TexasMD Anderson Cancer CenterGraduate School of Biomedical Sciences in Partial Fulfillment of the Requirements for the Degree of MASTER OF SCIENCE by James Leroy Neihart, B.S. Houston, Texas Date of Graduation August, 2013 Acknowledgments I would like to acknowledge my advisory committee members, chair David Followill, Ph.D., Peter Balter, Ph.D, Narayan Sahoo, Ph.D., Kenneth Hess, Ph.D., Paige Summers M.S. and, for their time and effort contributed to this project. I would additionally like to thank the faculty and staff at the PTC-H and the RPC who assisted in many aspects of this project. Falk Pӧnisch, Ph.D. for his breath hold proton therapy treatment expertise, Matt Palmer and Jaques Bluett for proton dosimetry assistance, Matt Kerr for verification plan assistance, Carrie Amador, Nadia Hernandez, Trang Nguyen, Andrea Molineu, Lynda McDonald for TLD and film dosimetry assistance. Finally, I would like to thank my wife and family for their support and encouragement during my research and studies. Development and implementation of a dynamic heterogeneous proton equivalent anthropomorphic thorax phantom for the assessment of scanned proton beam therapy By: James Leroy Neihart, B.S. Chair of Advisory Committee: David Followill, Ph.D Proton therapy has been gaining ground recently in radiation oncology. To date, the most successful utilization of proton therapy is in head and neck cases as well as prostate cases. These tumor locations do not suffer from the resulting difficulties of treatment delivery as a result of respiratory motion. Lung tumors require either breath hold or motion tracking, neither of which have been assessed with an end-to-end phantom for proton treatments. Currently, the RPC does not have a dynamic thoracic phantom for proton therapy procedure assessment. Additionally, such a phantom could be an excellent means of assessing quality assurance of the procedures of proton therapy centers wishing to participate in clinical trials. An eventual goal of this phantom is to have a means of evaluating and auditing institutions for the ability to start clinical trials utilizing proton therapy procedures for lung cancers. Therefore, the hypothesis of this study is that a dynamic anthropomorphic thoracic phantom can be created to evaluate end-to-end proton therapy treatment procedures for lung cancer to assure agreement between the measured and calculated dose within 5% / 5 mm with a reproducibility of 2%. Multiple materials were assessed for thoracic heterogeneity equivalency. The phantom was designed from the materials found to be in greatest agreement. The phantom was treated in an end-to-end treatment four times, which included simulation, treatment planning and treatment delivery. Each treatment plan was delivered three times to assess reproducibility. The dose measured within the phantom was compared to that of the treatment plan. The hypothesis was fully supported for three of the treatment plans, but failed the reproducibility requirement for the most aggressive treatment plan.

Volumetric Modulated Arc Therapy Evaluation with the Radiological Physics Center Head and Neck Phantom

Validation of treatment plan quality and dose calculation accuracy is essential for new radiotherapy techniques, including volumetric modulated arc therapy (VMAT). VMAT delivers intensity modulated radiotherapy treatments while simultaneously rotating the gantry, adding an additional level of complexity to both the dose calculation and delivery of VMAT treatments compared to static gantry IMRT. The purpose of this project was to compare two VMAT systems, Elekta VMAT and Varian RapidArc, to the current standard of care, IMRT, in terms of both treatment plan quality and dosimetric delivery accuracy using the Radiological Physics Center (RPC) head and neck (H&N) phantom. Clinically relevant treatment plans were created for the phantom using typical prescription and dose constraints for Elekta VMAT (planned with Pinnacle3 Smart Arc) and RapidArc and IMRT (both planned with Eclipse). The treatment plans were evaluated to determine if they were clinically comparable using several dosimetric criteria, including ability to meet dose objectives, hot spots, conformity index, and homogeneity index. The planned treatments were delivered to the phantom and absolute doses and relative dose distributions were measured with thermoluminescent dosimeters (TLDs) and radiochromic film, respectively. The measured and calculated doses of each treatment were compared to determine if they were clinically acceptable based upon RPC criteria of ±7% dose difference and 4 mm distance-to-agreement. Gamma analysis was used to assess dosimetric accuracy, as well. All treatment plans were able to meet the dosimetric objectives set by the RPC and had similar hot spots in the normal tissue. The Elekta VMAT plan was more homogenous but less conformal than the RapidArc and IMRT plans. When comparing the measured and calculated doses, all plans met the RPC ±7%/4 mm criteria. The percent of points passing the gamma analysis for each treatment delivery was acceptable. Treatment plan quality of the Elekta VMAT, RapidArc and IMRT treatments were comparable for consistent dose prescriptions and constraints. Additionally, the dosimetric accuracy of the Elekta VMAT and RapidArc treatments was verified to be within acceptable tolerances.

MODIFICATION AND IMPLEMENTATION OF THE RPC HETEROGENEOUS THORAX PHANTOM FOR VERIFICATION OF PROTON THERAPY TREATMENT PROCEDURES

The Radiological Physics Center (RPC) provides heterogeneous phantoms that are used to evaluate radiation treatment procedures as part of a comprehensive quality assurance program for institutions participating in clinical trials. It was hypothesized that the existing RPC heterogeneous thorax phantom can be modified to assess lung tumor proton beam therapy procedures involving patient simulation, treatment planning, and treatment delivery, and could confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 5%. The Hounsfield Units (HU) for lung equivalent materials (balsa wood and cork) was measured using a CT scanner. The relative linear stopping power (RLSP) of these materials was measured. The linear energy transfer (LET) of Gafchromic EBT2 film was analyzed utilizing parallel and perpendicular orientations in a water tank and compared to ion chamber readings. Both parallel and perpendicular orientations displayed a quenching effect underperforming the ion chamber, with the parallel orientation showing an average 31 % difference and the perpendicular showing an average of 15% difference. Two treatment plans were created that delivered the prescribed dose to the target volume, while achieving low entrance doses. Both treatment plans were designed using smeared compensators and expanded apertures, as would be utilized for a patient in the clinic. Plan 1a contained two beams that were set to orthogonal angles and a zero degree couch kick. Plan 1b utilized two beams set to 10 and 80 degrees with a 15 degree couch kick. EBT2 film and TLD were inserted and the phantom was irradiated 3 times for each plan. Both plans passed the criteria for the TLD measurements where the TLD values were within 7% of the dose calculated by Eclipse. Utilizing the 5%/3mm criteria, the 3 trial average of overall pass rate was 71% for Plan 1a. The 3 trial average for the overall pass rate was 76% for Plan 1b. The trials were then analyzed using RPC conventional lung treatment guidelines set forth by the RTOG: 5%/5mm, and an overall pass rate of 85%. Utilizing these criteria, only Plan 1b passed for all 3 trials, with an average overall pass rate of 89%.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

RANGE ADAPTIVE PROTON THERAPY FOR PROSTATE CANCER

Purpose: The rapid distal falloff of a proton beam allows for sparing of normal tissues distal to the target. However proton beams that aim directly towards critical structures are avoided due to concerns of range uncertainties, such as CT number conversion and anatomy variations. We propose to eliminate range uncertainty and enable prostate treatment with a single anterior beam by detecting the proton’s range at the prostate-rectal interface and adaptively adjusting the range in vivo and in real-time. Materials and Methods: A prototype device, consisting of an endorectal liquid scintillation detector and dual-inverted Lucite wedges for range compensation, was designed to test the feasibility and accuracy of the technique. Liquid scintillation filled volume was fitted with optical fiber and placed inside the rectum of an anthropomorphic pelvic phantom. Photodiode-generated current signal was generated as a function of proton beam distal depth, and the spatial resolution of this technique was calculated by relating the variance in detecting proton spills to its maximum penetration depth. The relative water-equivalent thickness of the wedges was measured in a water phantom and prospectively tested to determine the accuracy of range corrections. Treatment simulation studies were performed to test the potential dosimetric benefit in sparing the rectum. Results: The spatial resolution of the detector in phantom measurement was 0.5 mm. The precision of the range correction was 0.04 mm. The residual margin to ensure CTV coverage was 1.1 mm. The composite distal margin for 95% treatment confidence was 2.4 mm. Planning studies based on a previously estimated 2mm margin (90% treatment confidence) for 27 patients showed a rectal sparing up to 51% at 70 Gy and 57% at 40 Gy relative to IMRT and bilateral proton treatment. Conclusion: We demonstrated the feasibility of our design. Use of this technique allows for proton treatment using a single anterior beam, significantly reducing the rectal dose.

Implementation of an Anthropomorphic Phantom for the Evaluation of Proton Therapy Treatment Procedures

With an increasing number of institutions offering proton therapy, the number of multi-institutional clinical trials involving proton therapy will also increase in the coming years. The Radiological Physics Center monitors sites involved in clinical trials through the use of site visits and remote auditing with thermoluminescent dosimeters (TLD) and mailable anthropomorphic phantoms. Currently, there are no heterogeneous phantoms that have been commissioned to evaluate proton therapy. It was hypothesized that an anthropomorphic pelvis phantom can be designed to audit treatment procedures (patient simulation, treatment planning and treatment delivery) at proton facilities to confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 3%. A pelvis phantom originally designed for use with photon treatments was retrofitted for use in proton therapy. The relative stopping power (SP) of each phantom material was measured. Hounsfield Units (HU) for each phantom material were measured with a CT scanner and compared to the relative stopping power calibration curve. The tissue equivalency for each material was calculated. Two proton treatment plans were created; one which did not correct for material SP differences (Plan 1) and one plan which did correct for SP differences (Plan 2). Film and TLD were loaded into the phantom and the phantom was irradiated 3 times per plan. The measured values were compared to the HU-SP calibration curve and it was found that the stopping powers for the materials could be underestimated by 5-10%. Plan 1 passed the criteria for the TLD and film margins with reproducibility under 3% between the 3 trials. Plan 2 failed because the right-left film dose profile average displacement was -9.0 mm on the left side and 6.0 mm on the right side. Plan 2 was intended to improve the agreements and instead introduced large displacements along the path of the beam. Plan 2 more closely represented the actual phantom composition with corrected stopping powers and should have shown an agreement between the measured and calculated dose within 5%/3mm. The hypothesis was rejected and the pelvis phantom was found to be not suitable to evaluate proton therapy treatment procedures.

The Development and Implementation of an Anthropomorphic Head Phantom for the Assessment of Proton Therapy Treatment Procedures

Proton therapy has become an increasingly more common method of radiation therapy, with the dose sparing to distal tissue making it an appealing option, particularly for treatment of brain tumors. This study sought to develop a head phantom for the Radiological Physics Center (RPC), the first to be used for credentialing of institutions wishing to participate in clinical trials involving brain tumor treatment of proton therapy. It was hypothesized that a head phantom could be created for the evaluation of proton therapy treatment procedures (treatment simulation, planning, and delivery) to assure agreement between the measured dose and calculated dose within ±5%/3mm with a reproducibility of ±3%. The relative stopping power (RSP) and Hounsfield Units (HU) were measured for potential phantom materials and a human skull was cast in tissue-equivalent Alderson material (RLSP 1.00, HU 16) with anatomical airways and a cylindrical hole for imaging and dosimetry inserts drilled into the phantom material. Two treatment plans, proton passive scattering and proton spot scanning, were created. Thermoluminescent dosimeters (TLDs) and film were loaded into the phantom dosimetry insert. Each treatment plan was delivered three separate times. Each treatment plan passed our 5%/3mm criteria, with a reproducibility of ±3%. The hypothesis was accepted and the phantom was found to be suitable for remote audits of proton therapy treatment facilities.

Monte Carlo and analytical dose calculations for ocular proton therapy

Uveal melanoma is a rare but life-threatening form of ocular cancer. Contemporary treatment techniques include proton therapy, which enables conservation of the eye and its useful vision. Dose to the proximal structures is widely believed to play a role in treatment side effects, therefore, reliable dose estimates are required for properly evaluating the therapeutic value and complication risk of treatment plans. Unfortunately, current simplistic dose calculation algorithms can result in errors of up to 30% in the proximal region. In addition, they lack predictive methods for absolute dose per monitor unit (D/MU) values. ^ To facilitate more accurate dose predictions, a Monte Carlo model of an ocular proton nozzle was created and benchmarked against measured dose profiles to within ±3% or ±0.5 mm and D/MU values to within ±3%. The benchmarked Monte Carlo model was used to develop and validate a new broad beam dose algorithm that included the influence of edgescattered protons on the cross-field intensity profile, the effect of energy straggling in the distal portion of poly-energetic beams, and the proton fluence loss as a function of residual range. Generally, the analytical algorithm predicted relative dose distributions that were within ±3% or ±0.5 mm and absolute D/MU values that were within ±3% of Monte Carlo calculations. Slightly larger dose differences were observed at depths less than 7 mm, an effect attributed to the dose contributions of edge-scattered protons. Additional comparisons of Monte Carlo and broad beam dose predictions were made in a detailed eye model developed in this work, with generally similar findings. ^ Monte Carlo was shown to be an excellent predictor of the measured dose profiles and D/MU values and a valuable tool for developing and validating a broad beam dose algorithm for ocular proton therapy. The more detailed physics modeling by the Monte Carlo and broad beam dose algorithms represent an improvement in the accuracy of relative dose predictions over current techniques, and they provide absolute dose predictions. It is anticipated these improvements can be used to develop treatment strategies that reduce the incidence or severity of treatment complications by sparing normal tissue. ^

Proton therapy versus intensity modulated x-ray therapy in the treatment of prostate cancer: Estimating secondary cancer risks

External beam radiation therapy is used to treat nearly half of the more than 200,000 new cases of prostate cancer diagnosed in the United States each year. During a radiation therapy treatment, healthy tissues in the path of the therapeutic beam are exposed to high doses. In addition, the whole body is exposed to a low-dose bath of unwanted scatter radiation from the pelvis and leakage radiation from the treatment unit. As a result, survivors of radiation therapy for prostate cancer face an elevated risk of developing a radiogenic second cancer. Recently, proton therapy has been shown to reduce the dose delivered by the therapeutic beam to normal tissues during treatment compared to intensity modulated x-ray therapy (IMXT, the current standard of care). However, the magnitude of stray radiation doses from proton therapy, and their impact on this incidence of radiogenic second cancers, was not known. ^ The risk of a radiogenic second cancer following proton therapy for prostate cancer relative to IMXT was determined for 3 patients of large, median, and small anatomical stature. Doses delivered to healthy tissues from the therapeutic beam were obtained from treatment planning system calculations. Stray doses from IMXT were taken from the literature, while stray doses from proton therapy were simulated using a Monte Carlo model of a passive scattering treatment unit and an anthropomorphic phantom. Baseline risk models were taken from the Biological Effects of Ionizing Radiation VII report. A sensitivity analysis was conducted to characterize the uncertainty of risk calculations to uncertainties in the risk model, the relative biological effectiveness (RBE) of neutrons for carcinogenesis, and inter-patient anatomical variations. ^ The risk projections revealed that proton therapy carries a lower risk for radiogenic second cancer incidence following prostate irradiation compared to IMXT. The sensitivity analysis revealed that the results of the risk analysis depended only weakly on uncertainties in the risk model and inter-patient variations. Second cancer risks were sensitive to changes in the RBE of neutrons. However, the findings of the study were qualitatively consistent for all patient sizes and risk models considered, and for all neutron RBE values less than 100. ^

DEVELOPMENT OF A BEAM-SPECIFIC PLANNING TARGET VOLUME AND A ROBUST PLAN ANALYSIS TOOL FOR PROTON THERAPY

Proton therapy is growing increasingly popular due to its superior dose characteristics compared to conventional photon therapy. Protons travel a finite range in the patient body and stop, thereby delivering no dose beyond their range. However, because the range of a proton beam is heavily dependent on the tissue density along its beam path, uncertainties in patient setup position and inherent range calculation can degrade thedose distribution significantly. Despite these challenges that are unique to proton therapy, current management of the uncertainties during treatment planning of proton therapy has been similar to that of conventional photon therapy. The goal of this dissertation research was to develop a treatment planning method and a planevaluation method that address proton-specific issues regarding setup and range uncertainties. Treatment plan designing method adapted to proton therapy: Currently, for proton therapy using a scanning beam delivery system, setup uncertainties are largely accounted for by geometrically expanding a clinical target volume (CTV) to a planning target volume (PTV). However, a PTV alone cannot adequately account for range uncertainties coupled to misaligned patient anatomy in the beam path since it does not account for the change in tissue density. In order to remedy this problem, we proposed a beam-specific PTV (bsPTV) that accounts for the change in tissue density along the beam path due to the uncertainties. Our proposed method was successfully implemented, and its superiority over the conventional PTV was shown through a controlled experiment.. Furthermore, we have shown that the bsPTV concept can be incorporated into beam angle optimization for better target coverage and normal tissue sparing for a selected lung cancer patient. Treatment plan evaluation method adapted to proton therapy: The dose-volume histogram of the clinical target volume (CTV) or any other volumes of interest at the time of planning does not represent the most probable dosimetric outcome of a given plan as it does not include the uncertainties mentioned earlier. Currently, the PTV is used as a surrogate of the CTV’s worst case scenario for target dose estimation. However, because proton dose distributions are subject to change under these uncertainties, the validity of the PTV analysis method is questionable. In order to remedy this problem, we proposed the use of statistical parameters to quantify uncertainties on both the dose-volume histogram and dose distribution directly. The robust plan analysis tool was successfully implemented to compute both the expectation value and its standard deviation of dosimetric parameters of a treatment plan under the uncertainties. For 15 lung cancer patients, the proposed method was used to quantify the dosimetric difference between the nominal situation and its expected value under the uncertainties.

Dynamic Chemical Shift Imaging for Image-Guided Thermal Therapy

Magnetic resonance temperature imaging (MRTI) is recognized as a noninvasive means to provide temperature imaging for guidance in thermal therapies. The most common method of estimating temperature changes in the body using MR is by measuring the water proton resonant frequency (PRF) shift. Calculation of the complex phase difference (CPD) is the method of choice for measuring the PRF indirectly since it facilitates temperature mapping with high spatiotemporal resolution. Chemical shift imaging (CSI) techniques can provide the PRF directly with high sensitivity to temperature changes while minimizing artifacts commonly seen in CPD techniques. However, CSI techniques are currently limited by poor spatiotemporal resolution. This research intends to develop and validate a CSI-based MRTI technique with intentional spectral undersampling which allows relaxed parameters to improve spatiotemporal resolution. An algorithm based on autoregressive moving average (ARMA) modeling is developed and validated to help overcome limitations of Fourier-based analysis allowing highly accurate and precise PRF estimates. From the determined acquisition parameters and ARMA modeling, robust maps of temperature using the k-means algorithm are generated and validated in laser treatments in ex vivo tissue. The use of non-PRF based measurements provided by the technique is also investigated to aid in the validation of thermal damage predicted by an Arrhenius rate dose model.

The Role of Cell Sterilization in Population Based Studies of Radiogenic Second Cancers Following Radiation Therapy

Advances in radiotherapy have generated increased interest in comparative studies of treatment techniques and their effectiveness. In this respect, pediatric patients are of specific interest because of their sensitivity to radiation induced second cancers. However, due to the rarity of childhood cancers and the long latency of second cancers, large sample sizes are unavailable for the epidemiological study of contemporary radiotherapy treatments. Additionally, when specific treatments are considered, such as proton therapy, sample sizes are further reduced due to the rareness of such treatments. We propose a method to improve statistical power in micro clinical trials. Specifically, we use a more biologically relevant quantity, cancer equivalent dose (DCE), to estimate risk instead of mean absorbed dose (DMA). Our objective was to demonstrate that when DCE is used fewer subjects are needed for clinical trials. Thus, we compared the impact of DCE vs. DMA on sample size in a virtual clinical trial that estimated risk for second cancer (SC) in the thyroid following craniospinal irradiation (CSI) of pediatric patients using protons vs. photons. Dose reconstruction, risk models, and statistical analysis were used to evaluate SC risk from therapeutic and stray radiation from CSI for 18 patients. Absorbed dose was calculated in two ways: with (1) traditional DMA and (2) with DCE. DCE and DMA values were used to estimate relative risk of SC incidence (RRCE and RRMA, respectively) after proton vs. photon CSI. Ratios of RR for proton vs. photon CSI (RRRCE and RRRMA) were then used in comparative estimations of sample size to determine the minimal number of patients needed to maintain 80% statistical power when using DCE vs. DMA. For all patients, we found that protons substantially reduced the risk of developing a second thyroid cancer when compared to photon therapy. Mean RRR values were 0.052±0.014 and 0.087±0.021 for RRRMA and RRRCE, respectively. However, we did not find that use of DCE reduced the number of patents needed for acceptable statistical power (i.e, 80%). In fact, when considerations were made for RRR values that met equipoise requirements and the need for descriptive statistics, the minimum number of patients needed for a micro-clinical trial increased from 17 using DMA to 37 using DCE. Subsequent analyses revealed that for our sample, the most influential factor in determining variations in sample size was the experimental standard deviation of estimates for RRR across the patient sample. Additionally, because the relative uncertainty in dose from proton CSI was so much larger (on the order of 2000 times larger) than the other uncertainty terms, it dominated the uncertainty in RRR. Thus, we found that use of corrections for cell sterilization, in the form of DCE, may be an important and underappreciated consideration in the design of clinical trials and radio-epidemiological studies. In addition, the accurate application of cell sterilization to thyroid dose was sensitive to variations in absorbed dose, especially for proton CSI, which may stem from errors in patient positioning, range calculation, and other aspects of treatment planning and delivery.

Fast magnetic resonance temperature imaging for focused ultrasound thermal therapy

The current standard for temperature sensitive imaging using magnetic resonance (MR) is 2-D, spoiled, fast gradient-echo (fGRE) phase-difference imaging exploiting temperature dependent changes in the proton resonance frequency (PRF). The echo-time (TE) for optimal sensitivity is larger than the typical repetition time (TR) of an fGRE sequence. Since TE must be less than TR in the fGRE sequence, this limits the technique's achievable sensitivity, spatial, and temporal resolution. This adversely affects both accuracy and volume coverage of the measurements. Accurate measurement of the rapid temperature changes associated with pulsed thermal therapies, such as high-intensity focused ultrasound (FUS), at optimal temperature sensitivity requires faster acquisition times than those currently available. ^ Use of fast MR acquisition strategies, such as interleaved echo-planar and spiral imaging, can provide the necessary increase in temporal performance and sensitivity while maintaining adequate signal-to-noise and in-plane spatial resolution. This research explored the adaptation and optimization of several fast MR acquisition methods for thermal monitoring of pulsed FUS thermal therapy. Temperature sensitivity, phase-difference noise and phase-difference to phase-difference-to noise ratio for the different pulse sequences were evaluated under varying imaging parameters in an agar gel phantom to establish optimal sequence parameters for temperature monitoring. The temperature sensitivity coefficient of the gel phantom was measured, allowing quantitative temperature extrapolations. ^ Optimized fast sequences were compared based on the ability to accurately monitor temperature changes at the focus of a high-intensity focused ultrasound unit, volume coverage, and contrast-to-noise ratio in the temperature maps. Operating parameters, which minimize complex phase-difference measurement errors introduced by use of the fast-imaging methods, were established. ^