Ex) Article Title, Author, Keywords
Ex) Article Title, Author, Keywords
Progress in Medical Physics 2020; 31(4): 135-144
Published online December 31, 2020
https://doi.org/10.14316/pmp.2020.31.4.135
Copyright © Korean Society of Medical Physics.
Moonhee Lee1 , Sunghwan Ahn2 , Wonjoong Cheon1* , Youngyih Han1,2
Correspondence to:Youngyih Han
(youngyih@skku.edu)
Tel: 82-2-3410-2604
Fax: 82-2-3410-2619
*Present address: Department of Radiation Oncology, National Cancer Center (NCC), Goyang 10408, Korea
Purpose: Gafchromic films for proton dosimetry are dependent on linear energy transfers (LETs), resulting in dose underestimation for high LETs. Despite efforts to resolve this problem for singleenergy beams, there remains a need to do so for multi-energy beams. Here, a bimolecular reaction model was applied to correct the under-response of spread-out Bragg peaks (SOBPs).
Methods: For depth-dose measurements, a Gafchromic EBT3 film was positioned in water perpendicular to the ground. The gantry was rotated at 15° to avoid disturbances in the beam path. A set of films was exposed to a uniformly scanned 112-MeV pristine proton beam with six different dose intensities, ranging from 0.373 to 4.865 Gy, at a 2-cm depth. Another set of films was irradiated with SOBPs with maximum energies of 110, 150, and 190 MeV having modulation widths of 5.39, 4.27, and 5.34 cm, respectively. The correction function was obtained using 150.8-MeV SOBP data. The LET of the SOBP was then analytically calculated. Finally, the model was validated for a uniform cubic dose distribution and compared with multilayered ionization chamber data.
Results: The dose error in the plateau region was within 4% when normalized with the maximum dose. The discrepancy of the range was <1 mm for all measured energies. The highest errors occurred at 70 MeV owing to the steep gradient with the narrowest Bragg peak.
Conclusions: With bimolecular model-based correction, an EBT3 film can be used to accurately verify the depth dose of scanned proton beams and could potentially be used to evaluate the depth-dose distribution for patient plans.
KeywordsProton beam, Spread-out Bragg peak, Depth dose, Linear energy transfer, EBT3 film
In proton therapy using a dynamically scanned beam, the desired depth-dose distribution is composed of stacked beam layers with unique energies, and each layer is dynamically scanned with a pencil beam that is magnetically deflected [1]. Despite having advantages over the use of passive scattered beams [2], the verification of depth-dose distribution using an ionization chamber, which is a point-by-point detector, is extremely inefficient and laborious. Devices with multiple ionization chambers, such as multilayer ionization chambers (MLICs), enable a shorter measurement time owing to their short setup time and rapid estimation of depth dose. Therefore, MLICs have been used in many facilities to obtain depth-dose data. However, MLICs are expensive devices and have a detector spacing of 2 mm, which may be insufficient to provide the required spatial resolution for a low-energy proton beam, owing to the high-gradient distal depth dose. Moreover, MLICs only allow for measurements of fields smaller than 10×10 cm, rendering them unsuitable for measuring the depth dose of large field sizes. Therefore, a backup measurement system for depth-dose profiles that complements the limitations of MLICs is required.
Gafchromic EBT3 film (Ashland, Covington, GA, USA) possesses a high spatial resolution down to 25 µm, which is particularly beneficial for the high-dose-gradient region of proton beams. By placing the film along the direction of the dynamically delivered beam, the depth-dose distribution can be determined with a single irradiation, rendering it an extremely time-efficient tool compared with point-by-point dosimeters, which require a number of irradiations. However, the energy dependence of Gafchromic EBT3 film is still a critical issue in proton dosimetry, as proton beams undergo a significant change in energy and linear energy transfer (LET) along their paths, particularly in the Bragg peak region, where protons lose most of their energy and then stop, immediately after intensively depositing energy [3]. The energy dependence of the film leads to an underdose response in the Bragg peak region [4-6]. Therefore, determining an optimal method for correcting the underdose effect is essential to allow for the use of Gafchromic EBT films in the depth-dose measurement of proton beams.
The protocol for the calibration and dosimetry of EBT films provided by the manufacturer states that the calibration or characteristic curves are to be achieved by relating the film responses and the corresponding dose measurement results from a well-calibrated exposure, meaning it is necessary to apply a fitting function to the plot. However, a one-dimensional (1D) calibration function used to determine the absorbed dose based on film response cannot account for the implicit physics involved regarding the energy dependence of the film in the proton beam. Thus, conventional calibration methods for proton beams have resulted in underestimations of doses in the distal area.
Several researchers have attempted to resolve this underdosing issue [6-12]. In particular, numerical or mathematical correction methods that do not involve physical representation of the LET effect have been widely suggested. These methods complicate the correction of complex beam compositions, and therefore, previous studies have mainly focused on single-energy beams or a particular range of energies. Zhao and Das [7] reported polynomial fitting function, which is only based on a pristine Bragge peak, as a correction method for EBT3 films and evaluated the function by applying it to eight different energies ranging from 76 to 186 MeV. Park et al. [9] attempted to correct spread-out Bragg peaks (SOBPs) by separately correcting the data of each of the single-energy beams composing the SOBP beam. Subsequently, they combined these with a weighting function based on pristine Bragg peaks, thereby requiring a measurement for each single-energy beam with several sheets of films and the respective corrections [9]. Gambarini et al. [10] presented a method that uses the ratio of the percentage depth dose (PDD), measured using an ionization chamber, to the EBT3 film as a correction factor. They used the method on a multi-energy beam ranging from 62.73 to 81.56 MeV, which comprises the low-energy range in proton therapy [10]. Fiorini et al. [11,12] suggested a correction method for SOBP beams based on Monte Carlo simulations using a beam quality correction factor. The beam quality factor was successfully used to rectify underestimated doses for SOBPs. However, it was benchmarked using a 29-MeV modulated beam that is only applicable to low-energy applications, such as eye treatments. Perles et al. [8] suggested a bimolecular model as a method for film response-to-dose calibrations using an LET correction-incorporated system, as the radiation-induced chemical reaction of radiochromic films occurred between two adjacent monomers within the same crystal lattice. This theoretically derived formalism is based on the law of mass action, which accounts for the reaction rate of materials. The bimolecular model is valid within a 5% accuracy for the depth dose of an 85.6-MeV beam [8].
For the bimolecular model to be clinically applicable, further studies are required to extend the depth-dose measurements to multi-energy beams for a diverse range of energies. In this study, a method was thus developed for multi-energy beams with a larger energy range (i.e., 70–190 MeV), with a focus on clinical applications. To achieve this, an SOBP was formed by a number of different energy layers. Overall, the aim of this study was to devise a correction method for the sum of the beam layers.
The bimolecular model has been proposed as a correction method for the energy dependence of EBT2 films [8] and is generally presented as a formal derivation of the Schumann plate equation. The Schumann plate equation was originally derived for the detection of ultraviolet radiation using photographic plates [13] and is expressed as shown in Eq. (1). Subsequently, Candler [14,15] used the law of mass action to derive the Schumann equation by modeling a diatomic reaction. Regarding the chemical reaction to a radiation dose, at least two adjacent monomers must provoke a chemical response, or polymerization, which is modeled as a bimolecular chemical reaction. Perles et al. [8] modified the bimolecular model to account for the LET effect by making the reaction cross-section constant a function of the LET. In this model, the Schumann plate equation is used to represent the relationship between the dose and film response, with consideration of LET dependence.
It is noteworthy that the parameters in the bimolecular model are represented differently from those in the original Schumann plate equation. In Eq. (1),
Film measurements were performed using dynamic line-scanned beams comprising single and multiple energies to form an SOBP using the Sumitomo Proton Therapy System (Sumitomo Heavy Industries, Tokyo, Japan) at the Proton Therapy Center of the Samsung Medical Center (SMC, Seoul, Korea). To measure the pristine Bragg peaks, beams of energy ranging from 70 to 190 MeV, as well as a 112-MeV beam with six different dose intensities ranging from 0.373 to 4.865 Gy, were used. Furthermore, SOBPs were produced with three different maximum energies with modulation widths of 5.39, 4.27, and 5.34 cm. The field size of all irradiated beams was 10×10 cm. For continuous line-scanned beams, the dose intensity was determined by the scan speed (mph) and dose rate (MU/s). Once the range and modulation width were determined, an SOBP was created by consecutively scanning each layer with different energies and dose weights of the SOBP. For beam irradiations, the scan speed was changed according to the dose weight, whereas the dose rate was fixed for all energies. The irradiated beams are listed in Table 1.
Table 1 Specifications of irradiated proton beams
Specification | Pristine Bragg peak | SOBP | ||||||
---|---|---|---|---|---|---|---|---|
Snout Position (mm) | 400 | 400 | ||||||
Collimator (cm×cm) | 10×10 | 10×10 | ||||||
Energy/max. Energy of SOBP (MeV) | 70 | 92 | 112 | 150 | 190 | 110 | 150 | 190 |
Modulation width of SOBP (cm) | - | 5.39 | 4.27 | 5.34 |
Gafchromic EBT3 film (Ashland ISP Advanced Materials, Wayne, NJ, USA) was used to obtain the depth-dose distribution. All measurements were performed using film sheets from the same batch (Lot no. 05181501). To obtain the depth-dose profile, the film was positioned perpendicular to the ground and fixed by an in-house-designed film holder, as shown in Fig. 1. The in-house-designed phantom was composed of an acrylic framework and bolts. The films could be inserted easily into the holder, even in water, and the bolts were tightened to ensure that the EBT3 film was fixed firmly. However, if the bolts were fixed too tightly, it was difficult to insert the film.
Because the film was set parallel to the beam axis in water, perturbations occurred in the proton beam range. To achieve an optimum range accuracy, it is essential to avoid perturbations. First, the film was completely immersed in water to prevent perturbations caused by the surface tension of the water with the film [16]. Fig. 1c shows that the distance between the upper film edge and water surface was 1 cm. Next, the gantry was rotated 15° with respect to the film to prevent the disturbance of the proton beam path by the film. When the film was in the plane composed of the beam axis, perturbation occurred owing to the EBT3 film's lack of water equivalence for proton beams and the interruption of the beam path by different materials. These factors caused inaccuracies in the range verification owing to dosimetric perturbation. To avoid the aforementioned issue, the gantry can be rotated by a small amount, such as 3° or 5° [7,8,16]. However, protons penetrate through tissues in a relatively straight manner [17], and therefore, rotating the gantry by the maximum angle would yield more accurate measurements, as most of the beams would travel in water rather than in the EBT3 film. An angle of 15° was determined to be sufficiently large while ensuring that measurements were not affected at the edge by the penumbra of the 10×10 cm field.
Two factors could compensate for the perturbations: the gap between the film edge and water surface and oblique beam irradiation. A geometrical correction for the depth distance should be applied to verify the range of each proton beam with the minimum error, as follows:
Here, θ is the correction angle, which is the angle between the beam axis and the film;
The films were scanned twice for background subtraction, once before and once after irradiation, at the same location and position on the scanner bed. The post-irradiation scan was performed 48 hours after irradiation using an Epson 11000XL flatbed scanner (Seiko Epson Co., Tokyo, Japan). Post-irradiation film analyses, including background subtraction and dose conversion with LET correction, were performed using MATLAB version 9.0, R2016a (MathWorks Inc., Natick, MA, USA).
For a read-out process, a characteristic curve must be determined. Because the film response depends on the energy and intensity of the proton beam, the formalism used to convert the film response to the absorbed dose is dependent on both variables. First, the 1D relation between film response and absorbed dose for a fixed LET value was determined, and the LET dependence was subsequently considered.
A PPC05 parallel plate-type ionization chamber (IBA dosimetry, Schwarzenbruck, Germany) was used to measure the absolute point dose in the water at a specified depth. Measurements were performed for 112-MeV beams with six different dose intensities at a depth of 2 cm, which is the plateau region of the pristine Bragg peak. The absolute dose could be determined from the direct relationship between film response and absolute dose at a certain depth, namely, at a fixed LET, and Fig. 2 shows the scanned image of the film used to obtain the characteristic curve. By fitting the six plots using Eq. (1),
Wilkens and Oelfke [18] have presented a model for the analytical calculation of LET and validated it via Monte Carlo simulations. An analytical calculation of LET was adopted here owing to the fast calculation speed. Although the model of Wilkens and Oelfke [18] does not consider the effect of secondary particles, secondary particles are primarily produced by non-elastic nuclear interactions, which are most significant in the plateau region and negligible in the Bragg peak region [19,20]. Wilkens and Oelfke [18] describe the LET as the mean value of the stopping power. In this study, particle fluence appeared at a particular point
The SOBP is composed of each of the pristine Bragg peaks along with their respective dose weights. Likewise, the LET of the SOBP was calculated as the weighted sum of each Bragg peak, as expressed in Eq. (3)-2, where
Once the constant values of
The correction process was performed for an automated and optimized beam planned by a TPS, unlike the manually planned SOBPs presented in Table 1. One of the treatment plans was exported from the TPS in RayStation version 6 (6.2.0.7; RaySearch Laboratories AB, Stockholm, Sweden) to evaluate the application of the LET correction method to the TPS data. A dose with a uniform cubic distribution was irradiated in water, which is a uniform SOBP plane, with a square beam measuring 10×10 cm and comprising 34 energy layers ranging from 81.6 to 155.2 MeV. As shown in Fig. 4, the dose distribution was calculated using the TPS and then compared with the corrected film data. Each layer had a different line spacing and dose distribution, determined automatically from the TPS. The relative dose weights of each energy layer were exported from the treatment plan report to calculate the LET.
The film responses for the six different dose intensities of the 112-MeV beam and their corresponding doses measured using an ionization chamber were plotted (blue inverted triangle) and fitted (red solid line) with the curve of Eq. (1) (Fig. 5). At a 2-cm depth, the LET was 1.54 keV/µm, and the corresponding
Table 2 Fifth-degree Gaussian curve fitting parameters, as presented in Eq. (4)
xi | |||||
---|---|---|---|---|---|
2.95×105 | 1.78×105 | 3.12×105 | –1.53×105 | 4.43×105 | |
7.17 | 6.26 | 5.31 | 3.4 | 1.69 | |
0.9 | 1.04 | 1.91 | 0.18 | 5.09 |
Once
Table 3 Range (R80) comparison of single-energy, multi-energy, and TPS-planned beams between the EBT3 film and reference measurements or calculations
Energy | Single-energy | 70 MeV | 92 MeV | 112 MeV | 150 MeV | 170 MeV | 190 MeV |
---|---|---|---|---|---|---|---|
R80 (cm) | EBT3 | 3.79 | 6.38 | 9.22 | 15.66 | 19.54 | 23.63 |
Reference | 3.77 | 6.39 | 9.21 | 15.61 | 19.45 | 23.54 | |
Energy | Multi-energy | SOBP 110 MeV | SOBP 150 MeV | SOBP 190 MeV | TPS cubic plan | ||
R80 (cm) | EBT3 | 8.78 | 15.62 | 23.54 | 15.62 | ||
Reference | 8.85 | 15.54 | 23.45 | 15.53 |
The use of the Gafchromic EBT3 film for proton dosimetry has led to significant issues. In particular, the energy or LET dependence of the film has resulted in the underestimation of measured doses at low energies or for high-LET regions. As the bimolecular model is based on chemical reaction laws, which are affected by physical phenomena (known as the LET effect of protons), the performance of previously developed correction methods is adequate when the LET along the depth of the irradiated beam is known.
In this study, the bimolecular model was adapted and successfully extended for application to clinically meaningful multi-energy beams or SOBPs. The LET of a multienergy beam was calculated using an analytical method that considered the dose weight and LET of each single energy beam constituting the SOBP. The model was validated for a clinically important range of energies, from 70 to 190 MeV. The depth-dose distribution of the EBT3 film had an error within 4%, and it demonstrated submillimeter accuracy in range verification when compared with MLIC measurements and TPS calculations. The submillimeter accuracy of this method also enables range measurements of unknown beam energies, meaning that, conversely, an unknown beam energy can be calculated by measuring the range. The observed error of range verification contained uncertainties caused by the non-water characteristics of the EBT3 film material, and the effect of secondary particles on LET calculation was disregarded. Hence, the EBT3 film can be used with acceptable uncertainty using the investigated correction method.
Measuring the depth dose of a dynamically scanned proton beam using an ionization chamber is challenging owing to the extremely time-consuming dosimetry procedure. By constructing a generalized correction method for energy dependence, the EBT3 film becomes an excellent candidate for the depth-dose measurements of scanned beams owing to its efficient measurement procedure, which requires only one beam irradiation, and the high spatial resolution that can be obtained. A high spatial resolution provides excellent usability for low-energy beams that have steep distal fall-offs of the SOBP, such as eye treatment beams. A depth-wise dosimetric quality assurance (QA) is feasible for clinical beams using the EBT3 film. Hence, this could be a promising tool for patient-related QA. However, it still needs to be verified with a non-uniform depth-dose distribution by computing the LET at each location or depth. Calculating the correction parameter for the depth-dose distribution requires the dose weight of each energy along the depth, which can be obtained by analyzing the log file from the TPS. Therefore, further studies regarding nonuniform beams derived from TPS will warrant application to clinical cases.
This research was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2019M2A2B4096537 and 2019R1F1A1062775).
The authors have nothing to disclose.
All relevant data are within the paper and its Supporting Information files.
Conceptualization: Youngyih Han. Data curation: Moonhee Lee. Funding acquisition: Youngyih Han. Investigation: Moonhee Lee, Wonjoong Cheon, Sunghwan Ahn. Methodology: Moonhee Lee, Wonjoong Cheon. Project administration: Moonhee Lee, Youngyih Han. Resources: Youngyih Han. Supervision: Youngyih Han. Validation: Moonhee Lee. Visualization: Moonhee Lee. Writing–Original draft: Moonhee Lee. Writing–review & editing: Moonhee Lee, Younghyih Han.
Progress in Medical Physics 2020; 31(4): 135-144
Published online December 31, 2020 https://doi.org/10.14316/pmp.2020.31.4.135
Copyright © Korean Society of Medical Physics.
Moonhee Lee1 , Sunghwan Ahn2 , Wonjoong Cheon1* , Youngyih Han1,2
1Samsung Advanced Institute for Health Sciences & Technology (SAIHST), Sungkyunkwan University, 2Department of Radiation Oncology, Samsung Medical Center (SMC), Seoul, Korea
Correspondence to:Youngyih Han
(youngyih@skku.edu)
Tel: 82-2-3410-2604
Fax: 82-2-3410-2619
*Present address: Department of Radiation Oncology, National Cancer Center (NCC), Goyang 10408, Korea
Purpose: Gafchromic films for proton dosimetry are dependent on linear energy transfers (LETs), resulting in dose underestimation for high LETs. Despite efforts to resolve this problem for singleenergy beams, there remains a need to do so for multi-energy beams. Here, a bimolecular reaction model was applied to correct the under-response of spread-out Bragg peaks (SOBPs).
Methods: For depth-dose measurements, a Gafchromic EBT3 film was positioned in water perpendicular to the ground. The gantry was rotated at 15° to avoid disturbances in the beam path. A set of films was exposed to a uniformly scanned 112-MeV pristine proton beam with six different dose intensities, ranging from 0.373 to 4.865 Gy, at a 2-cm depth. Another set of films was irradiated with SOBPs with maximum energies of 110, 150, and 190 MeV having modulation widths of 5.39, 4.27, and 5.34 cm, respectively. The correction function was obtained using 150.8-MeV SOBP data. The LET of the SOBP was then analytically calculated. Finally, the model was validated for a uniform cubic dose distribution and compared with multilayered ionization chamber data.
Results: The dose error in the plateau region was within 4% when normalized with the maximum dose. The discrepancy of the range was <1 mm for all measured energies. The highest errors occurred at 70 MeV owing to the steep gradient with the narrowest Bragg peak.
Conclusions: With bimolecular model-based correction, an EBT3 film can be used to accurately verify the depth dose of scanned proton beams and could potentially be used to evaluate the depth-dose distribution for patient plans.
Keywords: Proton beam, Spread-out Bragg peak, Depth dose, Linear energy transfer, EBT3 film
In proton therapy using a dynamically scanned beam, the desired depth-dose distribution is composed of stacked beam layers with unique energies, and each layer is dynamically scanned with a pencil beam that is magnetically deflected [1]. Despite having advantages over the use of passive scattered beams [2], the verification of depth-dose distribution using an ionization chamber, which is a point-by-point detector, is extremely inefficient and laborious. Devices with multiple ionization chambers, such as multilayer ionization chambers (MLICs), enable a shorter measurement time owing to their short setup time and rapid estimation of depth dose. Therefore, MLICs have been used in many facilities to obtain depth-dose data. However, MLICs are expensive devices and have a detector spacing of 2 mm, which may be insufficient to provide the required spatial resolution for a low-energy proton beam, owing to the high-gradient distal depth dose. Moreover, MLICs only allow for measurements of fields smaller than 10×10 cm, rendering them unsuitable for measuring the depth dose of large field sizes. Therefore, a backup measurement system for depth-dose profiles that complements the limitations of MLICs is required.
Gafchromic EBT3 film (Ashland, Covington, GA, USA) possesses a high spatial resolution down to 25 µm, which is particularly beneficial for the high-dose-gradient region of proton beams. By placing the film along the direction of the dynamically delivered beam, the depth-dose distribution can be determined with a single irradiation, rendering it an extremely time-efficient tool compared with point-by-point dosimeters, which require a number of irradiations. However, the energy dependence of Gafchromic EBT3 film is still a critical issue in proton dosimetry, as proton beams undergo a significant change in energy and linear energy transfer (LET) along their paths, particularly in the Bragg peak region, where protons lose most of their energy and then stop, immediately after intensively depositing energy [3]. The energy dependence of the film leads to an underdose response in the Bragg peak region [4-6]. Therefore, determining an optimal method for correcting the underdose effect is essential to allow for the use of Gafchromic EBT films in the depth-dose measurement of proton beams.
The protocol for the calibration and dosimetry of EBT films provided by the manufacturer states that the calibration or characteristic curves are to be achieved by relating the film responses and the corresponding dose measurement results from a well-calibrated exposure, meaning it is necessary to apply a fitting function to the plot. However, a one-dimensional (1D) calibration function used to determine the absorbed dose based on film response cannot account for the implicit physics involved regarding the energy dependence of the film in the proton beam. Thus, conventional calibration methods for proton beams have resulted in underestimations of doses in the distal area.
Several researchers have attempted to resolve this underdosing issue [6-12]. In particular, numerical or mathematical correction methods that do not involve physical representation of the LET effect have been widely suggested. These methods complicate the correction of complex beam compositions, and therefore, previous studies have mainly focused on single-energy beams or a particular range of energies. Zhao and Das [7] reported polynomial fitting function, which is only based on a pristine Bragge peak, as a correction method for EBT3 films and evaluated the function by applying it to eight different energies ranging from 76 to 186 MeV. Park et al. [9] attempted to correct spread-out Bragg peaks (SOBPs) by separately correcting the data of each of the single-energy beams composing the SOBP beam. Subsequently, they combined these with a weighting function based on pristine Bragg peaks, thereby requiring a measurement for each single-energy beam with several sheets of films and the respective corrections [9]. Gambarini et al. [10] presented a method that uses the ratio of the percentage depth dose (PDD), measured using an ionization chamber, to the EBT3 film as a correction factor. They used the method on a multi-energy beam ranging from 62.73 to 81.56 MeV, which comprises the low-energy range in proton therapy [10]. Fiorini et al. [11,12] suggested a correction method for SOBP beams based on Monte Carlo simulations using a beam quality correction factor. The beam quality factor was successfully used to rectify underestimated doses for SOBPs. However, it was benchmarked using a 29-MeV modulated beam that is only applicable to low-energy applications, such as eye treatments. Perles et al. [8] suggested a bimolecular model as a method for film response-to-dose calibrations using an LET correction-incorporated system, as the radiation-induced chemical reaction of radiochromic films occurred between two adjacent monomers within the same crystal lattice. This theoretically derived formalism is based on the law of mass action, which accounts for the reaction rate of materials. The bimolecular model is valid within a 5% accuracy for the depth dose of an 85.6-MeV beam [8].
For the bimolecular model to be clinically applicable, further studies are required to extend the depth-dose measurements to multi-energy beams for a diverse range of energies. In this study, a method was thus developed for multi-energy beams with a larger energy range (i.e., 70–190 MeV), with a focus on clinical applications. To achieve this, an SOBP was formed by a number of different energy layers. Overall, the aim of this study was to devise a correction method for the sum of the beam layers.
The bimolecular model has been proposed as a correction method for the energy dependence of EBT2 films [8] and is generally presented as a formal derivation of the Schumann plate equation. The Schumann plate equation was originally derived for the detection of ultraviolet radiation using photographic plates [13] and is expressed as shown in Eq. (1). Subsequently, Candler [14,15] used the law of mass action to derive the Schumann equation by modeling a diatomic reaction. Regarding the chemical reaction to a radiation dose, at least two adjacent monomers must provoke a chemical response, or polymerization, which is modeled as a bimolecular chemical reaction. Perles et al. [8] modified the bimolecular model to account for the LET effect by making the reaction cross-section constant a function of the LET. In this model, the Schumann plate equation is used to represent the relationship between the dose and film response, with consideration of LET dependence.
It is noteworthy that the parameters in the bimolecular model are represented differently from those in the original Schumann plate equation. In Eq. (1),
Film measurements were performed using dynamic line-scanned beams comprising single and multiple energies to form an SOBP using the Sumitomo Proton Therapy System (Sumitomo Heavy Industries, Tokyo, Japan) at the Proton Therapy Center of the Samsung Medical Center (SMC, Seoul, Korea). To measure the pristine Bragg peaks, beams of energy ranging from 70 to 190 MeV, as well as a 112-MeV beam with six different dose intensities ranging from 0.373 to 4.865 Gy, were used. Furthermore, SOBPs were produced with three different maximum energies with modulation widths of 5.39, 4.27, and 5.34 cm. The field size of all irradiated beams was 10×10 cm. For continuous line-scanned beams, the dose intensity was determined by the scan speed (mph) and dose rate (MU/s). Once the range and modulation width were determined, an SOBP was created by consecutively scanning each layer with different energies and dose weights of the SOBP. For beam irradiations, the scan speed was changed according to the dose weight, whereas the dose rate was fixed for all energies. The irradiated beams are listed in Table 1.
Table 1 . Specifications of irradiated proton beams.
Specification | Pristine Bragg peak | SOBP | ||||||
---|---|---|---|---|---|---|---|---|
Snout Position (mm) | 400 | 400 | ||||||
Collimator (cm×cm) | 10×10 | 10×10 | ||||||
Energy/max. Energy of SOBP (MeV) | 70 | 92 | 112 | 150 | 190 | 110 | 150 | 190 |
Modulation width of SOBP (cm) | - | 5.39 | 4.27 | 5.34 |
Gafchromic EBT3 film (Ashland ISP Advanced Materials, Wayne, NJ, USA) was used to obtain the depth-dose distribution. All measurements were performed using film sheets from the same batch (Lot no. 05181501). To obtain the depth-dose profile, the film was positioned perpendicular to the ground and fixed by an in-house-designed film holder, as shown in Fig. 1. The in-house-designed phantom was composed of an acrylic framework and bolts. The films could be inserted easily into the holder, even in water, and the bolts were tightened to ensure that the EBT3 film was fixed firmly. However, if the bolts were fixed too tightly, it was difficult to insert the film.
Because the film was set parallel to the beam axis in water, perturbations occurred in the proton beam range. To achieve an optimum range accuracy, it is essential to avoid perturbations. First, the film was completely immersed in water to prevent perturbations caused by the surface tension of the water with the film [16]. Fig. 1c shows that the distance between the upper film edge and water surface was 1 cm. Next, the gantry was rotated 15° with respect to the film to prevent the disturbance of the proton beam path by the film. When the film was in the plane composed of the beam axis, perturbation occurred owing to the EBT3 film's lack of water equivalence for proton beams and the interruption of the beam path by different materials. These factors caused inaccuracies in the range verification owing to dosimetric perturbation. To avoid the aforementioned issue, the gantry can be rotated by a small amount, such as 3° or 5° [7,8,16]. However, protons penetrate through tissues in a relatively straight manner [17], and therefore, rotating the gantry by the maximum angle would yield more accurate measurements, as most of the beams would travel in water rather than in the EBT3 film. An angle of 15° was determined to be sufficiently large while ensuring that measurements were not affected at the edge by the penumbra of the 10×10 cm field.
Two factors could compensate for the perturbations: the gap between the film edge and water surface and oblique beam irradiation. A geometrical correction for the depth distance should be applied to verify the range of each proton beam with the minimum error, as follows:
Here, θ is the correction angle, which is the angle between the beam axis and the film;
The films were scanned twice for background subtraction, once before and once after irradiation, at the same location and position on the scanner bed. The post-irradiation scan was performed 48 hours after irradiation using an Epson 11000XL flatbed scanner (Seiko Epson Co., Tokyo, Japan). Post-irradiation film analyses, including background subtraction and dose conversion with LET correction, were performed using MATLAB version 9.0, R2016a (MathWorks Inc., Natick, MA, USA).
For a read-out process, a characteristic curve must be determined. Because the film response depends on the energy and intensity of the proton beam, the formalism used to convert the film response to the absorbed dose is dependent on both variables. First, the 1D relation between film response and absorbed dose for a fixed LET value was determined, and the LET dependence was subsequently considered.
A PPC05 parallel plate-type ionization chamber (IBA dosimetry, Schwarzenbruck, Germany) was used to measure the absolute point dose in the water at a specified depth. Measurements were performed for 112-MeV beams with six different dose intensities at a depth of 2 cm, which is the plateau region of the pristine Bragg peak. The absolute dose could be determined from the direct relationship between film response and absolute dose at a certain depth, namely, at a fixed LET, and Fig. 2 shows the scanned image of the film used to obtain the characteristic curve. By fitting the six plots using Eq. (1),
Wilkens and Oelfke [18] have presented a model for the analytical calculation of LET and validated it via Monte Carlo simulations. An analytical calculation of LET was adopted here owing to the fast calculation speed. Although the model of Wilkens and Oelfke [18] does not consider the effect of secondary particles, secondary particles are primarily produced by non-elastic nuclear interactions, which are most significant in the plateau region and negligible in the Bragg peak region [19,20]. Wilkens and Oelfke [18] describe the LET as the mean value of the stopping power. In this study, particle fluence appeared at a particular point
The SOBP is composed of each of the pristine Bragg peaks along with their respective dose weights. Likewise, the LET of the SOBP was calculated as the weighted sum of each Bragg peak, as expressed in Eq. (3)-2, where
Once the constant values of
The correction process was performed for an automated and optimized beam planned by a TPS, unlike the manually planned SOBPs presented in Table 1. One of the treatment plans was exported from the TPS in RayStation version 6 (6.2.0.7; RaySearch Laboratories AB, Stockholm, Sweden) to evaluate the application of the LET correction method to the TPS data. A dose with a uniform cubic distribution was irradiated in water, which is a uniform SOBP plane, with a square beam measuring 10×10 cm and comprising 34 energy layers ranging from 81.6 to 155.2 MeV. As shown in Fig. 4, the dose distribution was calculated using the TPS and then compared with the corrected film data. Each layer had a different line spacing and dose distribution, determined automatically from the TPS. The relative dose weights of each energy layer were exported from the treatment plan report to calculate the LET.
The film responses for the six different dose intensities of the 112-MeV beam and their corresponding doses measured using an ionization chamber were plotted (blue inverted triangle) and fitted (red solid line) with the curve of Eq. (1) (Fig. 5). At a 2-cm depth, the LET was 1.54 keV/µm, and the corresponding
Table 2 . Fifth-degree Gaussian curve fitting parameters, as presented in Eq. (4).
xi | |||||
---|---|---|---|---|---|
2.95×105 | 1.78×105 | 3.12×105 | –1.53×105 | 4.43×105 | |
7.17 | 6.26 | 5.31 | 3.4 | 1.69 | |
0.9 | 1.04 | 1.91 | 0.18 | 5.09 |
Once
Table 3 . Range (R80) comparison of single-energy, multi-energy, and TPS-planned beams between the EBT3 film and reference measurements or calculations.
Energy | Single-energy | 70 MeV | 92 MeV | 112 MeV | 150 MeV | 170 MeV | 190 MeV |
---|---|---|---|---|---|---|---|
R80 (cm) | EBT3 | 3.79 | 6.38 | 9.22 | 15.66 | 19.54 | 23.63 |
Reference | 3.77 | 6.39 | 9.21 | 15.61 | 19.45 | 23.54 | |
Energy | Multi-energy | SOBP 110 MeV | SOBP 150 MeV | SOBP 190 MeV | TPS cubic plan | ||
R80 (cm) | EBT3 | 8.78 | 15.62 | 23.54 | 15.62 | ||
Reference | 8.85 | 15.54 | 23.45 | 15.53 |
The use of the Gafchromic EBT3 film for proton dosimetry has led to significant issues. In particular, the energy or LET dependence of the film has resulted in the underestimation of measured doses at low energies or for high-LET regions. As the bimolecular model is based on chemical reaction laws, which are affected by physical phenomena (known as the LET effect of protons), the performance of previously developed correction methods is adequate when the LET along the depth of the irradiated beam is known.
In this study, the bimolecular model was adapted and successfully extended for application to clinically meaningful multi-energy beams or SOBPs. The LET of a multienergy beam was calculated using an analytical method that considered the dose weight and LET of each single energy beam constituting the SOBP. The model was validated for a clinically important range of energies, from 70 to 190 MeV. The depth-dose distribution of the EBT3 film had an error within 4%, and it demonstrated submillimeter accuracy in range verification when compared with MLIC measurements and TPS calculations. The submillimeter accuracy of this method also enables range measurements of unknown beam energies, meaning that, conversely, an unknown beam energy can be calculated by measuring the range. The observed error of range verification contained uncertainties caused by the non-water characteristics of the EBT3 film material, and the effect of secondary particles on LET calculation was disregarded. Hence, the EBT3 film can be used with acceptable uncertainty using the investigated correction method.
Measuring the depth dose of a dynamically scanned proton beam using an ionization chamber is challenging owing to the extremely time-consuming dosimetry procedure. By constructing a generalized correction method for energy dependence, the EBT3 film becomes an excellent candidate for the depth-dose measurements of scanned beams owing to its efficient measurement procedure, which requires only one beam irradiation, and the high spatial resolution that can be obtained. A high spatial resolution provides excellent usability for low-energy beams that have steep distal fall-offs of the SOBP, such as eye treatment beams. A depth-wise dosimetric quality assurance (QA) is feasible for clinical beams using the EBT3 film. Hence, this could be a promising tool for patient-related QA. However, it still needs to be verified with a non-uniform depth-dose distribution by computing the LET at each location or depth. Calculating the correction parameter for the depth-dose distribution requires the dose weight of each energy along the depth, which can be obtained by analyzing the log file from the TPS. Therefore, further studies regarding nonuniform beams derived from TPS will warrant application to clinical cases.
This research was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2019M2A2B4096537 and 2019R1F1A1062775).
The authors have nothing to disclose.
All relevant data are within the paper and its Supporting Information files.
Conceptualization: Youngyih Han. Data curation: Moonhee Lee. Funding acquisition: Youngyih Han. Investigation: Moonhee Lee, Wonjoong Cheon, Sunghwan Ahn. Methodology: Moonhee Lee, Wonjoong Cheon. Project administration: Moonhee Lee, Youngyih Han. Resources: Youngyih Han. Supervision: Youngyih Han. Validation: Moonhee Lee. Visualization: Moonhee Lee. Writing–Original draft: Moonhee Lee. Writing–review & editing: Moonhee Lee, Younghyih Han.
Table 1 Specifications of irradiated proton beams
Specification | Pristine Bragg peak | SOBP | ||||||
---|---|---|---|---|---|---|---|---|
Snout Position (mm) | 400 | 400 | ||||||
Collimator (cm×cm) | 10×10 | 10×10 | ||||||
Energy/max. Energy of SOBP (MeV) | 70 | 92 | 112 | 150 | 190 | 110 | 150 | 190 |
Modulation width of SOBP (cm) | - | 5.39 | 4.27 | 5.34 |
Table 2 Fifth-degree Gaussian curve fitting parameters, as presented in Eq. (4)
xi | |||||
---|---|---|---|---|---|
2.95×105 | 1.78×105 | 3.12×105 | –1.53×105 | 4.43×105 | |
7.17 | 6.26 | 5.31 | 3.4 | 1.69 | |
0.9 | 1.04 | 1.91 | 0.18 | 5.09 |
Table 3 Range (R80) comparison of single-energy, multi-energy, and TPS-planned beams between the EBT3 film and reference measurements or calculations
Energy | Single-energy | 70 MeV | 92 MeV | 112 MeV | 150 MeV | 170 MeV | 190 MeV |
---|---|---|---|---|---|---|---|
R80 (cm) | EBT3 | 3.79 | 6.38 | 9.22 | 15.66 | 19.54 | 23.63 |
Reference | 3.77 | 6.39 | 9.21 | 15.61 | 19.45 | 23.54 | |
Energy | Multi-energy | SOBP 110 MeV | SOBP 150 MeV | SOBP 190 MeV | TPS cubic plan | ||
R80 (cm) | EBT3 | 8.78 | 15.62 | 23.54 | 15.62 | ||
Reference | 8.85 | 15.54 | 23.45 | 15.53 |
pISSN 2508-4445
eISSN 2508-4453
Formerly ISSN 1226-5829
Frequency: Quarterly