Search

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

The development of new biological imaging technologies offers the opportunity to further individualize radiotherapy. Biologically conformal radiation therapy (BCRT) implies the use of the spatial distribution of one or more radiobiological parameters to guide the IMRT dose prescription. Our aim was to implement BCRT in an algorithmic segmentation-based planning approach. A biology-based segmentation tool was developed to generate initial beam segments that reflect the biological signal intensity pattern. The weights and shapes of the initial segments are optimized by means of an objective function that minimizes the root mean square deviation between the actual and intended dose values within the PTV. As proof of principle, [(18)F]FDG-PET-guided BCRT plans for two different levels of dose escalation were created for an oropharyngeal cancer patient. Both plans proved to be dosimetrically feasible without violating the planning constraints for the expanded spinal cord and the contralateral parotid gland as organs at risk. The obtained biological conformity was better for the first (2.5 Gy per fraction) than for the second (3 Gy per fraction) dose escalation level.

A new method is presented to decouple the parameters of the incident e(-) beam hitting the target of the linear accelerator, which consists essentially in optimizing the agreement between measurements and calculations when the difference filter, which is an additional filter inserted in the linac head to obtain uniform lateral dose-profile curves for the high energy photon beam, and flattening filter are removed from the beam path. This leads to lateral dose-profile curves, which depend only on the mean energy of the incident electron beam, since the effect of the radial intensity distribution of the incident e- beam is negligible when both filters are absent. The location of the primary collimator and the thickness and density of the target are not considered as adjustable parameters, since a satisfactory working Monte Carlo model is obtained for the low energy photon beam (6 MV) of the linac using the same target and primary collimator. This method was applied to conclude that the mean energy of the incident e- beam for the high energy photon beam (18 MV) of our Elekta SLi Plus linac is equal to 14.9 MeV. After optimizing the mean energy, the modelling of the filters, in accordance with the information provided by the manufacturer, can be verified by positioning only one filter in the linac head while the other is removed. It is also demonstrated that the parameter setting for Bremsstrahlung angular sampling in BEAMnrc ('Simple' using the leading term of the Koch and Motz equation or 'KM' using the full equation) leads to different dose-profile curves for the same incident electron energy for the studied 18 MV beam. It is therefore important to perform the calculations in 'KM' mode. Note that both filters are not physically removed from the linac head. All filters remain present in the linac head and are only rotated out of the beam. This makes the described method applicable for practical usage since no recommissioning process is required.

Monte Carlo based treatment planning of two different patient groups treated with step-and-shoot IMRT (head-and-neck and lung treatments) with different CT resolutions and scoring methods is performed to determine the effect of geometrical and scoring voxel sizes on DVHs and calculation times. Dose scoring is performed in two different ways: directly into geometrical voxels (or in a number of grouped geometrical voxels) or into scoring voxels defined by a separate scoring grid superimposed on the geometrical grid. For the head-and-neck cancer patients, more than 2% difference is noted in the right optical nerve when using voxel dimensions of 4 x 4 x 4 mm3 compared to the reference calculation with 1 x 1 x 2 mm3 voxel dimensions. For the lung cancer patients, 2% difference is noted in the spinal cord when using voxel dimensions of 4 x 4 x 10 mm3 compared to the 1 x 1 x 5 mm3 calculation. An independent scoring grid introduces several advantages. In cases where a relatively high geometrical resolution is required and where the scoring resolution is less important, the number of scoring voxels can be limited while maintaining a high geometrical resolution. This can be achieved either by grouping several geometrical voxels together into scoring voxels or by superimposing a separate scoring grid of spherical voxels with a user-defined radius on the geometrical grid. For the studied lung cancer cases, both methods produce accurate results and introduce a speed increase by a factor of 10-36. In cases where a low geometrical resolution is allowed, but where a high scoring resolution is required, superimposing a separate scoring grid on the geometrical grid allows a reduction in geometrical voxels while maintaining a high scoring resolution. For the studied head-and-neck cancer cases, calculations performed with a geometrical resolution of 2 x 2 x 2 mm3 and a separate scoring grid containing spherical scoring voxels with a radius of 2 mm produce accurate results and introduce a speed increase by a factor of 13. The scoring grid provides an additional degree of freedom for limiting calculation time and memory requirements by selecting optimized scoring and geometrical voxel dimensions in an independent way.

DOSSCORE is an accelerated version of DOSXYZnrc that allows photons to cross voxel boundaries of the same medium and utilizes a separate scoring grid superimposed on the geometrical grid. Two different stepping algorithms, the hownear method and the scaling method are implemented in DOSSCORE. The hownear method allows particles to travel larger distances in homogeneous regions where there is no interest in the dose deposition of these particles, whilst the scaling method utilizes a stepping algorithm in which particles are only slowed down by the boundaries of the geometrical voxels and not by the boundaries of the scoring voxels. For CT-based phantoms, only photon ray tracing is applied, which results in a rather modest speed gain of factor 1.2 compared to DOSXYZnrc. The hownear method and scaling method do not increase the speed for CT-based phantoms, but only for homogeneous phantoms and phantoms with a limited number of small heterogeneities. In cases where a small number of scoring voxels are needed, the hownear method performs better than the scaling method, whilst the opposite is true for cases when many scoring voxels are needed. The photon transport is accelerated by almost a factor of 2 for all phantoms (homogeneous, heterogeneous with much homogeneity and CT-based phantoms) compared to DOSXYZnrc. For a small number of scoring voxels, the hownear method is up to a factor of 2.6 and 1.9 faster than DOSXYZnrc for homogeneous and heterogeneous phantoms in the case of photon beams. For an electron beam, a speed gain of factor 2.4 is obtained. For a full scoring grid like the one used in DOSXYZnrc, the scaling method is up to a factor of 2.2 and 1.7 faster than DOSXYZnrc for homogeneous and heterogeneous phantoms in the case of photon beams. For an electron beam, a speed gain of factor 2 is obtained. A speed increase without biasing the results is very relevant. The use of two separate grids, the more efficient stepping algorithms and the accelerated photon transport can be applied to every EGS-based or other Monte Carlo code.

A new component module (CM), called MLCE, has been implemented in the BEAM program. The CM takes into account the particular 'tongue-and-groove' design of the Elekta multi-leaf collimator (MLC) and the air gap between the leaves. The model was validated by two series of measurements and simulations. The first benchmarking series focuses on the interleaf leakage and the intraleaf transmission. The measurement showed a total transmission through the MLC of 1.42% of the open field dose. Two Monte Carlo (MC) simulations were made, the first with the new CM MLCE (inclusive of air gap) and the second with the CM MLCQ (exclusive of air gap), which is available in the BEAM distribution. When the air gap between the leaves was determined by varying the parameters of the leaf geometry within tolerance limits on the technical drawing, the total measured transmission of 1.42% was well reproduced by the CM MLCE. In contrast, MC simulations with MLCQ showed that the transmission through the MLC calculated without the interleaf leakage is only 44% of the total transmitted radiation. The relevance of the detailed MLC modelling was demonstrated also by studying the 'adjacent' tongue-and-groove effect, where two adjacent (not opposing) leaves are complementary, opened or closed. The two complementary leaf settings were simulated both with the CM MLCE and MLCQ. A comparison with measurements was made. In regions covered by two or more leaves, the dose increased by 14% for two leaves and by 40% for more than two leaves when the interleaf leakage was included in the transmission. The tongue-and-groove effect was perfectly reproduced by the MLCE module.

In recent years, a change has been proposed from air kerma based reference dosimetry to absorbed dose based reference dosimetry for all radiotherapy beams of ionizing radiation. In this paper, a dosimetry study is presented in which absorbed dose based dosimetry using recently developed formalisms was compared with air kerma based dosimetry using older formalisms. Three ionization chambers of each of three different types were calibrated in terms of absorbed dose to water and air kerma and sent to five hospitals. There, reference dosimetry with all the chambers was performed in a total of eight high-energy clinical photon beams. The selected chamber types were the NE2571, the PTW-30004 and the Wellhöfer-FC65G (previously Wellhöfer-IC70). Having a graphite wall, they exhibit a stable volume and the presence of an aluminium electrode ensures the robustness of these chambers. The data were analysed with the most important recommendations for clinical dosimetry: IAEA TRS-398, AAPM TG-51, IAEA TRS-277, NCS report-2 (presently recommended in Belgium) and AAPM TG-21. The necessary conversion factors were taken from those protocols, or calculated using the data in the different protocols if data for a chamber type are lacking. Polarity corrections were within 0.1% for all chambers in all beams. Recombination corrections were consistent with theoretical predictions, did not vary within a chamber type and only slightly between different chamber types. The maximum chamber-to-chamber variations of the dose obtained with the different formalisms within the same chamber type were between 0.2% and 0.6% for the NE2571, between 0.2% and 0.6% for the PTW-30004 and 0.1% and 0.3% for the Wellhöfer-FC65G for the different beams. The absorbed dose results for the NE2571 and Wellhöfer-FC65G chambers were in good agreement for all beams and all formalisms. The PTW-30004 chambers gave a small but systematically higher result compared to the result for the NE2571 chambers (on the average 0.1% for IAEA TRS-277, 0.3% for NCS report-2 and AAPM TG-21 and 0.4% for IAEA TRS-398 and AAPM TG-51). Within the air kerma based protocols, the results obtained with the TG-21 protocol were 0.4-0.8% higher mainly due to the differences in the data used. Both absorbed dose to water based formalisms resulted in consistent values within 0.3%. The change from old to new formalisms is discussed together with the traceability of calibration factors obtained at the primary absorbed dose and air kerma standards in the reference beams (60Co). For the particular situation in Belgium (calibrations at the Laboratory for Standard Dosimetry of Ghent) the change amounts to 0.1-0.6%. This is similar to the magnitude of the change determined in other countries.

Current dosimetry protocols for clinical protons using ionization chambers do not take into account ionization chamber-dependent perturbation correction factors. In the present investigation, the relative response of 17 cylindrical ionization chambers was evaluated at three proton beam qualities: at two points in a modulated beam and one point in a non-modulated beam, all with an incident energy of 75 MeV. Thirteen of the ionization chambers had a Farmer-type geometry but consisted of different combinations of wall and central electrode materials. All ionization chambers were calibrated in terms of air kerma as well as in terms of absorbed dose to water in a 60Co beam. The relative response of the ionization chambers was compared with results of Monte Carlo simulations of proton and secondary electron transport in the phantom and the ionization chamber geometry. The results of the measurements for cylindrical ionization chambers show relative perturbation effects that are limited to 0.5-1%, resulting in perturbation correction factors that are larger than unity compared with an NE2571 ionization chamber. The experimental relative wall and total perturbation correction factors agree with Monte Carlo calculated values, indicating that the differences between the responses of different ionization chambers are due to secondary electron effects. This conclusion is supported by the comparison of our results with those from other investigators after re-analysis of their data. The central electrode perturbation correction factor for an aluminium electrode in a Farmer-type geometry was found to be unity within the experimental uncertainties.

An in-phantom calibration technique for 192Ir sources used for endovascular brachytherapy is presented. Three different source lengths were investigated. The calibration was performed in a solid phantom using a Farmer-type ionization chamber at source to detector distances ranging from 1 cm to 5 cm. The dosimetry protocol for medium-energy x-rays extended with a volume-averaging correction factor was used to convert the chamber reading to dose to water. The air kerma strength of the sources was determined as well. EGS4 Monte Carlo calculations were performed to determine the depth dose distribution at distances ranging from 0.6 mm to 10 cm from the source centre. In this way we were able to convert the absolute dose rate at 1 cm distance to the reference point chosen at 2 mm distance. The Monte Carlo results were confirmed by radiochromic film measurements, performed with a double-exposure technique. The dwell times to deliver a dose of 14 Gy at the reference point were determined and compared with results given by the source supplier (CORDIS). They determined the dwell times from a Sievert integration technique based on the source activity. The results from both methods agreed to within 2% for the 12 sources that were evaluated. A Visual Basic routine that superimposes dose distributions, based on the Monte Carlo calculations and the in-phantom calibration, onto intravascular ultrasound images is presented. This routine can be used as an online treatment planning program.

Dose to water (Dw) determination in clinical high-energy photon beams with ionization chambers calibrated in terms of absorbed dose to water has been proposed as an alternative to ionization chamber dosimetry based on air kerma calibrations. Dw in the clinical beam is derived using a kappaQ factor that scales the absorbed dose calibration factor in the reference beam to the absorbed dose calibration factor in the user beam. In the present study kappaQ values were determined for the NE2571 chamber in a 5 MV and a 10 MV high-energy photon beam generated at the 15 MeV high-intensity electron linac of the University of Gent. A set of three NE2571 chambers was calibrated relative to the Gent sealed water calorimeter both in 60Co and in the linac beam at a depth of 5 cm and a source to detector distance of 100 cm. Two high-purity chemical water systems were used in the detection vessel of the calorimeter, H2-saturated and Ar-saturated pure water, which are both supposed to give a zero heat defect. TPR20(10) and %dd(10) have been evaluated as beam quality specifiers. Simulations using the BEAM/DOSXYZ Monte Carlo system were performed to evaluate potential corrections on the measured beam qualities. The average kappaQ values measured for the three NE2571 chambers in the 5 MV and 10 MV photon beams are 0.995 +/- 0.005 and 0.979 +/- 0.005 respectively. For the three chambers used, the maximum deviation of individual kappaQ values is 0.2%. The measured beam quality specifiers %dd(10) and TPR20(10) are 67.0 and 0.705 for the 5 MV beam and 75.0 and 0.759 for the 10 MV beam. Although our beam design is very different from those used by other investigators for the measurement of kappaQ values, the agreement with their results is satisfactory showing a slightly better agreement when %dd(10) is used as the beam quality specifier.

An ionometric calibration procedure for 192Ir PDR brachytherapy sources in terms of dose rate to water is presented. The calibration of the source is performed directly in a water phantom at short distances (1.0, 2.5 and 5.0 cm) using an NE2571 Farmer type ion chamber. To convert the measured air-kerma rate in water to dose rate to water a conversion factor (CF) was calculated by adapting the medium-energy x-ray dosimetry protocol for a point source geometry (diverging beam). The obtained CF was verified using two different methods. Firstly, the CF was calculated by Monte Carlo simulations, where the source-ionization chamber geometry was modelled accurately. In a second method, a combination of Monte Carlo simulations and measurements of the air-kerma rate in water (at 1.0, 2.5 and 5.0 cm distance) and in air (1 m distance) was used to determine the CF. The obtained CFs were also compared with conversion factors calculated with the adapted dosimetry protocol for high-energy photons introduced by Tölli. All calculations were done for a Gammamed PDR 192Ir source-NE2571 chamber geometry. The conversion factors obtained with the four different methods agree to within 1% at the three distances of interest. We obtained the following values (medium-energy x-ray protocol): CF(1 cm) = 1.458; CF(2.5 cm) = 1.162; CF(5.0 cm) = 1.112 (1 sigma = 0.7% for the three distances of interest). The obtained results were checked with TLD measurements. The values of the specific dose rate constant and the radial dose function calculated in this work are in accordance with the literature data.