Moving target irradiation with fast rescanning and gating in particle therapy.

PURPOSE: In moving target irradiation with pencil beam scanning, the interplay effect between the target motion and the scanned beam is a problem because this effect causes over or under dosage in the target volume. To overcome this, we have studied rescanning using a gating technique. METHODS: A simulation and experimental study was carried out. In the experiment, we used the fast scanning system developed at the HIMAC to verify the validity of phase controlled rescanning method, in which the time for rescanning irradiation of each slice is matched to the gating duration. RESULTS: Simulation and experimental results showed that controlling the scan speed to match the respiration cycle with rescans can deliver the blurred dose distribution. In the comparison between the static measurements and the moving measurements with the phase controlled rescanning method, the dose difference was less than 2% for pinpoint chambers in the target volume. CONCLUSIONS: The simulation and experimental results demonstrated that the phase controlled rescanning method makes it possible to deliver the dose distribution close to the expected one. As an experimental result for 3D irradiation, it was estimated that blurring by the probability density function was not only for a lateral distribution, but also for a distal distribution, even in the lateral rescanning.

Performance of the NIRS fast scanning system for heavy-ion radiotherapy.

PURPOSE: A project to construct a new treatment facility, as an extension of the existing HIMAC facility, has been initiated for the further development of carbon-ion therapy at NIRS. This new treatment facility is equipped with a 3D irradiation system with pencil-beam scanning. The challenge of this project is to realize treatment of a moving target by scanning irradiation. To achieve fast rescanning within an acceptable irradiation time, the authors developed a fast scanning system. METHODS: In order to verify the validity of the design and to demonstrate the performance of the fast scanning prior to use in the new treatment facility, a new scanning-irradiation system was developed and installed into the existing HIMAC physics-experiment course. The authors made strong efforts to develop (1) the fast scanning magnet and its power supply, (2) the high-speed control system, and (3) the beam monitoring. The performance of the system including 3D dose conformation was tested by using the carbon beam from the HIMAC accelerator. RESULTS: The performance of the fast scanning system was verified by beam tests. Precision of the scanned beam position was less than +/-0.5 mm. By cooperating with the planning software, the authors verified the homogeneity of the delivered field within +/-3% for the 3D delivery. This system took only 20 s to deliver the physical dose of 1 Gy to a spherical target having a diameter of 60 mm with eight rescans. In this test, the average of the spot-staying time was considerably reduced to 154 micros, while the minimum staying time was 30 micros. CONCLUSIONS: As a result of this study, the authors verified that the new scanning delivery system can produce an accurate 3D dose distribution for the target volume in combination with the planning software.

The basic study of a bi-material range compensator for improving dose uniformity for proton therapy.

A range compensator (abbreviated as a RC hereafter) is used to form a conformal dose distribution for heavy-charged-particle therapy. However, it induces distortion of the dose distribution. The induced inhomogeneity may result in a calibration error of a monitor unit (MU) assigned to a transmission ionization chamber. By using a bi-material RC made from a low-Z material and a high-Z material instead of the regular RC, the dose inhomogeneity has been obviously reduced by equalizing the lateral dose distributions formed by pencil beams traversing elements of the RC with different base thicknesses at the same water-equivalent depth. We designed and manufactured a 4 x 4 matrix-shaped single-material RC and a bi-material RC with the same range losses at corresponding elements of the RCs. The bi-material RC is made from chemical wood (the main chemical component is an ABS resin) as a low-Z material and from brass as a high-Z material. Sixteen segments of the RC are designed so that the range-loss differences of the adjacent segments of the RC range from 0 to 50 mm in steps of 5 mm. We measured dose distributions in water formed by a 160 MeV proton beam traversing the single-material RC or the bi-material RC, using the HIMAC biology beam port. Large dips and bumps were observed in the dose distribution formed by the use of the single-material RC; the dose uniformity has been significantly improved in the target region by the use of the bi-material RC. The improvement has been obtained at the expense of blurring lateral penumbra. For clinical application of this method to a patient with large density inhomogeneity, a simple modification method of the original calculation model has been given.

Estimation of late rectal normal tissue complication probability parameters in carbon ion therapy for prostate cancer.

PURPOSE: The aim of this study was to estimate normal tissue complication probability (NTCP) parameters for late rectal complications after carbon ion radiotherapy (C-ion RT) for prostate cancer. METHODS AND MATERIALS: A total of 163 patients were used to derive NTCP parameters. These patients were treated with relative biological effectiveness (RBE)-weighted dose ranging from 57.6 Gy (RBE) up to 72 Gy (RBE) and included in dose escalation trials. The Lyman-Kutcher-Burman (LKB) model was used and the model parameters were fit to the relation between dose and complication observed after C-ion RT. RESULTS: The resulting NTCP parameters were the volume effect parameter; n=0.035 (95% CI: 0.024-0.047), the steepness of the NTCP curve; m=0.10 (0.084-0.13), the tolerance dose associated with 50% probability of complication; TD50=63.6 Gy (RBE) (61.8-65.4 Gy (RBE)) for Grade⩾1, n=0.012 (0.0050-0.023), m=0.046 (0.033-0.062), TD50=69.1 Gy (RBE) (67.6-70.9 Gy (RBE)) for Grade⩾2. CONCLUSION: A new set of rectal NTCP parameters in C-ion RT was determined. The rather small n values suggest that the rectum was consistent with being strictly serial organ. The new derived parameter values facilitate estimation of rectal NTCP in C-ion RT.

Influence of nuclear interactions in polyethylene range compensators for carbon-ion radiotherapy.

PURPOSE: A recent study revealed that polyethylene (PE) would cause extra carbon-ion attenuation per range shift by 0.45%/cm due to compositional differences in nuclear interactions. The present study aims to assess the influence of PE range compensators on tumor dose in carbon-ion radiotherapy. METHODS: Carbon-ion radiation was modeled to be composed of primary carbon ions and secondary particles, for each of which the dose and the relative biological effectiveness (RBE) were estimated at a tumor depth in the middle of spread-out Bragg peak. Assuming exponential behavior for attenuation and yield of these components with depth, the PE effect on dose was calculated for clinical carbon-ion beams and was partly tested by experiment. The two-component model was integrated into a treatment-planning system and the PE effect was estimated in two clinical cases. RESULTS: The attenuation per range shift by PE was 0.1%-0.3%/cm in dose and 0.2%-0.4%/cm in RBE-weighted dose, depending on energy and range-modulation width. This translates into reduction of RBE-weighted dose by up to 3% in extreme cases. In the treatment-planning study, however, the effect on RBE-weighted dose to tumor was typically within 1% reduction. CONCLUSIONS: The extra attenuation of primary carbon ions in PE was partly compensated by increased secondary particles for tumor dose. In practical situations, the PE range compensators would normally cause only marginal errors as compared to intrinsic uncertainties in treatment planning, patient setup, beam delivery, and clinical response.

[Study of fast acquisition technique for heavy ion CT system based on the measurement of residual range distribution].

We proposed a new technique for the fast acquisition heavy ion CT (HICT) system based on the measurement of residual range distribution using an intensifying screen and charge coupled device camera. The previously used fast acquisition HICT system had poor electron density resolution. In the new technique, the range shifter thickness is varied over the required dynamic range in the spill of the heavy ion beam at each projection angle and the residual range distribution is determined by a series of acquisition data. We examined the image quality using the contrast noise ratio and the noise power spectrum, and estimated the electron density resolution, using a low-contrast phantom for measurement of electron density resolution. The image quality of the new technique was superior to that of the previous fast acquisition HICT system. Furthermore, the relative electron density resolution was 0.011, which represented an improvement of about 12-fold. Therefore we showed that the new technique was potentially useful in clinical use of HICT, including treatment and quality assurance of heavy ion therapy.

Improvement of spread-out Bragg peak flatness for a carbon-ion beam by the use of a ridge filter with a ripple filter.

We have developed a novel design method of ridge filters for carbon-ion therapy using a broad-beam delivery system to improve the flatness of a biologically effective dose in the spread-out Bragg peak (SOBP). So far, the flatness of the SOBP is limited to about ±5% for carbon beams since the weight control of component Bragg curves composing the SOBP is difficult. This difficulty arises from using a large number of ridge-bar steps (e.g. about 100 for a SOBP width of 60 mm) required to form the SOBP for the pristine Bragg curve with an extremely sharp distal falloff. Instead of using a single ridge filter, we introduce a ripple filter to broaden the Bragg peak so that the number of ridge-bar steps can be reduced to about 30 for SOBP with of 60 mm for the ridge filter designed for the broadened Bragg peak. Thus we can manufacture the ridge filter more accurately and then attain a better flatness of the SOBP due to well-controlled weights of the component Bragg curves. We placed the ripple filter on the same frame of the ridge filter and arranged the direction of the ripple-filter-bar array perpendicular to that of the ridge-filter-bar array. We applied this method to a 290 MeV u(-1) carbon-ion beam in Heavy Ion Medical Accelerator in Chiba and verified the effectiveness by measurements.

Microdosimetric approach to NIRS-defined biological dose measurement for carbon-ion treatment beam.

The RBE-weighted absorbed dose, called "biological dose," has been routinely used for carbon-ion treatment planning in Japan to formulate dose prescriptions for treatment protocols. This paper presents a microdosimetric approach to measuring the biological dose, which was redefined to be derived from microdosimetric quantities measured by a tissue-equivalent proportional counter (TEPC). The TEPC was calibrated in (60)Co gamma rays to assure a traceability of the TEPC measurement to Japanese standards and to eliminate the discrepancies among matching counters. The absorbed doses measured by the TEPC were reasonably coincident with those measured by a reference ionization chamber. The RBE value was calculated from the microdosimetric spectrum on the basis of the microdosimetric kinetic model. The biological doses obtained by the TEPC were compared with those prescribed in the carbon-ion treatment planning system. We found that it was reasonable for the measured biological doses to decrease with depth around the rear SOBP region because of beam divergence, scattering effect, and fragmentation reaction. These results demonstrate that the TEPC can be an effective tool to assure the radiation quality in carbon-ion radiotherapy.

Preliminary calculation of RBE-weighted dose distribution for cerebral radionecrosis in carbon-ion treatment planning.

Cerebral radionecrosis is a significant side effect in radiotherapy for brain cancer. The purpose of this study is to calculate the relative biological effectiveness (RBE) of carbon-ion beams on brain cells and to show RBE-weighted dose distributions for cerebral radionecrosis speculation in a carbon-ion treatment planning system. The RBE value of the radionecrosis for the carbon-ion beam is calculated by the modified microdosimetric kinetic model on the assumption of a typical clinical α/ß ratio of 2 Gy for cerebral radionecrosis in X-rays. This calculation method for the RBE-weighted dose is built into the treatment planning system for the carbon-ion radiotherapy. The RBE-weighted dose distributions are calculated on computed tomography (CT) images of four patients who had been treated by carbon-ion radiotherapy for astrocytoma (WHO grade 2) and who suffered from necrosis around the target areas. The necrotic areas were detected by brain scans via magnetic resonance imaging (MRI) after the treatment irradiation. The detected necrotic areas are easily found near high RBE-weighted dose regions. The visual comparison between the RBE-weighted dose distribution and the necrosis region indicates that the RBE-weighted dose distribution will be helpful information for the prediction of radionecrosis areas after carbon-ion radiotherapy.

Improved dose-calculation accuracy in proton treatment planning using a simplified Monte Carlo method verified with three-dimensional measurements in an anthropomorphic phantom.

Treatment planning for proton tumor therapy requires a fast and accurate dose-calculation method. We have implemented a simplified Monte Carlo (SMC) method in the treatment planning system of the National Cancer Center Hospital East for the double-scattering beam delivery scheme. The SMC method takes into account the scattering effect in materials more accurately than the pencil beam algorithm by tracking individual proton paths. We confirmed that the SMC method reproduced measured dose distributions in a heterogeneous slab phantom better than the pencil beam method. When applied to a complex anthropomorphic phantom, the SMC method reproduced the measured dose distribution well, satisfying an accuracy tolerance of 3 mm and 3% in the gamma index analysis. The SMC method required approximately 30 min to complete the calculation over a target volume of 500 cc, much less than the time required for the full Monte Carlo calculation. The SMC method is a candidate for a practical calculation technique with sufficient accuracy for clinical application.

Analytical design method of optimum ridge filters for wobbled and collimated proton beams.

A novel design method of ridge filters (RFs) has been developed for general proton beam lines which use a single-radius beam wobbling method. It can be applied to beam lines that transport both protons and carbon ions which are about three times longer than regular beam lines dedicated to protons. We designed an RF with an SOBP (spread-out Bragg peak) width of 60 mm in water for the 160-MeV proton beam of the HIMAC (Heavy Ion Medical Accelerator in Chiba) biology beam line using an existing model of the RF. Yet we observed a slope in the SOBP region when we used the RF. To elucidate the source of the slope, we have developed a new calculation model taking into account the geometry of the RF and a beam-limiting device. The source for the slope was found to be the large scattering effect of protons in the RF and beam restriction by a ring collimator (aperture diameter: 160 mm) placed just before the RF. When both fluence reduction by the scattering effect of protons in the RF and the beam-collimation effect are taken into account, proper RFs can be designed universally for a given beam line arrangement using the single-radius beam-wobbling method from the start without any trial-and-error process. This will serve to reduce the commissioning time of newly designed beam delivery systems.