High Dose Gradient Regions Research Articles

Improving 3D dose prediction for breast radiotherapy using novel glowing masks and gradient-weighted loss functions.

The quality of treatment plans for breast cancer can vary greatly. This variation could be reduced by using dose prediction to automate treatment planning. Our work investigates novel methods for training deep-learning models that are capable of producing high-quality dose predictions for breast cancer treatment planning. The goal of this work was to compare the performance impact of two novel techniques for deep learning dose prediction models for tangent field treatments for breast cancer. The first technique, a "glowing" mask algorithm, encodes the distance from a contour into each voxel in a mask. The second, a gradient-weighted mean squared error (MSE) loss function, emphasizes the error in high-dose gradient regions in the predicted image. Four 3D U-Net deep learning models were trained using the planning CT and contours of the heart, lung, and tumor bed as inputs. The dataset consisted of 305 treatment plans split into 213/46/46 training/validation/test sets using a 70/15/15% split. We compared the impact of novel "glowing" anatomical mask inputs and a novel gradient-weighted MSE loss function to their standard counterparts, binary anatomical masks, and MSE loss, using an ablation study methodology. To assess performance, we examined the mean error and mean absolute error (ME/MAE) in dose across all within-body voxels, the error in mean dose to heart, ipsilateral lung, and tumor bed, dice similarity coefficient (DSC) across isodose volumes defined by 0%-100% prescribed dose thresholds, and gamma analysis (3%/3mm). The combination of novel glowing masks and gradient weighted loss function yielded the best-performing model in this study. This model resulted in a mean ME of 0.40%, MAE of 2.70%, an error in mean dose to heart and lung of -0.10 and 0.01Gy, and an error in mean dose to the tumor bed of -0.01%. The median DSC at 50/95/100% isodose levels were 0.91/0.87/0.82. The mean 3D gamma pass rate (3%/3mm) was 93%. This study found the combination of novel anatomical mask inputs and loss function for dose prediction resulted in superior performance to their standard counterparts. These results have important implications for the field of radiotherapy dose prediction, as the methods used here can be easily incorporated into many other dose prediction models for other treatment sites. Additionally, this dose prediction model for breast radiotherapy has sufficient performance to be used in an automated planning pipeline for tangent field radiotherapy and has the major benefit of not requiring a PTV for accurate dose prediction.

Development and first implementation of a novel multi-modality cardiac motion and dosimetry phantom for radiotherapy applications.

Cardiac applications in radiation therapy are rapidly expanding including magnetic resonance guided radiation therapy (MRgRT) for real-time gating for targeting and avoidance near the heart or treating ventricular tachycardia (VT). This work describes the development and implementation of a novel multi-modality and magnetic resonance (MR)-compatible cardiac phantom. The patient-informed 3D model was derived from manual contouring of a contrast-enhanced Coronary Computed Tomography Angiography scan, exported as a Stereolithography model, then post-processed to simulate female heart with an average volume. The model was 3D-printed using Elastic50A to provide MR contrast to water background. Two rigid acrylic modules containing cardiac structures were designed and assembled, retrofitting to an MR-safe programmable motor to supply cardiac and respiratory motion in superior-inferior directions. One module contained a cavity for an ion chamber (IC), and the other was equipped with multiple interchangeable cavities for plastic scintillation detectors (PSDs). Images were acquired on a 0.35 T MR-linac for validation of phantom geometry, motion, and simulated online treatment planning and delivery. Three motion profiles were prescribed: patient-derived cardiac (sine waveform, 4.3mm peak-to-peak, 60 beats/min), respiratory (cos4 waveform, 30mm peak-to-peak, 12 breaths/min), and a superposition of cardiac (sine waveform, 4mm peak-to-peak, 70 beats/min) and respiratory (cos4 waveform, 24mm peak-to-peak, 12 breaths/min). The amplitude of the motion profiles was evaluated from sagittal cine images at eight frames/s with a resolution of 2.4mm × 2.4mm. Gated dosimetry experiments were performed using the two module configurations for calculating dose relative to stationary. A CT-based VT treatment plan was delivered twice under cone-beam CT guidance and cumulative stationary doses to multi-point PSDs were evaluated. No artifacts were observed on any images acquired during phantom operation. Phantom excursions measured 49.3±25.8%/66.9±14.0%, 97.0±2.2%/96.4±1.7%, and 90.4±4.8%/89.3±3.5% of prescription for cardiac, respiratory, and cardio-respiratory motion profiles for the 2-chamber (PSD) and 12-substructure (IC) phantom modules respectively. In the gated experiments, the cumulative dose was<2% from expected using the IC module. Real-time dose measured for the PSDs at 10Hz acquisition rate demonstrated the ability to detect the dosimetric consequences of cardiac, respiratory, and cardio-respiratory motion when sampling of different locations during a single delivery, and the stability of our phantom dosimetric results over repeated cycles for the high dose and high gradient regions. For the VT delivery, high dose PSD was<1% from expected (5-6 cGy deviation of 5.9Gy/fraction) and high gradient/low dose regions had deviations<3.6% (6.3 cGy less than expected 1.73Gy/fraction). A novel multi-modality modular heart phantom was designed, constructed, and used for gated radiotherapy experiments on a 0.35 T MR-linac. Our phantom was capable of mimicking cardiac, cardio-respiratory, and respiratory motion while performing dosimetric evaluations of gated procedures using IC and PSD configurations. Time-resolved PSDs with small sensitive volumes appear promising for low-amplitude/high-frequency motion and multi-point data acquisition for advanced dosimetric capabilities. Illustrating VT planning and delivery further expands our phantom to address the unmet needs of cardiac applications in radiotherapy.

Practical experience with initial quality assurance of a high dose rate brachytherapy grid-based Boltzmann solver algorithm.

The purpose of this study was to validate the use of a model-based dose calculation algorithm (MBDCA), Acuros BV, for high dose rate brachytherapy treatment planning for a community-based hospital with a Bravos afterloader. Based on published AAPM recommendations, this work details a practical approach for community-based clinics to complete initial validation of Acuros BV, in order to add a MBDCA to a TG-43 based brachytherapy treatment planning program. Source dimensions and materials used in Acuros BV and TG-43 source models were compared to the physical source. TG-186 testing was completed with standardized test cases externally calculated with Monte Carlo compared to locally calculated with Acuros BV. Point doses calculated using TG-43 were compared to those calculated with Acuros BV in water at various dose grid settings. Secondary dose check software was used to evaluate dose distributions resembling clinical patient plans, both in water and on CT datasets representative of patient anatomy. The major source of discrepancy of source models was the length of modeled steel cable. TG-186 testing showed that the largest differences between Monte Carlo and Acuros BV dose distributions were located along the source axis for cases calculated in water, as well as located in regions of high dose gradients and within the applicator for the case calculated with a generic shielded applicator. An audit of point doses calculated with both TG-43 and Acuros BV in water found that dose grid settings significantly affected agreement. Secondary dose check software indicated that Acuros BV functioned satisfactorily, and a 5% threshold was adopted for secondary dose checks on gynecologic plans. This validation process indicated that Acuros BV met expected standards and affirmed its suitability for integration into this clinical practice's brachytherapy treatment planning.

Comparison of single catheter versus dual catheter-based EBT3 film calibration for the Ir-192 beam energy

Purpose. Films and TLDs have been the common choices for passive in-vivo dose measurement in radiotherapy. In the brachytherapy applications, it is very difficult to report and verify the dose at multiple localized high dose gradient regions and also the dose to organ at risk. This study was carried out to introduce a new and accurate calibration method for GafChromic EBT3 films irradiated using Ir-192 photon energy from miniature High Dose Rate (HDR) Brachytherapy source. Materials and methods. Film holder made of Styrofoam was used to hold the EBT3 film at its center. It was placed inside the mini water phantom and the films were irradiated by Ir-192 source of microSelectron HDR afterloading brachytherapy system. Two different setups: Single catheter-based film exposure and dual catheter-based film exposure were compared. The films scanned on a flatbed scanner were analysed in three different color channels: red, green, and blue using Image J software. The dose calibration graphs were generated using the third-order polynomial equations fitted on the data points from two different methods of calibration procedure. Maximum and mean dose difference between TPS calculated and measured was analyzed. Results. The measured dose difference from the TPS calculated doses were evaluated for the three groups of dose ranges (low, medium and high). In the high dose range, standard uncertainty of dose difference are ±2.3%, ±2.9%, and ±2.4% respectively for the red, green, and blue color channel when the TPS calculated dose was compared with single catheter based film calibration equation. Whereas it is observed as 1.3%, 1.4% and 3.1% for the red, green, and blue color channels respectively when compared with the dual catheter based film calibration equation. A test film was exposed to a TPS calculated dose of 666 cGy to validate the calibration equations, single catheter based film calibration equation estimated the dose difference as −9.2%, −7.8% and −3.6% respectively in the red, green, and blue color channels whereas the same were observed as 0.1%, 0.2% and 6.1% respectively when dual catheter based film calibration equation was applied. Conclusion. Source miniature size, reproducible positioning of the film and catheter system inside water medium are the major challenges in the film calibration with Ir-192 beam. To overcome these situations dual catheter-based film calibration was found more accurate and reproducible as compare to the single catheter based film calibration.

Protocol determination for OSL in vivo measurements of absorbed dose in the oral mucosa in oral cancer patients: A pilot study

Radiation-induced oral mucositis (RIOM) is one of the major oral complications caused by radiotherapy (RT), and can lead a negative impact on patients’ quality of life. Considering that there is no specific radiation therapy dose constraint for oral mucosa, the need to study the correlation between absorbed dose and effect in this area is crucial. This pilot study aims to report in vivo OSL dosimetry protocol to investigate the dose distribution in the oral mucosa in patients with oral cancer using Intraoral customised Stent (IOS). For measurements, aluminium oxide nanoDot dosimeters and a microStar ii OSL reader from Landauer Inc. were used. The OSL dosimetry system was calibrated with a 6 MV photon beam using an Elekta Synergy linear accelerator and solid water phantoms following the recommendations of the AAPM TG 191 protocol. The dose distribution in the oral mucosa using seven OSL dosimeters fixed on the intraoral stent was evaluated and compared with the dose distribution available from the RT planning for the patient. Our experimental results showed the expected linearity response of the nanoDot dosimeters from 20 to 200 cGy. The repeatability of readings was better than 1.0% in all commissioning measurements. The uncertainty budget of the commissioning analysis resulted in an overall type An uncertainty of ∼3.4% at the 1-sigma level for dose calculation, compatible with the TG 191 protocol. The dose distribution in the oral mucosa was assessed and dose deviations up to 35% can be found in high dose gradient regions of the treatment planning. This type of deviation may impact on the development and progression of RIOM, as analysed in this protocol. It is feasible to establish this procedure in vivo dosimetry for a patient undergoing RT using a multidisciplinary approach. Additionally, intraoral customised stents can decrease radiation doses in the oral mucosa, and the understanding of the dose distribution in the oral mucosa may help to minimise damages from radiation.

EVALUATION OF PATIENT-SPECIFIC IMRT QUALITY ASSURANCE AND POINT DOSE MEASUREMENT FOR COMPLEX HEAD AND NECK AND BRAIN CANCER USING GAFCHROMIC EBT3 FILM.

Patient-specific intensity-modulated radiation therapy (IMRT) quality assurance (QA) is essential for complex radiotherapy treatment as it involves complex intensity modulation and high-dose gradient regions. IMRT QA was performed by point dose verification and two-dimensional (2D) dose distribution measurement using gamma method. Calibrated External Beam Therapy 3 (EBT3) film was used for point dose and pre-treatment verification of 10 IMRT plans, five complex Head and Neck (HN) and five brain cases. The gamma passing rate (GPR) was evaluated for 3%/3mm gamma criteria and compared with 2D array. Isocentre dose was measured for all 10 IMRT plans on EBT3 film. Percentage deviation of point dose measurement from TPS calculated was found 0.4% for brain cases and 2.9% for HN cases. The GPR for 3%/3mm criteria was obtained higher than 95% for brain and HN cases. Results suggest that film dosimetry is also a reliable verification system for patient-specific IMRT QA as the 2D array.

Evaluation and Comparison of Dosimetric Characteristics of Semiflex®3D and Microdiamond in Relative Dosimetry under 6 and 15 MV Photon Beams in Small Fields

Background: In modern radiotherapy techniques, the frequently small and non-uniformed fields can increase treatment efficiency due to their highly conformal dose distribution. Particular features including lack of Lateral Charge Particle Equilibrium (LCPE) lead to detectors with high resolution since any error in obtained dosimetric data could cause patient mistreatments.Objective: This study aims to evaluate and compare two small detectors (Semiflex®3D and microdiamond) dosimetric characteristics in small field relative dosimetry.Material and Methods: In this experimental study, the dosimetric properties of Semiflex®3D and microdiamond were assessed under 6 and 15 MV photon beams. The linearity and stability of the detector’s response and dose rate were measured. Square-field sizes ranging from 0.6×0.6 - 5×5 cm2 were used for obtaining percentage depth dose curves (PDDs) and in-plane profiles. The angular and temperature dependence of both detectors’ responses were also studied.Results: The detector response shows good stability, no deviation from linearity, and low dose rate dependence (≤1.6%). PDDs and in-plan profiles of both detectors are in good agreement and no significant difference was observed except for the high dose gradient regions (P-value≤0.017). Both detectors demonstrated low angular dependence (<0.3%) with temperature dependence lower than 1% for both detectors.Conclusion: The results indicate both investigated detectors were well performed in small field relative dosimetry and for measuring penumbra, it is better to use microdiamond detector.

Development of a new Hybrid Virtual Source Model to simulate Elekta Synergy MLCi2 linac

Monte Carlo (MC) simulations of medical linear accelerators applying the phase space approach provide excellent results and substantially reduce execution time. This approach becomes impractical because a new complete MC simulation is needed for any new generation of linac heads. To overcome this disadvantage, the investigation exhibits an innovative hybrid Virtual Source Model (VSM) based on a theoretical reconstruction model utilizing the unchangeable parameters of linacs, taking as input the distribution of dose profiles and energy spectra weight. This VSM involves a couple of principal components, virtual source of photons and scattered particle sources from the Virtual Flattening Filter (VFF) and the asymmetric jaws. Extensive calculations and comparisons between the experimental dose distributions and the hybrid VSM are performed to evaluate VSM accuracy. The TOpas and Gate MC codes are utilized to achieve dose computational speedup and memory management optimization. For validation, several regular squared fields of varying sizes ranging from 3 × 3 to 20×20cm2 are considered, as well as two complex IMRT plans based on head&neck and pelvic treatments. Examining the VSM accuracy in the penumbra and the high dose gradient regions, applying gamma index (2%/2mm) criterion confirmed satisfactory results. A two and three dimensional γ index examination demonstrated an excellent agreement concerning the (3%/3mm) criteria for IMRT treatments. Furthermore, the speedup factor reached is 290 when applying this hybrid VSM compared to the entire simulation. Consequently, and according to the results achieved in this task, the hybrid VSM proposed is considered suitable for treatment planning systems (TPS) calculation purposes.

On measuring the 3D dose distribution for notched and circular Ru-106 plaque shapes through Gafchromic film dosimetry approach

Ruthenium-106 (106Ru) plaques are specific applicators dedicated for intraocular brachytherapy. Small tumor volume and presence of vital organs in its close vicinity as well as appearance of very high dose gradient regions around the employed plaques cause unique challenges in plaque-based brachytherapy. Therefore, access to 3D dose distribution for each brachytherapy plaque is quite helpful to decide about tumor volume coverage and overexposure of surrounding sensitive organs. Current study aims to measure the 3D dose distribution for two commercially available 106Ru plaques using film dosimetry.COB (a notched-shape plaque) and CGD (a fully circular plaque) 106Ru plaques were used and dose distributions in different planes were measured by simultaneously employing a 3D printed eyeball phantom and Gafchromic EBT3 films. Finally, 3D dose distributions were acquired by interpolating the measured 2D dose data at different planes.The results showed that employing the COB plaque tilts the isodose curves toward the non-notched edge of the plaque, while no isodose tilting was observed for CGD plaque. The obtained 3D dose distributions were different for studied plaques. Non-uniform 3D dose distributions were found for both plaques which can be attributed to the plaque design and reference point position on plaque central axis.Finally, it can be concluded that the obtained 3D dose distributions can specify the optimal tumor shape which can be adequately covered by brachytherapy plaque with minimizing the dose non-uniformity and hot point formation inside the treated volume. Obtained 3D dose distribution for each plaque can be employed for patient treatment planning.

Biologically consistent dose accumulation using daily patient imaging

BackgroundThis work addresses a basic inconsistency in the way dose is accumulated in radiotherapy when predicting the biological effect based on the linear quadratic model (LQM). To overcome this inconsistency, we introduce and evaluate the concept of the total biological dose, bEQDd.MethodsDaily computed tomography imaging of nine patients treated for prostate carcinoma with intensity-modulated radiotherapy was used to compute the delivered deformed dose on the basis of deformable image registration (DIR). We compared conventional dose accumulation (DA) with the newly introduced bEQDd, a new method of accumulating biological dose that considers each fraction dose and tissue radiobiology. We investigated the impact of the applied fractionation scheme (conventional/hypofractionated), uncertainties induced by the DIR and by the assigned α/β-value.ResultsbEQDd was systematically higher than the conventionally accumulated dose with difference hot spots of 3.3–4.9 Gy detected in six out of nine patients in regions of high dose gradient in the bladder and rectum. For hypofractionation, differences are up to 8.4 Gy. The difference amplitude was found to be in a similar range to worst-case uncertainties induced by DIR and was higher than that induced by α/β. ConclusionUsing bEQDd for dose accumulation overcomes a potential systematic inaccuracy in biological effect prediction based on accumulated dose. Highest impact is found for serial-type late responding organs at risk in dose gradient regions and for hypofractionation. Although hot spot differences are in the order of several Gray, in dose-volume parameters there is little difference compared with using conventional or biological DA. However, when local dose information is used, e.g. dose surface maps, difference hot spots can potentially change outcomes of dose-response modelling and adaptive treatment strategies.

Derivative-based gamma index: a novel methodology for stringent patient-specific quality assurance in the stereotactic treatment planning of liver cancer

Objective: The development of a stringent derivative-based gamma (DBG) index for patient-specific QA in stereotactic radiotherapy treatment planning (SRTP) to account for the spatial change in dose. Methods: Twenty-five patients of liver SBRT were selected retrospectively for this study. Deliberately, two different kinds of treatment planning approaches were used for each patient. Firstly, the treatment plans were generated using a conventional treatment planning (CTP) approach in which the target was covered with a homogeneous dose along with the nominal dose fall-off around the treatment field. Subsequently, the other treatment plans were generated using an SRTP approach with the intent of heterogeneous dose within the target region along with a steeper dose gradient outside the treatment field as much as possible. For both kinds of treatment plans, two dimensional (2D) conventional gamma (CG) and DBG analysis were performed using the 2D ion chamber array and radiochromic film. Results: Difference in the DBG index was statistically significant whereas, for CG analysis, the difference in CG index was insignificant for both types of treatment plans (CTP and SRTP). A significant positive correlation was observed between the difference in the DBG index and the difference in HI for high gamma criteria. Conclusion: The DBG evaluation is found to be more rigorous, and sensitive to the only SRTP. The proposed method could be opted-in the routine clinical practice in addition to CG. Advances in knowledge: DBG is more sensitive to detect the spatial change of dose, especially in high dose gradient regions.

Error detection and classification in patient-specific IMRT QA with dual neural networks.

Despite being the standard metric in patient-specific quality assurance (QA) for intensity-modulated radiotherapy (IMRT), gamma analysis has two shortcomings: (a) it lacks sensitivity to small but clinically relevant errors (b) it does not provide efficient means to classify the error sources. The purpose of this work is to propose a dual neural network method to achieve simultaneous error detection and classification in patient-specific IMRT QA. For a pair of dose distributions, we extracted the dose difference histogram (DDH) for the low dose gradient region and two signed distance-to-agreement (sDTA) maps (one in x direction and one in y direction) for the high dose gradient region. An artificial neural network (ANN) and a convolutional neural network (CNN) were designed to analyze the DDH and the two sDTA maps, respectively. The ANN was trained to detect and classify six classes of dosimetric errors: incorrect multileaf collimator (MLC) transmission (±1%) and four types of monitor unit (MU) scaling errors (±1% and ±2%). The CNN was trained to detect and classify seven classes of spatial errors: incorrect effective source size, 1mm MLC leaf bank overtravel or undertravel, 2mm single MLC leaf overtravel or undertravel, and device misalignment errors (1mm in x- or y direction). An in-house planar dose calculation software was used to simulate measurements with errors and noise introduced. Both networks were trained and validated with 13 IMRT plans (totaling 88 fields). A fivefold cross-validation technique was used to evaluate their accuracy. Distinct features were found in the DDH and the sDTA maps. The ANN perfectly identified all four types of MU scaling errors and the specific accuracies for the classes of no error, MLC transmission increase, MLC transmission decrease were 98.9%, 96.6%, and 94.3%, respectively. For the CNN, the largest confusion occurred between the 1-mm-MLC bank overtravel class and the 1-mm-device alignment error in x-direction class, which brought the specific accuracies down to 90.9% and 92.0%, respectively. The specific accuracy for the 2-mm-single MLC leaf undertravel class was 93.2% as it misclassified 5.7% of the class as being error free (false negative). Otherwise, the specific accuracy was above 95%. The overall accuracies across the fivefold were 98.3±0.7% and 95.6%±1.5% for the ANN and the CNN, respectively. Both the DDH and the sDTA maps are suitable features for error classification in IMRT QA. The proposed dual neural network method achieved simultaneous error detection and classification with excellent accuracy. It could be used in complement with the gamma analysis to potentially shift the IMRT QA paradigm from passive pass/fail analysis to active error detection and root cause identification.

Statistical approach to the selection of the tolerances for distance to agreement improves the quality control of the dose delivery in radiotherapy

In the study, a local approach to setting reference tolerance values for the distance-to-agreement (DTA) component of the gamma index is proposed. The reference tolerance values are calculated in simulations, following a dose delivery model presented in a previous work. An analytical model for determining the quantiles of DTA distribution is also proposed and verified. It is shown that the distributions of DTA values normalized with either quantiles or standard deviation of DTA distributions are universal over analyzed plans and points within a single plan. This enables statistically sound inference about the quality of dose delivery. In particular, based on the normalized distributions the comparison of planned and delivered doses can be formulated within the framework of statistical inference as a problem of multiple statistical testing. For every evaluated point P of a plan, one may formulate and test a null hypothesis that there is no delivery error against an alternative hypothesis that there is a delivery error in P. It is also shown that the proposed approach is more sensitive than the current standard approach to shift errors in high dose gradient regions.

Evaluation of Delta4DVH Anatomy in 3D Patient-Specific IMRT Quality Assurance.

Purpose:To evaluate the performance of Delta4DVH Anatomy in patient-specific intensity-modulated radiotherapy quality assurance.Materials and Methods:Dose comparisons were performed between Anatomy doses calculated with treatment plan dose measured modification and pencil beam algorithms, treatment planning system doses, film doses, and ion chamber measured doses in homogeneous and inhomogeneous geometries. The sensitivity of Anatomy doses to machine errors and output calibration errors was also investigated.Results:For a Volumetric Modulated Arc Therapy (VMAT) plan evaluated on the Delta4 geometry, the conventional gamma passing rate was 99.6%. For a water-equivalent slab geometry, good agreements were found between dose profiles in film, treatment planning system, and Anatomy treatment plan dose measured modification and pencil beam calculations. Gamma passing rate for Anatomy treatment plan dose measured modification and pencil beam doses versus treatment planning system doses was 100%. However, gamma passing rate dropped to 97.2% and 96% for treatment plan dose measured modification and pencil beam calculations in inhomogeneous head & neck phantom, respectively. For the 10 patients’ quality assurance plans, good agreements were found between ion chamber measured doses and the planned ones (deviation: 0.09% ± 1.17%). The averaged gamma passing rate for conventional and Anatomy treatment plan dose measured modification and pencil beam gamma analyses in Delta4 geometry was 99.6% ± 0.89%, 98.54% ± 1.60%, and 98.95% ± 1.27%, respectively, higher than averaged gamma passing rate of 97.75% ± 1.23% and 93.04% ± 2.69% for treatment plan dose measured modification and pencil beam in patients’ geometries, respectively. Anatomy treatment plan dose measured modification dose profiles agreed well with those in treatment planning system for both Delta4 and patients’ geometries, while pencil beam doses demonstrated substantial disagreement in patients’ geometries when compared to treatment planning system doses. Both treatment planning system doses are sensitive to multileaf collimator and monitor unit (MU) errors for high and medium dose metrics but not sensitive to the gantry and collimator rotation error smaller than 3°.Conclusions:The new Delta4DVH Anatomy with treatment plan dose measured modification algorithm is a useful tool for the anatomy-based patient-specific quality assurance. Cautions should be taken when using pencil beam algorithm due to its limitations in handling heterogeneity and in high-dose gradient regions.

Analysis of central axis and build up region dose calculated by different dose calculation algorithms used for radiation treatment planning

Introduction: Computerized Treatment planning system (TPS) is an integral part of radiation treatment procedure. The main part of treatment planning system is a calculation algorithm incorporated in the Treatment Planning System. Various algorithms used in the TPS have capabilities and limitation on specific situations. Objective: The objective of the work is to analyse the depth dose distribution generated by different dose calculation algorithms used in TPS in regions of low and high dose gradient. Materials and Methods: Four different dose calculation algorithms used for treatment planning such as Modified Clarkson, Sector integration, Fast Fourier Transform, Convolution- Superposition algorithm were used to generate the percentage depth dose distribution in a water phantom up to 25 cm depth. The values are compared with experimental results using Radiation Field Analyser. The PDD values also generated at build up region and the values were compared with experimental results using a Parallel plate ionization chamber. Results: The results show that the dose calculated by different algorithms shows difference at build up region and at larger depths.

Secondary monitor unit calculations for VMAT using parallelized Monte Carlo simulations

We have developed a fast and accurate in‐house Monte Carlo (MC) secondary monitor unit (MU) check method, based on the EGSnrc system, for independent verification of volumetric modulated arc therapy (VMAT) treatment planning system dose calculations, in accordance with TG‐114 recommendations. For a VMAT treatment plan created for a Varian Trilogy linac, DICOM information was exported from Eclipse. An open‐source platform was used to generate input files for dose calculations using the EGSnrc framework. The full VMAT plan simulation employed 107 histories, and was parallelized to run on a computer cluster. The resulting 3ddose matrices were converted to the DICOM format using CERR and imported into Eclipse. The method was evaluated using 35 clinical VMAT plans of various treatment sites. For each plan, the doses calculated with the MC approach at four three‐dimensional reference points were compared to the corresponding Eclipse calculations, as well as calculations performed using the clinical software package, MUCheck. Each MC arc simulation of 107 particles required 13–25 min of total time, including processing and calculation. The average discrepancies in calculated dose values between the MC method and Eclipse were 2.03% (compared to 3.43% for MUCheck) for prostate cases, 2.45% (3.22% for MUCheck) for head and neck cases, 1.7% (5.51% for MUCheck) for brain cases, and 2.84% (5.64% for MUCheck) for miscellaneous cases. Of 276 comparisons, 201 showed greater agreement between the treatment planning system and MC vs MUCheck. The largest discrepancies between MC and MUCheck were found in regions of high dose gradients and heterogeneous densities. By parallelizing the calculations, point‐dose accuracies of 2‐7%, sufficient for clinical secondary checks, can be achieved in a reasonable amount of time. As computer clusters and/or cloud computing become more widespread, this method will be useful in most clinical setups.

OA087] IVD methodology in apbi brachytherapy

Purpose To establish a methodology for the comparison of the dose measured by dosimeters used in in vivo dosimetry (IVD) of 192Ir brachytherapy (BT) treatments with the dose planned by treatment planning system (TPS) with two different algorithms. Methods IVD was performed on 9 accelerated partial breast irradiation (APBI) patients. The dosimeter used was a micromosfet placed on the skin surface in 3 positions marked with radiopaque markers: a position near the tumor site, and two other positions located at the beginning and end of the treatment tubes. For the calculation of the dose planned by TPS, and taking into account a possible positional error of the dosimeter due to the markers artifact in the planning CT image, we created a volume of interest (VOI) expanded by 2 mm around these radiopaque markers, calculated minimum and maximum doses in that VOI, and compared the dosimeter data within that dose range. Treatment plans were performed using the traditional TG43 method and a commercially model-based dose calculation algorithm (MBDCA) following TG186 recommendations. Results Regarding the planned dose range calculated with TG43, 4 of the total 27 measurements performed fell outside this range. The minimum and maximum deviation found was 2 and 16%. Regarding the TG186 planned dose range, 3 of the 27 measurements fell outside, with a minimum and maximum deviation of 4 and 23%. Conclusions The dosimeter deviation from the planned dose may have been influenced by a variation of its position relative to the radiopaque marker, an important factor in regions of high dose gradients, and by the fact that TPS may not include heterogeneities in its dose calculation, an important factor in the tissue-air barrier. Nevertheless, it happened in a maximum of 4 out of 27 treatments, thus highlighting the good performance of these dosimeters in BT IVD within our methodology of using a planned dose range. As future work, we intend to observe the influence of a bolus material on the skin in the dosimeter behavior and acquire a deeper knowledge of this new algorithm and its applications in a more accurate dose calculation.

A mathematical deconvolution formulation for superficial dose distribution measurement by Cerenkov light dosimetry.

Cerenkov photons are created by high-energy radiation beams used for radiation therapy. In this study, we developed a Cerenkov light dosimetry technique to obtain a two-dimensional dose distribution in a superficial region of medium from the images of Cerenkov photons by using a deconvolution method. An integral equationwas derived to represent the Cerenkov photon image acquired by a camera for a given incident high-energy photon beam by using convolution kernels. Subsequently, an equationrelating the planar dose at a depth to a Cerenkov photon image using the well-known relationship between the incident beam fluence and the dose distribution in a medium was obtained. The final equationcontained a convolution kernel called the Cerenkov dose scatter function (CDSF). The CDSF function was obtained by deconvolving the Cerenkov scatter function (CSF) with the dose scatter function (DSF). The GAMOS (Geant4-based Architecture for Medicine-Oriented Simulations) Monte Carlo particle simulation software was used to obtain the CSF and DSF. The dose distribution was calculated from the Cerenkov photon intensity data using an iterative deconvolution method with the CDSF. The theoretical formulation was experimentally evaluated by using an optical phantom irradiated by high-energy photon beams. The intensity of the deconvolved Cerenkov photon image showed linear dependence on the dose rate and the photon beam energy. The relative intensity showed a field size dependence similar to the beam output factor. Deconvolved Cerenkov images showed improvement in dose profiles compared with the raw image data. In particular, the deconvolution significantly improved the agreement in the high dose gradient region, such as in the penumbra. Deconvolution with a single iteration was found to provide the most accurate solution of the dose. Two-dimensional dose distributions of the deconvolved Cerenkov images agreed well with the reference distributions for both square fields and a multileaf collimator (MLC) defined, irregularly shaped field. The proposed technique improved the accuracy of the Cerenkov photon dosimetry in the penumbra region. The results of this study showed initial validation of the deconvolution method for beam profile measurements in a homogeneous media. The new formulation accounted for the physical processes of Cerenkov photon transport in the medium more accurately than previously published methods.

Validation of a virtual source model of medical linac for Monte Carlo dose calculation using multi-threaded Geant4

The use of phase space in medical linear accelerator Monte Carlo (MC) simulations significantly improves the execution time and leads to results comparable to those obtained from full calculations. The classical representation of phase space stores directly the information of millions of particles, producing bulky files. This paper presents a virtual source model (VSM) based on a reconstruction algorithm, taking as input a compressed file of roughly 800 kb derived from phase space data freely available in the International Atomic Energy Agency (IAEA) database. This VSM includes two main components; primary and scattered particle sources, with a specific reconstruction method developed for each. Energy spectra and other relevant variables were extracted from IAEA phase space and stored in the input description data file for both sources. The VSM was validated for three photon beams: Elekta Precise 6 MV/10 MV and a Varian TrueBeam 6 MV. Extensive calculations in water and comparisons between dose distributions of the VSM and IAEA phase space were performed to estimate the VSM precision. The Geant4 MC toolkit in multi-threaded mode (Geant4-[mt]) was used for fast dose calculations and optimized memory use. Four field configurations were chosen for dose calculation validation to test field size and symmetry effects, , , and for squared fields, and for an asymmetric rectangular field. Good agreement in terms of formalism, for 3%/3 mm and 2%/3 mm criteria, for each evaluated radiation field and photon beam was obtained within a computation time of 60 h on a single WorkStation for a 3 mm voxel matrix. Analyzing the VSM’s precision in high dose gradient regions, using the distance to agreement concept (DTA), showed also satisfactory results. In all investigated cases, the mean DTA was less than 1 mm in build-up and penumbra regions. In regards to calculation efficiency, the event processing speed is six times faster using Geant4-[mt] compared to sequential Geant4, when running the same simulation code for both. The developed VSM for 6 MV/10 MV beams widely used, is a general concept easy to adapt in order to reconstruct comparable beam qualities for various linac configurations, facilitating its integration for MC treatment planning purposes.

A 3D correction method for predicting the readings of a PinPoint chamber on the CyberKnife® M6™ machine

The use of small fields in radiation therapy techniques has increased substantially in particular in stereotactic radiosurgery (SRS) and stereotactic body radiation therapy (SBRT). However, as field size reduces further still, the response of the detector changes more rapidly with field size, and the effects of measurement uncertainties become increasingly significant due to the lack of lateral charged particle equilibrium, spectral changes as a function of field size, detector choice, and subsequent perturbations of the charged particle fluence. This work presents a novel 3D dose volume-to-point correction method to predict the readings of a 0.015 cc PinPoint chamber (PTW 31014) for both small static-fields and composite-field dosimetry formed by fixed cones on the CyberKnife® M6™ machine.A 3D correction matrix is introduced to link the 3D dose distribution to the response of the PinPoint chamber in water. The parameters of the correction matrix are determined by modeling its 3D dose response in circular fields created using the 12 fixed cones (5 mm–60 mm) on a CyberKnife® M6™ machine. A penalized least-square optimization problem is defined by fitting the calculated detector reading to the experimental measurement data to generate the optimal correction matrix; the simulated annealing algorithm is used to solve the inverse optimization problem. All the experimental measurements are acquired for every 2 mm chamber shift in the horizontal planes for each field size. The 3D dose distributions for the measurements are calculated using the Monte Carlo calculation with the MultiPlan® treatment planning system (Accuray Inc., Sunnyvale, CA, USA). The performance evaluation of the 3D conversion matrix is carried out by comparing the predictions of the output factors (OFs), off-axis ratios (OARs) and percentage depth dose (PDD) data to the experimental measurement data. The discrepancy of the measurement and the prediction data for composite fields is also performed for clinical SRS plans.The optimization algorithm used for generating the optimal correction factors is stable, and the resulting correction factors were smooth in the spatial domain. The measurement and prediction of OFs agree closely with percentage differences of less than 1.9% for all the 12 cones. The discrepancies between the prediction and the measurement PDD readings at 50 mm and 80 mm depth are 1.7% and 1.9%, respectively. The percentage differences of OARs between measurement and prediction data are less than 2% in the low dose gradient region, and 2%/1 mm discrepancies are observed within the high dose gradient regions. The differences between the measurement and prediction data for all the CyberKnife based SRS plans are less than 1%.These results demonstrate the existence and efficiency of the novel 3D correction method for small field dosimetry. The 3D correction matrix links the 3D dose distribution and the reading of the PinPoint chamber. The comparison between the predicted reading and the measurement data for static small fields (OFs, OARs and PDDs) yield discrepancies within 2% for low dose gradient regions and 2%/1 mm for high dose gradient regions; the discrepancies between the predicted and the measurement data are less than 1% for all the SRS plans. The 3D correction method provides an access to evaluate the clinical measurement data and can be applied to non-standard composite fields intensity modulated radiation therapy point dose verification.