Sofie Ceberg

RapidArc treatment verification in 3D using polymer gel dosimetry and Monte Carlo simulation

The aim of this study was to verify the advanced inhomogeneous dose distribution produced by a volumetric arc therapy technique (RapidArc) using 3D gel measurements and Monte Carlo (MC) simulations. The TPS (treatment planning system)-calculated dose distribution was compared with gel measurements and MC simulations, thus investigating any discrepancy between the planned dose delivery and the actual delivery. Additionally, the reproducibility of the delivery was investigated using repeated gel measurements. A prostate treatment plan was delivered to a 1.3 liter nPAG gel phantom using one single arc rotation and a target dose of 3.3 Gy. Magnetic resonance imaging of the gel was carried out using a 1.5 T scanner. The MC dose distributions were calculated using the VIMC-Arc code. The relative absorbed dose differences were calculated voxel-by-voxel, within the volume enclosed by the 90% isodose surface (VOI(90)), for the TPS versus gel and TPS versus MC. The differences between the verification methods, MC versus gel, and between two repeated gel measurements were investigated in the same way. For all volume comparisons, the mean value was within 1% and the standard deviation of the differences was within 2.5% (1SD). A 3D gamma analysis between the dose matrices were carried out using gamma criteria 3%/3 mm and 5%/5 mm (% dose difference and mm distance to agreement) within the volume enclosed by the 50% isodose surface (VOI(50)) and the 90% isodose surface (VOI(90)), respectively. All comparisons resulted in very high pass rates. More than 95% of the TPS points were within 3%/3 mm of both the gel measurement and MC simulation, both inside VOI(50) and VOI(90). Additionally, the repeated gel measurements showed excellent consistency, indicating reproducible delivery. Using MC simulations and gel measurements, this verification study successfully demonstrated that the RapidArc plan was both accurately calculated and delivered as planned.

Verification of dynamic radiotherapy: the potential for 3D dosimetry under respiratory-like motion using polymer gel.

Following the implementation of advanced treatment procedures in radiotherapy, there is a need for dynamic dose verification in 3D. Gel dosimetry could potentially be used for such measurements. However, recently published data show that certain types of gels have a dose rate and fractionation dependence. The aim of this study was to investigate the feasibility of using a polymer gel dosimeter for dose verification of dynamic radiotherapy. To investigate the influence of dose rate dependence during respiratory-like motion in and out of the beam, a respiration robot together with two types of gel systems (normoxic methacrylic acid gel (nMAG) and normoxic polyacrylamide gel (nPAG)) were used. Reference measurements were obtained using a linear diode array (LDA). Expected results, if there was no influence of the dose rate variation, were calculated by convolving the static irradiated gel data with the motion function controlling the robot. To investigate the fractionation dependence, the gels were irradiated using gated and ungated deliveries. Magnetic resonance imaging was used to evaluate the absorbed dose response of the gel. The measured gel data coincided well with the LDA data. Also, the calculated data agreed well with the measured dynamic gel data, i.e. no dose rate dependence due to motion was observed. The difference in the R2 response for the gels receiving ungated and gated, i.e. fractionated, deliveries was less than 1% for the nPAG and 4% for the nMAG, for absorbed doses up to 2 Gy. The maximum difference was 1.2% for the nPAG and 9% for the nMAG, which occurred at the highest given dose (4 Gy). The investigated gels were found to be feasible detectors for dose measurements under respiratory-like motion. For dose verification of dynamic RT involving gated delivery, e.g. breathing-adapted radiotherapy, relative absorbed dose evaluation should be used in order to minimize the effects of fractionated irradiation.

Comparison of doses and NTCP to risk organs with enhanced inspiration gating and free breathing for left-sided breast cancer radiotherapy using the AAA algorithm

BackgroundThe purpose of this study was to investigate the potential dose reduction to the heart, left anterior descending (LAD) coronary artery and the ipsilateral lung for patients treated with tangential and locoregional radiotherapy for left-sided breast cancer with enhanced inspiration gating (EIG) compared to free breathing (FB) using the AAA algorithm. The radiobiological implication of such dose sparing was also investigated.MethodsThirty-two patients, who received tangential or locoregional adjuvant radiotherapy with EIG for left-sided breast cancer, were retrospectively enrolled in this study. Each patient was CT-scanned during FB and EIG. Similar treatment plans, with comparable target coverage, were created in the two CT-sets using the AAA algorithm. Further, the probability of radiation induced cardiac mortality and pneumonitis were calculated using NTCP models.ResultsFor tangential treatment, the median V25Gy for the heart and LAD was decreased for EIG from 2.2% to 0.2% and 40.2% to 0.1% (p < 0.001), respectively, whereas there was no significant difference in V20Gy for the ipsilateral lung (p = 0.109). For locoregional treatment, the median V25Gy for the heart and LAD was decreased for EIG from 3.3% to 0.2% and 51.4% to 5.1% (p < 0.001), respectively, and the median ipsilateral lung V20Gy decreased from 27.0% for FB to 21.5% (p = 0.020) for EIG. The median excess cardiac mortality probability decreased from 0.49% for FB to 0.02% for EIG (p < 0.001) for tangential treatment and from 0.75% to 0.02% (p < 0.001) for locoregional treatment. There was no significant difference in risk of radiation pneumonitis for tangential treatment (p = 0.179) whereas it decreased for locoregional treatment from 6.82% for FB to 3.17% for EIG (p = 0.004).ConclusionsIn this study the AAA algorithm was used for dose calculation to the heart, LAD and left lung when comparing the EIG and FB techniques for tangential and locoregional radiotherapy of breast cancer patients. The results support the dose and NTCP reductions reported in previous studies where dose calculations were performed using the pencil beam algorithm.

Technical evaluation of a laser-based optical surface scanning system for prospective and retrospective breathing adapted computed tomography

Abstract Background. For breathing adapted radiotherapy, the same motion monitoring system can be used for imaging and triggering of the accelerator. Purpose. To evaluate a new technique for prospective gated computed tomography (CT) and four-dimensional CT (4DCT) using a laser based surface scanning system (Sentinel™, C-RAD, Uppsala, Sweden). The system was compared to the AZ-733V respiratory gating system (Anzai Medical, Tokyo, Japan) and the Real-Time Position Management System (RPM™) (Varian Medical Systems, Palo Alto, CA, USA). Material and methods. Temporal accuracy was evaluated using a moving phantom programmed to move a platform along trajectories following a sin6(ωt) function with amplitudes from 6 to 20 mm and periods from 2 to 5 s during 120 s while the motion was recorded. The recorded data was Fourier transformed and the peak area at the fundamental and harmonic frequencies compared to data generated using the same sinusoidal function. For verification of the 4DCT reconstruction process, the phantom was programmed to move along a sinusoidal trajectory. Ten phase series were reconstructed. The distance from the couch to the platform was measured in each image. By fitting the function sin(ωt-φ) to the values measured in the images corresponding to each slice, the phase of each image was verified. Results and conclusion. In the recorded data, the peak area at the fundamental frequency covered on average 104 ± 4%, 102 ± 4% and 91 ± 27% of the peak area in the generated data for the Sentinel™, RPM™ and AZ-733V systems, respectively. All systems managed to resolve both harmonic frequencies. The second experiment showed that all images were sorted into the correct series using breathing data recorded by each system. The systems generated very similar results, however, it is preferable to use the same system both for imaging and treatment.

Dosimetric effects of intrafractional isocenter variation during deep inspiration breath‐hold for breast cancer patients using surface‐guided radiotherapy

Abstract The aim of this study was to investigate potential dose reductions to the heart, left anterior descending coronary artery (LAD), and ipsilateral lung for left‐sided breast cancer using visually guided deep inspiration breath‐hold (DIBH) with the optical surface scanning system Catalyst™, and how these potential dosimetric benefits are affected by intrafractional motion in between breath holds. For both DIBH and free breathing (FB), treatment plans were created for 20 tangential and 20 locoregional left‐sided breast cancer patients. During DIBH treatment, beam‐on was triggered by a region of interest on the xiphoid process using a 3 mm gating window. Using a novel nonrigid algorithm, the Catalyst™ system allows for simultaneous real‐time tracking of the isocenter position, which was used to calculate the intrafractional DIBH isocenter reproducibility. The 50% and 90% cumulative probabilities and maximum values of the intrafractional DIBH isocenter reproducibility were calculated and to obtain the dosimetric effect isocenter shifts corresponding to these values were performed in the treatment planning system. For both tangential and locoregional treatment, the dose to the heart, LAD and ipsilateral lung was significantly reduced for DIBH compared to FB. The intrafractional DIBH isocenter reproducibility was very good for the majority of the treatment sessions, with median values of approximately 1 mm in all three translational directions. However, for a few treatment sessions, intrafractional DIBH isocenter reproducibility of up to 5 mm was observed, which resulted in large dosimetric effects on the target volume and organs at risk. Hence, it is of importance to set tolerance levels on the intrafractional isocenter motion and not only perform DIBH based on the xiphoid process.

Dose integration and dose rate characteristics of a NiPAM polymer gel MRI dosimeter system

The normoxic polymer gel dosimeter based on N-isopropyl acrylamide (NiPAM) is a promising full 3D-dosimeter with high spatial resolution and near tissue equivalency. NiPAM gel samples were irradiated to different doses using a linear accelerator. The absorbed dose was evaluated using MRI and statistical significance of the analysed data was calculated. The analysis was carried out using an in-house developed software. It was found that the gel dosimeter responded linearly to the absorbed dose. The gel exhibited a dose rate dependence, as well as a dependence on the sequential beam irradiation scheme. A higher dose rate, as well as a higher dose per sequential beam, resulted in a lower dose response.

Evaluation of breathing interplay effects during VMAT by using 3D gel measurements

Respiratory motion during dynamic radiotherapy may affect the absorbed dose distribution both by dose-reducing smoothing and by more complicated interplay effects. In this study we present a novel method to determine the relative importance of these two effects. For the two dynamic deliveries studied in this work, the expected target dose reduction due to the smoothing effect was estimated by measurements convolved by the motion function. Remaining absorbed dose differences were attributed to interplay effects between the motion of the gel phantom and the movement of the modulating MLC leaves during modulated arc radiotherapy. The total dosimetric effect due to breathing motion and dynamic MLC motion during VMAT delivery resulted in an average of about 4% target dose reduction. Comparing with only the smoothing effect, the average difference was decreased to around 1%, and the remaining distribution was attributed to interplay effects. Although the interplay effects were small compared to the smoothing effect, the standard deviations of 1.4-2.3% (1SD) were larger than the narrow distribution for repeated stationary measurement with a standard deviation between 0.5-0.9% (1SD). (Less)

Modelling the dynamic dose response of an nMAG polymer gel dosimeter

Gel dosimetry measures the absorbed radiation dose with high spatial resolution in 3D. However, recently published data show that the response of metacrylic-based polymer gels depends on the segmented delivery pattern, which could potentially be a considerable disadvantage for measurements of modern dynamic radiotherapy techniques. The aim of this study is to design a dynamic compartment model for the response of a gel dosimeter, exposed to an arbitrary irradiation pattern (segmented delivery and intensity modulation), in order to evaluate the associated effects on absorbed dose measurements. The model is based on the separation of the protons affecting the magnetic resonance signal (i.e. the R2 value) into six compartments, described by a set of differential equations. The model is used to calculate R2 values for a number of different segmented delivery patterns between 0-4 Gy over 1-33 fractions. Very good agreement is found between calculated and measured R2 values, with an average difference of 0.3 ± 1.1% (1 SD). The model is also used to predict the behaviour of a gel dosimeter exposed to irradiation according to typical IMRT, VMAT and respiratory gating scenarios. The calculated R2 values are approximately independent of the segmented delivery, given that the same total dose is delivered during the same total time. It is concluded that this study helps to improve the theoretical understanding of the dependence of metacrylic-based polymer gel response to segmented radiation delivery.

Tumor-tracking radiotherapy of moving targets; verification using 3D polymer gel, 2D ion-chamber array and biplanar diode array

The aim of this study was to carry out a dosimetric verification of a dynamic multileaf collimator (DMLC)-based tumor-tracking delivery during respiratory-like motion. The advantage of tumor-tracking radiation delivery is the ability to allow a tighter margin around the target by continuously following and adapting the dose delivery to its motion. However, there are geometric and dosimetric uncertainties associated with beam delivery system constraints and output variations, and several investigations have to be accomplished before a clinical integration of this tracking technique. Two types of delivery were investigated in this study I) a single beam perpendicular to a target with a one dimensional motion parallel to the MLC moving direction, and II) an intensity modulated arc delivery (RapidArc®) with a target motion diagonal to the MLC moving direction. The feasibility study (I) was made using an 2D ionisation chamber array and a true 3D polymer gel. The arc delivery (II) was verified using polymer gel and a biplanar diode array. Good agreement in absorbed dose was found between delivery to a static target and to a moving target with DMLC tracking using all three detector systems. However, due to the limited spatial resolution of the 2D array a detailed comparison was not possible. The RapidArc® plan delivery was successfully verified using the biplanar diode array and true 3D polymer gel, and both detector systems could verify that the DMLC-based tumor-tracking delivery system has a very good ability to account for respiratory target motion.

Motion induced interplay effects for VMAT radiotherapy

The purpose of this study was to develop a method to simulate breathing motion induced interplay effects for volumetric modulated arc therapy (VMAT), to verify the proposed method with measurements, and to use the method to investigate how interplay effects vary with different patient- and machine specific parameters. VMAT treatment plans were created on a virtual phantom in a treatment planning system (TPS). Interplay effects were simulated by dividing each plan into smaller sub-arcs using an in-house developed software and shifting the isocenter for each sub-arc to simulate a sin6 breathing motion in the superior-inferior direction. The simulations were performed for both flattening-filter (FF) and flattening-filter free (FFF) plans and for different breathing amplitudes, period times, initial breathing phases, dose levels, plan complexities, CTV sizes, and collimator angles. The resulting sub-arcs were calculated in the TPS, generating a dose distribution including the effects of motion. The interplay effects were separated from dose blurring and the relative dose differences to 2% and 98% of the CTV volume (ΔD98% and ΔD2%) were calculated. To verify the simulation method, measurements were carried out, both static and during motion, using a quasi-3D phantom and a motion platform. The results of the verification measurements during motion were comparable to the results of the static measurements. Considerable interplay effects were observed for individual fractions, with the minimum ΔD98% and maximum ΔD2% being  -16.7% and 16.2%, respectively. The extent of interplay effects was larger for FFF compared to FF and generally increased for higher breathing amplitudes, larger period times, lower dose levels, and more complex treatment plans. Also, the interplay effects varied considerably with the initial breathing phase, and larger variations were observed for smaller CTV sizes. In conclusion, a method to simulate motion induced interplay effects was developed and verified with measurements, which allowed for a large number of treatment scenarios to be investigated. The simulations showed large interplay effects for individual fractions and that the extent of interplay effects varied with the breathing pattern, FFF/FF, dose level, CTV size, collimator angle, and the complexity of the treatment plan.