M. A. Duch

Comparison of dose calculation algorithms in phantoms with lung equivalent heterogeneities under conditions of lateral electronic disequilibrium

An extensive set of benchmark measurement of PDDs and beam profiles was performed in a heterogeneous layer phantom, including a lung equivalent heterogeneity, by means of several detectors and compared against the predicted dose values by different calculation algorithms in two treatment planning systems. PDDs were measured with TLDs, plane parallel and cylindrical ionization chambers and beam profiles with films. Additionally, Monte Carlo simulations by means of the PENELOPE code were performed. Four different field sizes (10 x 10, 5 x 5, 2 x 2, and 1 x 1 cm2) and two lung equivalent materials (CIRS, p(w)e=0.195 and St. Bartholomew Hospital, London, p(w)e=0.244-0.322) were studied. The performance of four correction-based algorithms and one based on convolution-superposition was analyzed. The correction-based algorithms were the Batho, the Modified Batho, and the Equivalent TAR implemented in the Cadplan (Varian) treatment planning system and the TMS Pencil Beam from the Helax-TMS (Nucletron) treatment planning system. The convolution-superposition algorithm was the Collapsed Cone implemented in the Helax-TMS. The only studied calculation methods that correlated successfully with the measured values with a 2% average inside all media were the Collapsed Cone and the Monte Carlo simulation. The biggest difference between the predicted and the delivered dose in the beam axis was found for the EqTAR algorithm inside the CIRS lung equivalent material in a 2 x 2 cm2 18 MV x-ray beam. In these conditions, average and maximum difference against the TLD measurements were 32% and 39%, respectively. In the water equivalent part of the phantom every algorithm correctly predicted the dose (within 2%) everywhere except very close to the interfaces where differences up to 24% were found for 2 x 2 cm2 18 MV photon beams. Consistent values were found between the reference detector (ionization chamber in water and TLD in lung) and Monte Carlo simulations, yielding minimal differences (0.4%+/-1.2%). The penumbra broadening effect in low density media was not predicted by any of the correction-based algorithms, and the only one that matched the experimental values and the Monte Carlo simulations within the estimated uncertainties was the Collapsed Cone Algorithm.

SBRT of lung tumours: Monte Carlo simulation with PENELOPE of dose distributions including respiratory motion and comparison with different treatment planning systems.

The purpose of this work was to simulate with the Monte Carlo (MC) code PENELOPE the dose distribution in lung tumours including breathing motion in stereotactic body radiation therapy (SBRT). Two phantoms were modelled to simulate a pentagonal cross section with chestwall (unit density), lung (density 0.3 g cm(-3)) and two spherical tumours (unit density) of diameters respectively of 2 cm and 5 cm. The phase-space files (PSF) of four different SBRT field sizes of 6 MV from a Varian accelerator were calculated and used as beam sources to obtain both dose profiles and dose-volume histograms (DVHs) in different volumes of interest. Dose distributions were simulated for five beams impinging on the phantom. The simulations were conducted both for the static case and including the influence of respiratory motion. To reproduce the effect of breathing motion different simulations were performed keeping the beam fixed and displacing the phantom geometry in chosen positions in the cranial and caudal and left-right directions. The final result was obtained by combining the different position with two motion patterns. The MC results were compared with those obtained with three commercial treatment planning systems (TPSs), two based on the pencil beam (PB) algorithm, the TMS-HELAX (Nucletron, Sweden) and Eclipse (Varian Medical System, Palo Alto, CA), and one based on the collapsed cone algorithm (CC), Pinnacle(3) (Philips). Some calculations were also carried out with the analytical anisotropic algorithm (AAA) in the Eclipse system. All calculations with the TPSs were performed without simulated breathing motion, according to clinical practice. In order to compare all the TPSs and MC an absolute dose calibration in Gy/MU was performed. The analysis shows that the dose (Gy/MU) in the central part of the gross tumour volume (GTV) is calculated for both tumour sizes with an accuracy of 2-3% with PB and CC algorithms, compared to MC. At the periphery of the GTV the TPSs overestimate the dose up to 10%, while in the lung tissue close to the GTV PB algorithms overestimate the dose and the CC underestimates it. When clinically relevant breathing motions are included in the MC simulations, the static calculations with the TPSs still give a relatively accurate estimate of the dose in the GTV. On the other hand, the dose at the periphery of the GTV is overestimated, compared to the static case.

Dose distributions in SBRT of lung tumors: Comparison between two different treatment planning algorithms and Monte-Carlo simulation including breathing motions

In stereotactic body radiotherapy (SBRT) of lung tumors, dosimetric problems arise from: 1) the limited accuracy in the dose calculation algorithms in treatment planning systems, and 2) the motions with the respiration of the tumor during treatment. Longitudinal dose distributions have been calculated with Monte Carlo simulation (MC), a pencil beam algorithm (PB) and a collapsed cone algorithm (CC) for two spherical lung tumors (2 cm and 5 cm diameter) in lung tissue, in a phantom situation. Respiratory motions were included by a convolution method, which was validated. In the static situation, the PB significantly overestimates the dose, relative to MC, while the CC gives a relatively accurate estimate. Four different respiratory motion patterns were included in the dose calculation with the MC. A “narrowing” of the longitudinal dose profile of up to 20 mm (at about 90% dose level) is seen relative the static dose profile calculated with the PB.

Comparison of dose calculation algorithms in slab phantoms with cortical bone equivalent heterogeneities

To evaluate the dose values predicted by several calculation algorithms in two treatment planning systems, Monte Carlo (MC) simulations and measurements by means of various detectors were performed in heterogeneous layer phantoms with water- and bone-equivalent materials. Percentage depth doses (PDDs) were measured with thermoluminescent dosimeters (TLDs), metal-oxide semiconductor field-effect transistors (MOSFETs), plane parallel and cylindrical ionization chambers, and beam profiles with films. The MC code used for the simulations was the PENELOPE code. Three different field sizes (10 x 10, 5 x 5, and 2 x 2 cm2) were studied in two phantom configurations and a bone equivalent material. These two phantom configurations contained heterogeneities of 5 and 2 cm of bone, respectively. We analyzed the performance of four correction-based algorithms and one based on convolution superposition. The correction-based algorithms were the Batho, the Modified Batho, the Equivalent TAR implemented in the Cadplan (Varian) treatment planning system (TPS), and the Helax-TMS Pencil Beam from the Helax-TMS (Nucletron) TPS. The convolution-superposition algorithm was the Collapsed Cone implemented in the Helax-TMS. All the correction-based calculation algorithms underestimated the dose inside the bone-equivalent material for 18 MV compared to MC simulations. The maximum underestimation, in terms of root-mean-square (RMS), was about 15% for the Helax-TMS Pencil Beam (Helax-TMS PB) for a 2 x 2 cm2 field inside the bone-equivalent material. In contrast, the Collapsed Cone algorithm yielded values around 3%. A more complex behavior was found for 6 MV where the Collapsed Cone performed less well, overestimating the dose inside the heterogeneity in 3%-5%. The rebuildup in the interface bone-water and the penumbra shrinking in high-density media were not predicted by any of the calculation algorithms except the Collapsed Cone, and only the MC simulations matched the experimental values within the estimated uncertainties. The TLD and MOSFET detectors were suitable for dose measurement inside bone-equivalent materials, while parallel ionization chambers, applying the same calibration and correction factors as in water, systematically underestimated dose by 3%-5%.

Monte Carlo simulation of MOSFET detectors for high-energy photon beams using the PENELOPE code.

The aim of this work was the Monte Carlo (MC) simulation of the response of commercially available dosimeters based on metal oxide semiconductor field effect transistors (MOSFETs) for radiotherapeutic photon beams using the PENELOPE code. The studied Thomson&Nielsen TN-502-RD MOSFETs have a very small sensitive area of 0.04 mm(2) and a thickness of 0.5 microm which is placed on a flat kapton base and covered by a rounded layer of black epoxy resin. The influence of different metallic and Plastic water build-up caps, together with the orientation of the detector have been investigated for the specific application of MOSFET detectors for entrance in vivo dosimetry. Additionally, the energy dependence of MOSFET detectors for different high-energy photon beams (with energy >1.25 MeV) has been calculated. Calculations were carried out for simulated 6 MV and 18 MV x-ray beams generated by a Varian Clinac 1800 linear accelerator, a Co-60 photon beam from a Theratron 780 unit, and monoenergetic photon beams ranging from 2 MeV to 10 MeV. The results of the validation of the simulated photon beams show that the average difference between MC results and reference data is negligible, within 0.3%. MC simulated results of the effect of the build-up caps on the MOSFET response are in good agreement with experimental measurements, within the uncertainties. In particular, for the 18 MV photon beam the response of the detectors under a tungsten cap is 48% higher than for a 2 cm Plastic water cap and approximately 26% higher when a brass cap is used. This effect is demonstrated to be caused by positron production in the build-up caps of higher atomic number. This work also shows that the MOSFET detectors produce a higher signal when their rounded side is facing the beam (up to 6%) and that there is a significant variation (up to 50%) in the response of the MOSFET for photon energies in the studied energy range. All the results have shown that the PENELOPE code system can successfully reproduce the response of a detector with such a small active area.

Measurements of eye lens doses in interventional cardiology using OSL and electronic dosemeters

The purpose of this paper is to test the appropriateness of OSL and electronic dosemeters to estimate eye lens doses at interventional cardiology environment. Using TLD as reference detectors, personal dose equivalent was measured in phantoms and during clinical procedures. For phantom measurements, OSL dose values resulted in an average difference of -15 % vs. TLD. Tests carried out with other electronic dosemeters revealed differences up to ±20 % versus TLD. With dosemeters positioned outside the goggles and when TLD doses were >20 μSv, the average difference OSL vs. TLD was -9 %. Eye lens doses of almost 700 μSv per procedure were measured in two cases out of a sample of 33 measurements in individual clinical procedures, thus showing the risk of high exposure to the lenses of the eye when protection rules are not followed. The differences found between OSL and TLD are acceptable for the purpose and range of doses measured in the survey.

Influence of dosemeter position for the assessment of eye lens dose during interventional cardiology

The equivalent dose limit for the eye lens for occupational exposure recommended by the ICRP has been reduced to 20 mSv y(-1) averaged over defined periods of 5 y, with no single year exceeding 50 mSv. The compliance with this new requirement could not be easy in some workplace such as interventional radiology and cardiology. The aim of this study is to evaluate different possible approaches in order to have a good estimate of the eye lens dose during interventional procedures. Measurements were performed with an X-ray system Philips Allura FD-10, using a PMMA phantom to simulate the patient scattered radiation and a Rando phantom to simulate the cardiologist. Thermoluminescence (TL) whole-body and TL eye lens dosemeters together with Philips DoseAware active dosemeters were located on different positions of the Rando phantom to estimate the eye lens dose in typical cardiology procedures. The results show that, for the studied conditions, any of the analysed dosemeter positions are suitable for eye lens dose assessment. However, the centre of the thyroid collar and the left ear position provide a better estimate. Furthermore, in practice, improper use of the ceiling-suspended screen can produce partial protection of some parts of the body, and thus large differences between the measured doses and the actual exposure of the eye could arise if the dosemeter is not situated close to the eye.

An algorithm to assess the need for clinical Monte Carlo dose calculation for small proton therapy fields based on quantification of tissue heterogeneity

PURPOSE In proton therapy, complex density heterogeneities within the beam path constitute a challenge to dose calculation algorithms. This might question the reliability of dose distributions predicted by treatment planning systems based on analytical dose calculation. For cases in which substantial dose errors are expected, resorting to Monte Carlo dose calculations might be essential to ensure a successful treatment outcome and therefore the benefit is worth a presumably long computation time. The aim of this study was to define an indicator for the accuracy of dose delivery based on analytical dose calculations in treatment planning systems for small proton therapy fields to identify those patients for which Monte Carlo dose calculation is warranted. METHODS Fourteen patients treated at our facility with small passively scattered proton beams (apertures diameters below 7 cm) were selected. Plans were generated in the XiO treatment planning system in combination with a pencil beam algorithm developed at the Massachusetts General Hospital and compared to Monte Carlo dose calculations. Differences in the dose to the 50% of the gross tumor volume (D50, GTV) were assessed in a field-by-field basis. A simple and fast methodology was developed to quantify the inhomogeneity of the tissue traversed by a single small proton beam using a heterogeneity index (HI)-a concept presented by Plugfelder et al. [Med. Phys. 34, 1506-1513 (2007)] for scanned proton beams. Finally, the potential correlation between the error made by the pencil beam based treatment planning algorithm for each field and the level of tissue heterogeneity traversed by the proton beam given by the HI was evaluated. RESULTS Discrepancies up to 5.4% were found in D50 for single fields, although dose differences were within clinical tolerance levels (<3%) when combining all of the fields involved in the treatment. The discrepancies found for each field exhibited a strong correlation to their associated HI-values (Spearmans ρ=0.8, p<0.0001); the higher the level of tissue inhomogeneities for a particular field, the larger the error by the analytical algorithm. With the established correlation a threshold for HI can be set by choosing a tolerance level of 2-3%-commonly accepted in radiotherapy. CONCLUSIONS The HI is a good indicator for the accuracy of proton field delivery in terms of GTV prescription dose coverage when small fields are delivered. Each HI-value was obtained from the CT image in less than 3 min on a computer with 2 GHz CPU allowing implementation of this methodology in clinical routine. For HI-values exceeding the threshold, either a change in beam direction (if feasible) or a recalculation of the dose with Monte Carlo would be highly recommended.

EYE LENS DOSE IN INTERVENTIONAL CARDIOLOGY

The ICRP has recently recommended reducing the occupational exposure dose limit for the lens of the eye to 20 mSv y(-1), averaged over a period of 5 y, with no year exceeding 50 mSv, instead of the current 150 mSv y(-1). This reduction will have important implications for interventional cardiology and radiology (IC/IR) personnel. In this work, lens dose received by a staff working in IC is studied in order to determine whether eye lens dose monitoring or/and additional radiological protection measures are required. Eye lens dose exposure was monitored in 10 physicians and 6 nurses. The major IC procedures performed were coronary angiography and percutaneous transluminal coronary angioplasty. The personnel were provided with two thermoluminescent dosemeters (TLDs): one calibrated in terms of Hp(3) located close to the left ear of the operator and a whole-body dosemeter calibrated in terms of Hp(10) and Hp(0.07) positioned on the lead apron. The estimated annual eye lens dose for physicians ranged between 8 and 60 mSv, for a workload of 200 procedures y(-1). Lower doses were collected for nurses, with estimated annual Hp(3) between 2 and 4 mSv y(-1). It was observed that for nurses the Hp(0.07) measurement on the lead apron is a good estimate of eye lens dose. This is not the case for physicians, where the influence of both the position and use of protective devices such as the ceiling shield is very important and produces large differences among doses both at the eyes and on the thorax. For physicians, a good correlation between Hp(3) and dose area product is shown.

AIR KERMA TO Hp(3) CONVERSION COEFFICIENTS FOR IEC 61267 RQR X-RAY RADIATION QUALITIES: APPLICATION TO DOSE MONITORING OF THE LENS OF THE EYE IN MEDICAL DIAGNOSTICS

Recent studies highlight the fact that the new eye lens dose limit can be exceeded in interventional radiology procedures and that eye lens monitoring could be required for these workers. The recommended operational quantity for monitoring of eye lens exposure is the personal dose equivalent at 3 mm depth Hp(3) (ICRU 51). However, there are no available conversion coefficients in international standards, while in the literature coefficients have only been calculated for monoenergetic beams and for ISO 4037-1 X-ray qualities. The aim of this article is to provide air kerma to Hp(3) conversion coefficients for a cylindrical phantom made of ICRU-4 elements tissue-equivalent material for RQR radiation qualities (IEC-61267) from 40 to 120 kV and for angles of incidence from 0 to 180°, which are characteristic of medical workplace. Analytic calculations using interpolation techniques and Monte Carlo modelling have been compared.