Taisuke Takayanagi

Study on the use of electron-tracking Compton gamma-ray camera to monitor the therapeutic proton dose distribution in real time

Radiation therapy with proton and heavy-ion beams has been better established lately and the patient throughput is increasing. Although the therapy beam is controlled with high accuracy, it is difficult to know the location of distal dose falloff in the body. If real-time monitoring of the location is realized, the treatment quality would be improved. We have developed an electron-tracking Compton camera (ETCC) for real-time monitoring on the proton therapy. Our ETCC has a wide energy dynamic range of 200-1300 keV and a wide field of view. Therefore, ETCC has a potential as a quality assurance tool for proton therapy. We simulated and conducted an experiment with a 155 MeV proton beam and a water phantom. We succeeded in imaging a Bragg peak with prompt gamma rays.

Biological effect of dose distortion by fiducial markers in spot-scanning proton therapy with a limited number of fields: A simulation study

PURPOSE In accurate proton spot-scanning therapy, continuous target tracking by fluoroscopic x ray during irradiation is beneficial not only for respiratory moving tumors of lung and liver but also for relatively stationary tumors of prostate. Implanted gold markers have been used with great effect for positioning the target volume by a fluoroscopy, especially for the cases of liver and prostate with the targets surrounded by water-equivalent tissues. However, recent studies have revealed that gold markers can cause a significant underdose in proton therapy. This paper focuses on prostate cancer and explores the possibility that multiple-field irradiation improves the underdose effect by markers on tumor-control probability (TCP). METHODS A Monte Carlo simulation was performed to evaluate the dose distortion effect. A spherical gold marker was placed at several characteristic points in a water phantom. The markers were with two different diameters of 2 and 1.5 mm, both visible on fluoroscopy. Three beam arrangements of single-field uniform dose (SFUD) were examined: one lateral field, two opposite lateral fields, and three fields (two opposite lateral fields + anterior field). The relative biological effectiveness (RBE) was set to 1.1 and a dose of 74 Gy (RBE) was delivered to the target of a typical prostate size in 37 fractions. The ratios of TCP to that without the marker (TCP(r)) were compared with the parameters of the marker sizes, number of fields, and marker positions. To take into account the dependence of biological parameters in TCP model, α∕β values of 1.5, 3, and 10 Gy (RBE) were considered. RESULTS It was found that the marker of 1.5 mm diameter does not affect the TCPs with all α∕β values when two or more fields are used. On the other hand, if the marker diameter is 2 mm, more than two irradiation fields are required to suppress the decrease in TCP from TCP(r) by less than 3%. This is especially true when multiple (two or three) markers are used for alignment of a patient. CONCLUSIONS It is recommended that 1.5-mm markers be used to avoid the reduction of TCP as well as to spare the surrounding critical organs, as long as the markers are visible on x-ray fluoroscopy. When 2-mm markers are implanted, more than two fields should be used and the markers should not be placed close to the distal edge of any of the beams.

Evaluation of the influence of double and triple Gaussian proton kernel models on accuracy of dose calculations for spot scanning technique.

PURPOSE The main purpose in this study was to present the results of beam modeling and how the authors systematically investigated the influence of double and triple Gaussian proton kernel models on the accuracy of dose calculations for spot scanning technique. METHODS The accuracy of calculations was important for treatment planning software (TPS) because the energy, spot position, and absolute dose had to be determined by TPS for the spot scanning technique. The dose distribution was calculated by convolving in-air fluence with the dose kernel. The dose kernel was the in-water 3D dose distribution of an infinitesimal pencil beam and consisted of an integral depth dose (IDD) and a lateral distribution. Accurate modeling of the low-dose region was important for spot scanning technique because the dose distribution was formed by cumulating hundreds or thousands of delivered beams. The authors employed a double Gaussian function as the in-air fluence model of an individual beam. Double and triple Gaussian kernel models were also prepared for comparison. The parameters of the kernel lateral model were derived by fitting a simulated in-water lateral dose profile induced by an infinitesimal proton beam, whose emittance was zero, at various depths using Monte Carlo (MC) simulation. The fitted parameters were interpolated as a function of depth in water and stored as a separate look-up table. These stored parameters for each energy and depth in water were acquired from the look-up table when incorporating them into the TPS. The modeling process for the in-air fluence and IDD was based on the method proposed in the literature. These were derived using MC simulation and measured data. The authors compared the measured and calculated absolute doses at the center of the spread-out Bragg peak (SOBP) under various volumetric irradiation conditions to systematically investigate the influence of the two types of kernel models on the dose calculations. RESULTS The authors investigated the difference between double and triple Gaussian kernel models. The authors found that the difference between the two studied kernel models appeared at mid-depths and the accuracy of predicting the double Gaussian model deteriorated at the low-dose bump that appeared at mid-depths. When the authors employed the double Gaussian kernel model, the accuracy of calculations for the absolute dose at the center of the SOBP varied with irradiation conditions and the maximum difference was 3.4%. In contrast, the results obtained from calculations with the triple Gaussian kernel model indicated good agreement with the measurements within ±1.1%, regardless of the irradiation conditions. CONCLUSIONS The difference between the results obtained with the two types of studied kernel models was distinct in the high energy region. The accuracy of calculations with the double Gaussian kernel model varied with the field size and SOBP width because the accuracy of prediction with the double Gaussian model was insufficient at the low-dose bump. The evaluation was only qualitative under limited volumetric irradiation conditions. Further accumulation of measured data would be needed to quantitatively comprehend what influence the double and triple Gaussian kernel models had on the accuracy of dose calculations.

Reduction of the number of stacking layers in proton uniform scanning

Uniform scanning with a relatively large beam size can improve beam utilization efficiency more than conventional irradiation methods using scatterers and can achieve a large-field, long-range and large spread-out Bragg peak (SOBP). The SOBP is obtained by energy stacking in uniform scanning, but its disadvantage is that the number of stacking layers is large, especially in the low-energy region, because the Bragg peak of the pristine beam is very sharp. We applied a mini-ridge filter to broaden the pristine Bragg peak up to a stacked layer thickness of 1 or 2 cm in order to decrease the number of stacking layers. The number of stacking layers can be reduced to 20% or less than that in the case of pristine beam stacking. Although the distal falloff of the SOBP is deteriorated by applying the mini-ridge filter, we can improve the distal falloff to that of pristine beam stacking by introducing the distal filter to the irradiation of the most distal layer. Uniform scanning in combination with mini-ridge filter use can more than double the beam utilization efficiency over that of passive irradiation techniques.

Development and evaluation of a short-range applicator for treating superficial moving tumors with respiratory-gated spot-scanning proton therapy using real-time image guidance.

Treatment of superficial tumors that move with respiration (e.g. lung tumors) using spot-scanning proton therapy (SSPT) is a high-priority research area. The recently developed real-time image-gated proton beam therapy (RGPT) system has proven to be useful for treating moving tumors deep inside the liver. However, when treating superficial tumors, the protons range is small and so is the sizes of range straggling, making the Bragg-peaks extremely sharp compared to those located in deep-seated tumors. The extreme sharpness of Bragg-peaks is not always beneficial because it necessitates a large number of energy layers to make a spread-out Bragg-peak, resulting in long treatment times, and is vulnerable to motion-induced dose deterioration. We have investigated a method to treat superficial moving tumors in the lung by the development of an applicator compatible with the RGPT system. A mini-ridge filter (MRF) was developed to broaden the pristine Bragg-peak and, accordingly, decrease the number of required energy layers to obtain homogeneous irradiation. The applicator position was designed so that the fiducial markers trajectory can be monitored by fluoroscopy during proton beam-delivery. The treatment plans for three lung cancer patients were made using the applicator, and four-dimensional (4D) dose calculations for the RGPT were performed using patient respiratory motion data. The effect of the MRF on the dose distributions and treatment time was evaluated. With the MRF, the number of energy layers was decreased to less than half of that needed without it, whereas the target volume coverage values (D99%, D95%, D50%, D2%) changed by less than 1% of the prescribed dose. Almost no dose distortion was observed after the 4D dose calculation, whereas the treatment time decreased by 26%-37%. Therefore, we conclude that the developed applicator compatible with RGPT is useful to solve the issue in the treatment of superficial moving tumors with SSPT.

Design of a ridge filter structure based on the analysis of dose distributions.

Dose distributions distorted by a periodic structure, such as a ridge filter, are analytically investigated. Based on the beam optics, the fluence distributions of scanned beams passing through the ridge filter are traced. It is shown that the periodic lateral dose distribution blurred by multiple Coulomb scattering can be expressed by a sum of cosine functions through Fourier transform. The result shows that the dose homogeneity decreases exponentially as the period of the structure becomes longer. This analysis is applied to the example case of a mini-ridge filter. The mini-ridge filter is designed to broaden sharp Bragg peaks for an energy-stacking irradiation method. The dose distributions depend on the period of the ridge filter structure and the angular straggling at the ridge filter position. Several cases are prepared where the period and angular straggling are supposed to be probable values. In these cases, the lateral distributions obtained by the analytical method are compared to Monte Carlo simulation results. Both distributions show good agreement with each other within 1%, which means that this analysis allows estimation of the dose distribution downstream of the ridge filter quantitatively. The appropriate period of grooves and scatterer width can be determined which ensures sufficient homogeneity.

TH-C-BRD-04: Beam Modeling and Validation with Triple and Double Gaussian Dose Kernel for Spot Scanning Proton Beams

PURPOSE To present the validity of our beam modeling with double and triple Gaussian dose kernels for spot scanning proton beams in Nagoya Proton Therapy Center. This study investigates the conformance between the measurements and calculation results in absolute dose with two types of beam kernel. METHODS A dose kernel is one of the important input data required for the treatment planning software. The dose kernel is the 3D dose distribution of an infinitesimal pencil beam of protons in water and consists of integral depth doses and lateral distributions. We have adopted double and triple Gaussian model as lateral distribution in order to take account of the large angle scattering due to nuclear reaction by fitting simulated inwater lateral dose profile for needle proton beam at various depths. The fitted parameters were interpolated as a function of depth in water and were stored as a separate look-up table for the each beam energy. The process of beam modeling is based on the method of MDACC [X.R.Zhu 2013]. RESULTS From the comparison results between the absolute doses calculated by double Gaussian model and those measured at the center of SOBP, the difference is increased up to 3.5% in the high-energy region because the large angle scattering due to nuclear reaction is not sufficiently considered at intermediate depths in the double Gaussian model. In case of employing triple Gaussian dose kernels, the measured absolute dose at the center of SOBP agrees with calculation within ±1% regardless of the SOBP width and maximum range. CONCLUSION We have demonstrated the beam modeling results of dose distribution employing double and triple Gaussian dose kernel. Treatment planning system with the triple Gaussian dose kernel has been successfully verified and applied to the patient treatment with a spot scanning technique in Nagoya Proton Therapy Center.

A simplified Monte Carlo algorithm considering large‐angle scattering for fast and accurate calculation of proton dose

Abstract Purpose The purpose of this study is to improve dose calculation accuracy of the simplified Monte Carlo (SMC) algorithm in the low‐dose region. Because conventional SMC algorithms calculate particle scattering in consideration of multiple Coulomb scattering (MCS) only, they approximate lateral dose profiles by a single Gaussian function. However, it is well known that the low‐dose region spreads away from the beam axis, and it has been pointed out that modeling of the low‐dose region is important to calculated dose accurately. Methods A SMC algorithm, which is named modified SMC and considers not only MCS but also large angle scattering resembling hadron elastic scattering, was developed. In the modified SMC algorithm, the particle fluence varies in the longitudinal direction because the large‐angle scattering decreases residual range of particles in accordance with their scattering angle and tracking of the particles with large scattering angle is terminated at a short distance downstream from the scattering. Therefore, modified integrated depth dose (m‐IDD) tables, which are converted from measured IDD in consideration of the fluence loss, are used to calculate dose. Results In the case of a 1‐liter cubic target, the calculation accuracy was improved in comparison with that of a conventional algorithm, and the modified algorithm results agreed well with Geant4‐based simulation results; namely, 98.8% of the points satisfied the 2% dose/2 mm distance‐to‐agreement (DTA) criterion. The calculation time of the modified SMC algorithm was 1972 s in the case of 4.4 × 108 particles and 16‐threading operation of an Intel Xeon E5‐2643 (3.3‐GHz clock). Conclusions An SMC algorithm that can reproduce a laterally widespread low‐dose region was developed. According to the comparison with a Geant4‐based simulation, it was concluded that the modified SMC algorithm is useful for calculating dose of proton radiotherapy.

Dual ring multilayer ionization chamber and theory-based correction technique for scanning proton therapy

PURPOSE To develop a multilayer ionization chamber (MLIC) and a correction technique that suppresses differences between the MLIC and water phantom measurements in order to achieve fast and accurate depth dose measurements in pencil beam scanning proton therapy. METHODS The authors distinguish between a calibration procedure and an additional correction: 1-the calibration for variations in the air gap thickness and the electrometer gains is addressed without involving measurements in water; 2-the correction is addressed to suppress the difference between depth dose profiles in water and in the MLIC materials due to the nuclear interaction cross sections by a semiempirical model tuned by using measurements in water. In the correction technique, raw MLIC data are obtained for each energy layer and integrated after multiplying them by the correction factor because the correction factor depends on incident energy. The MLIC described here has been designed especially for pencil beam scanning proton therapy. This MLIC is called a dual ring multilayer ionization chamber (DRMLIC). The shape of the electrodes allows the DRMLIC to measure both the percentage depth dose (PDD) and integrated depth dose (IDD) because ionization electrons are collected from inner and outer air gaps independently. RESULTS IDDs for which the beam energies were 71.6, 120.6, 159, 180.6, and 221.4 MeV were measured and compared with water phantom results. Furthermore, the measured PDDs along the central axis of the proton field with a nominal field size of 10 × 10 cm(2) were compared. The spread out Bragg peak was 20 cm for fields with a range of 30.6 and 3 cm for fields with a range of 6.9 cm. The IDDs measured with the DRMLIC using the correction technique were consistent with those that of the water phantom; except for the beam energy of 71.6 MeV, all of the points satisfied the 1% dose/1 mm distance to agreement criterion of the gamma index. The 71.6 MeV depth dose profile showed slight differences in the shallow region, but 94.5% of the points satisfied the 1%/1 mm criterion. The 90% ranges, defined at the 90% dose position in distal fall off, were in good agreement with those in the water phantom, and the range differences from the water phantom were less than ±0.3 mm. The PDDs measured with the DRMLIC were also consistent with those that of the water phantom; 97% of the points passed the 1%/1 mm criterion. CONCLUSIONS It was demonstrated that the new correction technique suppresses the difference between the depth dose profiles obtained with the MLIC and those obtained from a water phantom, and a DRMLIC enabling fast measurements of both IDD and PDD was developed. The IDDs and PDDs measured with the DRMLIC and using the correction technique were in good agreement with those that of the water phantom, and it was concluded that the correction technique and DRMLIC are useful for depth dose profile measurements in pencil beam scanning proton therapy.

SU‐E‐T‐561: Development of Depth Dose Measurement Technique Using the Multilayer Ionization Chamber for Spot Scanning Method

PURPOSE To develop a measurement technique which suppresses the difference between profiles obtained with a multilayer ionization chamber (MLIC) and with a water phantom. METHODS The developed technique multiplies the raw MLIC data by a correction factor that depends on the initial beam range and water equivalent depth. The correction factor is derived based on a Bragg curve calculation formula considering range straggling and fluence loss caused by nuclear reactions. Furthermore, the correction factor is adjusted based on several integrated depth doses measured with a water phantom and the MLIC. The measured depth dose profiles along the central axis of the proton field with a nominal field size of 10 by 10 cm were compared between the MLIC using the new technique and the water phantom. The spread out Bragg peak was 20 cm for fields with a range of 30.6 cm and 6.9 cm. Raw MLIC data were obtained with each energy layer, and integrated after multiplying by the correction factor. The measurements were performed by a spot scanning nozzle at Nagoya Proton Therapy Center, Japan. RESULTS The profile measured with the MLIC using the new technique is consistent with that of the water phantom. Moreover, 97% of the points passed the 1% dose /1mm distance agreement criterion of the gamma index. CONCLUSION We have demonstrated that the new technique suppresses the difference between profiles obtained with the MLIC and with the water phantom. It was concluded that this technique is useful for depth dose measurement in proton spot scanning method.