Experimental validation of absorbed dose-to-medium calculation algorithms in heterogeneous media.

Objective.The aim of this work was to determine heterogeneous correction factorshQclin,Qreffclin,frefdetm,wto validate absorbed dose-to-mediumDm,Qclinm,fclincalculation algorithms from detector readings. The impact of detector orientation perpendicular and parallel to the beam central axis on the correction factors was also investigated.Approach.ThehQclin,Qreffclin,frefdetm,wfactors were calculated for four types of detectors (PTW PinPoint T31016, PTW microDiamond T60019, PTW microSilicon T60023 and EBT3 film) placed in different media (cortical bone, lung, adipose tissue, Teflon and RW3) for the 6 MV energy beam with a 10 × 10 cm2field size. These corrections were then applied to the detector measurements performed at different depths in heterogeneous phantoms.Main results.ThehQclin,Qreffclin,frefdetm,wfactors mainly depended on the media and slightly on the type of detector. Considering all detectors, the largest corrections were found in high-density media with values ranging from 0.911 to 0.934 in cortical bone. For comparison, the corrections in other media were closer to unity with values from 0.966 (lung and RW3) to 0.991 (adipose tissue). Except for the PinPoint T31016, detector orientation-dependence was observed especially in high-density media. A good agreement (≤1.5%) was found betweenDm,Qclinm,fclincalculations and the detector readings corrected with thehQclin,Qreffclin,frefdetm,wfactor for all studied heterogeneous phantoms.Significance.This paper could serve as an initial guideline for medical physicists involved in the validation of the advanced type-b dose calculation algorithms reportingDm,Qclinm,fclin. To our knowledge, this is the first study to assess the impact of the orientation of different detectors in heterogeneous media. The orientation dependence of the detector response observed in water may not reflect what is observed in heterogeneous media, especially in high-density media. The knowledge of thehQclin,Qreffclin,frefdetm,wfactors becomes mandatory for accurate interpretation of detector readings and comparisons withDm,Qclinm,fclincalculations.

Universal evaluation of MLC models in treatment planning systems based on a common set of dynamic tests.

PURPOSE: To demonstrate the feasibility of characterising MLCs and MLC models implemented in TPSs using a common set of dynamic beams. MATERIALS AND METHODS: A set of tests containing synchronous (SG) and asynchronous sweeping gaps (aSG) was distributed among twenty-five participating centres. Doses were measured with a Farmer-type ion chamber and computed in TPSs, which provided a dosimetric characterisation of the leaf tip, tongue-and-groove, and MLC transmission of each MLC, as well as an assessment of the MLC model in each TPS. Five MLC types and four TPSs were evaluated, covering the most frequent combinations used in radiotherapy departments. RESULTS: Measured differences within each MLC type were minimal, while large differences were found between MLC models implemented in clinical TPSs. This resulted in some concerning discrepancies, especially for the HD120 and Agility MLCs, for which differences between measured and calculated doses for some MLC-TPS combinations exceeded 10%. These large differences were particularly evident for small gap sizes (5 and 10 mm), as well as for larger gaps in the presence of tongue-and-groove effects. A much better agreement was found for the Millennium120 and Halcyon MLCs, differences being within ± 5% and ± 2.5%, respectively. CONCLUSIONS: The feasibility of using a common set of tests to assess MLC models in TPSs was demonstrated. Measurements within MLC types were very similar, but TPS dose calculations showed large variations. Standardisation of the MLC configuration in TPSs is necessary. The proposed procedure can be readily applied in radiotherapy departments and can be a valuable tool in IMRT and credentialing audits.

Practical and technical key challenges in head and neck adaptive radiotherapy: The GORTEC point of view.

Anatomical variations occur during head and neck (H&N) radiotherapy (RT) treatment. These variations may result in underdosage to the target volume or overdosage to the organ at risk. Replanning during the treatment course can be triggered to overcome this issue. Due to technological, methodological and clinical evolutions, tools for adaptive RT (ART) are becoming increasingly sophisticated. The aim of this paper is to give an overview of the key steps of an H&N ART workflow and tools from the point of view of a group of French-speaking medical physicists and physicians (from GORTEC). Focuses are made on image registration, segmentation, estimation of the delivered dose of the day, workflow and quality assurance for an implementation of H&N offline and online ART. Practical recommendations are given to assist physicians and medical physicists in a clinical workflow.

Field output correction factors and electron fluence perturbation of the microSilicon and microSilicon X detectors.

Objective.The aim of this study was to determine field output correction factorskQclin,Qreffclin,frefand electron fluence perturbation for new PTW unshielded microSilicon and shielded microSilicon X detectors.Approach.kQclin,Qreffclin,freffactors were calculated for 6 and 10 MV with and without flattening filter beams delivered by a TrueBeam STx. Correction factors were determined for field sizes ranging from 0.5 × 0.5 cm2to 3 × 3 cm2using both experimental and numerical methods. To better understand the underlying physics of their response, total electron (+positron) fluence spectra were scored in the sensitive volume considering the various component-dependent perturbations.Main results.The microSilicon and microSilicon X detectors can be used down to the smallest studied field size by applying corrections factors fulfilling the tolerance of 5% recommended by the IAEA TRS483. Electron fluence perturbation in both microSilicon detectors was greater than that in water but to a lesser extent than their predecessors. The main contribution of the overall perturbation of the detectors comes from the materials surrounding their sensitive volume, especially the epoxy in the case of unshielded diodes and the shielding for shielded diodes. This work demonstrated that the decrease in the density of the epoxy for the microSilicon led to a decrease in the electron fluence perturbation.Significance.A real improvement was observed regarding the design of the microSilicon and microSilicon X detectors compared to their predecessors.

Dosimetric leaf gap and leaf trailing effect in a double-stacked multileaf collimator.

PURPOSE: To investigate (i) the dosimetric leaf gap (DLG) and the effect of the "trailing distance" between leaves from different multileaf collimator (MLC) layers in Halcyon systems and (ii) the ability of the currently available treatment planning systems (TPSs) to approximate this effect. METHODS: DICOM plans with transmission beams and sweeping gap tests were created in Python for measuring the DLG for each MLC layer independently and for both layers combined. In clinical Halcyon plans both MLC layers are interchangeably used and leaves from different layers are offset, thus forming a trailing pattern. To characterize the impact of such configuration, new tests called "trailing sweeping gaps" were designed and created where the leaves from one layer follow the leaves from the other layer at a fixed "trailing distance" t between the tips. Measurements were carried out on five Halcyons SX2 from different institutions and calculations from both the Eclipse and RayStation TPSs were compared with measurements. RESULTS: The dose accumulated during a sweeping gap delivery progressively increased with the trailing distance t . We call this "the trailing effect." It is most pronounced for t between 0 and 5 mm, although some changes were obtained up to 20 mm. The dose variation was independent of the gap size. The measured DLG values also increased with t up to 20 mm, again with the steepest variation between 0 and 5 mm. Measured DLG values were negative at t  = 0 (the leaves from both layers at the same position) but changed sign for t  ≥ 1 mm, in line with the positive DLG sign usually observed with single-layer rounded-end MLCs. The Eclipse TPS does not explicitly model the leaf tip and, as a consequence, could not predict the dose reduction due to the trailing effect. This resulted in dose discrepancies up to +10% and -8% for the 5 mm sweeping gap and up to ±5% for the 10 mm one depending on the distance t . RayStation implements a simple model of the leaf tip that was able to approximate the trailing effect and improved the agreement with measured doses. In particular, with a prototype version of RayStation that assigned a higher transmission at the leaf tip the agreement with measured doses was within ±3% even for the 5 mm gap. The five Halcyon systems behaved very similarly but differences in the DLG around 0.2 mm were found across different treatment units and between MLC layers from the same system. The DLG for the proximal layer was consistently higher than for the distal layer, with differences ranging between 0.10 mm and 0.24 mm. CONCLUSIONS: The trailing distance between the leaves from different layers substantially affected the doses delivered by sweeping gaps and the measured DLG values. Stacked MLCs introduce a new level of complexity in TPSs, which ideally need to implement an explicit model of the leaf tip in order to reproduce the trailing effect. Dynamic tests called "trailing sweeping gaps" were designed that are useful for characterizing and commissioning dual-layer MLC systems.

On the need for tuning the dosimetric leaf gap for stereotactic treatment plans in the Eclipse treatment planning system.

The dosimetric leaf gap (DLG) and tongue-and-groove (T&G) effects are critical aspects in the modeling of multileaf collimators (MLC) in the treatment planning system (TPS). In this study, we investigated the dosimetric impact of limitations associated with the T&G modeling in stereotactic plans and its relationship with the need for tuning the DLG in the Eclipse TPS. Measurements were carried out using Varian TrueBeam STx systems from two different institutions. Test fields presenting MLC patterns with several MLC gap sizes (meanGap) and different amounts of T&G effect (TGi) were first evaluated. Secondly, dynamic conformal arc (DCA) and volumetric modulated arc therapy (VMAT) deliveries of stereotactic cases were analyzed in terms of meanGap and TGi. Two DLG values were used in the TPS: the measured DLG (DLGmeas ) and an optimal DLG (DLGopt ). Measured and calculated doses were compared according to dose differences and gamma passing rates (GPR) with strict local gamma criteria of 2%/2 mm. The discrepancies were analyzed for DLGmeas and DLGopt , and their relationships with both TGi and meanGap were investigated. DCA arcs involved significantly lower TGi and larger meanGap than VMAT arcs (P < 0.0001). By using DLGmeas in the TPS, the dose discrepancies increased as TGi increased and meanGap decreased for both test fields and clinical plans. Dose discrepancies dramatically increased with the ratio TGi/meanGap. Adjusting the DLG value was then required to achieve acceptable calculations and configuring the TPS with DLGopt led to an excellent agreement with median GPRs (2%/2 mm) > 99% for both institutions. We also showed that DLGopt could be obtained from the results of the test fields. We demonstrated that the need for tuning the DLG is due to the limitations of the T&G modeling in the Eclipse TPS. A set of sweeping gap tests modified to incorporate T&G effects can be used to determine the optimal DLG value.

A new method for modelling the tongue-and-groove in treatment planning systems.

Commercial TPSs typically model the tongue-and-groove (TG) by extending the projections of the leaf sides by a certain constant width. However, this model may produce discrepancies of as much as 7%-10% in the calculated average doses, especially for the High Definition multi-leaf collimator (MLC) (Hernandez et al 2017 Phys. Med. Biol. 62 6688-707). The purpose of the present study is to introduce and validate a new method for modelling the TG that uses a non constant TG width. We provide the theoretical background and a detailed methodology to determine the optimal shape of this TG width from measurements and we fit an empirical function to the TG width that depended on two parameters [Formula: see text] and [Formula: see text]. Parameter [Formula: see text] represents the TG width and [Formula: see text] introduces a curvature correction in the width near the leaf tip end. The new TG model was implemented in MATLAB and when the curvature correction was zero ([Formula: see text]) it caused the same discrepancies as the constant width model used by the Eclipse TPS. On the other hand, when the experimentally determined [Formula: see text] was used the new model's calculations were in close agreement with measurements, with all differences in average doses [Formula: see text]1%. Additionally, film dosimetry was used to successfully validate the potential of the new TG model to recreate the fine spatial details associated to TG effects. We also showed that the parameters [Formula: see text], [Formula: see text] depend solely on the MLC design by evaluating three different linear accelerators for each MLC model considered, namely Varian's High Definition and Millennium120 MLCs. In conclusion, a new method was presented that greatly improves the TG modelling. The present method can be easily implemented in commercial TPSs and has the potential to further increase their accuracy, especially for MLCs with rounded leaf ends.

Dosimetrical and radiobiological approach to manage the dosimetric shift in the transition of dose calculation algorithm in radiation oncology: how to improve high quality treatment and avoid unexpected outcomes?

BACKGROUND: For a given prescribed dose of radiotherapy, with the successive generations of dose calculation algorithms, more monitor units (MUs) are generally needed. This is due to the implementation of successive improvements in dose calculation: better heterogeneity correction and more accurate estimation of secondary electron transport contribution. More recently, there is the possibility to report the dose-to-medium, physically more accurate compared to the dose-to-water as the reference one. This last point is a recent concern and the main focus of this study. METHODS: In this paper, we propose steps for a general analysis procedure to estimate the dosimetric alterations, and the potential clinical changes, between a reference algorithm and a new one. This includes dosimetric parameters, gamma index, radiobiology indices based on equivalent uniform dose concept and statistics with bootstrap simulation. Finally, we provide a general recommendation on the clinical use of new algorithms regarding the dose prescription or dose limits to the organs at risks. RESULTS: The dosimetrical and radiobiological data showed a significant effect, which might exceed 5-10%, of the calculation method on the dose the distribution and clinical outcomes for lung cancer patients. Wilcoxon signed rank paired comparisons indicated that the delivered dose in MUs was significantly increased (> 2%) using more advanced dose calculation methods as compared to the reference one. CONCLUSION: This paper illustrates and explains the use of dosimetrical, radiobiologcal and statistical tests for dosimetric comparisons in radiotherapy. The change of dose calculation algorithm may induce a dosimetric shift, which has to be evaluated by the physicists and the oncologists. This includes the impact on tumor control and on the risk of toxicity based on normal tissue dose constraints. In fact, the alteration in dose distribution makes it hard to keep exactly the same tumor control probability along with the same normal tissue complication probability.

Monte Carlo simulation of TrueBeam flattening-filter-free beams using varian phase-space files: comparison with experimental data.

PURPOSE: Phase-space files for Monte Carlo simulation of the Varian TrueBeam beams have been made available by Varian. The aim of this study is to evaluate the accuracy of the distributed phase-space files for flattening filter free (FFF) beams, against experimental measurements from ten TrueBeam Linacs. METHODS: The phase-space files have been used as input in PRIMO, a recently released Monte Carlo program based on the PENELOPE code. Simulations of 6 and 10 MV FFF were computed in a virtual water phantom for field sizes 3 × 3, 6 × 6, and 10 × 10 cm(2) using 1 × 1 × 1 mm(3) voxels and for 20 × 20 and 40 × 40 cm(2) with 2 × 2 × 2 mm(3) voxels. The particles contained in the initial phase-space files were transported downstream to a plane just above the phantom surface, where a subsequent phase-space file was tallied. Particles were transported downstream this second phase-space file to the water phantom. Experimental data consisted of depth doses and profiles at five different depths acquired at SSD = 100 cm (seven datasets) and SSD = 90 cm (three datasets). Simulations and experimental data were compared in terms of dose difference. Gamma analysis was also performed using 1%, 1 mm and 2%, 2 mm criteria of dose-difference and distance-to-agreement, respectively. Additionally, the parameters characterizing the dose profiles of unflattened beams were evaluated for both measurements and simulations. RESULTS: Analysis of depth dose curves showed that dose differences increased with increasing field size and depth; this effect might be partly motivated due to an underestimation of the primary beam energy used to compute the phase-space files. Average dose differences reached 1% for the largest field size. Lateral profiles presented dose differences well within 1% for fields up to 20 × 20 cm(2), while the discrepancy increased toward 2% in the 40 × 40 cm(2) cases. Gamma analysis resulted in an agreement of 100% when a 2%, 2 mm criterion was used, with the only exception of the 40 × 40 cm(2) field (∼95% agreement). With the more stringent criteria of 1%, 1 mm, the agreement reduced to almost 95% for field sizes up to 10 × 10 cm(2), worse for larger fields. Unflatness and slope FFF-specific parameters are in line with the possible energy underestimation of the simulated results relative to experimental data. CONCLUSIONS: The agreement between Monte Carlo simulations and experimental data proved that the evaluated Varian phase-space files for FFF beams from TrueBeam can be used as radiation sources for accurate Monte Carlo dose estimation, especially for field sizes up to 10 × 10 cm(2), that is the range of field sizes mostly used in combination to the FFF, high dose rate beams.