Combined clinical and research training in medical physics in a multi-institutional setting: 13-year experience of Harvard Medical Physics Residency Program.

PURPOSE: This manuscript describes the structure, management and outcomes of a multi-institutional clinical and research medical physics residency program (Harvard Medical Physics Residency Program, or HMPRP) to provide potentially useful information to the centers considering a multi-institutional approach for their training programs. METHODS: Data from the program documents and public records was used to describe HMPRP and obtain statistics about participating faculty, enrolled residents, and graduates. Challenges associated with forming and managing a multi-institutional program and developed solutions for effective coordination between several clinical centers are described. RESULTS: HMPRP was formed in 2009 and was accredited by the Commission on Accreditation of Medical Physics Education Programs (CAMPEP) in 2011. It is a 3-year therapy program, with a dedicated year of research and the 2 years of clinical training at three academic hospitals. A CAMPEP-accredited Certificate Program is embedded in HMPRP to allow enrolled residents to complete a formal didactic training in medical physics if necessary. The clinical training covers the material required by CAMPEP. In addition, training in protons, CyberKnife, MR-linac, and at network locations is included. The clinical training and academic record of the residents is outstanding. All graduates have found employment within clinical medical physics, mostly at large academic centers and graduates had a 100% pass rate at the oral American Board of Radiology exams. On average, three manuscripts per resident are published during residency, and multiple abstracts are presented at conferences. CONCLUSIONS: A multi-institutional medical physics residency program can be successfully formed and managed. With a collaborative administrative structure, the program creates an environment for high-quality clinical training of the residents and high productivity in research. The main advantage of such program is access to a wide variety of resources. The main challenge is creating a structure for efficient management of multiple resources at different locations. This report may provide valuable information to centers considering starting a multi-institutional residency program.

Practice Patterns of Stereotactic Radiotherapy in Pediatrics: Results From an International Pediatric Research Consortium.

PURPOSE/OBJECTIVES: There is little consensus regarding the application of stereotactic radiotherapy (SRT) in pediatrics. We evaluated patterns of pediatric SRT practice through an international research consortium. MATERIALS AND METHODS: Eight international institutions with pediatric expertise completed a 124-item survey evaluating patterns of SRT use for patients 21 years old and younger. Frequencies of SRT use and median margins applied with and without SRT were evaluated. RESULTS: Across institutions, 75% reported utilizing SRT in pediatrics. SRT was used in 22% of brain, 18% of spine, 16% of other bone, 16% of head and neck, and <1% of abdomen/pelvis, lung, and liver cases across sites. Of the hypofractionated SRT cases, 42% were delivered with definitive intent. Median gross tumor volume to planning target volume margins for SRT versus non-SRT plans were 0.2 versus 1.4 cm for brain, 0.3 versus 1.5 cm for spine/other bone, 0.3 versus 2.0 cm for abdomen/pelvis, 0.7 versus 1.5 cm for head and neck, 0.5 versus 1.7 cm for lung, and 0.5 versus 2.0 cm for liver sites. CONCLUSIONS: SRT is commonly utilized in pediatrics across a range of treatment sites. Margins used for SRT were substantially smaller than for non-SRT planning, highlighting the utility of this approach in reducing treatment volumes.

Practice patterns of palliative radiation therapy in pediatric oncology patients in an international pediatric research consortium.

BACKGROUND/OBJECTIVES: The practice of palliative radiation therapy (RT) is based on extrapolation from adult literature. We evaluated patterns of pediatric palliative RT to describe regimens used to identify opportunity for future pediatric-specific clinical trials. DESIGN/METHODS: Six international institutions with pediatric expertise completed a 122-item survey evaluating patterns of palliative RT for patients ≤21 years old from 2010 to 2015. Two institutions use proton RT. Palliative RT was defined as treatment with the goal of symptom control or prevention of immediate life-threatening progression. RESULTS: Of 3,225 pediatric patients, 365 (11%) were treated with palliative intent to a total of 427 disease sites. Anesthesia was required in 10% of patients. Treatment was delivered to metastatic disease in 54% of patients. Histologies included neuroblastoma (30%), osteosarcoma (18%), leukemia/lymphoma (12%), rhabdomyosarcoma (12%), medulloblastoma/ependymoma (12%), Ewing sarcoma (8%), and other (8%). Indications included pain (43%), intracranial symptoms (23%), respiratory compromise (14%), cord compression (8%), and abdominal distention (6%). Sites included nonspine bone (35%), brain (16% primary tumors, 6% metastases), abdomen/pelvis (15%), spine (12%), head/neck (9%), and lung/mediastinum (5%). Re-irradiation comprised 16% of cases. Techniques employed three-dimensional conformal RT (41%), intensity-modulated RT (23%), conventional RT (26%), stereotactic body RT (6%), protons (1%), electrons (1%), and other (2%). The most common physician-reported barrier to consideration of palliative RT was the concern about treatment toxicity (83%). CONCLUSION: There is significant diversity of practice in pediatric palliative RT. Combined with ongoing research characterizing treatment response and toxicity, these data will inform the design of forthcoming clinical trials to establish effective regimens and minimize treatment toxicity for this patient population.

Geometric and dosimetric uncertainties in intracranial stereotatctic treatments for multiple nonisocentric lesions.

The purpose of this study was to determine the effects of geometric uncertainties of patient position on treatments of multiple nonisocentric intracranial lesions. The average distance between lesions in patients with multiple targets was determined by a retrospective survey of patients with multiple lesions. Retrospective patient imaging data from fractionated stereotactic patients were used to calculate interfractional and intrafractional patient position uncertainty. Three different immobilization devices were included in the positioning study. The interfractional and intrafractional patient positioning error data were used to calculate the geometric offset of a lesion located at varying distances from the mechanical isocenter for treatments of multiple lesions with a single arc, assuming that no intrafractional position correction is employed during an arc rotation. Dosimetric effects were studied using two representative lesions of two sizes, 6 mm and 13 mm maximum dimensions, and prescribed to 20 Gy and 18 Gy, respectively. Distances between lesions ranged from < 10 mm to 150 mm, which would correspond to a range of isocenter to lesion separations of < 10 mm to 75 mm, assuming an isocenter located at the geometric mean. In the presence of a full six degree of freedom patient correction system, the effects of the intrafractional patient positioning uncertainties were less than 1.8 mm (3.6mm) for 1σ (2σ) deviations for lesion spacing up to 75 mm assuming a quadratic summation of 1σ and 2σ. Without the benefit of a six DOF correction device, only correcting for three translations, the effects of the intrafractional patient positioning uncertainties were within 3.1 mm (7.2 mm) for 1σ (2σ) deviations for distances up to 75 mm. 1σ and 2σ deviations along all six axes were observed in 3.6% and 0.3%, respectively, of 974 fractions analyzed. Dosimetric effects for 2 mm and 4 mm offsets were most significant for the small lesion with minimum dose (Dmin) decreasing from 20 Gy to 13.6 Gy and 5.7 Gy and volume receiving the prescription (V20Gy) reducing from 100% to 57% and 16%, respectively. The dosimetric effects on the larger lesion were less pronounced with Dmin reducing from 18 Gy to 17.5 Gy and 14.2 Gy, and V18Gy reducing from 100% to 98.3% and 85.4%, for 2 mm and 4 mm offsets, respectively. In the 1σ scenario (3.6% of patients) angular uncertainties in patient positioning can introduce 1.0 mm shifts in the location of the lesion position at distances of 75 mm, compared to an isocentric treatment even with a full six DOF correction. Without the ability to correct angular positioning errors, a lesion positioned 75 mm away from the mechanical isocenter can be located in 3.6% of patients > 3.0 mm distant from the planned position. Dosimetric results depend upon the distance from isocenter and the size of the target. Single isocenter treatments for multiple lesions should be considered only when full six DOF corrections can be applied, the intrafractional immobilization precision is well quantified, and a PTV expansion is included for more distant lesions to account for unavoidable residual patient positioning uncertainties.

Feasibility of Using Distal Endpoints for In-room PET Range Verification of Proton Therapy.

In an effort to verify the dose delivery in proton therapy, Positron Emission Tomography (PET) scans have been employed to measure the distribution of ß+ radioactivity produced from nuclear reactions of the protons with native nuclei. Because the dose and PET distributions are difficult to compare directly, the range verification is currently carried out by comparing measured and Monte Carlo (MC) simulation predicted PET distributions. In order to reduce the reliance on MC, simulated PET (simPET) and dose distal endpoints were compared to explore the feasibility of using distal endpoints for in-room PET range verification. MC simulations were generated for six head and neck patients with corrections for radiological decay, biological washout, and PET resolution. One-dimensional profiles of the dose and simPET were examined along the direction of the beam and covering the cross section of the beam. The chosen endpoints of the simPET (x-intercept of the linear fit to the distal falloff) and planned dose (20-50% of maximum dose) correspond to where most of the protons are below the threshold energy for the nuclear reactions. The difference in endpoint range between the distal surfaces of the dose and MC-PET were compared and the spread of range differences were assessed. Among the six patients, the mean difference between MC-PET and dose depth was found to be -1.6 mm to +0.5 mm between patients, with a standard deviation of 1.1 to 4.0 mm across the individual beams. In clinical practice, regions with deviations beyond the safety margin need to be examined more closely and can potentially lead to adjustments to the treatment plan.

Evaluating the effect of setup uncertainty reduction and adaptation to geometric changes on normal tissue complication probability using online adaptive head and neck intensity modulated proton therapy.

Objective. To evaluate the impact of setup uncertainty reduction (SUR) and adaptation to geometrical changes (AGC) on normal tissue complication probability (NTCP) when using online adaptive head and neck intensity modulated proton therapy (IMPT).Approach.A cohort of ten retrospective head and neck cancer patients with daily scatter corrected cone-beam CT (CBCT) was studied. For each patient, two IMPT treatment plans were created: one with a 3 mm setup uncertainty robustness setting and one with no explicit setup robustness. Both plans were recalculated on the daily CBCT considering three scenarios: the robust plan without adaptation, the non-robust plan without adaptation and the non-robust plan with daily online adaptation. Online-adaptation was simulated using an in-house developed workflow based on GPU-accelerated Monte Carlo dose calculation and partial spot-intensity re-optimization. Dose distributions associated with each scenario were accumulated on the planning CT, where NTCP models for six toxicities were applied. NTCP values from each scenario were intercompared to quantify the reduction in toxicity risk induced by SUR alone, AGC alone and SUR and AGC combined. Finally, a decision tree was implemented to assess the clinical significance of the toxicity reduction associated with each mechanism.Main results. For most patients, clinically meaningful NTCP reductions were only achieved when SUR and AGC were performed together. In these conditions, total reductions in NTCP of up to 30.48 pp were obtained, with noticeable NTCP reductions for aspiration, dysphagia and xerostomia (mean reductions of 8.25, 5.42 and 5.12 pp respectively). While SUR had a generally larger impact than AGC on NTCP reductions, SUR alone did not induce clinically meaningful toxicity reductions in any patient, compared to only one for AGC alone.SignificanceOnline adaptive head and neck proton therapy can only yield clinically significant reductions in the risk of long-term side effects when combining the benefits of SUR and AGC.

Large anatomical changes in head-and-neck cancers - A dosimetric comparison of online and offline adaptive proton therapy.

Purpose: This work evaluates an online adaptive (OA) workflow for head-and-neck (H&N) intensity-modulated proton therapy (IMPT) and compares it with full offline replanning (FOR) in patients with large anatomical changes. Methods: IMPT treatment plans are created retrospectively for a cohort of eight H&N cancer patients that previously required replanning during the course of treatment due to large anatomical changes. Daily cone-beam CTs (CBCT) are acquired and corrected for scatter, resulting in 253 analyzed fractions. To simulate the FOR workflow, nominal plans are created on the planning-CT and delivered until a repeated-CT is acquired; at this point, a new plan is created on the repeated-CT. To simulate the OA workflow, nominal plans are created on the planning-CT and adapted at each fraction using a simple beamlet weight-tuning technique. Dose distributions are calculated on the CBCTs with Monte Carlo for both delivery methods. The total treatment dose is accumulated on the planning-CT. Results: Daily OA improved target coverage compared to FOR despite using smaller target margins. In the high-risk CTV, the median D98 degradation was 1.1 % and 2.1 % for OA and FOR, respectively. In the low-risk CTV, the same metrics yield 1.3 % and 5.2 % for OA and FOR, respectively. Smaller setup margins of OA reduced the dose to all OARs, which was most relevant for the parotid glands. Conclusion: Daily OA can maintain prescription doses and constraints over the course of fractionated treatment, even in cases of large anatomical changes, reducing the necessity for manual replanning in H&N IMPT.

Optically stimulated luminescence dosimeters for simultaneous measurement of point dose and dose-weighted LET in an adaptive proton therapy workflow.

Purpose: To demonstrate the suitability of optically stimulated luminescence detectors (OSLDs) for accurate simultaneous measurement of the absolute point dose and dose-weighted linear energy transfer (LETD) in an anthropomorphic phantom for experimental validation of daily adaptive proton therapy. Methods: A clinically realistic intensity-modulated proton therapy (IMPT) treatment plan was created based on a CT of an anthropomorphic head-and-neck phantom made of tissue-equivalent material. The IMPT plan was optimized with three fields to deliver a uniform dose to the target volume covering the OSLDs. Different scenarios representing inter-fractional anatomical changes were created by modifying the phantom. An online adaptive proton therapy workflow was used to recover the daily dose distribution and account for the applied geometry changes. To validate the adaptive workflow, measurements were performed by irradiating Al2O3:C OSLDs inside the phantom. In addition to the measurements, retrospective Monte Carlo simulations were performed to compare the absolute dose and dose-averaged LET (LETD) delivered to the OSLDs. Results: The online adaptive proton therapy workflow was shown to recover significant degradation in dose conformity resulting from large anatomical and positioning deviations from the reference plan. The Monte Carlo simulations were in close agreement with the OSLD measurements, with an average relative error of 1.4% for doses and 3.2% for LETD. The use of OSLDs for LET determination allowed for a correction for the ionization quenched response. Conclusion: The OSLDs appear to be an excellent detector for simultaneously assessing dose and LET distributions in proton irradiation of an anthropomorphic phantom. The OSLDs can be cut to almost any size and shape, making them ideal for in-phantom measurements to probe the radiation quality and dose in a predefined region of interest. Although we have presented the results obtained in the experimental validation of an adaptive proton therapy workflow, the same approach can be generalized and used for a variety of clinical innovations and workflow developments that require accurate assessment of point dose and/or average LET.

Influence of imaging source and panel position uncertainties on the accuracy of 2D∕3D image registration of cranial images.

PURPOSE: To determine the effects of imager source and panel positioning uncertainties on the accuracy of dual intensity-based 2D∕3D image registration of cranial images. METHODS: An open source 2D∕3D image registration algorithm has been developed for registration of two orthogonal x-rays to a 3D volumetric image. The initialization files of the algorithm allow for nine degrees of freedom system calibration including x, y, z positions of the source and panel, and three rotational degrees of freedom of the panel about each of the three translational axes. A baseline system calibration was established and a baseline 2D∕3D registration between two orthogonal x-rays and the volumetric image was determined. The calibration file was manipulated to insert errors into each of the nine calibration variables of both imager geometries. Rigid six degrees of freedom registrations were iterated for each panel or source positional error over a range of predetermined calibration errors to determine the resulting error in the registration versus the baseline registration due to the manipulated error of the panel or source calibration. RESULTS: Panel and source translational errors orthogonal to the imager∕panel axis introduced the greatest errors in the registration accuracy (4.0 mm geometric error results in up to 2.7 mm registration error). Panel rotation about the imaging direction also resulted in errors of the registration (2.0° geometric error results in up to 1.7° registration error). Differences in magnification and panel tilt and roll, i.e., source and∕or panel translation along the imaging direction and panel rotations about the orthogonal axes had minimal effects on the registration accuracy (below 0.3 mm and 0.2° registration error). CONCLUSIONS: While five of the nine imaging system variables were found to have a considerable effect on 2D∕3D registration accuracy of cranial images, the other four variables showed minimal effects. Vendors typically provide simplified calibration procedures which aim to remove encountered geometric uncertainties by accounting for two panel translations. This study shows that at least the five relevant positional variables should be separately calibrated, if accurate alignment is required for 2D∕3D registration.

Radiation recall myositis in pediatric Ewing sarcoma.

Radiation recall is a rare and poorly understood phenomenon, characterized by an acute inflammatory reaction within the previously irradiated area, triggered by a precipitating systemic agent. This reaction typically affects the skin, and radiation recall myositis in the absence of cutaneous involvement has rarely been described in the literature. In this report, we present two cases of radiation recall in pediatric Ewing sarcoma patients receiving successive proton radiotherapy and chemotherapy, with magnetic resonance imaging (MRI) of muscle edema within the prior radiation fields.

Immobilization precision of a modified GTC frame.

The purpose of this study was to evaluate and quantify the interfraction reproducibility and intrafraction immobilization precision of a modified GTC frame. The error of the patient alignment and imaging systems were measured using a cranial skull phantom, with simulated, predetermined shifts. The kV setup images were acquired with a room-mounted set of kV sources and panels. Calculated translations and rotations provided by the computer alignment software relying upon three implanted fiducials were compared to the known shifts, and the accuracy of the imaging and positioning systems was calculated. Orthogonal kV setup images for 45 proton SRT patients and 1002 fractions (average 22.3 fractions/patient) were analyzed for interfraction and intrafraction immobilization precision using a modified GTC frame. The modified frame employs a radiotransparent carbon cup and molded pillow to allow for more treatment angles from posterior directions for cranial lesions. Patients and the phantom were aligned with three 1.5 mm stainless steel fiducials implanted into the skull. The accuracy and variance of the patient positioning and imaging systems were measured to be 0.10 ± 0.06 mm, with the maximum uncertainty of rotation being ±0.07°. 957 pairs of interfraction image sets and 974 intrafraction image sets were analyzed. 3D translations and rotations were recorded. The 3D vector interfraction setup reproducibility was 0.13 mm ± 1.8 mm for translations and the largest uncertainty of ± 1.07º for rotations. The intrafraction immobilization efficacy was 0.19 mm ± 0.66 mm for translations and the largest uncertainty of ± 0.50º for rotations. The modified GTC frame provides reproducible setup and effective intrafraction immobilization, while allowing for the complete range of entrance angles from the posterior direction.

Using four-dimensional computed tomography images to optimize the internal target volume when using volume-modulated arc therapy to treat moving targets.

In this work we used 4D dose calculations, which include the effects of shape deformations, to investigate an alternative approach to creating the ITV. We hypothesized that instead of needing images from all the breathing phases in the 4D CT dataset to create the outer envelope used for treatment planning, it is possible to exclude images from the phases closest to the inhale phase. We used 4D CT images from 10 patients with lung cancer. For each patient, we drew a gross tumor volume on the exhale-phase image and propagated this to the images from other phases in the 4D CT dataset using commercial image registration software. We created four different ITVs using the N phases closest to the exhale phase (where N = 10, 8, 7, 6). For each ITV contour, we created a volume-modulated arc therapy plan on the exhale-phase CT and normalized it so that the prescribed dose covered at least 95% of the ITV. Each plan was applied to CT images from each CT phase (phases 1-10), and the calculated doses were then mapped to the exhale phase using deformable registration. The effect of the motion was quantified using the dose to 95% of the target on the exhale phase (D95) and tumor control probability. For the three-dimensional and 4D dose calculations of the plan where N = 10, differences in the D95 value varied from 3% to 14%, with an average difference of 7%. For 9 of the 10 patients, the reduction in D95 was less than 5% if eight phases were used to create the ITV. For three of the 10 patients, the reduction in the D95 was less than 5% if seven phases were used to create the ITV. We were unsuccessful in creating a general rule that could be used to create the ITV. Some reduction (8/10 phases) was possible for most, but not all, of the patients, and the ITV reduction was small.

Pencil beam scanning dose calibration at reduced source-to-axis distance.

BACKGROUND: Pencil beam scanning (PBS) monitoring chambers use an ionization control signal, monitor units (MUs), or gigaprotons (Gp) to irradiate a pencil beam and normalize dose calculations. The nozzle deflects the beam from the nozzle axis by an angle subtended at the source-to-axis distance (τ) from the isocenter. If the angle is not correctly considered in calibrations or calculations, it can lead to systematic errors. PURPOSE: Aspects to consider for machines of various τs are fourfold. First, for the machine, there is a pathlength change of proton tracks in the monitor chamber. Second, for measurements, a uniform-square irradiation over a plane, with constant Gp per spot, does not deliver uniform dose in a measurement plane. Third, for Monte Carlo (MC) simulations, Gp (and not MU) is proportional to simulating a number of protons. Fourth, for pencil beam algorithms (PBA), MU or Gp may be used for pencil beam weight, but usage needs to be consistent. Another consideration is the beam shape change from circular to oval in the projection onto voxels. METHODS: Coordinate systems for PBS delivery are described. RESULTS: Users of intermediate-τ machines, corresponding to the onset of 1% pathlength corrections within the scanned field size, must not assume that MUs are proportional to the number of particles in MC simulations, and the PBA may need pathlength corrections. For a field size of 24 × 24 cm2 , intermediate-τ machines correspond to 59 cm ≤ τ < 120 cm. For a field size of 40 × 40 cm2 , intermediate-τ machines correspond to 98 cm ≤ τ < 200 cm. Small-τ machines correspond to τ < 59 and 98 cm at these field sizes, respectively, which require corrections in projecting the beam shape onto voxels. CONCLUSIONS: Identifying corrections due to the pencil beam angle and their onset are important for reducing the outer diameter of proton therapy gantries. The use of Gp (or the number of protons) meterset standardizes data interchange and helps to reduce systematic errors due to the angle of the beam.

MOQUI: an open-source GPU-based Monte Carlo code for proton dose calculation with efficient data structure.

Objective.Monte Carlo (MC) codes are increasingly used for accurate radiotherapy dose calculation. In proton therapy, the accuracy of the dose calculation algorithm is expected to have a more significant impact than in photon therapy due to the depth-dose characteristics of proton beams. However, MC simulations come at a considerable computational cost to achieve statistically sufficient accuracy. There have been efforts to improve computational efficiency while maintaining sufficient accuracy. Among those, parallelizing particle transportation using graphic processing units (GPU) achieved significant improvements. Contrary to the central processing unit, a GPU has limited memory capacity and is not expandable. It is therefore challenging to score quantities with large dimensions requiring extensive memory. The objective of this study is to develop an open-source GPU-based MC package capable of scoring those quantities.Approach.We employed a hash-table, one of the key-value pair data structures, to efficiently utilize the limited memory of the GPU and score the quantities requiring a large amount of memory. With the hash table, only voxels interacting with particles will occupy memory, and we can search the data efficiently to determine their address. The hash-table was integrated with a novel GPU-based MC code, moqui.Main results.The developed code was validated against an MC code widely used in proton therapy, TOPAS, with homogeneous and heterogeneous phantoms. We also compared the dose calculation results of clinical treatment plans. The developed code agreed with TOPAS within 2%, except for the fall-off and regions, and the gamma pass rates of the results were >99% for all cases with a 2 mm/2% criteria.Significance.We can score dose-influence matrix and dose-rate on a GPU for a 3-field H&N case with 10 GB of memory using moqui, which would require more than 100 GB of memory with the conventionally used array data structure.

Integrating Structure Propagation Uncertainties in the Optimization of Online Adaptive Proton Therapy Plans.

Currently, adaptive strategies require time- and resource-intensive manual structure corrections. This study compares different strategies: optimization without manual structure correction, adaptation with physician-drawn structures, and no adaptation. Strategies were compared for 16 patients with pancreas, liver, and head and neck (HN) cancer with 1-5 repeated images during treatment: 'reference adaptation', with structures drawn by a physician; 'single-DIR adaptation', using a single set of deformably propagated structures; 'multi-DIR adaptation', using robust planning with multiple deformed structure sets; 'conservative adaptation', using the intersection and union of all deformed structures; 'probabilistic adaptation', using the probability of a voxel belonging to the structure in the optimization weight; and 'no adaptation'. Plans were evaluated using reference structures and compared using a scoring system. The reference adaptation with physician-drawn structures performed best, and no adaptation performed the worst. For pancreas and liver patients, adaptation with a single DIR improved the plan quality over no adaptation. For HN patients, integrating structure uncertainties brought an additional benefit. If resources for manual structure corrections would prevent online adaptation, manual correction could be replaced by a fast 'plausibility check', and plans could be adapted with correction-free adaptation strategies. Including structure uncertainties in the optimization has the potential to make online adaptation more automatable.

Low-Dose Computed Tomography Scanning Protocols for Online Adaptive Proton Therapy of Head-and-Neck Cancers.

PURPOSE: To evaluate the suitability of low-dose CT protocols for online plan adaptation of head-and-neck patients. METHODS: We acquired CT scans of a head phantom with protocols corresponding to CT dose index volume CTDIvol in the range of 4.2-165.9 mGy. The highest value corresponds to the standard protocol used for CT simulations of 10 head-and-neck patients included in the study. The minimum value corresponds to the lowest achievable tube current of the GE Discovery RT scanner used for the study. For each patient and each low-dose protocol, the noise relative to the standard protocol, derived from phantom images, was applied to a virtual CT (vCT). The vCT was obtained from a daily CBCT scan corresponding to the fraction with the largest anatomical changes. We ran an established adaptive workflow twice for each low-dose protocol using a high-quality daily vCT and the corresponding low-dose synthetic vCT. For a relative comparison of the adaptation efficacy, two adapted plans were recalculated in the high-quality vCT and evaluated with the contours obtained through deformable registration of the planning CT. We also evaluated the accuracy of dose calculation in low-dose CT volumes using the standard CT protocol as reference. RESULTS: The maximum differences in D98 between low-dose protocols and the standard protocol for the high-risk and low-risk CTV were found to be 0.6% and 0.3%, respectively. The difference in OAR sparing was up to 3%. The Dice similarity coefficient between propagated contours obtained with low-dose and standard protocols was above 0.982. The mean 2%/2 mm gamma pass rate for the lowest-dose image, using the standard protocol as reference, was found to be 99.99%. CONCLUSION: The differences between low-dose protocols and the standard scanning protocol were marginal. Thus, low-dose CT protocols are suitable for online adaptive proton therapy of head-and-neck cancers. As such, considering scanning protocols used in our clinic, the imaging dose associated with online adaption of head-and-neck cancers treated with protons can be reduced by a factor of 40.

A fast double template convolution isocenter evaluation algorithm with subpixel accuracy.

PURPOSE: To design a fast Winston Lutz (fWL) algorithm for accurate analysis of radiation isocenter from images without edge detection or center of mass calculations. METHODS: An algorithm has been developed to implement the Winston Lutz test for mechanical/ radiation isocenter agreement using an electronic portal imaging device (EPID). The algorithm detects the position of the radiation shadow of a tungsten ball within a stereotactic cone. The fWL algorithm employs a double convolution to independently find the position of the sphere and cone centers. Subpixel estimation is used to achieve high accuracy. Results of the algorithm were compared to (1) a human observer with template guidance and (2) an edge detection/center of mass (edCOM) algorithm. Testing was performed with high resolution (0.05 mm/px, film) and low resolution (0.78 mm/px, EPID) image sets. RESULTS: Sphere and cone center relative positions were calculated with the fWL algorithm for high resolution test images with an accuracy of 0.002 +/- 0.061 mm compared to 0.042 +/- 0.294 mm for the human observer, and 0.003 +/- 0.038 mm for the edCOM algorithm. The fWL algorithm required 0.01 s per image compared to 5 s for the edCOM algorithm and 20 s for the human observer. For lower resolution images the fWL algorithm localized the centers with an accuracy of 0.083 +/- 0.12 mm compared to 0.03 +/- 0.5514 mm for the edCOM algorithm. CONCLUSIONS: A fast (subsecond) subpixel algorithm has been developed that can accurately determine the center locations of the ball and cone in Winston Lutz test images without edge detection or COM calculations.

Effect of respiratory trace shape on optimal treatment margin.

PURPOSE: To evaluate the effect of target trajectory shape on the optimal treatment margin. METHODS: Intensity-modulated radiation therapy and volumetric modulated arc therapy plans were created for three spherical targets (3, 5, and 7 cm diameter) simulated in exhalation phases, each with margins of 2, 4, 6, 8, and 10 mm to account for motion. The plans were delivered to a stationary 2D ion chamber array, and dose movies were obtained of the delivered doses. The dose movie frames were then displaced to simulate different respiratory traces. Five traces were used: sin2, sin4 sin6, and two patient traces. The optimal margin was defined as the margin for which the dose delivered to 95% of the target was closest to that obtained with no margin or motion. The equivalent uniform dose was also investigated as an alternative cost function. RESULTS: The optimal margin was always smaller than the peak-to-peak motion. When the respiratory trace spent less time in the inhale phases, the optimal margin was consistently smaller than when more time was spent in the inhale phases. The target size and treatment modality also affected the optimal margin. CONCLUSIONS: The necessary margin for targets that spend less time in the exhale phase (sin6) is 2-4 mm smaller than for targets that spend equal time in the inhale and exhale phases (sin).

Adaptive proton therapy.

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

CT-on-Rails Versus In-Room CBCT for Online Daily Adaptive Proton Therapy of Head-and-Neck Cancers.

PURPOSE: To compare the efficacy of CT-on-rails versus in-room CBCT for daily adaptive proton therapy. METHODS: We analyzed a cohort of ten head-and-neck patients with daily CBCT and corresponding virtual CT images. The necessity of moving the patient after a CT scan is the most significant difference in the adaptation workflow, leading to an increased treatment execution uncertainty σ. It is a combination of the isocenter-matching σi and random patient movements induced by the couch motion σm. The former is assumed to never exceed 1 mm. For the latter, we studied three different scenarios with σm = 1, 2, and 3 mm. Accordingly, to mimic the adaptation workflow with CT-on-rails, we introduced random offsets after Monte-Carlo-based adaptation but before delivery of the adapted plan. RESULTS: There were no significant differences in accumulated dose-volume histograms and dose distributions for σm = 1 and 2 mm. Offsets with σm = 3 mm resulted in underdosage to CTV and hot spots of considerable volume. CONCLUSION: Since σm typically does not exceed 2 mm for in-room CT, there is no clinically significant dosimetric difference between the two modalities for online adaptive therapy of head-and-neck patients. Therefore, in-room CT-on-rails can be considered a good alternative to CBCT for adaptive proton therapy.