Real-time dose reconstruction for wedged photon beams: a generalized procedure.

A practical and accurate generalized procedure to reconstruct the isocenter dose D(iso) for 3D conformal radiotherapy (3DCRT) has been developed for X-ray open beams supplied by linacs of different manufacturers and equipped with aSi electronic portal imaging devices (aSi EPIDs). This paper reports an extension of the method, to be applied at the wedged X-ray beams characterized by the wedge attenuation factor W(AF). Using water-equivalent solid phantoms (SPs) of different thicknesses, w, and photon square fields of sizes, L, the generalized midplane doses D(0)(W(AF), w/2,L) and generalized transit signals s(t)(0)(W(AF),w,L) by 38 beams of six different linacs were determined. The generalized data were fitted by surface equations and used together with the information of the 'record & verify' network of the centers. In this manner, for every beam, the D(iso) reconstruction was obtained in about 25 seconds after the treatment. To test the in vivo dosimetric procedure, six pelvic treatments that used conformed wedged beams were carried out with three linacs of different manufacturers. For every beam, the comparison between the reconstructed D(iso) and the D(iso,TPS) computed by the TPS, resulted in an acceptable tolerance level of ±5%, estimated for this kind of treatment. Generally the in vivo dosimetry methods that use EPIDs require: (i) a special effort for the dosimetric commissioning with SPs of different thicknesses, and (ii) extra time for the analysis of the EPID signals. The proposed procedure simplifies the commissioning step and supplies for Varian, Elekta, and Siemens linacs equipped with the aSi EPIDs a quasi-real time in vivo dosimetry for open and wedged 3DCRT fields.

Spatially fractionated radiotherapy (SFRT) targeting the hypoxic tumor segment for the intentional induction of non-targeted effects: An in silico study to exploit a new treatment paradigm.

INTRODUCTION: The possibility of intentionally triggering non targeted effects (NTEs) using spatially fractionated radiotherapy (SFRT) alone or combined with immunotherapy is an intriguing and fascinating area of research. Among different techniques for SFRT, stereotactic body radiotherapy targeting exclusively the central hypoxic segment of bulky tumors, (SBRT-PATHY) might trigger immunogenic cell death more efficiently. This in silico study aims to identify the best possible dosimetric trade-off for prescribing SFRT with volumetric modulated arc (VMAT) based stereotactic radiotherapy (SRT). MATERIAL AND METHODS: Eight spherical volumes defined "Gross Tumor Volumes" (GTVs) were generated with diameters of 3-10 cm (with incremental steps of 1 cm), simulating tumor lesions. The inner third part of each GTV (GTVcentral) was selected to simulate the central hypoxic area and a ring structure was derived around it to simulate the tumor periphery (GTVperipheral). Volumetric modulated arc radiation treatment (VMAT) plans were calculated to deliver a single fraction of 10 Gy to each GTVcentral with different dose prescription methods: target mean and isodose driven (40, 50, 60, 70, 80 and 90%).The volume of GTVperipheral receiving less than 2 Gy was recorded as dosimetric performance indicator. RESULTS: 56 possible dosimetric scenarios were analyzed. The largest percentage of GTVperipheral spared from the dose of 2 Gy was achieved with dose prescription methods to the 70% isodose line for lesions smaller than 6 cm (range 42.9-48.4%) and to the target mean for larger ones (range 52.9-64.5%). CONCLUSIONS: Optimizing the dose prescription method may reduce the dose to tumor periphery in VMAT-based SFRT, thus potentially sparing tumor infiltrating immune cells. The optimal method may vary according to the size of the lesion. This should be taken into account when designing prospective trials using SFRT.

A validation study of a dedicated software for an automated in vivo dosimetry control in radiotherapy.

In vivo dosimetry (IVD) is the last step of a radiotherapy quality control program aimed to ensure that the dose delivered is in agreement with that prescribed. IVD procedures based on single detectors are time-consuming and impossible to use for the modern radiotherapy techniques, based on static or kinetic beams (modulated in intensity fluence); this means that more efficient and practical methods are highly recommended. The practical method SOFTDISO, based on the use of electronic portal image device (EPID), provides two tests (i) the R ratio between the reconstructed and the planned isocenter doses to verify an agreement within 5% and (ii) the γ-analysis of the EPID images, to verify γ% ≥ 90% and γmean ≤ 0.4. This paper reports the results of 11,357 IVD tests carried out for 823 patients treated by three-dimensional conformal radiation therapy and volumetric modulated arc therapy techniques. In particular, the dose disagreements are reported distinguishing two kinds of causes, those of (i) class 1 that includes the errors due to inadequate quality controls and (ii) the class 2, due to patient morphological changes. About the tests out of tolerance, 6% were by VMAT and 21% by 3DCRT, but taking into account the only class 1 of errors, i.e., removing the causes of class 2, only 7% of patients examined presented at least one of the three mean indexes out of tolerance. The workload for IVD on 9 patients/day per linac is about 52 min/day but recently, a new automated SOFTDISO version has been implemented to reduce the time to about 34 min/day.

Portal dose measurements by a 2D array.

A 2D array (PTW, type 10024), equipped with 729 vented plane parallel ion-chambers, has been calibrated as a detector for the in vivo comparison between measured and predicted portal doses for head-neck tumors. The comparison of absolute portal doses measured to ones predicted by a commercial treatment planning system within the field of view of the CT scanner, can help the delivered dose verification during different treatment fractions, in particular when the patient's present weight loss. This paper reports the preliminary results of the comparison of the portal doses measured by a PTW 2D array during several radiotherapy fractions and the predicted portal doses for seven patients undergoing head-neck tumor radiotherapy. The gamma index analysis supplied an agreement of more than 95% of the dose-point P(gamma)>95% within acceptance criteria, in terms of dose difference, DeltaD(max), and distance-agreement, Deltad(max), equal to 5% and 4mm, respectively. After the third week, one patient showed a decrease of P(gamma) values due to the markedly reduced patient's thickness. Even if the spatial resolution of the 2D array was 1cm, there were two advantages in the use of this 2D array as a portal dose device for IMRT quality control. The first one was the use of a stable and efficient absolute dosimeter for in vivo verification, although its construction and behavior for other gantry angles need to be tested, and the second one was the time efficiency in verifying the correct dose delivery in several fractions of the therapy. This study presents acceptance criteria for the comparison of TPS-predicted portal dose images with in vivo 2D ion-chamber measurements for IMRT. In particular, portal dose measurements offer clues for additional studies as to which indicators can signal the need for replanning during treatment.

Step-and-Shoot IMRT by Siemens Beams: An EPID Dosimetry Verification During Treatment.

PURPOSE: This work reports the extension of a semiempirical method based on the correlation ratios to convert electronic portal imaging devices transit signals into in vivo doses for the step-and-shoot intensity-modulated radiotherapy Siemens beams. The dose reconstructed at the isocenter point Diso, compared to the planned dose, Diso,TPS, and a γ-analysis between 2-dimensional electronic portal imaging device images obtained day to day, seems to supply a practical method to verify the beam delivery reproducibility. METHOD: The electronic portal imaging device images were obtained by the superposition of many segment fields, and the algorithm for the Diso reconstruction for intensity-modulated radiotherapy step and shoot was formulated using a set of simulated intensity-modulated radiotherapy beams. Moreover, the in vivo dose-dedicated software was integrated with the record and verify system of the centers. RESULTS: Three radiotherapy centers applied the in vivo dose procedure at 30 clinical intensity-modulated radiotherapy treatments, each one obtained with 5 or 7 beams, and planned for patients undergoing radiotherapy for prostatic tumors. Each treatment beam was checked 5 times, obtaining 900 tests of the ratios R = Diso/Diso,TPS. The average R value was equal to 1.002 ± 0.056 (2 standard deviation), while the mean R value for each patient was well within 5%, once the causes of errors were removed. The γ-analysis of the electronic portal imaging device images, with 3% 3 mm acceptance criteria, showed 90% of the tests with Pγ < 1 ≥ 95% and γmean ≤ 0.5. The off-tolerance tests were found due to incorrect setup or presence of morphological changes. This preliminary experience shows the great utility of obtaining the in vivo dose results in quasi real time and close to the linac, where the radiotherapy staff may immediately spot possible causes of errors. The in vivo dose procedure presented here is one of the objectives of a project, for the development of practical in vivo dose procedures, financially supported by the Istituto Nazionale di Fisica Nucleare.

A preliminary dosimetric characterization of chemical vapor deposition diamond detector prototypes in photon and electron radiotherapy beams.

Three radiation detectors based on polycrystalline diamond films with different thickness and resistivity, obtained by microwave chemical vapor deposition, were tested to assess their suitability for relative dosimetry of photon and electron beams supplied by clinical linear accelerators. All samples showed a linear response as a function of the absorbed dose. The sensitivity per unit of detector sensitive volume spanned between 7 and 43 nC Gy(-1) mm(-3) with an applied electric field of 40 kV/cm. The dose rate dependence was evaluated following the Fowler theory and delta coefficient values between 0.95 and 1.00 were found for the three samples when polarized at 40 kV/cm. Percentage depth dose curves, output factors, and normalized dose profiles were determined for 6 and 10 MV photon beams and for 6 and 15 MeV electron beams. The results obtained with the diamond detectors were in good agreement with those obtained by reference detector measurements [all the data were within the experimental uncertainty of 1% (1sigma)].

PTW-diamond detector: dose rate and particle type dependence.

In this paper the suitability of a PTW natural diamond detector (DD) for relative and reference dosimetry of photon and electron beams, with dose per pulse between 0.068 mGy and 0.472 mGy, was studied and the results were compared with those obtained by a stereotactic silicon detector (SFD). The results show that, in the range of the examined dose per pulse the DD sensitivity changes up to 1.8% while the SFD sensitivity changes up to 4.5%. The fitting parameter, delta, used to correct the dose per pulse dependence of solid state detectors, was delta = 0.993 +/- 0.002 and delta = 1.025 +/- 0.002 for the diamond detector and for the silicon diode, respectively. The delta values were found to be independent of particle type of two conventional beams (a 10 MV x-ray beam and a 21 MeV electron beam). So if delta is determined for a radiotherapy beam, it can be used to correct relative dosimetry for other conventional radiotherapy beams. Moreover the diamond detector shows a calibration factor which is independent of beam quality and particle type, so an empirical dosimetric formalism is proposed here to obtain the reference dosimetry. This formalism is based on a dose-to-water calibration factor and on an empirical coefficient, that takes into account the reading dependence on the dose per pulse.

p-type silicon detector for brachytherapy dosimetry.

The sensitivity of a cylindrical p-type silicon detector was studied by means of air and water measurements using different photon beams. A lead filter cap around the diode was used to minimize the dependence of the detector response as a function of the brachytherapy photon energy. The radial dose distribution of a high-activity 192Ir source in a brachytherapy phantom was measured by means of the shielded diode and the agreement of these data with theoretical evaluations confirms the method used to compensate diode response in the intermediate energy range. The diode sensitivity was constant over a wide range of dose rates of clinical interest; this allowed one to have a small detector calibrated in terms of absorbed dose in a medium. Theoretical evaluations showed that a single shielding filter around the p-type diode is sufficient to obtain accurate dosimetry for 192Ir, 137Cs, and 60Co brachytherapy sources.

A correction method for diamond detector signal dependence with proton energy.

Small dosimeters as solid state detectors can be useful for the dosimetric characterization and periodic quality control of radiotherapy proton beams. The calibration of solid state detectors for proton beams is not a solved problem especially for ophthalmologic proton beams, where these detectors present a LET-dependent signal. In this work a PTW diamond detector has been selected because of its good signal reproducibility (0.3%) and stable response with accumulated dose. A method that takes into account the LET dependence of the diamond detector signal, at 62 MeV proton beam, is here proposed. In particular an empirical correction factor, kDD(Eo) (Rres), has been determined as a function of the residual range quality index, to correct the diamond detector signal for a proton beam of incident effective energy E0= 62 MeV. A dedicated software allows us to use the diamond detector as an on-line reference dosimeter, where an ionization chamber may be difficult to use, or for periodic quality control procedures. The article also reports a comparison between the signal dependence on proton energy of silicon, diamond, and radiochromic film detectors.

A study of quality of bremsstrahlung spectra reconstructed from transmission measurements.

A numerical method for the reconstruction of bremsstrahlung spectra has been applied to the analysis of simulated data. The method is found to be reliable in reconstructing x-ray spectra of maximum energy up to 10 MeV. Measurements of aluminum transmission data carried out for five linear accelerators have also been analyzed by this method to determine the fractional energy fluence. The values of the Spencer-Attix water/air stopping power ratio, SS.A. w,air, reported as a function of the ratios TPR20 10, are in good agreement with computed results. Quality index values were correlated to mean incident energies of x-ray beams used in radiotherapy. The experimental setup has also provided information regarding the softening effect on the off-axis beam, which can be used in clinical computer dosimetry to correct conventional zero field size tissue maximum ratio.

Radiochromic film dosimetry of a low energy proton beam.

In this work some dosimetric characteristics of MD-55-2 GafChromic films were studied in a low energy proton beam (21.5 MeV) directly in a water phantom. The nonlinearity of the optical density was quantified by a factor P(lin). A correction factor P(en), that accounts for optical density dependence on the energy, was empirically determined. The effects of detector thickness in depth dose measurements and of the film orientation with respect to beam direction were investigated. The results show that the MD-55-2 films provide dose measurements with the films positioned perpendicularly to the proton beam. A dosimetric formalizm is proposed to determine the dose to water at depth d, with films oriented perpendicularly to the beam axis. This formalism uses a calibration factor of the radiochromic film determined directly on the proton beam at a reference depth in water, and the P(lin) factor, that takes into account the nonlinearity of the calibration curve and the P(en) factor that, in turn takes into account the change of proton beam energy in water. The MD-55-2 films with their high spatial resolution and the quasiwater equivalent material are attractive, positioned perpendicularly along the beam axis, for the absolute dose determination of very small beam sizes and modulated proton beams.

Applications of the Italian protocol for the calibration of brachytherapy sources.

The Associazione Italianà di Fisica Biomedica (AIFB; Italian Association of Biomedical Physics) has adopted the Italian protocol for the calibration of brachytherapy sources. The AIFB protocol allows measurements of the reference air kerma rate, Kr, within 1.7% (1 sigma). To measure Kr the AIFB protocol has identified a direct and an indirect procedure. The direct procedure is based on the use of spherical or cylindrical ionization chambers as local reference dosimeters positioned along the transverse bisector axis of the source. Once the source is specified by a Kr value, this can be used to calibrate a field instrument, such as a well-type ionization chamber, for further source calibrations by means of an indirect procedure. This paper reports the results obtained by the Physics Laboratory of the Università Cattolica del S Cuore (PL-UCSC), in terms of Kr calibration of five types of source, 169Yb, 192Ir and 137Cs. The role of the Kr determination for a brachytherapy source has been underlined when a new source such as the 169Yb seed model X1267 has been proposed for clinical use. The Kr values for 137Cs spherical sources differed by 5% from the vendor's mean value. The five types of source calibrated in terms of Kr were used to obtain the calibration factor. NKrsource, of an HDR-1000 well-type ionization chamber.

Implantation guidelines for 169 Yb seed interstitial treatments.

The adequacy of an interstitial implant carried out with a new radioactive source, the 169 Yb seed model X1267, has been examined by computing volumetric indices based on dose-volume histograms. The comparison of these indices with the ones computed for 125I seed implantations shows that the use of ytterbium seeds presents an improvement of the dose homogeneity in interstitial implants. This is due to the significant build-up associated with 169 Yb photons that reduces the rapid dose fall-off with the distance from the source. Moreover, relative to 192Ir, the lower photon energy gives 169 Yb the advantage in clinical use of reduced radiation exposure (i) to health care workers, (ii) to relatives of treated patients and (iii) to healthy neighbouring tissues of the patients if appropriate thin shielding is used.

Ion-recombination correction factor ksat for spherical ion chambers irradiated by continuous photon beams.

The large range of reference air kerma rates of brachytherapy sources involves the use of large-volume ionization chambers. When such ionization chambers are used the ion-recombination correction factor ksat has to be determined. In this paper three spherical ion chambers with volumes ranging from 30 to 10(4) cm3 have been irradiated by photons of a 192Ir source to determine the ksat factors. The ionization currents of the ion chambers as a function of the applied voltage and the air kerma rate have been analysed to determine the contribution of the initial and general ion recombination. The ksat values for large-volume ionization chambers obtained by considering the general ion recombination as predominant (Almond's approach) are in disagreement with the results obtained using methods that consider both initial and general ion-recombination contributions (Niatel's approach). Such disagreement can reach 0.7% when high currents are measured for a high-activity source calibration in terms of reference air kerma rate. In this study a new 'two-voltage' method, independent of the voltage ratio given by a dosimetry system, is proposed for practical dosimetry of continuous x- and gamma-radiation beams. In the case where the Almond approach is utilized, the voltage ratio V1/V2 should be less than 2 instead of Almond's limit of V1/V2 < 5.

Alanine/EPR dosimetry in brachytherapy.

The paper reports the experimental procedure adopted to determine the absorbed dose rate in water per reference air kerma rate, D(Kr), (d, theta), along the transverse bisector axis of a 137Cs brachytherapy source. The dose rate measurements have been carried out at difference distances, d, from the source using alanine dosemeters in a water phantom. The reference air kerma rate, Kr, was determined adopting a 'direct procedure' that uses a spherical ionization chamber in air. The dose rate constant of the source examined was D(Kr) (1, pi/2) = 0.99 +/- 0.03 cGy h(-1) (microGy h(-1))(-1). The values of the radial dose function along the transverse axis, g(d), determined with an uncertainty of 3.4% (1sigma), were found to be in good agreement with the results reported in the literature. The uncertainty in dose rate value has been estimated as 2.8% (1sigma) for distances from the source up to 7 cm. Kr has been determined with 1.2% (1sigma) uncertainty. So D(r) (d, pi/2) values were determined with 3% (1sigma) uncertainty.

Experimental dosimetry of a 32P catheter-based endovascular brachytherapy source.

The experimental dosimetry in a water phantom of a 32P linear source, 20 mm in length, used for the brachytherapy of coronary vessels is reported. The source content activity, A, was determined by means of a calibrated well ion-chamber and the value was compared with the contained activity reported in the manufacturer's certification. In this field of brachytherapy dosimetry, radiochromic film supplies a high enough spatial resolution. A highly sensitive radiochromic film, that presents only one active layer, was used in this work for the source dosimetry in a water phantom. The radiochromic film was characterized by electron beams produced by a clinical linac. A Monte Carlo calculation of beta spectra in water at different distances along the source transverse bisector axis allowed to take into account the low dependence of film response from the electron beam energy. The adopted experimental set-up, with the source in its catheter positioned on the film plane inside the water phantom, supplies accurate dosimetric information. The measured dose rate to water per unit of source activity at reference distance, D(r0, theta0)/A, in units of cGy s(-1) GBq(-1), was in agreement with the value reported in the manufacturer's certification within the experimental uncertainty. The radial dose function, g(r), is in good agreement with the literature data. The anisotropy function F(r, theta) is also reported. The analysis of the dose profile obtained at 2 mm from the source longitudinal axis shows that the uniformity is within 10% along 75% of the 20 mm treatment length. The adopted experimental set-up seems to be adequate for the quality control procedure of the dose homogeneity distribution in the water medium.

A standard dosimetry procedure for 192Ir sources used for endovascular brachytherapy.

The experimental dosimetry of a high dose rate (HDR) 192Ir source used for the brachytherapy of peripheral vessels is reported. The direct determination of the reference air kerma rate Kr agrees, within the experimental uncertainty, with the results obtained by a well ionization chamber calibrated at the NIST and the manufacturer's certification. A highly sensitive (HS) radiochromic film (RCF), that presents only one active layer, was used for the source dosimetry in a water phantom. The adopted experimental set-up, with the source in its catheter positioned on the RCF plane, seems to have given better accuracy of the RCF optical density measurements. The agreement between the measurement of the dose rate constant DKr (10 mm, pi/2) and the literature data confirmed the coherence of the HS RCF calibration obtained by the kerma in air measurements. The RCF measurements supplied dosimetric information about the dose to water per reference air kerma rate D(r, theta)/Kr along the source transverse bisector axis, the radial dose function g(r) and the anisotropy function F(r, theta). The value D(2 mm, pi/2)/Kr = 22.4 +/- 1.2 cGy h(-1)/(microGy h(-1)) is supplied with a dose uncertainty that is essentially due to the indeterminacy of the source position in the catheter. The data of the radial and anisotropy functions have been compared with Monte Carlo determinations reported in the literature.

Dosimetry of 169 Yb seed model X1267.

Unlike previous brachytherapy sources a number of published studies have been addressed to the dosimetry of 169 Yb seeds, manufactured in several prototypes, before widespread clinical use has been made. Discrepancies seen in the dosimetry obtained for ytterbium seed prototypes appear to be related to inconsistency and non-reproducibility in the vendor's calibration procedure to determine contained activity. Av. The comparison of 169 Yb seed dosimetries demonstrates a need for more accurate implementation of calibration procedures to determine the air kerma rate for the definitive 169 Yb seed design. This paper reports an experimental procedure to determine the reference air kerma rate, Kr (mu Gy h-1), defined as the kerma rate at 1 m along the source transverse axis in free space for the new 169 Yb seed, model X1267. A mean value of the ratio Kr/Av = 1.53 mu Gy h-1 mCi-1 was obtained from determining the Kr value of eleven seeds. Since this ratio is only 3% less than the air kerma rate constant for the 169 Yb point source, (gamma delta)k = 1.58 mu Gy h-1 m2 mCi-1, this means that the Av is closer to an apparent activity than a contained activity, Ac. A Monte Carlo simulation to determine the ratio between reference air kerma rate and the contained activity gave Kr/Ac = 1.33 mu Gy h-1 mCi-1. For the dose rate constant in water we obtained DKr (1, pi/2) = 1.20 +/- 0.05 cGy h-1 (mu Gy h-1)-1, using calibrated thermoluminescent dosimeters (TLDs) and DKr (1, pi/2) = 1.21 +/- 0.03 cGy h-1 (mu Gy h-1)-1 by Monte Carlo simulation. TLDs were used both to determine the radial dose distribution along the seed transverse axis and to calibrate GAFChromic films to obtain the two-dimensional dose distribution around the seed.

The saturation loss for plane parallel ionization chambers at high dose per pulse values.

The use of plane parallel ionization chambers with electron beams with high dose per pulse entails dose uncertainties due to the overestimation of the ion recombination factor, k, up to 20% if conventional dosimetric protocols are used. In this work MD-55-2 radiochromic films have been used as reference dosimeters to obtain dose to water per pulse DGAF(w) values for three Novac7 (Hitesys) electron beams of E0 = 5.8 MeV. However, the beam calibration by MD-55-2 films is time consuming and the use of plane parallel chambers is fundamental for a periodic quality control procedure. Three plane parallel chambers have been used and the general formula for the k determination has been tested using the calibration doses, DGAF(w). In particular, consistent ion recombination factors ksat(V0) (with the ion chamber polarized at V0), that follow the Boag theory, have been estimated at different dose per pulse values for the three plane parallel ionization chambers. This means that at present any ion chamber needs a specific ksat (V0) determination by using a reference dosimeter for which the response is independent of the dose rate. An accurate determination of ksat(V0), using a reference quality beam, can be used to determine the dose to water per pulse for electron beams of different quality and geometrical configuration.

Characteristics of silicon and diamond detectors in a 60 MeV proton beam.

The calibration factor variation for a PTW natural diamond detector and a Scanditronix p-type stereotactic silicon diode (designed for use in photon beams) was studied in the 10-59 MeV range. Irradiations were performed in a water phantom with the 60 MeV ocular therapy beam at the CCO (UK). The diamond detector showed a sensitivity increase with energy, underestimating the dose by about 18% at the Bragg peak, by 7% at the centre and by 17% at the distal end of the SOBP region. The silicon diode did not show any significant sensitivity change with energy. However, a decrease in response of 24% was observed for an accumulated dose of 300 Gy.