Cone beam CT doses in radiotherapy patient positioning in Finland-prostate treatments.

Imaging parameters, frequencies and resulting patient organ doses in treatments of prostate cancer were assessed in Finnish radiotherapy centres. Based on a questionnaire to the clinics, Monte Carlo method was used to estimate organ doses in International Commission on Radiological Protection standard phantom for prostate, bladder, rectum and femoral head. The results show that doses from cone beam computed tomography imaging have reduced compared to earlier studies and are between 3.6 and 34.5 mGy per image for the above-mentioned organs and for normal sized patients. There still is room for further optimization of the patient exposure, as many centres use the default imaging parameters, and the length of the imaged region may not be optimal for the purpose.

Validation of HDR brachytherapy doses in the treatment of keloid scars using the egs_brachy Monte Carlo application.

Objective. Radiotherapy is a well-known alternative in the treatment of keloid scars to reduce the recurrence of scars. The purpose of this study was to investigate the feasibility and accuracy of dose delivered from a high-dose-rate (HDR) afterloaders in keloid scar brachytherapy using Monte Carlo (MC) simulations and measurements.Approach. Treatment doses and central axis dose profiles were measured using radiophotoluminescence dosimeters and radiochromic films, respectively, with two HDR afterloaders, both using an Ir-192 source, in a phantom made of solid water and polycarbonate sheets. The nominal treatment dose calculated by the AAPM Task Group No. 43 (TG-43) dose model was set to 8.5 Gy at a distance of 0.5 cm laterally from the middle of the source line located in a plastic applicator simulating a 15 cm long surgically removed scar treatment with 30 equally spaced (0.5 cm) source positions. The dose profiles were measured at three different distances from the applicator and the absolute doses at four points at different distances. MC simulations were performed using the egs_brachy, which is based on EGSnrc code system.Main results. The measured and simulated dose profiles match well, especially at 10.0 mm (difference <1%) and 15.0 mm depths (difference <4%), and with a small dose difference at 5.0 mm depth (difference <4%). Point dose measurements agreed well in the dose maximum area (difference <7%) with the simulated dose profiles, although the largest difference near the edge of the profile was <30%. The dose differences between the TG-43 dose model and the MC simulation were small (differences <4%).Significance. Simulated and measured dose levels at a depth of 0.5 cm showed that the nominal treatment dose can be achieved with the utilized setup. The measurement results of the absolute dose agree well with the corresponding simulation results.

Profeel-An open source dosimetry data visualization and analysis software.

BACKGROUND AND OBJECTIVES: This article presents Profeel, a Matlab (MathWorks Inc., MA) based open source dosimetry data visualization and analysis software. Profeel aims to support quality assurance, dosimetry and research in the field of radiotherapy by providing an environment to visualize, process and analyse measured and simulated dosimetry data from several data sources used in radiotherapy practice and research. METHODS: The processing and analysis tools are based on routinely used dosimetry analysis methods, such as gamma analysis, different data normalizations and data filtering. Additionally the Profeel performs an automatic 1 dimensional profile and percentage depth dose analysis in accordance with International Electrotechnical Commission definitions. All data can be operated by user created custom functions and lower dimensionality data can be extracted from volume doses and dose planes. RESULTS: Profeel supports data import in all 3 dimensions and offers an intuitive user interface to perform data visualization, processing and analysis between simulated and measured data. Profeel and its source code are distributed free of charge under the General Public Licence (GPL). CONCLUSIONS: Profeel has shown to be an agile tool for fulfilling various needs of several researchers and since Profeel is under constant development and is an open source project, community needs, issues and bug reports are taken into account in the development.

Properties of IBA Razor Nano Chamber in small-field radiation therapy using 6 MV FF, 6 MV FFF, and 10 MV FFF photon beams.

BACKGROUND: Small megavoltage photon fields are increasingly used in modern radiotherapy techniques such as stereotactic radiotherapy. Therefore, it is important to study the reliability of dosimetry in the small-field conditions. The IBA Razor Nano Chamber (Nano chamber) ionization chamber is particularly intended for small-field measurements. In this work, properties of the Nano chamber were studied with both measurements and Monte Carlo (MC) simulations. MATERIAL AND METHODS: The measurements and MC simulations were performed with 6 MV, 6 MV FFF and 10 MV FFF photon beams from the Varian TrueBeam linear accelerator. The source-to-surface distance was fixed at 100 cm. The measurements and MC simulations included profiles, percentage depth doses (PDD), and output factors (OF) in square jaw-collimated fields. The MC simulations were performed with the EGSnrc software system in a large water phantom. RESULTS: The measured profiles and PDDs obtained with the Nano chamber were compared against IBA Razor Diode, PTW microDiamond and the PTW Semiflex ionization chamber. These results indicate that the Nano chamber is a high-resolution detector and thus suitable for small field profile measurements down to field sizes 2 × 2 cm2 and appropriate for the PDD measurements. The field output correction factors kQclin, Qmsrfclin, fmsr and field OFs ΩQclin, Qmsrfclin, fmsr were determined according to TRS-483 protocol In the 6 MV FF and FFF beams, the determined correction factors kQclin, Qmsrfclin, fmsr were within 1.2% for the field sizes of 1 × 1 cm2-3 × 3 cm2 and the experimental and MC defined field output factors ΩQclin,Qmsrfclin,fmsr showed good agreement. CONCLUSION: The Nano chamber with its small cavity volume is a potential detector for the small-field dosimetry. In this study, the properties of this detector were characterized with measurements and MC simulations. The determined correction factors kQclin, Qmsrfclin, fmsr are novel results for the NC in the TrueBeam fields.

Dose-area product ratio in external small-beam radiotherapy: beam shape, size and energy dependencies in clinical photon beams.

In small-field radiotherapy (RT), a significant challenge is to define the amount of radiation dose absorbed in the patient where the quality of the beam has to be measured with high accuracy. The properties of a proposed new beam quality specifier, namely the dose-area-product ratio at 20 and 10 cm depths in water or DAPR20,10, were studied to yield more information on its feasibility over the conventional quality specifier tissue-phantom ratio or TPR20,10. The DAPR20,10may be measured with a large-area ionization chamber (LAC) instead of small volume chambers or semi-conductors where detector, beam and water phantom positioning and beam perturbations introduce uncertainties. The effects of beam shape, size and energy on the DAPR20,10were studied and it was shown that the DAPR20,10increases with increasing beam energy similarly to TPR20,10but in contrast exhibits a small beam size and shape dependence. The beam profile outside the beam limiting devices has been shown to have a large contribution to the DAPR20,10. There is potential in large area chambers to be used in DAPR measurement and its use in dosimetry of small-beam RT for beam quality measurements.

Monte Carlo investigation on the effect of air gap under bolus in post-mastectomy radiotherapy.

OBJECTIVES: To investigate the dosimetric effect of air gaps under bolus on skin dose for left-sided post-mastectomy radiotherapy with loco regional involvement. METHODS: Eight patients were planned retrospectively with volume modulated arc therapy (VMAT) and conventional static Field-in-Field (FinF) methods. Three different setups were applied for the 5-mm bolus over the chest wall having 0, 5 or 10 mm air gap under the bolus. The dose calculation was performed using Monte Carlo (MC) simulation. In addition, Analytic Anisotropic Algorithm (AAA) was used to demonstrate the differences observed in clinical setting. RESULTS: The investigated air gaps under the bolus had minimal effect on surface dose for FinF plans (relative difference ≤ 2.6%), whereas for VMAT plans the surface dose decreased 13.6% when compared to the case with no air gap. In both FinF and VMAT, the largest differences between AAA and MC were seen at the surface where AAA underestimated the dose by 1.5 Gy (p < 0.05) on average; while the dose in the target volume excluding the surface was relatively similar being on average 0.3 Gy (p > 0.05) larger with AAA than with MC calculations. CONCLUSIONS: The surface dose was significantly lower with VMAT technique than with FinF technique. Possible air gaps under the bolus reduced the surface dose significantly further for VMAT but not for FinF treatments, which may have clinical impact on recurrence rate. AAA was shown to underestimate the surface dose when compared to MC calculation.

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.

Measurement and properties of the dose-area product ratio in external small-beam radiotherapy.

In small-beam radiation therapy (RT) the measurement of the beam quality parameter, i.e. the tissue-phantom ratio or TPR20,10, using a conventional point detector is a challenge. To obtain reliable results, one has to consider potential sources of error, including volume averaging and adjustment of the point detector into the narrow beam. To overcome these challenges, a different type of beam quality parameter in small beams was studied, namely the dose-area product ratio, or DAPR20,10. With this method, the measurement of a dose-area product (DAP) using a large-area plane-parallel chamber (LAC) eliminates the uncertainties in detector positioning and volume averaging that are present when using a point detector. In this study, the properties of the DAPR20,10 of a cone-collimated 6 MV photon beam were investigated using Monte Carlo (MC) calculations and the obtained values were compared to measurements obtained using two LAC detectors, PTW Type 34073 and PTW Type 34070. In addition, the possibility of determining the DAP using EBT3 film and a Razor diode detector was studied. The determination of the DAPR20,10 value was found to be feasible in external small-beam radiotherapy using cone-collimated beams with diameters from 4-40 mm, based on the results of the two LACs, the MC calculations and the Razor diode. The measurements indicated a constant DAPR20,10 value for fields 20-40 mm in diameter, with a maximum relative change of 0.6%, but an increase of 7.0% for fields from 20-4 mm in diameter for the PTW Type 34070 chamber. Simulations and measurements showed an increase of DAPR20,10 with increasing LAC size or dose integral area for the studied 4-40 mm cone-collimated 6 MV photon beams. This has the consequence that there should be a reference to the size of the used LAC active area or the DAP integration area with the reported DAPR20,10 value.

Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range.

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4MeV to 16MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2mm with the first material set. Larger differences were observed for a large 4MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.

Gafchromic EBT3 film dosimetry in electron beams - energy dependence and improved film read-out.

For megavoltage photon radiation, the fundamental dosimetry characteristics of Gafchromic EBT3 film were determined in 60Co gamma ray beam with addition of experimental and Monte Carlo (MC)-simulated energy dependence of the film for 6 MV photon beam and 6 MeV, 9 MeV, 12 MeV, and 16 MeV electron beams in water phantom. For the film read-out, two phase correction of scanner sensitivity was applied: a matrix correction for scanning area and dose-dependent correction by iterative procedure. With these corrections, the uniformity of response can be improved to be within ± 50 pixel values (PVs). To improve the read-out accuracy, a procedure with flipped film orientations was established. With the method, scanner uniformity can be improved further and dust particles, scratches and/or dirt on scan-ner glass can be detected and eliminated. Responses from red and green channels were averaged for read-out, which decreased the effect of noise present in values from separate channels. Since the signal level with the blue channel is considerably lower than with other channels, the signal variation due to different perturbation effects increases the noise level so that the blue channel is not recommended to be used for dose determination. However, the blue channel can be used for the detection of emulsion thickness variations for film quality evaluations with unexposed films. With electron beams ranging from 6 MeV to 16 MeV and at reference measurement conditions in water, the energy dependence of the EBT3 film is uniform within 0.5%, with uncertainties close to 1.6% (k = 2). Including 6 MV photon beam and the electron beams mentioned, the energy dependence is within 1.1%. No notable differences were found between the experimental and MC-simulated responses, indicating negligible change in intrinsic energy dependence of the EBT3 film for 6 MV photon beam and 6 MeV-16 MeV electron beams. Based on the dosimetric characteristics of the EBT3 film, the read-out procedure established, the nearly uniform energy dependence found and the estimated uncertainties, the EBT3 film was concluded to be a suitable 2D dosimeter for measuring electron or mixed photon/electron dose distributions in water phantom. Uncertainties of 3.7% (k = 2) for absolute and 2.3% (k = 2) for relative dose were estimated.

Quantification of dose differences between two versions of Acuros XB algorithm compared to Monte Carlo simulations--the effect on clinical patient treatment planning.

A commercialized implementation of linear Boltzmann transport equation solver, the Acuros XB algorithm (AXB), represents a class of most advanced type 'c' photon radiotherapy dose calculation algorithms. The purpose of the study was to quantify the effects of the modifications implemented in the more recent version 11 of the AXB (AXB11) compared to the first commercial implementation, version 10 of the AXB (AXB10), in various anatomical regions in clinical treatment planning. Both versions of the AXB were part of Varian's Eclipse clinical treatment planning system and treatment plans for 10 patients were created using intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc radiotherapy (VMAT). The plans were first created with the AXB10 and then recalculated with the AXB11 and full Monte Carlo (MC) simulations. Considering the full MC simulations as reference, a DVH analysis for gross tumor and planning target volumes (GTV and PTV) and organs at risk was performed, and also 3D gamma agreement index (GAI) values within a 15% isodose region and for the PTV were determined. Although differences up to 12% in DVH analysis were seen between the MC simulations and the AXB, based on the results of this study no general conclusion can be drawn that the modifications made in the AXB11 compared to the AXB10 would imply that the dose calculation accuracy of the AXB10 would be inferior to the AXB11 in the clinical patient treatment planning. The only clear improvement with the AXB11 over the AXB10 is the dose calculation accuracy in air cavities. In general, no large deviations are present in the DVH analysis results between the two versions of the algorithm, and the results of 3D gamma analysis do not favor one or the other. Thus it may be concluded that the results of the comprehensive studies assessing the accuracy of the AXB10 may be extended to the AXB11.

The accuracy of Acuros XB algorithm for radiation beams traversing a metallic hip implant - comparison with measurements and Monte Carlo calculations.

In this study, the clinical benefit of the improved accuracy of the Acuros XB (AXB) algorithm, implemented in a commercial radiotherapy treatment planning system (TPS), Varian Eclipse, was demonstrated with beams traversing a high-Z material. This is also the first study assessing the accuracy of the AXB algorithm applying volumetric modulated arc therapy (VMAT) technique compared to full Monte Carlo (MC) simulations. In the first phase the AXB algorithm was benchmarked against point dosimetry, film dosimetry, and full MC calculation in a water-filled anthropometric phantom with a unilateral hip implant. Also the validity of the full MC calculation used as reference method was demonstrated. The dose calculations were performed both in original computed tomography (CT) dataset, which included artifacts, and in corrected CT dataset, where constant Hounsfield unit (HU) value assignment for all the materials was made. In the second phase, a clinical treatment plan was prepared for a prostate cancer patient with a unilateral hip implant. The plan applied a hybrid VMAT technique that included partial arcs that avoided passing through the implant and static beams traversing the implant. Ultimately, the AXB-calculated dose distribution was compared to the recalculation by the full MC simulation to assess the accuracy of the AXB algorithm in clinical setting. A recalculation with the anisotropic analytical algorithm (AAA) was also performed to quantify the benefit of the improved dose calculation accuracy of type 'c' algorithm (AXB) over type 'b' algorithm (AAA). The agreement between the AXB algorithm and the full MC model was very good inside and in the vicinity of the implant and elsewhere, which verifies the accuracy of the AXB algorithm for patient plans with beams traversing through high-Z material, whereas the AAA produced larger discrepancies.

Performance of dose calculation algorithms from three generations in lung SBRT: comparison with full Monte Carlo-based dose distributions.

The accuracy of dose calculation is a key challenge in stereotactic body radiotherapy (SBRT) of the lung. We have benchmarked three photon beam dose calculation algorithms--pencil beam convolution (PBC), anisotropic analytical algorithm (AAA), and Acuros XB (AXB)--implemented in a commercial treatment planning system (TPS), Varian Eclipse. Dose distributions from full Monte Carlo (MC) simulations were regarded as a reference. In the first stage, for four patients with central lung tumors, treatment plans using 3D conformal radiotherapy (CRT) technique applying 6 MV photon beams were made using the AXB algorithm, with planning criteria according to the Nordic SBRT study group. The plans were recalculated (with same number of monitor units (MUs) and identical field settings) using BEAMnrc and DOSXYZnrc MC codes. The MC-calculated dose distributions were compared to corresponding AXB-calculated dose distributions to assess the accuracy of the AXB algorithm, to which then other TPS algorithms were compared. In the second stage, treatment plans were made for ten patients with 3D CRT technique using both the PBC algorithm and the AAA. The plans were recalculated (with same number of MUs and identical field settings) with the AXB algorithm, then compared to original plans. Throughout the study, the comparisons were made as a function of the size of the planning target volume (PTV), using various dose-volume histogram (DVH) and other parameters to quantitatively assess the plan quality. In the first stage also, 3D gamma analyses with threshold criteria 3%/3mm and 2%/2 mm were applied. The AXB-calculated dose distributions showed relatively high level of agreement in the light of 3D gamma analysis and DVH comparison against the full MC simulation, especially with large PTVs, but, with smaller PTVs, larger discrepancies were found. Gamma agreement index (GAI) values between 95.5% and 99.6% for all the plans with the threshold criteria 3%/3 mm were achieved, but 2%/2 mm threshold criteria showed larger discrepancies. The TPS algorithm comparison results showed large dose discrepancies in the PTV mean dose (D50%), nearly 60%, for the PBC algorithm, and differences of nearly 20% for the AAA, occurring also in the small PTV size range. This work suggests the application of independent plan verification, when the AAA or the AXB algorithm are utilized in lung SBRT having PTVs smaller than 20-25 cc. The calculated data from this study can be used in converting the SBRT protocols based on type 'a' and/or type 'b' algorithms for the most recent generation type 'c' algorithms, such as the AXB algorithm.