Pelvic organ motion and dosimetric implications during horizontal patient rotation for prostate radiation therapy

PURPOSE

Gantry-free radiation therapy systems utilizing patient rotation would be simpler and more cost effective than the conventional gantry-based systems. Such a system could enable the expansion of radiation therapy to meet global demand and reduce capital costs. Recent advances in adaptive radiation therapy could potentially be applied to correct for gravitational deformation during horizontal patient rotation. This study aims to quantify the pelvic organ motion and the dosimetric implications of horizontal rotation for prostate intensity-modulated radiation therapy (IMRT) treatments.

METHODS

Eight human participants who previously received prostate radiation therapy were imaged in a clinical magnetic resonance imaging (MRI) scanner using a bespoke patient rotation system (PRS). The patients were imaged every 45 degrees during a full roll rotation (0-360 degrees). Whole pelvic bone, prostate, rectum, and bladder motion were compared to the supine position using dice similarity coefficient (DSC) and mean absolute surface distance (MASD). Prostate centroid motion was compared in the left-right (LR), superior-inferior (SI), and anterior-posterior (AP) direction prior to and following pelvic bone-guided rigid registration. Seven-field prostate IMRT treatment plans were generated for each patient rotation angles under three adaption scenarios: No plan adaption, rigid planning target volume (PTV)-guided alignment to the prostate, and plan re-optimization. Prostate, rectum, and bladder doses were compared for each adaption scenario.

RESULTS

Pelvic bone motion within the PRS of up to 53 mm relative to the supine position was observed for some participants. Internal organ motion was greatest at the 180-degree PRS couch angle (prone), with prostate centroid motion range < 2 mm LR, 0 mm to 14 mm SI, and -11 mm to 4 mm AP. Rotation with no adaption of the treatment plan resulted in an underdose to the PTV -- in some instances up to 75% (D95%: 78 ± 0.3 Gy at supine to 20 ± 15.0 Gy at the 225-degree PRS couch angle). Bladder dose was reduced during the rotation by up to 98% (V60 Gy: 15.0 ± 9.4% supine to 0.3 ± 0.5% at the 225-degree PRS couch angle). In some instances, the rectum dose increased during rotation (V60Gy: 20.0 ± 4.5% supine to 25.0 ± 15.0% at the 135-degree PRS couch angle). Rigid PTV-guided alignment resulted in PTV coverage which, though statistically lower (P < 0.05 for all D95% values), was within 1 Gy of the supine plans. Plan re-optimization resulted in a statistically equivalent PTV coverage compared to the supine plans (P > 0.05 for all D95% metrics and all within ±0.4 Gy). For both rigid PTV-guided alignment and plan re-optimization, rectum dose volume metrics were reduced compared to the supine position between the 90- and 225-degree PRS couch angles (P < 0.05). Bladder dose volume metrics were not impacted by rotation.

CONCLUSION

Pelvic bone and internal organ motion are present during patient rotation. Rigid PTV-guided alignment to the prostate will be a requirement if prostate IMRT is to be safely delivered using patient rotation. Plan re-optimization for each PRS couch angle to account for anatomical deformations further improves the PTV coverage.

An interdimensional correlation framework for real-time estimation of six degree of freedom target motion using a single x-ray imager during radiotherapy

Increasing evidence suggests that intrafraction tumour motion monitoring needs to include both 3D translations and 3D rotations. Presently, methods to estimate the rotation motion require the 3D translation of the target to be known first. However, ideally, translation and rotation should be estimated concurrently. We present the first method to directly estimate six-degree-of-freedom (6DoF) motion from the target's projection on a single rotating x-ray imager in real-time. This novel method is based on the linear correlations between the superior-inferior translations and the motion in the other five degrees-of-freedom. The accuracy of the method was evaluated in silico with 81 liver tumour motion traces from 19 patients with three implanted markers. The ground-truth motion was estimated using the current gold standard method where each marker's 3D position was first estimated using a Gaussian probability method, and the 6DoF motion was then estimated from the 3D positions using an iterative method. The 3D position of each marker was projected onto a gantry-mounted imager with an imaging rate of 11 Hz. After an initial 110° gantry rotation (200 images), a correlation model between the superior-inferior translations and the five other DoFs was built using a least square method. The correlation model was then updated after each subsequent frame to estimate 6DoF motion in real-time. The proposed algorithm had an accuracy (±precision) of  -0.03  ±  0.32 mm, -0.01  ±  0.13 mm and 0.03  ±  0.52 mm for translations in the left-right (LR), superior-inferior (SI) and anterior-posterior (AP) directions respectively; and, 0.07  ±  1.18°, 0.07  ±  1.00° and 0.06  ±  1.32° for rotations around the LR, SI and AP axes respectively on the dataset. The first method to directly estimate real-time 6DoF target motion from segmented marker positions on a 2D imager was devised. The algorithm was evaluated using 81 motion traces from 19 liver patients and was found to have sub-mm and sub-degree accuracy.

Radiotherapy beyond cancer: target localization in real-time MRI and treatment planning for cardiac radiosurgery

PURPOSE

Atrial fibrillation (AFib) is the most common cardiac arrhythmia that affects millions of patients world-wide. AFib is usually treated with minimally invasive, time consuming catheter ablation techniques. While recently noninvasive radiosurgery to the pulmonary vein antrum (PVA) in the left atrium has been proposed for AFib treatment, precise target location during treatment is challenging due to complex respiratory and cardiac motion. A MRI linear accelerator (MRI-Linac) could solve the problems of motion tracking and compensation using real-time image guidance. In this study, the authors quantified target motion ranges on cardiac magnetic resonance imaging (MRI) and analyzed the dosimetric benefits of margin reduction assuming real-time motion compensation was applied.

METHODS

For the imaging study, six human subjects underwent real-time cardiac MRI under free breathing. The target motion was analyzed retrospectively using a template matching algorithm. The planning study was conducted on a CT of an AFib patient with a centrally located esophagus undergoing catheter ablation, representing an ideal case for cardiac radiosurgery. The target definition was similar to the ablation lesions at the PVA created during catheter treatment. Safety margins of 0 mm (perfect tracking) to 8 mm (untracked respiratory motion) were added to the target, defining the planning target volume (PTV). For each margin, a 30 Gy single fraction IMRT plan was generated. Additionally, the influence of 1 and 3 T magnetic fields on the treatment beam delivery was simulated using Monte Carlo calculations to determine the dosimetric impact of MRI guidance for two different Linac positions.

RESULTS

Real-time cardiac MRI showed mean respiratory target motion of 10.2 mm (superior-inferior), 2.4 mm (anterior-posterior), and 2 mm (left-right). The planning study showed that increasing safety margins to encompass untracked respiratory motion leads to overlapping structures even in the ideal scenario, compromising either normal tissue dose constraints or PTV coverage. The magnetic field caused a slight increase in the PTV dose with the in-line MRI-Linac configuration.

CONCLUSIONS

The authors' results indicate that real-time tracking and motion compensation are mandatory for cardiac radiosurgery and MRI-guidance is feasible, opening the possibility of treating cardiac arrhythmia patients completely noninvasively.

Quality assurance for the clinical implementation of kilovoltage intrafraction monitoring for prostate cancer VMAT

PURPOSE

Kilovoltage intrafraction monitoring (KIM) is a real-time 3D tumor monitoring system for cancer radiotherapy. KIM uses the commonly available gantry-mounted x-ray imager as input, making this method potentially more widely available than dedicated real-time 3D tumor monitoring systems. KIM is being piloted in a clinical trial for prostate cancer patients treated with VMAT (NCT01742403). The purpose of this work was to develop clinical process and quality assurance (QA) practices for the clinical implementation of KIM.

METHODS

Informed by and adapting existing guideline documents from other real-time monitoring systems, KIM-specific QA practices were developed. The following five KIM-specific QA tests were included: (1) static localization accuracy, (2) dynamic localization accuracy, (3) treatment interruption accuracy, (4) latency measurement, and (5) clinical conditions accuracy. Tests (1)-(4) were performed using KIM to measure static and representative patient-derived prostate motion trajectories using a 3D programmable motion stage supporting an anthropomorphic phantom with implanted gold markers to represent the clinical treatment scenario. The threshold for system tolerable latency is <1 s. The tolerances for all other tests are that both the mean and standard deviation of the difference between the programmed trajectory and the measured data are <1 mm. The (5) clinical conditions accuracy test compared the KIM measured positions with those measured by kV/megavoltage (MV) triangulation from five treatment fractions acquired in a previous pilot study.

RESULTS

For the (1) static localization, (2) dynamic localization, and (3) treatment interruption accuracy tests, the mean and standard deviation of the difference are <1.0 mm. (4) The measured latency is 350 ms. (5) For the tests with previously acquired patient data, the mean and standard deviation of the difference between KIM and kV/MV triangulation are <1.0 mm.

CONCLUSIONS

Clinical process and QA practices for the safe clinical implementation of KIM, a novel real-time monitoring system using commonly available equipment, have been developed and implemented for prostate cancer VMAT.

Determining appropriate imaging parameters for kilovoltage intrafraction monitoring: an experimental phantom study

Kilovoltage intrafraction monitoring (KIM) utilises the kV imager during treatment for real-time tracking of prostate fiducial markers. However, its effectiveness relies on sufficient image quality for the fiducial tracking task. To guide the performance characterisation of KIM under different clinically relevant conditions, the effect of different kV parameters and patient size on image quality, and quantification of MV scatter from the patient to the kV detector panel were investigated in this study. Image quality was determined for a range of kV acquisition frame rates, kV exposure, MV dose rates and patient sizes. Two methods were used to determine image quality; the ratio of kV signal through the patient to the MV scatter from the patient incident on the kilovoltage detector, and the signal-to-noise ratio (SNR). The effect of patient size and frame rate on MV scatter was evaluated in a homogeneous CIRS pelvis phantom and marker segmentation was determined utilising the Rando phantom with embedded markers. MV scatter incident on the detector was shown to be dependent on patient thickness and frame rate. The segmentation code was shown to be successful for all frame rates above 3 Hz for the Rando phantom corresponding to a kV to MV ratio of 0.16 and an SNR of 1.67. For a maximum patient dimension less than 36.4 cm the conservative kV parameters of 5 Hz at 1 mAs can be used to reduce dose while retaining image quality, where the current baseline kV parameters of 10 Hz at 1 mAs is shown to be adequate for marker segmentation up to a patient dimension of 40 cm. In conclusion, the MV scatter component of image quality noise for KIM has been quantified. For most prostate patients, use of KIM with 10 Hz imaging at 1 mAs is adequate however image quality can be maintained and imaging dose reduced by altering existing acquisition parameters.

Estimation of effective imaging dose for kilovoltage intratreatment monitoring of the prostate position during cancer radiotherapy

Kilovoltage intratreatment monitoring (KIM) is a novel real-time localization modality where the tumor position is continuously measured during intensity modulated radiation therapy (IMRT) or intensity modulated arc therapy (IMAT) by a kilovoltage (kV) x-ray imager. Adding kV imaging during therapy adds radiation dose. The additional effective dose is quantified for prostate radiotherapy and compared to dose from other localization modalities. The software PCXMC 2.0 was used to calculate the effective dose delivered to a phantom as a function of imager angle and field size for a Varian On-Board Imager. The average angular effective dose was calculated for a field size of 6 cm × 6 cm. The average angular effective dose was used in calculations for different treatment scenarios. Treatment scenarios considered were treatment type and fractionation. For all treatment scenarios, (i.e. conventionally fractionated and stereotactic body radiotherapy (SBRT), IMRT and IMAT), the total KIM dose at 1 Hz ranged from 2-10 mSv. This imaging dose is less than the Navotek radioactive implant dose (64 mSv) and a standard SBRT cone beam computed tomography pretreatment scan dose (22 mSv) over an entire treatment regime. KIM delivers an acceptably low effective dose for daily use as a real-time image-guidance method for prostate radiotherapy.

Electron contamination modeling and skin dose in 6 MV longitudinal field MRIgRT: Impact of the MRI and MRI fringe field

PURPOSE

In recent times, longitudinal field MRI-linac systems have been proposed for 6 MV MRI-guided radiotherapy (MRIgRT). The magnetic field is parallel with the beam axis and so will alter the transport properties of any electron contamination particles. The purpose of this work is to provide a first investigation into the potential effects of the MR and fringe magnetic fields on the electron contamination as it is transported toward a phantom, in turn, providing an estimate of the expected patient skin dose changes in such a modality.

METHODS

Geant4 Monte Carlo simulations of a water phantom exposed to a 6 MV x-ray beam were performed. Longitudinal magnetic fields of strengths between 0 and 3 T were applied to a 30 × 30 × 20 cm(3) phantom. Surrounding the phantom there is a region where the magnetic field is at full MRI strength, consistent with clinical MRI systems. Beyond this the fringe magnetic field entering the collimation system is also modeled. The MRI-coil thickness, fringe field properties, and isocentric distance are varied and investigated. Beam field sizes of 5 × 5, 10 × 10, 15 × 15 and 20 × 20 cm(2) were simulated. Central axis dose, 2D virtual entry skin dose films, and 70 μm skin depth doses were calculated using high resolution scoring voxels.

RESULTS

In the presence of a longitudinal magnetic field, electron contamination from the linear accelerator is encouraged to travel almost directly toward the patient surface with minimal lateral spread. This results in a concentration of electron contamination within the x-ray beam outline. This concentration is particularly encouraged if the fringe field encompasses the collimation system. Skin dose increases of up to 1000% were observed for certain configurations and increases above Dmax were common. In nonmagnetically shielded cases, electron contamination generated from the jaw faces and air column is trapped and propagated almost directly to the phantom entry region, giving rise to intense dose hot spots inside the x-ray treatment field. These range up to 1000% or more of Dmax at the CAX, depending on field size, isocenter, and coil thickness. In the case of a fully magnetically shielded collimation system and the lowest MRI field of 0.25 T, the entry skin dose is expected to increase to at least 40%, 50%, 65%, and 80% of Dmax for 5 × 5, 10 × 10, 15 × 15, and 20 × 20 cm(2), respectively.

CONCLUSIONS

Electron contamination from the linac head and air column may cause considerable skin dose increases or hot spots at the beam central axis on the entry side of a phantom or patient in longitudinal field 6 MV MRIgRT. This depends heavily on the properties of the magnetic fringe field entering the linac beam collimation system. The skin dose increase is also related to the MRI-coil thickness, the fringe field, and the isocenter distance of the linac. The results of this work indicate that the properties of the MRI fringe field, electron contamination production, and transport must be considered carefully during the design stage of a longitudinal MRI-linac system.

Accuracy in the localization of thoracic and abdominal tumors using respiratory displacement, velocity, and phase

Current four dimensional (4D) computed tomography (CT) reconstruction techniques are retrospectively created based on either the phase or displacement of the respiratory signal. Both techniques have known limitations which cause clinically significant motion artifacts in 4D CT images. These artifacts, which appear as undefined or irregular boundaries in the 4D CT images, cause systematic errors in patient contouring and dose calculations. The aim of this work was to evaluate the reproducibility of tumor position as a function of displacement, phase, and velocity of the respiratory signal, respectively, in order to determine the optimum parameter or combination of parameters to use in order to minimize artifacts in 4D CT images or to accurately deliver radiation to relevant structures during treatment.

METHOD AND MATERIALS

Estimated tumor centroid position and respiratory signal data were acquired with the Cyberknife Synchrony system for 26 thoracic radiotherapy patients (52 fractions). A reference respiratory cycle was calculated for each patient. Displacement, phase, and velocity of ten data points were calculated from this reference respiratory cycle, where each point represents an image bin. The corresponding tumor position was then sorted into these image bins if the phase, displacement, simultaneous displacement and phase, or simultaneous displacement and velocity of the respiratory signal were within tolerances of 0.5 mm for displacement and 0.5 mm/s for velocity, respectively, from the corresponding data of the reference cycle for each image bin.

RESULTS

The mean of the standard deviations of tumor positions over all bins and all fractions for the superior-inferior direction were 2.13 +/- 1.01 mm for phase sorting, 1.20 +/- 0.76 mm for displacement sorting, 1.20 +/- 0.71 mm for simultaneous displacement and phase sorting, and 1.10 +/- 0.71 mm for simultaneous displacement and velocity sorting, with maximum deviations of 43.0, 16.1, 15.5, and 14.1 mm for each scenario, respectively. The same trend was observed for the anterior-posterior and left-right directions. A linear dependence was observed between the mean of the standard deviations of tumor positions over all fractions as a function of the velocity of the respiratory signal at each bin for all the sorting scenarios. A substantially larger gradient for the phase sorting scenario, compared to the other scenarios, suggests that tumor localization will become increasingly less accurate as the velocity of the tumor increases during a breathing cycle, e.g., if the amplitude of motion increases while the period of the respiratory cycle stays constant or during mid inhale or exhale phases of the respiratory cycle.

CONCLUSION

This study illustrates that position of a tumor can be determined more accurately if displacement and velocity are used simultaneously as sorting parameters for 4D CT images or during treatment. A real-time displacement and velocity based 4D CT image sorting method may therefore produce fewer and smaller artifacts in 4D CT images than current retrospective sorting methods.

Prospective displacement and velocity-based cine 4D CT

Four dimensional (4D) computed tomography (CT) image sorting is currently a retrospective procedure. Mismatches in displacement and/or phase of a patient's respiratory signal, corresponding with two dimensional images taken at subsequent couch positions, become visible as artifacts in reconstructed 4D CT images. These artifacts appear as undefined or irregular boundaries in the 4D CT images and cause systematic errors in patient contouring and dose calculations. In addition, the substantially higher dose required for 4D CT, compared with 3D CT, is of concern. To minimize these problems, we developed a prospective respiratory displacement and velocity based cine 4D CT (PDV CT) method to trigger image acquisition if the displacement and velocity of the respiratory signal occurred within predetermined tolerances simultaneously. The use of velocity avoids real-time phase estimation. Real-time image acquisition ensures data sufficiency, while avoiding the need for redundant data. This may potentially result in a lower dose to the patient. PDV CT was compared with retrospective 4D CT acquisition methods, using respiratory signals of 24 lung cancer patients (103 sessions) under free breathing conditions. Image acquisition was simulated for each of these sessions from the respiratory signal. The root mean square (RMS) of differences between displacements and velocities of the respiratory signal corresponding to subsequent images was calculated in order to evaluate the image-sorting accuracy of each method. Patient dose reductions of 22 to 50% were achieved during image acquisition depending on the model parameters chosen. The mean RMS differences over all sessions and image bins show that PDV CT produces similar results to retrospective displacement sorting overall, although improvements of the RMS difference up to 20% were achieved depending on the model parameters chosen. Velocity RMS differences improved between 30 and 45% when compared with retrospective phase sorting. The efficiency in acquisition compared with retrospective phase sorting varied from approximately 10% for displacement and velocity tolerances of 1 mm and 4 mm/s, respectively, to 80 to 93% for 4 mm and 4 mm/s. The lower variation in the displacement and velocity of the respiratory signal in each image bin indicates that PDV CT could be a valuable tool for reducing artifacts in 4D CT images and lowering patient dose, although the cost may be increased acquisition time.

Quantifying the accuracy of automated structure segmentation in 4D CT images using a deformable image registration algorithm

Four-dimensional (4D) radiotherapy is the explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy. One key component of 4D radiotherapy planning is the ability to automatically ("auto") create contours on all of the respiratory phase computed tomography (CT) datasets comprising a 4D CT scan, based on contours manually drawn on one CT image set from one phase. A tool that can be used to automatically propagate manually drawn contours to CT scans of other respiratory phases is deformable image registration. The purpose of the current study was to geometrically quantify the difference between automatically generated contours with manually drawn contours. Four-DCT data sets of 13 patients consisting of ten three-dimensional CT image sets acquired at different respiratory phases were used for this study. Tumor and normal tissue structures [gross tumor volume (GTV), esophagus, right lung, left lung, heart and cord] were manually drawn on each respiratory phase of each patient. Large deformable diffeomorphic image registration was performed to map each CT set from the peak-inhale respiration phase to the CT image sets corresponding with subsequent respiration phases. The calculated displacement vector fields were used to deform contours automatically drawn on the inhale phase to the other respiratory phase CT image sets. The code was interfaced to a treatment planning system to view the resulting images and to obtain the volumetric, displacement, and surface congruence information; 692 automatically generated structures were compared with 692 manually drawn structures. The auto- and manual methods showed similar trends, with a smaller difference observed between the GTVs than other structures. The auto-contoured structures agree with the manually drawn structures, especially in the case of the GTV, to within published interobserver variations. For the GTV, fractional volumes agree to within 0.2+/-0.1, center of mass displacements agree to within 0.5+/-1.5 mm, and agreement of surface congruence is 0.0+/-1.1 mm. The surface congruence between automatic and manual contours for the GTV, heart, left lung, right lung and esophagus was less than 5 mm in 99%, 94%, 94%, 91% and 89%, respectively. Careful assessment of the performance of automatic algorithms is needed in the presence of 4D CT artifacts.

A strategy to correct for intrafraction target translation in conformal prostate radiotherapy: simulation results

A strategy is proposed in which intrafraction internal target translation is corrected for by repositioning the multileaf collimator position aperture to conform to the new target pose in the beam projection, and the beam monitor units are adjusted to account for the change in the geometric relationship between the target and the beam. The purpose of this study was to investigate the dosimetric stability of the prostate and critical structures in the presence of internal target translation using the dynamic compensation strategy. Twenty-five previously treated prostate cancer patients were replanned using a four-field conformal technique to deliver 72 Gy to 95% of the planning target volume (PTV). Internal translation was introduced by displacing the prostate PTV (no rotation or deformation was considered). Thirty-six randomly selected isotropic displacements of magnitude 0.5, 1.0, 1.5 and 2.0 cm were sampled for each patient, for a total of 3600 errors. Due to their anatomic relation to the prostate, the rectum and bladder contours were also moved with the same magnitude and direction as the prostate. The dynamic compensation strategy was used to correct each of these errors by conforming the beam apertures to the new target pose and adjusting the monitor units using inverse-square and off-axis factor corrections. The dynamic compensation strategy plans were then compared to the original treatment plans via dose-volume histogram (DVH) analysis. Changes of more than 5% of the prescription dose (3.6 Gy) were deemed clinically significant. Compared to the original treatment plans, the dynamic compensation strategy produced small discrepancies in isodose distributions and DVH analyses for all structures considered apart from the femoral heads. These differences increased with the magnitude of the internal motion. Coverage of the PTV was excellent: D5, D95, and Dmean were not increased or decreased by more than 5% of the prescription dose for any of the 3600 simulated internal motion shifts. Dose increases for adjacent organs at risk were rare. D33 of the rectum and D20 of the bladder were increased by more than 5% of the prescription dose in 9 and 1 instances of the 3600 sampled internal motion shifts, respectively. Dmean of the right femoral head increased by more than 5% of the prescription dose in 651 (18%) internal motion shifts, predominantly due to the projection of the lateral beams through the femoral head for anterior prostate motion. However, D2 was not increased by more than 5% for any of the internal motion shifts. These data demonstrate the robustness of the proposed dynamic compensation strategy for correction of internal motion in conformal prostate radiotherapy, with minimal deviation from the original treatment plans even for errors exceeding those commonly encountered in the clinic. The compensation strategy could be performed automatically with appropriate enhancements to available delivery software.

Optimum parameters in a model for tumour control probability, including interpatient heterogeneity: evaluation of the log-normal distribution

The heterogeneity of human tumour radiation response is well known. Researchers have used the normal distribution to describe interpatient tumour radiosensitivity. However, many natural phenomena show a log-normal distribution. Log-normal distributions are common when mean values are low, variances are large and values cannot be negative. These conditions apply to radiosensitivity. The aim of this work was to evaluate the log-normal distribution to predict clinical tumour control probability (TCP) data and to compare the results with the homogeneous (delta-function with single alpha-value) and normal distributions. The clinically derived TCP data for four tumour types-melanoma, breast, squamous cell carcinoma and nodes-were used to fit the TCP models. Three forms of interpatient tumour radiosensitivity were considered: the log-normal, normal and delta-function. The free parameters in the models were the radiosensitivity mean, standard deviation and clonogenic cell density. The evaluation metric was the deviance of the maximum likelihood estimation of the fit of the TCP calculated using the predicted parameters to the clinical data. We conclude that (1) the log-normal and normal distributions of interpatient tumour radiosensitivity heterogeneity more closely describe clinical TCP data than a single radiosensitivity value and (2) the log-normal distribution has some theoretical and practical advantages over the normal distribution. Further work is needed to test these models on higher quality clinical outcome datasets.

Investigation of patient, tumour and treatment variables affecting residual motion for respiratory-gated radiotherapy

Respiratory gating can reduce the apparent respiratory motion during imaging and treatment; however, residual motion within the gating window remains. Respiratory training can improve respiratory reproducibility and, therefore, the efficacy of respiratory-gated radiotherapy. This study was conducted to determine whether residual motion during respiratory gating is affected by patient, tumour or treatment characteristics. The specific aims of this study were to: (1) identify significant characteristics affecting residual motion, (2) investigate time trends of residual motion over a period of days (inter-session) and (3) investigate time trends of residual motion within the same day (intra-session). Twenty-four lung cancer patients were enrolled in an Institutional Review Board (IRB)-approved protocol. For approximately five sessions, 331 four-minute, respiratory motion traces were acquired with free breathing, audio instructions and audio-visual biofeedback for each patient. The residual motion was quantified by the standard deviation of the displacement within the gating window. The generalized linear model was used to obtain coefficients for each variable within the model and to evaluate the clinical and statistical significance. The statistical significance was determined by a p-value<0.05, while effect sizes of 0.1 cm (one standard deviation) were considered clinically significant. This data analysis was applied to patient, tumour and treatment variables. Inter- and intra-session variations were also investigated. The only variable that was significant for both inhale- and exhale-based gating was disease type. In addition, visual-training displacement, breathing type and Karnofsky performance status (KPS) values were significant for inhale-based gating, and dose-per-fraction was significant for exhale-based gating. Temporal respiratory variations within and between sessions were observed for individual patients. However inter- and intra-session analyses did not show significant time trends on average for any of the variables considered. The lack of significant time trends within and between sessions indicates that on average (1) there is no significant learning period for breathing training, (2) the patients did not experience training-related fatigue and (3) the margin component to account for residual motion during gated radiotherapy appears to remain constant throughout the treatment.

Clinical comparison of head and neck and prostate IMRT plans using absorbed dose to medium and absorbed dose to water

Conventional photon radiation therapy dose-calculation algorithms typically compute and report the absorbed dose to water (D(w)). Monte Carlo (MC) dose-calculation algorithms, however, generally compute and report the absorbed dose to the material (D(m)). As MC-calculation algorithms are being introduced into routine clinical usage, the question as to whether there is a clinically significant difference between D(w) and D(m) remains. The goal of the current study is to assess the differences between dose-volume indices for D(m) and D(w) MC-calculated IMRT plans. Ten head-and-neck (H&N) and ten prostate cancer patients were selected for this study. MC calculations were performed using an EGS4-based system. Converting D(m) to D(w) for MC-based calculations was accomplished as a post-MC calculation process. D(w) and D(m) results for target and critical structures were evaluated using the dose-volume-based indices. For H&N IMRT plans, systematic differences between dose-volume indices computed with D(w) and D(m) were up to 2.9% for the PTV prescription dose (D(98)), up to 5.8% for maximum (D(2)) dose to the PTV and up to 2.7% for the critical structure dose indices. For prostate IMRT plans, the systematic differences between D(w)- and D(m)-based computed indices were up to 3.5% for the prescription dose (D(98)) to the PTVs, up to 2.0% for the maximum (D(2)) dose to the PTVs and up to 8% for the femoral heads due to their higher water/bone mass stopping power ratio. This study showed that converting D(m) to D(w) in MC-calculated IMRT treatment plans introduces a systematic error in target and critical structure DVHs. In some cases, this systematic error may reach up to 5.8% for H&N and 8.0% for prostate cases when the hard-bone-containing structures such as femoral heads are present. Ignoring differences between D(m) and D(w) will result in systematic dose errors ranging from 0% to 8%.

The application of the sinusoidal model to lung cancer patient respiratory motion

Accurate modeling of the respiratory cycle is important to account for the effect of organ motion on dose calculation for lung cancer patients. The aim of this study is to evaluate the accuracy of a respiratory model for lung cancer patients. Lujan et al. [Med. Phys. 26(5), 715-720 (1999)] proposed a model, which became widely used, to describe organ motion due to respiration. This model assumes that the parameters do not vary between and within breathing cycles. In this study, first, the correlation of respiratory motion traces with the model f(t) as a function of the parameter n (n = 1, 2, 3) was undertaken for each breathing cycle from 331 four-minute respiratory traces acquired from 24 lung cancer patients using three breathing types: free breathing, audio instruction, and audio-visual biofeedback. Because cos2 and cos4 had similar correlation coefficients, and cos2 and cos1 have a trigonometric relationship, for simplicity, the cos1 value was consequently used for further analysis in which the variations in mean position (z0), amplitude of motion (b) and period (tau) with and without biofeedback or instructions were investigated. For all breathing types, the parameter values, mean position (z0), amplitude of motion (b), and period (tau) exhibited significant cycle-to-cycle variations. Audio-visual biofeedback showed the least variations for all three parameters (z0, b, and tau). It was found that mean position (z0) could be approximated with a normal distribution, and the amplitude of motion (b) and period (tau) could be approximated with log normal distributions. The overall probability density function (pdf) of f(t) for each of the three breathing types was fitted with three models: normal, bimodal, and the pdf of a simple harmonic oscillator. It was found that the normal and the bimodal models represented the overall respiratory motion pdfs with correlation values from 0.95 to 0.99, whereas the range of the simple harmonic oscillator pdf correlation values was 0.71 to 0.81. This study demonstrates that the pdfs of mean position (z0), amplitude of motion (b), and period (tau) can be used for sampling to obtain more realistic respiratory traces. The overall standard deviations of respiratory motion were 0.48, 0.57, and 0.55 cm for free breathing, audio instruction, and audio-visual biofeedback, respectively.

Radiological properties of normoxic polymer gel dosimeters

The radiological properties of the normoxic polymer gel dosimeters MAGIC, MAGAS, and MAGAT [methacrylic and ascorbic acid in gelatin initiated by copper; methacrylic acid gelatine gel with ascorbic acid; and methacrylic acid gelatine and tetrakis (hydroxymethyl) phosphonium chloride, respectively] have been investigated. The radiological water equivalence was determined by comparing the polymer gel macroscopic photon and electron interaction cross sections over the energy range from 10 keV to 20 MeV and by Monte Carlo modeling of depth doses. Normoxic polymer gel dosimeters have a high gelatine and monomer concentration and therefore mass density (kg m(-3)) up to 3.8% higher than water. This results in differences between the cross-section ratios of the normoxic polymer gels and water of up to 3% for the attenuation, energy absorption, and collision stopping power coefficient ratios through the Compton dominant energy range. The mass cross-section ratios were within 2% of water except for the mass attenuation and energy absorption coefficients ratios, which showed differences with water of up to 6% for energies less than 100 keV. Monte Carlo modeling was undertaken for the polymer gel dosimeters to model the electron and photon transport resulting from a 6 MV photon beam. The absolute percentage differences between gel and water were within 1% and the relative percentage differences were within 3.5%. The results show that the MAGAT gel formulation is the most radiological water equivalent of the normoxic polymer gel dosimeters investigated due to its lower mass density measurement compared with MAGAS and MAGIC gels.

Determination of maximum leaf velocity and acceleration of a dynamic multileaf collimator: Implications for 4D radiotherapy

The dynamic multileaf collimator (MLC) can be used for four-dimensional (4D), or tumor tracking radiotherapy. However, the leaf velocity and acceleration limitations become a crucial factor as the MLC leaves need to respond in near real time to the incoming respiration signal. The aims of this paper are to measure maximum leaf velocity, acceleration, and deceleration to obtain the mechanical response times for the MLC, and determine whether the MLC is suitable for 4D radiotherapy. MLC leaf sequence files, requiring the leaves to reach maximum acceleration and velocity during motion, were written. The leaf positions were recorded every 50 ms, from which the maximum leaf velocity, acceleration, and deceleration were derived. The dependence on the velocity and acceleration of the following variables were studied: leaf banks, inner and outer leaves, MLC-MLC variations, gravity, friction, and the stability of measurements over time. Measurement results show that the two leaf banks of a MLC behave similarly, while the inner and outer leaves have significantly different maximum leaf velocities. The MLC-MLC variations and the dependence of gravity on maximum leaf velocity are statistically significant. The average maximum leaf velocity at the isocenter plane of the MLC ranged from 3.3 to 3.9 cm/s. The acceleration and deceleration at the isocenter plane of the MLC ranged from 50 to 69 cm/s2 and 46 to 52 cm/s2, respectively. Interleaf friction had a negligible effect on the results, and the MLC parameters remained stable with time. Equations of motion were derived to determine the ability of the MLC response to fluoroscopymeasured diaphragm motion. Given the present MLC mechanical characteristics, 4D radiotherapy is feasible for up to 97% of respiratory motion. For the largest respiratory motion velocities observed, beam delivery should be temporarily stopped (beam hold).

Predicting respiratory motion for four-dimensional radiotherapy

Adapting radiation delivery to respiratory motion is made possible through corrective action based on real-time feedback of target position during respiration. The advantage of this approach lies with its ability to allow tighter margins around the target while simultaneously following its motion. A significant hurdle to the successful implementation of real-time target-tracking-based radiation delivery is the existence of a finite time delay between the acquisition of target position and the mechanical response of the system to the change in position. Target motion during the time delay leads to a resultant lag in the system's response to a change in tumor position. Predicting target position in advance is one approach to ensure accurate delivery. The aim of this manuscript is to estimate the predictive ability of sinusoidal and adaptive filter-based prediction algorithms on multiple sessions of patient respiratory patterns. Respiratory motion information was obtained from recordings of diaphragm motion for five patients over 60 sessions. A prediction algorithm that employed both prediction models-the sinusoidal model and the adaptive filter model-was developed to estimate prediction accuracy over all the sessions. For each session, prediction error was computed for several time instants (response time) in the future (0-1.8 seconds at 0.2-second intervals), based on position data collected over several signal-history lengths (1-7 seconds at 1-second intervals). Based on patient data included in this study, the following observations are made. Qualitative comparison of predicted and actual position indicated a progressive increase in prediction error with an increase in response time. A signal-history length of 5 seconds was found to be the optimal signal history length for prediction using the sinusoidal model for all breathing training modalities. In terms of overall error in predicting respiratory motion, the adaptive filter model performed better than the sinusoidal model. With the adaptive filter, average prediction errors of less than 0.2 cm (1sigma) are possible for response times less than 0.4 seconds. In comparing prediction error with system latency error (no prediction), the adaptive filter model exhibited lesser prediction errors as compared to the sinusoidal model, especially for longer response time values (>0.4 seconds). At smaller response time values (<0.4 seconds), improvements in prediction error reduction are required for both predictive models in order to maximize gains in position accuracy due to prediction. Respiratory motion patterns are inherently complex in nature. While linear prediction-based prediction models perform satisfactorily for shorter response times, their prediction accuracy significantly deteriorates for longer response times. Successful implementation of real-time target-tracking-based radiotherapy requires response times less than 0.4 seconds or improved prediction algorithms.

Is the diaphragm motion probability density function normally distributed?

During radiotherapy treatment planning, the margins given to the clinical target volume to form the planning target volume accounts for internal motion and set-up error. Most margin formulas assume that the underlying distributions are independent and normal. Clinical data suggests that the set-up error probability density function (pdf) can be considered to have an approximately normal distribution. However, there is evidence that internal motion does not have a normal distribution. Thus, in general, a convolution of the two pdfs should be performed to determine the total geometric error. The goals of this article were to (1) determine if the internal motion pdf due to respiration can be characterized using a normal distribution, and (2) if not, determine if the total geometric uncertainty for combining internal motion and set-up error can be characterized by a normal distribution. Sixty fluoroscopy diaphragm motion data sets were obtained using three breathing training types: free breathing, audio instruction, and visual feedback. Diaphragm motion was used as a surrogate for liver and lung cancer motion. The data were analyzed with normality tests in the following groups: (1) single motion measurements, (2) combined motion measurements for each patient, and (3) combined motion measurements for all patients. Following this analysis, the diaphragm motion pdfs were convolved with a set-up error pdf, and the standard deviation of the set-up error pdf at which the total geometric error pdf became normal was determined. At set-up error standard deviation values of at least 0.27 and 0.1 cm for free breathing, 0.57 and 0.42 cm for audio instruction, and 0.55 and 0 cm for visual feedback, for single motion measurements and combined motion measurements for each patient, respectively, total geometric error pdfs became approximately normal. When the motion measurements for all the patients were combined, diaphragm motion pdfs were approximately normal for all feedback types. Therefore, for treatment planning purposes in the absence of individual patient measurements, the diaphragm motion pdf can be considered an approximately normal distribution. However, care should be taken when determining a margin based on individual patients measurements as the total geometric error will, in general, not be normally distributed.

On the use of EPID-based implanted marker tracking for 4D radiotherapy

Four-dimensional (4D) radiotherapy delivery to dynamically moving tumors requires a real-time signal of the tumor position as a function of time so that the radiation beam can continuously track the tumor during the respiration cycle. The aim of this study was to develop and evaluate an electronic portal imaging device (EPID)-based marker-tracking system that can be used for real-time tumor targeting, or 4D radiotherapy. Three gold cylinders, 3 mm in length and 1 mm in diameter, were implanted in a dynamic lung phantom. The phantom range of motion was 4 cm with a 3-s "breathing" period. EPID image acquisition parameters were modified, allowing image acquisition in 0.1 s. Images of the stationary and moving phantom were acquired. Software was developed to segment automatically the marker positions from the EPID images. Images acquired in 0.1 s displayed higher noise and a lower signal-noise ratio than those obtained using regular (> 1 s) acquisition settings. However, the markers were still clearly visible on the 0.1-s images. The motion of the phantom blurred the images of the markers and further reduced the signal-noise ratio, though they could still be successfully segmented from the images in 10-30 ms of computation time. The positions of gold markers placed in the lung phantom were detected successfully, even for phantom velocities substantially higher than those observed for typical lung tumors. This study shows that using EPID-based marker tracking for 4D radiotherapy is feasible, however, changes in linear accelerator technology and EPID-based image acquisition as well as patient studies are required before this method can be implemented clinically.

Monte Carlo as a four-dimensional radiotherapy treatment-planning tool to account for respiratory motion

Four-dimensional (4D) radiotherapy is the explicit inclusion of the temporal changes in anatomy during the imaging, planning and delivery of radiotherapy. Temporal anatomic changes can occur for many reasons, though the focus of the current investigation was respiration motion for lung tumours. The aims of the current research were first to develop a 4D Monte Carlo methodology and second to apply this technique to an existing 4D treatment plan. A 4D CT scan consisting of a series of 3D CT image sets acquired at different respiratory phases was used. Deformable image registration was performed to map each CT set from the end-inhale respiration phase to the CT image sets corresponding with subsequent respiration phases. This deformable registration allowed the contours drawn on the end-inhale CT to be automatically drawn on the other respiratory phase CT image sets. A treatment plan was created on the end-inhale CT image set and then automatically created on each of the 3D CT image sets corresponding with subsequent respiration phases, based on the beam arrangement and dose prescription in the end-inhale plan. Dose calculation using Monte Carlo was simultaneously performed on each of the N (=8) 3D image sets with 1/N fewer particles per calculation than for a 3D plan. The dose distribution from each respiratory phase CT image set was mapped back to the end-inhale CT image set for analysis. This use of deformable image registration to merge all the statistically noisy dose distributions back onto one CT image set effectively yielded a 4D Monte Carlo calculation with a statistical uncertainty equivalent to a 3D calculation, with a similar calculation time for the 3D and 4D methods. Monte Carlo as a dose calculation tool for 4D radiotherapy planning has two advantages: (1) higher accuracy for calculation in electronic disequilibrium conditions, such as those encountered during lung radiotherapy, and (2) if deformable image registration is used, the calculation time for Monte Carlo is independent of the number of 3D CT image sets constituting a 4D CT, unlike other algorithms for which the calculation time scales linearly with the number of 3D CT image sets constituting a 4D CT.

Acquiring 4D thoracic CT scans using a multislice helical method

Respiratory motion degrades anatomic position reproducibility during imaging, necessitates larger margins during radiotherapy planning and causes errors during radiation delivery. Computed tomography (CT) scans acquired synchronously with the respiratory signal can be used to reconstruct 4D CT scans, which can be employed for 4D treatment planning to explicitly account for respiratory motion. The aim of this research was to develop, test and clinically implement a method to acquire 4D thoracic CT scans using a multislice helical method. A commercial position-monitoring system used for respiratory-gated radiotherapy was interfaced with a third generation multislice scanner. 4D cardiac reconstruction methods were modified to allow 4D thoracic CT acquisition. The technique was tested on a phantom under different conditions: stationary, periodic motion and non-periodic motion. 4D CT was also implemented for a lung cancer patient with audio-visual breathing coaching. For all cases, 4D CT images were successfully acquired from eight discrete breathing phases, however, some limitations of the system in terms of respiration reproducibility and breathing period relative to scanner settings were evident. Lung mass for the 4D CT patient scan was reproducible to within 2.1% over the eight phases, though the lung volume changed by 20% between end inspiration and end expiration (870 cm3). 4D CT can be used for 4D radiotherapy, respiration-gated radiotherapy, 'slow' CT acquisition and tumour motion studies.

Quantifying the predictability of diaphragm motion during respiration with a noninvasive external marker

The aim of this work was to quantify the ability to predict intrafraction diaphragm motion from an external respiration signal during a course of radiotherapy. The data obtained included diaphragm motion traces from 63 fluoroscopic lung procedures for 5 patients, acquired simultaneously with respiratory motion signals (an infrared camera-based system was used to track abdominal wall motion). During these sessions, the patients were asked to breathe either (i) without instruction, (ii) with audio prompting, or (iii) using visual feedback. A statistical general linear model was formulated to describe the relationship between the respiration signal and diaphragm motion over all sessions and for all breathing training types. The model parameters derived from the first session for each patient were then used to predict the diaphragm motion for subsequent sessions based on the respiration signal. Quantification of the difference between the predicted and actual motion during each session determined our ability to predict diaphragm motion during a course of radiotherapy. This measure of diaphragm motion was also used to estimate clinical target volume (CTV) to planning target volume (PTV) margins for conventional, gated, and proposed four-dimensional (4D) radiotherapy. Results from statistical analysis indicated a strong linear relationship between the respiration signal and diaphragm motion (p<0.001) over all sessions, irrespective of session number (p=0.98) and breathing training type (p=0.19). Using model parameters obtained from the first session, diaphragm motion was predicted in subsequent sessions to within 0.1 cm (1 sigma) for gated and 4D radiotherapy. Assuming a 0.4 cm setup error, superior-inferior CTV-PTV margins of 1.1 cm for conventional radiotherapy could be reduced to 0.8 cm for gated and 4D radiotherapy. The diaphragm motion is strongly correlated with the respiration signal obtained from the abdominal wall. This correlation can be used to predict diaphragm motion, based on the respiration signal, to within 0.1 cm (1 sigma) over a course of radiotherapy.

Quantifying the effect of intrafraction motion during breast IMRT planning and dose delivery

Respiratory motion during intensity modulated radiation therapy (IMRT) causes two types of problems. First, the clinical target volume (CTV) to planning target volume (PTV) margin needed to account for respiratory motion means that the lung and heart dose is higher than would occur in the absence of such motion. Second, because respiratory motion is not synchronized with multileaf collimator (MLC) motion, the delivered dose is not the same as the planned dose. The aims of this work were to evaluate these problems to determine (a) the effects of respiratory motion and setup error during breast IMRT treatment planning, (b) the effects of the interplay between respiratory motion and multileaf collimator (MLC) motion during breast IMRT delivery, and (c) the potential benefits of breast IMRT using breath-hold, respiratory gated, and 4D techniques. Seven early stage breast cancer patient data sets were planned for IMRT delivered with a dynamic MLC (DMLC). For each patient case, eight IMRT plans with varying respiratory motion magnitudes and setup errors (and hence CTV to PTV margins) were created. The effects of respiratory motion and setup error on the treatment plan were determined by comparing the eight dose distributions. For each fraction of these plans, the effect of the interplay between respiratory motion and MLC motion during IMRT delivery was simulated by superimposing the respiratory trace on the planned DMLC leaf motion, facilitating comparisons between the planned and expected dose distributions. When considering respiratory motion in the CTV-PTV expansion during breast IMRT planning, our results show that PTV dose heterogeneity increases with respiratory motion. Lung and heart doses also increase with respiratory motion. Due to the interplay between respiratory motion and MLC motion during IMRT delivery, the planned and expected dose distributions differ. This difference increases with respiratory motion. The expected dose varies from fraction to fraction. However, for the seven patients studied and respiratory trace used, for no breathing, shallow breathing, and normal breathing, there were no statistically significant differences between the planned and expected dose distributions. Thus, for breast IMRT, intrafraction motion degrades treatment plans predominantly by the necessary addition of a larger CTV to PTV margin than would be required in the absence of such motion. This motion can be limited by breath-hold, respiratory gated, or 4D techniques.

Acquiring a four-dimensional computed tomography dataset using an external respiratory signal

Four-dimensional (4D) methods strive to achieve highly conformal radiotherapy, particularly for lung and breast tumours, in the presence of respiratory-induced motion of tumours and normal tissues. Four-dimensional radiotherapy accounts for respiratory motion during imaging, planning and radiation delivery, and requires a 4D CT image in which the internal anatomy motion as a function of the respiratory cycle can be quantified. The aims of our research were (a) to develop a method to acquire 4D CT images from a spiral CT scan using an external respiratory signal and (b) to examine the potential utility of 4D CT imaging. A commercially available respiratory motion monitoring system provided an 'external' tracking signal of the patient's breathing. Simultaneous recording of a TTL 'X-Ray ON' signal from the CT scanner indicated the start time of CT image acquisition, thus facilitating time stamping of all subsequent images. An over-sampled spiral CT scan was acquired using a pitch of 0.5 and scanner rotation time of 1.5 s. Each image from such a scan was sorted into an image bin that corresponded with the phase of the respiratory cycle in which the image was acquired. The complete set of such image bins accumulated over a respiratory cycle constitutes a 4D CT dataset. Four-dimensional CT datasets of a mechanical oscillator phantom and a patient undergoing lung radiotherapy were acquired. Motion artefacts were significantly reduced in the images in the 4D CT dataset compared to the three-dimensional (3D) images, for which respiratory motion was not accounted. Accounting for respiratory motion using 4D CT imaging is feasible and yields images with less distortion than 3D images. 4D images also contain respiratory motion information not available in a 3D CT image.

A fluence-convolution method to calculate radiation therapy dose distributions that incorporate random set-up error

The International Commission on Radiation Units and Measurements Report 62 (ICRU 1999) introduced the concept of expanding the clinical target volume (CTV) to form the planning target volume by a two-step process. The first step is adding a clinically definable internal margin, which produces an internal target volume that accounts for the size, shape and position of the CTV in relation to anatomical reference points. The second is the use of a set-up margin (SM) that incorporates the uncertainties of patient beam positioning, i.e. systematic and random set-up errors. We propose to replace the random set-up error component of the SM by explicitly incorporating the random set-up error into the dose-calculation model by convolving the incident photon beam fluence with a Gaussian set-up error kernel. This fluence-convolution method was implemented into a Monte Carlo (MC) based treatment-planning system. Also implemented for comparison purposes was a dose-matrix-convolution algorithm similar to that described by Leong (1987 Phys. Med. Biol. 32 327-34). Fluence and dose-matrix-convolution agree in homogeneous media. However, for the heterogeneous phantom calculations, discrepancies of up to 5% in the dose profiles were observed with a 0.4 cm set-up error value. Fluence-convolution mimics reality more closely, as dose perturbations at interfaces are correctly predicted (Wang et al 1999 Med. Phys. 26 2626-34, Sauer 1995 Med. Phys. 22 1685-90). Fluence-convolution effectively decouples the treatment beams from the patient. and more closely resembles the reality of particle fluence distributions for many individual beam-patient set-ups. However, dose-matrix-convolution reduces the random statistical noise in MC calculations. Fluence-convolution can easily be applied to convolution/superposition based dose-calculation algorithms.

Potential radiotherapy improvements with respiratory gating

Gating is a relatively new and potentially useful therapeutic addition to external beam radiotherapy applied to regions affected by intra-fraction motion. The impact was of gating on treatment margins, image artifacts, and volume and positional accuracy was investigated by CT imaging of sinusoidally moving spheres. The motion of the spheres simulates target motion. During the CT imaging of dynamically moving spheres, gating reproduced the static volume to within 1%, whereas errors of over 20% were observed where gating was not used. Using a theoretical analysis of margins, gating alone or in combination with an electronic portal imaging device may allow a 2-11 mm reduction in the CTV to PTV margin, and thus less healthy tissue need be irradiated. Gating may allow a reduction of treatment margins, an improvement in image quality, and an improvement in positional and volumetric accuracy of the gross tumor volume.

A Monte Carlo study of radiation transport through multileaf collimators

Due to the significant increase in the number of monitor units used to deliver a dynamic IMRT treatment, the total MLC leakage (transmission plus scatter) can exceed 10% of the maximum in-field dose. To avoid dosimetric errors, this leakage must be accurately accounted for in the dose calculation and conversion of optimized intensity patterns to MLC trajectories used for treatment delivery. In this study, we characterized the leaf end transmission and leakage radiation for Varian 80- and 120-leaf MLCs using Monte Carlo simulations. The complex geometry of the MLC, including the rounded leaf end, leaf edges (tongue-and-groove and offset notch), mounting slots, and holes was modeled using MCNP4b. Studies were undertaken to determine the leakage as a function of field size, components of the leakage, electron contamination, beam hardening and leaf tip effects. The leakage radiation with the MLC configured to fully block the field was determined. Dose for 6 and 18 MV beams was calculated at 5 cm depth in a water phantom located at 95 cm SSD, and normalized to the dose for an open field. Dose components were scored separately for radiation transmitted through and scattered from the MLC. For the 80-leaf MLC at 6 MV, the average leakage dose is 1.6%, 1.7%, 1.8%, and 1.9% for 5 x 5, 10 x 10, 15 x 15, and 20 x 20cm2 fields, respectively. For the 120-leaf MLC at 6 MV, the average leakage dose is 1.6%, 1.6%, 1.7%, and 1.9% for the same field sizes. Measured leakage values for the 120-leaf MLC agreed with calculated values to within 0.1% of the open field dose. The increased leakage with field size is attributed to MLC scattered radiation. The fractional electron contamination for a blocked MLC field is greater than that for an open field. The MLC attenuation significantly affects the photon spectrum, resulting in an increase in percent depth dose at 6 MV, however, little effect is observed at 18 MV. Both phantom scatter and the finite source size contribute to the leaf tip profile observed in phantom. The results of this paper can be applied to fluence-to-trajectory and trajectory-to-fluence calculations for IMRT.

Determining parameters for respiration-gated radiotherapy

Respiration-gated radiotherapy for tumor sites affected by respiratory motion will potentially improve radiotherapy outcomes by allowing reduced treatment margins leading to decreased complication rates and/or increased tumor control. Furthermore, for intensity-modulated radiotherapy (IMRT), respiratory gating will minimize the hot and cold spot artifacts in dose distributions that may occur as a result of interplay between respiratory motion and leaf motion. Most implementations of respiration gating rely on the real time knowledge of the relative position of the internal anatomy being treated with respect to that of an external marker. A method to determine the amplitude of motion and account for any difference in phase between the internal tumor motion and external marker motion has been developed. Treating patients using gating requires several clinical decisions, such as whether to gate during inhale or exhale, whether to use phase or amplitude tracking of the respiratory signal, and by how much the intrafraction tumor motion can be decreased at the cost of increased delivery time. These parameters may change from patient to patient. A method has been developed to provide the data necessary to make decisions as to the CTV to PTV margins to apply to a gated treatment plan.

Monte Carlo dose calculations for dynamic IMRT treatments

Dose calculations for intensity modulated radiation therapy (IMRT) face new challenges due to the complex leaf geometry and time dependent nature of the delivery. A fast method of particle transport through a dynamic multileaf collimator (MLC) geometry that accounts for photon attenuation and first-scattered Compton photon production has been incorporated into an existing Monte Carlo code used for patient dose calculations. Dosimetric agreement between calculation and measurement for two photon energies and MLC types is within experimental error for the sliding window tests. For a patient IMRT field, the Monte Carlo calculations are closer to measured dose than similar superposition or pencil beam calculations.

Motion adaptive x-ray therapy: a feasibility study

Intrafraction motion caused by breathing requires increased treatment margins for chest and abdominal radiotherapy and may lead to 'motion artefacts' in dose distributions during intensity modulated radiotherapy (IMRT). Technologies such as gated radiotherapy may significantly increase the treatment time, while breath-hold techniques may be poorly tolerated by pulmonarily compromised patients. A solution that allows reduced margins and dose distribution artefacts, without compromising delivery time, is to synchronously follow the target motion by adapting the x-ray beam using a dynamic multileaf collimator (MLC), i.e. motion adaptive x-ray therapy, or MAX-T for short. Though the target is moving with time, in the MAX-T beam view the target is static. The MAX-T method superimposes the target motion due to respiration onto the beam originally planned for delivery. Thus during beam delivery the beam is dynamically changing position with respect to the isocentre using a dynamic MLC, the leaf positions of which are dependent upon the target position. Synchronization of the MLC motion and target motion occurs using respiration gated radiotherapy equipment. The concept and feasibility of MAX-T and the capability of the treatment machine to deliver such a treatment were investigated by performing measurements for uniform and IMRT fields using a mechanical sinusoidal oscillator to simulate target motion. Target dose measurements obtained using MAX-T for a moving target were found to be equivalent to those delivered to a static target by a static beam.

A method for determining multileaf collimator transmission and scatter for dynamic intensity modulated radiotherapy

The main purpose of this work is to demonstrate a practical means of determining the leaf transmission and scatter characteristics of a multileaf collimator (MLC) pertinent to the commissioning of dynamic intensity modulated radiotherapy, especially for the sweeping window technique. The data are necessary for the conversion of intensity distributions produced by intensity-modulated radiotherapy optimization systems into trajectories of MLC leaves for dynamic delivery. Measurements are described for two, tungsten alloy MLCs: a Mark II 80-leaf MLC on a Varian 2100C accelerator and a Millenium 120-leaf MLC on a Varian 2100EX accelerator. MLC leakage was measured by film for a series of field sizes. Measured MLC leakage was 1.68% for a 10 x 10 cm2 field for both 6 and 18 MV for the 80-leaf MLC. For the 6 MV field, the 1.68% leakage consisted of 1.48% direct transmission and 0.20% leaf scatter. Direct transmission through the 80-leaf MLC, including the rounded leaf tip, was calculated analytically taking into account the detailed leaf geometry and a Monte Carlo-generated energy spectrum of the accelerator. The integrated fluence under the leaf tip was equivalent to an inward shift of 0.06 cm of a hypothetical leaf with a flat, focused tip. Monte Carlo calculations of the dose to phantom beyond a closed 80-leaf MLC showed excellent agreement with the analytic results. The transmission depends on the density of the MLC alloy, which may differ among individual MLCs. Thus, it is important to measure the transmission of any particular MLC. Calculated doses for a series of uniform fields produced by dynamic sweeping windows of various widths agree with measurements within 2%.

Converting absorbed dose to medium to absorbed dose to water for Monte Carlo based photon beam dose calculations

Current clinical experience in radiation therapy is based upon dose computations that report the absorbed dose to water, even though the patient is not made of water but of many different types of tissue. While Monte Carlo dose calculation algorithms have the potential for higher dose accuracy, they usually transport particles in and compute the absorbed dose to the patient media such as soft tissue, lung or bone. Therefore, for dose calculation algorithm comparisons, or to report dose to water or tissue contained within a bone matrix for example, a method to convert dose to the medium to dose to water is required. This conversion has been developed here by applying Bragg-Gray cavity theory. The dose ratio for 6 and 18 MV photon beams was determined by computing the average stopping power ratio for the primary electron spectrum in the transport media. For soft tissue, the difference between dose to medium and dose to water is approximately 1.0%, while for cortical bone the dose difference exceeds 10%. The variation in the dose ratio as a function of depth and position in the field indicates that for photon beams a single correction factor can be used for each particular material throughout the field for a given photon beam energy. The only exception to this would be for the clinically non-relevant dose to air. Pre-computed energy spectra for 60Co to 24 MV are used to compute the dose ratios for these photon beams and to determine an effective energy for evaluation of the dose ratio.