Dose-response dependencies of OSL dosimeters in conventional linacs and 1.5T MR-linacs: an experimental and Monte Carlo study.

Objective. This work focuses on the optically stimulated luminescence dosimetry (OSLD) dose-response characterization, with emphasis on 1.5T MR-Linacs.Approach. Throughout this study, the nanoDots OSLDs (Landauer, USA) were considered. In groups of three, the mean OSLD response was measured in a conventional linac and an MR-Linac under various irradiation conditions to investigate (i) dose-response linearity with and without the 1.5T magnetic field, (ii) signal fading rate and its dependencies, (iii) beam quality, detector orientation and dose rate dependencies in a conventional linac, (iii) potential MR imaging related effects on OSLD response and (iv) detector orientation dependence in an MR-Linac. Monte Carlo calculations were performed to further quantify angular dependence after rotating the detector around its central axis parallel to the magnetic field, and determine the magnetic field correction factors,kB,Q,for all cardinal detector orientations.Main results. OSLD dose-response supralinearity in an MR-Linac setting was found to agree within uncertainties with the corresponding one in a conventional linac, for the axial detector orientation investigated. Signal fading rate does not depend on irradiation conditions for the range of 3-30 d considered. OSLD angular (orientation) dependence is more pronounced under the presence of a magnetic field. OSLDs irradiated with and without real-time T2w MR imaging enabled during irradiation yielded the same response within uncertainties.kB,Qvalues were determined for all three cardinal orientations. Corrections needed reached up to 6.4%. However, if OSLDs are calibrated in the axial orientation and then irradiated in an MR-Linac placed again in the axial orientation (perpendicular to the magnetic field), then simulations suggest thatkB,Qcan be considered unity within uncertainties, irrespective of the incident beam angle.Significance. This work contributes towards OSLD dose-response characterization and relevant correction factors availability. OSLDs are suitable for QA checks in MR-based beam gating applications andin vivodosimetry in MR-Linacs.

Dosimetric accuracy of the Convolution algorithm for Leksell Gamma Plan radiosurgery treatment planning: Evaluation in the presence of clinically relevant inhomogeneities.

PURPOSE: The Leksell Gamma Plan Convolution algorithm (LGP-Convolution) has not been widely adopted. This mainly stems from the higher calculated beam-on times relative to the standard ray tracing-based LGP-TMR10 dose calculation algorithm. This study aims to evaluate the accuracy of the LGP-Convolution in scenarios where the treated lesions are in the vicinity of or encompassed by bone and/or air inhomogeneities. METHODS: The solid water dosimetry phantom provided by the vendor was modified to include bone and air inhomogeneities. Two treatment planning scenarios were investigated involving a single shot and multiple shots, respectively. Treatment planning and dose prescription were performed using the LGP-Convolution algorithm. Triple channel film dosimetry was performed using GafChromic EBT3 films calibrated in terms of absorbed dose to water in a 60 Co beam. Monte Carlo (MC) simulation dosimetry was also performed in the inhomogeneous experimental geometry using the EGSnrc MC platform and a previously validated sector-based phase-space source model. MC simulations were also employed to determine correction factors required for converting EBT3 measurements at points within the bone and air inhomogeneities from dose-to-water values to the corresponding dose to medium values. RESULTS AND CONCLUSIONS: EBT3 dose to medium correction factors ranged with field size (4, 8, or 16 mm) within 0.941-0.946 for bone and 0.745-0.749 for air inhomogeneities. An excellent agreement was found between the LGP-Convolution calculations with corresponding EBT3 and MC dose to medium results at all measurement points, except those located inside the air inhomogeneity. The latter is of no clinical importance and excluding them yielded gamma index passing rates of nearly 100% for 3% local dose difference and 1 mm distance-to-agreement criteria. The excellent agreement observed between LGP-Convolution calculations and film as well as MC results of dose to medium indicates that the latter is the quantity reported by the LGP-Convolution.

Dosimetry in 1.5 T MR-Linacs: Monte Carlo determination of magnetic field correction factors and investigation of the air gap effect.

BACKGROUND: In Magnetic Resonance-Linac (MR-Linac) dosimetry formalisms, a new correction factor, kB,Q , has been introduced to account for corresponding changes to detector readings under the beam quality, Q, and the presence of magnetic field, B. PURPOSE: This study aims to develop and implement a Monte Carlo (MC)-based framework for the determination of kB,Q correction factors for a series of ionization chambers utilized for dosimetry protocols and dosimetric quality assurance checks in clinical 1.5 T MR-Linacs. Their dependencies on irradiation setup conditions are also investigated. Moreover, to evaluate the suitability of solid phantoms for dosimetry checks and end-to-end tests, changes to the detector readings due to the presence of small asymmetrical air gaps around the detector's tip are quantified. METHODS: Phase space files for three irradiation fields of the ELEKTA Unity 1.5 T/7 MV flattening-filter-free MR-Linac were provided by the manufacturer and used as source models throughout this study. Twelve ionization chambers (three farmer-type and nine small-cavity detectors, from three manufacturers) were modeled (including their dead volume) using the EGSnrc MC code package. kB,Q values were calculated for the 10 × 10 cm2 irradiation field and for four cardinal orientations of the detectors' axes with respect to the 1.5 T magnetic field. Potential dependencies of kB,Q values with respect to field size, depth, and phantom material were investigated by performing additional simulations. Changes to the detectors' readings due to the presence of small asymmetrical air gaps (0.1 up to 1 mm) around the chambers' sensitive volume in an RW3 solid phantom were quantified for three small-cavity chambers and two orientations. RESULTS: For both parallel (to the magnetic field) orientations, kB,Q values were found close to unity. The maximum correction needed was 1.1%. For each detector studied, the kB,Q values calculated for the two parallel orientations agreed within uncertainties. Larger corrections (up to 5%) were calculated when the detectors were oriented perpendicularly to the magnetic field. Results were compared with corresponding ones found in the literature, wherever available. No considerable dependence of kB,Q with respect to field size (down to 3 × 3 cm2 ), depth, or phantom material was noticed, for the detectors investigated. As compared to the perpendicular one, in the parallel to the magnetic field orientation, the air gap effect is minimized but is still considerable even for the smallest air gap considered (0.1 mm). CONCLUSION: For the 10 × 10 cm2 field, magnetic field correction factors for 12 ionization chambers and four orientations were determined. For each detector, the kB,Q value may be also applied for dosimetry procedures under different irradiation parameters provided that the orientation is taken into account. Moreover, if solid phantoms are used, even the smallest asymmetrical air gap may still bias small-cavity chamber response. This work substantially expands the availability and applicability of kB,Q correction factors that are detector- and orientation-specific, enabling more options in MR-Linac dosimetry checks, end-to-end tests, and quality assurance protocols.

Assessment of Radiation-Induced Bladder and Bowel Cancer Risks after Conventionally and Hypo-Fractionated Radiotherapy for the Preoperative Management of Rectal Carcinoma.

Preoperative management of rectal carcinoma can be performed by employing either conventionally or hypo-fractionated Radiotherapy (CFRT or HFRT, respectively), delivered by Intensity Modulated Radiotherapy (IMRT) or Volumetric Modulated Arc Therapy (VMAT) plans, employing 6 MV or 10 MV photon beams. This study aims to dosimetrically and radiobiologically compare all available approaches, with emphasis on the risk of radiation-induced second cancer to the bladder and bowel. Computed Tomography (CT) scans and relevant radiotherapy contours from 16 patients were anonymized and analyzed retrospectively. For each case, CFRT of 25 × 2 Gy and HFRT of 5 × 5 Gy were both considered. IMRT and VMAT plans using 6 MV and 10 MV photons were prepared. Plan optimization was performed, considering all clinically used plan quality indices and dose-volume constraints for the critical organs. Resulting dose distributions were analyzed and compared. Moreover, the Lifetime Attributable Risk (LAR) for developing radiation-induced bladder and bowel malignancies were assessed using a non-linear mechanistic model, assuming patient ages at treatment of 45, 50, 55 and 60 years. All 128 plans created were clinically acceptable. Risk of second bladder cancer reached 0.26% for HFRT (5 × 5 Gy) and 0.19% for CFRT (25 × 2 Gy) at the age of 45. Systematically higher risks were calculated for HFRT (5 × 5 Gy) as compared to CFRT (25 × 2 Gy), with 6 MV photons resulting in slightly increased LAR, as well. Similar or equal bowel cancer risks were calculated for all techniques and patient ages investigated (range 0.05-0.14%). This work contributes towards radiotherapy treatment protocol selection criteria for the preoperative irradiation of rectal carcinoma. However, more studies are needed to establish the associated radiation-induced risk of each RT protocol.

Geometric distortion assessment in 3T MR images used for treatment planning in cranial Stereotactic Radiosurgery and Radiotherapy.

Magnetic Resonance images (MRIs) are employed in brain Stereotactic Radiosurgery and Radiotherapy (SRS/SRT) for target and/or critical organ localization and delineation. However, MRIs are inherently distorted, which also impacts the accuracy of the Magnetic Resonance Imaging/Computed Tomography (MRI/CT) co-registration process. In this phantom-based study, geometric distortion is assessed in 3T T2-weighted images (T2WIs), while the efficacy of an MRI distortion correction technique is also evaluated. A homogeneous polymer gel-filled phantom was CT-imaged before being irradiated with 26 4-mm Gamma Knife shots at predefined locations (reference control points). The irradiated phantom was MRI-scanned at 3T, implementing a T2-weighted protocol suitable for SRS/SRT treatment planning. The centers of mass of all shots were identified in the 3D image space by implementing an iterative localization algorithm and served as the evaluated control points for MRI distortion detection. MRIs and CT images were spatially co-registered using a mutual information algorithm. The inverse transformation matrix was applied to the reference control points and compared with the corresponding MRI-identified ones to evaluate the overall spatial accuracy of the MRI/CT dataset. The mean image distortion correction technique was implemented, and resulting MRI-corrected control points were compared against the corresponding reference ones. For the scanning parameters used, increased MRI distortion (>1mm) was detected at areas distant from the MRI isocenter (>5cm), while median radial distortion was 0.76mm. Detected offsets were slightly higher for the MRI/CT dataset (0.92mm median distortion). The mean image distortion correction improves geometric accuracy, but residual distortion cannot be considered negligible (0.51mm median distortion). For all three datasets studied, a statistically significant positive correlation between detected spatial offsets and their distance from the MRI isocenter was revealed. This work contributes towards the wider adoption of 3T imaging in SRS/SRT treatment planning. The presented methodology can be employed in commissioning and quality assurance programmes of corresponding treatment workflows.

Monte Carlo-Based Radiobiological Investigation of the Most Optimal Ion Beam Forming SOBP for Particle Therapy.

Proton (p) and carbon (C) ion beams are in clinical use for cancer treatment, although other particles such as He, Be, and B ions have more recently gained attention. Identification of the most optimal ion beam for radiotherapy is a challenging task involving, among others, radiobiological characterization of a beam, which is depth-, energy-, and cell type- dependent. This study uses the FLUKA and MCDS Monte Carlo codes in order to estimate the relative biological effectiveness (RBE) for several ions of potential clinical interest such as p, 4He, 7Li, 10Be, 10B, and 12C forming a spread-out Bragg peak (SOBP). More specifically, an energy spectrum of the projectiles corresponding to a 5-cm SOBP at a depth of 8 cm was used. All secondary particles produced by the projectiles were considered and RBE was determined based on radiation-induced Double Strand Breaks (DSBs), as calculated by MCDS. In an attempt to identify the most optimal ion beam, using the latter data, biological optimization was performed and the obtained depth-dose distributions were inter-compared. The results showed that 12C ions are more effective inside the SOBP region, which comes at the expense of higher dose values at the tail (i.e., after the SOBP). In contrast, p beams exhibit a higher DSOPB/DEntrance ratio, if physical doses are considered. By performing a biological optimization in order to obtain a homogeneous biological dose (i.e., dose × RBE) in the SOBP, the corresponding advantages of p and 12C ions are moderated. 7Li ions conveniently combine a considerably lower dose tail and a DSOPB/DEntrance ratio similar to 12C. This work contributes towards identification of the most optimal ion beam for cancer therapy. The overall results of this work suggest that 7Li ions are of potential interest, although more studies are needed to demonstrate the relevant advantages. Future work will focus on studying more complex beam configurations.

Assessment of sequence dependent geometric distortion in contrast-enhanced MR images employed in stereotactic radiosurgery treatment planning.

This work focuses on MR-related sequence dependent geometric distortions, which are associated with B 0 inhomogeneity and patient-induced distortion (susceptibility differences and chemical shift effects), in MR images used in stereotactic radiosurgery (SRS) applications. Emphasis is put on characterizing distortion at target brain areas identified by gadolinium diethylenetriamine pentaacetic acid (Gd-DTPA) paramagnetic contrast agent uptake. A custom-made phantom for distortion detection was modified to accommodate two small cylindrical inserts, simulating small brain targets. The inserts were filled with Gd-DTPA solutions of various concentrations (0-20 mM). The phantom was scanned at 1.5 T unit using both the reversed read gradient polarity (to determine the overall distortion as reflected by the inserts centroid offset) and the field mapping (to determine B 0 inhomogeneity related distortion in the vicinity of the inserts) techniques. Post-Gd patient images involving a total of 10 brain metastases/targets were also studied using a similar methodology. For the specific imaging conditions, contrast agent presence was found to evidently affect phantom insert position, with centroid offset extending up to 0.068 mm mM-1 (0.208 ppm mM-1). The Gd-DTPA induced distortion in patient images was of the order of 0.5 mm for the MRI protocol used, in agreement with the phantom results. Total localization uncertainty of metastases-targets in patient images ranged from 0.35 mm to 0.87 mm, depending on target location, with an average value of 0.54 mm (2.24 ppm). This relative wide range of target localization uncertainty results from the fact that the B 0 inhomogeneity distortion vector in a specific location may add to or partly counterbalance Gd-DTPA induced distortion, thus increasing or decreasing, respectively, the total sequence dependent distortion. Although relatively small, the sequence dependent distortion in Gd-DTPA enhanced brain images can be easily taken into account for SRS treatment planning and target definition purposes by carefully inspecting both the forward and reversed polarity series.

On the use of a novel Ferrous Xylenol-orange gelatin dosimeter for HDR brachytherapy commissioning and quality assurance testing.

PURPOSE: To evaluate a commercially available Ferrous-Xylenol Orange-Gel (FXG) dosimeter (TrueView™) coupled with Optical-Computed Tomography (OCT) read out, for 3D dose verification in an Ir-192 superficial brachytherapy application. METHODS: Two identical polyethylene containers filled with gel from the same batch were used. One was irradiated with an 18 MeV electron field to examine the dose-response linearity and obtain a calibration curve. A flap surface applicator was attached to the other to simulate treatment of a skin lesion. The dose distribution in the experimental set up was calculated with the TG-43 and the model based dose calculation (MBCA) algorithms of a commercial treatment planning system (TPS), as well as Monte Carlo (MC) simulation using the MCNP code. Measured and calculated dose distributions were spatially registered and compared. RESULTS: Apart from a region close to the container's neck, where gel measurements exhibited an over-response relative to MC calculations (probably due to stray light perturbation), an excellent agreement was observed between measurements and simulations. More than 97% of points within the 10% isodose line (80 cGy) met the gamma index criteria established from uncertainty analysis (5%/2 mm). The corresponding passing rates for the comparison of experiment to calculations using the TG-43 and MBDCA options of the TPS were 57% and 92%, respectively. CONCLUSION: TrueView™ is suitable for the quality assurance of demanding radiotherapy applications. Experimental results of this work confirm the advantage of the studied MBDCA over TG-43, expected from the improved account of scatter radiation in the treatment geometry.

Time resolved dose rate distributions in brachytherapy.

PURPOSE: To investigate the biological significance of introducing time-resolved dose rate distributions (TR-DRD) in brachytherapy. MATERIALS AND METHODS: The treatment plan of a head and neck patient treated with pulsed-dose-rate (PDR) brachytherapy was considered. The TR-DRD was calculated on the basis of a Monte Carlo generated single source dose rate matrix taking into account the dose rate per source dwell position. Biologically Effective Dose (BED) was obtained considering either the mean dose rate per pulse (analytical method) or the TR-DRD (numerical method). Corresponding Tumor Control Probabilities (TCP) were calculated and compared for various PDR schemes and repair half-times from the literature. The dose of the biologically equivalent high-dose-rate (HDR) treatment schedule was also evaluated. RESULTS: The analytical method presents an overall BED underestimation (up to 2%) relative to TR-DRD results. This is associated with an analytical-based TCP underestimation which increases with dose/pulse, pulse duration and period time and decreases with total dose. The half-time of repair seems to have the largest impact on the TCP calculations, with significant differences (up to 39.1%) corresponding to the shorter repair half-times. Regarding the equivalent HDR treatment schedule, the analytical method resulted to a HDR isoeffective dose underestimation lower than 2.2% and thus does not warrant any change in the derivation of the equivalent HDR scheme. CONCLUSION: TR-DRD data should be taken into account for PDR biological effectiveness estimations, especially for short tissue repair half-times. This does not appear however to influence dose prescription of the equivalent HDR treatment schedule for mobile tongue carcinoma.

On the experimental validation of model-based dose calculation algorithms for 192Ir HDR brachytherapy treatment planning.

There is an acknowledged need for the design and implementation of physical phantoms appropriate for the experimental validation of model-based dose calculation algorithms (MBDCA) introduced recently in 192Ir brachytherapy treatment planning systems (TPS), and this work investigates whether it can be met. A PMMA phantom was prepared to accommodate material inhomogeneities (air and Teflon), four plastic brachytherapy catheters, as well as 84 LiF TLD dosimeters (MTS-100M 1 × 1 × 1 mm3 microcubes), two radiochromic films (Gafchromic EBT3) and a plastic 3D dosimeter (PRESAGE). An irradiation plan consisting of 53 source dwell positions was prepared on phantom CT images using a commercially available TPS and taking into account the calibration dose range of each detector. Irradiation was performed using an 192Ir high dose rate (HDR) source. Dose to medium in medium, [Formula: see text], was calculated using the MBDCA option of the same TPS as well as Monte Carlo (MC) simulation with the MCNP code and a benchmarked methodology. Measured and calculated dose distributions were spatially registered and compared. The total standard (k = 1) spatial uncertainties for TLD, film and PRESAGE were: 0.71, 1.58 and 2.55 mm. Corresponding percentage total dosimetric uncertainties were: 5.4-6.4, 2.5-6.4 and 4.85, owing mainly to the absorbed dose sensitivity correction and the relative energy dependence correction (position dependent) for TLD, the film sensitivity calibration (dose dependent) and the dependencies of PRESAGE sensitivity. Results imply a LiF over-response due to a relative intrinsic energy dependence between 192Ir and megavoltage calibration energies, and a dose rate dependence of PRESAGE sensitivity at low dose rates (<1 Gy min-1). Calculations were experimentally validated within uncertainties except for MBDCA results for points in the phantom periphery and dose levels <20%. Experimental MBDCA validation is laborious, yet feasible. Further work is required for the full characterization of dosimeter response for 192Ir and the reduction of experimental uncertainties.

On the impact of improved dosimetric accuracy on head and neck high dose rate brachytherapy.

PURPOSE: To study the effect of finite patient dimensions and tissue heterogeneities in head and neck high dose rate brachytherapy. METHODS AND MATERIALS: The current practice of TG-43 dosimetry was compared to patient specific dosimetry obtained using Monte Carlo simulation for a sample of 22 patient plans. The dose distributions were compared in terms of percentage dose differences as well as differences in dose volume histogram and radiobiological indices for the target and organs at risk (mandible, parotids, skin, and spinal cord). RESULTS: Noticeable percentage differences exist between TG-43 and patient specific dosimetry, mainly at low dose points. Expressed as fractions of the planning aim dose, percentage differences are within 2% with a general TG-43 overestimation except for the spine. These differences are consistent resulting in statistically significant differences of dose volume histogram and radiobiology indices. Absolute differences of these indices are however small to warrant clinical importance in terms of tumor control or complication probabilities. CONCLUSIONS: The introduction of dosimetry methods characterized by improved accuracy is a valuable advancement. It does not appear however to influence dose prescription or call for amendment of clinical recommendations for the mobile tongue, base of tongue, and floor of mouth patient cohort of this study.

A user-oriented procedure for the commissioning and quality assurance testing of treatment planning system dosimetry in high-dose-rate brachytherapy.

PURPOSE: To develop a user-oriented procedure for testing treatment planning system (TPS) dosimetry in high-dose-rate brachytherapy, with particular focus to TPSs using model-based dose calculation algorithms (MBDCAs). METHODS AND MATERIALS: Identical plans were prepared for three computational models using two commercially available systems and the same (192)Ir source. Reference dose distributions were obtained for each plan using the MCNP v.6.1 Monte Carlo (MC) simulation code with input files prepared via automatic parsing of plan information using a custom software tool. The same tool was used for the comparison of reference dose distributions with corresponding MBDCA exports. RESULTS: The single source test case yielded differences due to the MBDCA spatial discretization settings. These affect points at relatively increased distance from the source, and they are abated in test cases with multiple source dwells. Differences beyond MC Type A uncertainty were also observed very close to the source(s), close to the test geometry boundaries, and within heterogeneities. Both MBDCAs studied were found equivalent to MC within 5 cm from the target volume for a clinical breast brachytherapy test case. These are in agreement with previous findings of MBDCA benchmarking in the literature. CONCLUSIONS: The data and the tools presented in this work, that are freely available via the web, can serve as a benchmark for advanced clinical users developing their own tests, a complete commissioning procedure for new adopters of currently available TPSs using MBDCAs, a quality assurance testing tool for future updates of already installed TPSs, or as an admission prerequisite in multicentric clinical trials.

BrachyGuide: a brachytherapy-dedicated DICOM RT viewer and interface to Monte Carlo simulation software.

This work presents BrachyGuide, a brachytherapy-dedicated software tool for the automatic preparation of input files for Monte Carlo simulation from treatment plans exported in DICOM RT format, and results of calculations performed for its benchmarking. Three plans were prepared using two computational models, the image series of a water sphere and a multicatheter breast brachytherapy patient, for each of two commercially available treatment planning systems: BrachyVision and Oncentra Brachy. One plan involved a single source dwell position of an 192Ir HDR source (VS2000 or mHDR-v2) at the center of the water sphere using the TG43 algorithm, and the other two corresponded to the TG43 and advanced dose calculation algorithm for the multicatheter breast brachytherapy patient. Monte Carlo input files were prepared using BrachyGuide and simulations were performed with MCNP v.6.1. For the TG43 patient plans, the Monte Carlo computational model was manually edited in the prepared input files to resemble TG43 dosimetry assumptions. Hence all DICOM RT dose exports were equivalent to corresponding simulation results and their comparison was used for benchmarking the use of BrachyGuide. Monte Carlo simulation results and corresponding DICOM RT dose exports agree within type A uncertainties in the majority of points in the computational models. Treatment planning system, algorithm, and source specific differences greater than type A uncertainties were also observed, but these were explained by treatment planning system-related issues and other sources of type B uncertainty. These differences have to be taken into account in commissioning procedures of brachytherapy dosimetry algorithms. BrachyGuide is accurate and effective for use in the preparation of commissioning tests for new brachytherapy dosimetry algorithms as a user-oriented commissioning tool and the expedition of retrospective patient cohort studies of dosimetry planning.