Frederik Crop

The influence of small field sizes, penumbra, spot size and measurement depth on perturbation factors for microionization chambers

The purpose of this study was the investigation of perturbation factors for microionization chambers in small field dosimetry and the influence of penumbra for different spot sizes. To this purpose, correlated sampling was implemented in the EGSnrc Monte Carlo (MC) user code cavity: CScavity. CScavity was first benchmarked against results in the literature for an NE2571 chamber. An efficiency increase of 17 was attained for the calculation of a realistic chamber perturbation factor in a water phantom. Calculations have been performed for microionization chambers of type PinPoint 31006 and PinPoint 31016 in full BEAMnrc linac simulations. Investigating the physical backgrounds of the differences for these small field settings, perturbation factors have been split up into (1) central electrode perturbation, (2) wall perturbation, (3) air-to-water perturbation (chamber volume air-to-water) and (4) water volume perturbation (water chamber volume to 1 mm(3) voxel). The influence of different spot sizes, position in penumbra, measuring depth and detector geometry on these perturbation factors has been investigated, in a 0.8 x 0.8 cm(2) field setting. p(cel) for the PP31006 steel electrode shows a variation of up to 1% in the lateral position, but only 0.4% for the PP31016 with an Al electrode. The air-to-water perturbation in the optimal scanning direction for both profiles and depth is most influenced by the radiation field, and only to a small extent the chamber geometry. The PP31016 geometry (shorter, larger radius) requires less total perturbation within the central axis of the field, but results in slightly larger variations off axis in the optimal scanning direction. Smaller spot sizes (0.6 mm FWHM) and sharper penumbras, compared to larger spot sizes (2 mm FWHM), result in larger perturbation starting in the penumbra. The longer geometries of the PP31006/14/15 exhibit in the non-optimal scanning direction large variations in total perturbation (p(tot) 1.201(4) (0.6 mm spot, 3 mm off axis, type A MC uncertainty) to 0.803(4) (5 mm off axis)) mainly due to volume perturbation. Therefore in IMRT settings, when the detector is not always in the optimal scanning direction, the PP31016 geometry requires less extreme perturbation (max p(tot) 1.130(3)) and shows less variation. However, these results suggest that small variations in positioning, spot size or MLC result in large differences in perturbation factors. Therefore even these 0.016 cm(3) ionization chambers are limited in their use for a field setting of 0.8 x 0.8 cm(2), as used in this investigation.

Monte Carlo modeling of the ModuLeaf miniature MLC for small field dosimetry and quality assurance of the clinical treatment planning system

The purpose of this investigation was the verification of both the measured data and quality of the implementation of the add-on ModuLeaf miniature multileaf collimator (ML mMLC) into the clinical treatment planning system for conformal stereotactic radiosurgery treatment. To this end the treatment head with ML mMLC was modeled in the BEAMnrc Monte Carlo (MC) code. The 6 MV photon beams used in the setup were first benchmarked with a set of measurements. A total ML mMLC transmission of 1.13% of the 10 x 10 cm2 open field dose was measured and reproduced with the BEAMnrc/DOSXYZnrc code. Correspondence between calculated and measured output factors (OFs) was within 2%. Correspondence between MC and measured profiles was within 2% dose and 2 mm distance, only for the smallest 0.5 x 0.5 cm2 field the results were within 3% dose. In the next step, the MC model was compared with Gafchromic film measurements and Pinnacle(3) 7.4 f (convolution superposition algorithm) calculated dose distributions, using a gamma evaluation comparison, for a multi-beam patient setup delivered to a Lucytrade mark phantom. The gamma evaluation of the MC versus Gafchromic film resulted in 3.4% of points not fulfilling gamma <or= 1 for a 2%/2 mm criterion, the Pinnacle(3) 7.4 f versus Gafchromic results 3.8% and Pinnacle versus MC less than 1%. For specific patients with lesions of 8 cc and 0.2 cc, Monte Carlo and Pinnacle simulations of the plans were performed and compared using DVH evaluation. DVHs corresponded within 2% dose and 2% volume.

Use of a liquid ionization chamber for stereotactic radiotherapy dosimetry

Liquid ionization chambers (LICs) offer an interesting tool in the field of small beam dosimetry, allowing better spatial resolution and reduced perturbation effects. However, some aspects remain to be addressed, such as the higher recombination and the effects from the materials of the detector. Our aim was to investigate these issues and their impact. The first step was the evaluation of the recombination effects. Measurements were performed at different SSDs to vary the dose per pulse, and the collection efficiency was obtained. The BEAMnrc code was then used to model the Cyberknife head. Finally, the liquid ionization chamber itself was modelled using the EGSnrc-based code Cavity allowing the evaluation of the influence of the volume and the chamber materials. The liquid ionization charge collection efficiency is approximately 0.98 at 1.5 mGy pulse(-1), the highest dose per pulse that we have measured. Its impact on the accuracy of output factors is less than half a per cent. The detector modelling showed a significant contribution from the graphite electrode, up to 6% for the 5 mm collimator. The dependence of the average electronic mass collision stopping power of iso-octane with beam collimation is negligible and thus has no influence on output factor measurements. Finally, the volume effect reaches 5% for the small 5 mm collimator and becomes much smaller (<0.5%) for diameters above 10 mm. LICs can effectively be used for small beam relative dosimetry as long as adequate correction factors are applied, especially for the electrode and volume effects.

Improving dose calculations on tomotherapy MVCT images.

The purpose of this investigation was the creation of a new protocol allowing more precise dose calculations on megavoltage CT (MVCT) images for tomotherapy, both for adaptive and StatRT planning. Daily MVCT images offer, next to positioning purposes, the possibility for daily dose check and adaptive planning. Dose calculations use the image value to density table (IVDT) to calculate physical densities from Hounsfield Units (HUs). These measured HUs change over time, leading to a dose calculation error. We noticed dose calculation variations due to IVDT changes of: 0.2% dose during a day, up to 1.6% dose from long‐term variations, and up to 1.5% dose due to technical interventions. An analysis was performed applying the general methodology of a calibration problem. A model HU=bρc‐1020 was obtained using a weighted least squares inverse prediction method (HU as function of density) taking into account the heteroscedasticity. The b parameter is the major variable and depends also on the dose rate (DR). We demonstrate the correction for DR variations and the constance of the c parameter. Instead of scanning the whole tissue characterization phantom daily, we propose a simplified daily protocol: (a) morning airscan‐like procedure with only two inserts on the table (defining b and thus the IVDT curve), (b) DR variations throughout the day can be corrected for using the DR model. A patient‐specific protocol for which two inserts next to the patient are scanned could also be used, but results in equal uncertainties and is less practical. Therefore we recommend the morning procedure with dose rate variation correction. Applying the proposed transformations and the model, the correct IVDT of the moment can be reconstructed, with a simple measurement in the morning, and corrected with DR changes during the day. This corresponds with a linear mapping in time of the proposed IVDT function. The dosimetric variation is hereby reduced from up to 3% to 0.4 % for the tested pelvic and head‐and‐neck cases. In practice, several IVDT curves corresponding to “b” values can be entered. The correct IVDT curve of that moment can then be chosen from the list. Instead of the two high‐density inserts on table, any calibrated single density phantom could be used in order to create the IVDT curve of the day, but it should have a larger size than the current inserts. PACS number: 87.55.Gh

Measurement Techniques, Modeling Strategies and Pitfalls to Avoid when Implementing a Mini MLC in a Non Dedicated Planning System*

Background and Purpose:Ghent University Hospital investigated the feasibility of the Pinnacle® system for planning intracranial stereotactic treatments. The aim was to perform precise dose computation using the collapsed cone engine for treatment delivery with the Moduleaf mini-MLC mounted on an Elekta accelerator.Material and Methods:The Moduleaf® was commissioned using dose rate corrected data recorded by a diamond detector and using data measured by cylindrical chambers each limited to restricted field sizes.Results:Automatic modeling resulted in clinical relevant dose errors up to 10%. Using manual modeling in Pinnacle®, for clinical applicable fields a 2%/2 mm agreement between modeled data and measurements was obtained.Conclusion:The overall accuracy of the collapsed cone algorithm is within tolerances for single fraction stereotactic treatments.Hintergrund und Ziel:Am Universitätsklinikum Gent wurde die Anwendung eines Pinnacle®-Planungssystems für die intrakranielle stereotaktische Bestrahlung untersucht. Das Ziel bestand darin, für die Bestrahlung mit einem Moduleaf-Mini-MLC, der an einem Elekta-Beschleuniger® befestigt ist, mit dem „Collapsed cone“-Rechenalgorithmus eine präzise Dosisberechnung zu erstellen.Material und Methoden:Hierzu wurden Messungen mit einem Diamantdetektor durchgeführt, die bezüglich der Dosisleistung korrigiert wurden, sowie Messungen mit einer zylindrischen Ionisationskammer für kleine Felder.Ergebnisse:Die automatische Modellierung führte zu Fehlern von bis zu 10%, was bereits klinisch relevant ist. Mit der manuellen Modellierung in Pinnacle war es möglich, zwischen modellierten Daten und Messungen eine Übereinstimmung von 2%/2 mm zu erreichen.Schlussfolgerung:Die Gesamtgenauigkeit des „Collapsed-cone“-Algorithmus liegt innerhalb der Toleranzgrenzen für stereotaktische Einzeldosisbehandlungen.

On the calibration process of film dosimetry: OLS inverse regression versus WLS inverse prediction.

The purpose of this study was both putting forward a statistically correct model for film calibration and the optimization of this process. A reliable calibration is needed in order to perform accurate reference dosimetry with radiographic (Gafchromic) film. Sometimes, an ordinary least squares simple linear (in the parameters) regression is applied to the dose-optical-density (OD) curve with the dose as a function of OD (inverse regression) or sometimes OD as a function of dose (inverse prediction). The application of a simple linear regression fit is an invalid method because heteroscedasticity of the data is not taken into account. This could lead to erroneous results originating from the calibration process itself and thus to a lower accuracy. In this work, we compare the ordinary least squares (OLS) inverse regression method with the correct weighted least squares (WLS) inverse prediction method to create calibration curves. We found that the OLS inverse regression method could lead to a prediction bias of up to 7.3 cGy at 300 cGy and total prediction errors of 3% or more for Gafchromic EBT film. Application of the WLS inverse prediction method resulted in a maximum prediction bias of 1.4 cGy and total prediction errors below 2% in a 0-400 cGy range. We developed a Monte-Carlo-based process to optimize calibrations, depending on the needs of the experiment. This type of thorough analysis can lead to a higher accuracy for film dosimetry.

Characterization of Recombination Effects in a Liquid Ionization Chamber Used for the Dosimetry of a Radiosurgical Accelerator

Most modern radiation therapy devices allow the use of very small fields, either through beamlets in Intensity-Modulated Radiation Therapy (IMRT) or via stereotactic radiotherapy where positioning accuracy allows delivering very high doses per fraction in a small volume of the patient. Dosimetric measurements on medical accelerators are conventionally realized using air-filled ionization chambers. However, in small beams these are subject to nonnegligible perturbation effects. This study focuses on liquid ionization chambers, which offer advantages in terms of spatial resolution and low fluence perturbation. Ion recombination effects are investigated for the microLion detector (PTW) used with the Cyberknife system (Accuray). The method consists of performing a series of water tank measurements at different source-surface distances, and applying corrections to the liquid detector readings based on simultaneous gaseous detector measurements. This approach facilitates isolating the recombination effects arising from the high density of the liquid sensitive medium and obtaining correction factors to apply to the detector readings. The main difficulty resides in achieving a sufficient level of accuracy in the setup to be able to detect small changes in the chamber response.

On the conversion of dose to bone to dose to water in radiotherapy treatment planning systems

Background and purpose Conversion factors between dose to medium (Dm,m) and dose to water (Dw,w) provided by treatment planning systems that model the patient as water with variable electron density are currently based on stopping power ratios. In the current paper it will be illustrated that this conversion method is not correct. Materials and methods Monte Carlo calculations were performed in a phantom consisting of a 2 cm bone layer surrounded by water. Dw,w was obtained by modelling the bone layer as water with the electron density of bone. Conversion factors between Dw,w and Dm,m were obtained and compared to stopping power ratios and ratios of mass-energy absorption coefficients in regions of electronic equilibrium and interfaces. Calculations were performed for 6 MV and 20 MV photon beams. Results In the region of electronic equilibrium the stopping power ratio of water to bone (1.11) largely overestimates the conversion obtained using the Monte Carlo calculations (1.06). In that region the MC dose conversion corresponds to the ratio of mass energy absorption coefficients. Near the water to bone interface, the MC ratio cannot be determined from stopping powers or mass energy absorption coefficients. Conclusion Stopping power ratios cannot be used for conversion from Dm,m to Dw,w provided by treatment planning systems that model the patient as water with variable electron density, either in regions of electronic equilibrium or near interfaces. In regions of electronic equilibrium mass energy absorption coefficient ratios should be used. Conversions at interfaces require detailed MC calculations.