Abstract: We report our experience with extended usage of range compensators in heavy-ion radiotherapy with broad beams to lighten the management task of the beam data library, which is a collection of the standard beams to be referenced in treatment planning. Partly due to interference between lateral spreading and range shifting, as many as hundreds of beam entries may be necessary to cover all the possible clinical situations. We have introduced downstream secondary range shifting with a range compensator to reduce the interference and consequently to simplify the library. In our case, 30% reduction in beam entries is achieved without significantly degrading the beam quality nor increasing the material consumption by more than 3%, which is experimentally verified with carbon-ion beams or statistically estimated from the clinical records.
Abstract: The concept of internal target volume (ITV) is highly significant in radiotherapy for the lung, an organ which is hampered by organ motion. To date, different methods to obtain the ITV have been published and are therefore available. To define ITV, we developed a new method by adapting a time filter to the four-dimensional CT scan technique (4DCT) which is projection-data processing (4D projection data maximum attenuation (4DPM)), and compared it with reconstructed image processing (4D image maximum intensity projection (4DIM)) using a phantom and clinical evaluations. 4DIM and 4DPM captured accurate maximum intensity volume (MIV), that is tumour encompassing volume, easily. Although 4DIM increased the CT number 1.8 times higher than 4DPM, 4DPM provided the original tumour CT number for MIV via a reconstruction algorithm. In the patient with lung fibrosis honeycomb, the MIV with 4DIM is 0.7 cm larger than that for cine imaging in the cranio-caudal direction. 4DPM therefore provided an accurate MIV independent of patient characteristics and reconstruction conditions. These findings indicate the usefulness of 4DPM in determining ITV in radiotherapy.
Abstract: In carbon therapy, doses at center of spread-out Bragg peaks depend on field size. For a small field of 5 x 5 cm2, the central dose reduces to 96% of the central dose for the open field in case of 400 MeV/n carbon beam. Assuming the broad beam injected to the water phantom is made up of many pencil beams, the transverse dose distribution can be reconstructed by summing the dose distribution of the pencil beams. We estimated dose profiles of this pencil beam through measurements of dose distributions of broad uniform beams blocked half of the irradiation fields. The dose at a distance of a few cm from the edge of the irradiation field reaches up to a few percent of the central dose. From radiation quality measurements of this penumbra, the large-angle scattered particles were found to be secondary fragments which have lower LET than primary carbon beams. Carbon ions break up in beam modifying devices or in water phantom through nuclear interaction with target nuclei. The angular distributions of these fragmented nuclei are much broader than those of primary carbon particles. The transverse dose distribution of the pencil beam can be approximated by a function of the three-Gaussian form. For a simplest case of mono-energetic beam, contributions of the Gaussian components which have large mean deviations become larger as the depth in the water phantom increases.
Abstract: The clinical dose distributions of therapeutic carbon beams, currently used at NIRS HIMAC, are based on in-vitro Human Salivary Gland (HSG) cell survival response and clinical experience from fast neutron radiotherapy. Moderate radiosensitivity of HSG cells is expected to be a typical response of tumours to carbon beams. At first, the biological dose distribution is designed so as to cause a flat biological effect on HSG cells in the spread-out Bragg peak (SOBP) region. Then, the entire biological dose distribution is evenly raised in order to attain a RBE (relative biological effectiveness) = 3.0 at a depth where dose-averaged LET (linear energy transfer) is 80 keV/mum. At that point, biological experiments have shown that carbon ions can be expected to have a biological effect identical to fast neutrons, which showed a clinical RBE of 3.0 for fast neutron radiotherapy at NIRS. The resulting clinical dose distribution in this approximation is not dependent on dose level, tumour type or fractionation scheme and thus reduces the unknown parameters in the analysis of the clinical results. The width SOBP and the clinical / physical dose at the center of SOBP specify the dose distribution. The clinical results analysed in terms of TCP were found to show good agreement with the expected RBE value at higher TCP levels. The TCP analysis method was applied for the prospective dose estimation of hypofractionation.
Abstract: Clinical trials of carbon radiotherapy started at HIMAC in 1994 using three treatment rooms and four beam ports, two horizontal and two vertical. The broad beam method was adopted to make a three-dimensionally uniform field at an isocenter. A spot beam extracted from an accelerator was laterally spread out by using a pair of wobbler magnets and a scatterer. A bar ridge filter modulated the beam energy to obtain the spread out Bragg peak (SOBP). The SOBP was designed to be flat in terms of the biological dose based on the consideration that the field consisted of various beams with different LET. Finally, the field of 20 cm in diameter with +/- 2.5% uniformity was formed at the isocenter. The width of the maximum SOBP was 15 cm. When treating the lung or liver, organs that move due to breathing, the beam was irradiated only during the expiration period in a respiration-gated irradiation method. This reduced the treatment margin of the moving target. In order to prevent normal tissues adjacent to the target volume from irradiation by an unwanted dose, a layer-stacking method was developed. In this method, thin SOBP layers which have different ranges were piled up step by step from the distal end to the entrance of the target volume. At the same time, a multi-leaf collimator was used to change the aperture shape to match the shape of each layer to the cross-sectional shape of the target. This method has been applied to rather large volume cancers including bone and soft-tissue cancers. Only a few serious problems in the irradiation systems have been encountered since the beginning of the clinical trials. Overall the systems have been working stably and reliably.
Abstract: Although slow-rotation CT scanning (slow-scan CT: SSCT) has been used for radiation therapy planning, based on the rationale that the average duration of the human respiratory cycle is 4 s, a number of physical and quantitative questions require answering before it can be adopted for clinical use. This study was performed to evaluate SSCT physically in comparison with other scan methods, including respiratory-gated CT (RGCT), and to develop procedures to improve treatment accuracy. Evaluation items were geometrical accuracy, volume accuracy, water equivalent length and dose distribution using the 256-detector row CT with three scan methods. Fast-scan CT (FSCT) was defined as obtaining all respiratory phases in cine scan mode at 1.0 s per rotation. FSCT-ave was the averaged FSCT images in all respiratory phases, obtained by reconstructing short time intervals. SSCT has been defined as scanning with slow gantry rotation to capture the whole respiratory cycle in one rotation. RGCT was scanned at the most stable point in the respiratory cycle, which provides the same image as that by FSCT at the most stable point. Results showed that all evaluation items were dependent on motion characteristics. The findings of this study indicate that 3D planning based solely on SSCT under free breathing may result in underdosing of the target volume and increase toxicity to surrounding normal tissues. Of the three methods, RGCT showed the best ability to significantly increase the accuracy of dose distribution, and provided more information to minimize the margins. FSCT-ave is a satisfactory radiotherapy planning alternative if RGCT is not available.
Abstract: In the absence of a predictor of beam output in proton therapy using a broad beam, the beam output is obtained for individual treatments by calibrating the beam monitors. The calibration is carried out under conditions similar to the treatment conditions but with a phantom instead of the patient. However, the dose in the phantom a priori differs from that in the patient. In order to deliver the accurate dose, a correction factor has been introduced to correct the difference. This correction factor is referred to as a scatter factor in an analogy with photon therapy, and is defined as the ratio of the dose at the prescription point in the patient to the dose at the calibration point in the phantom. Under the calibration conditions at Hyogo Ion Beam Medical Center (HIBMC), the range compensator and the collimator, which are usually required in proton therapy with a broad beam, are not used. Therefore the scatter factor includes the effects of the devices as well as the difference between the dose in the patient and that in the phantom. We have developed an estimator using a dose calculation based on the pencil beam algorithm and implemented it in a treatment planning system (TPS) for clinical use. This estimator estimates the scatter factor by calculating the ratio of the doses under the same conditions in the TPS. In order to evaluate the performance of the estimator, demonstrations were carried out for cases with measurable outcomes using a gantry nozzle at HIBMC. We observed 2-3% differences between the measurements and the estimations. These differences were considered to result from the limitations of the dose calculation algorithm in modelling the beam and the patient.
Abstract: Treatment planning of heavy-ion radiotherapy involves predictive calculation of not only the physical dose but also the biological dose in a patient body. The biological dose is defined as the product of the physical dose and the relative biological effectiveness (RBE). In carbon-ion radiotherapy at National Institute of Radiological Sciences, the RBE value has been defined as the ratio of the 10% survival dose of 200 kVp x-rays to that of the radiation of interest for in vitro human salivary gland tumour cells. In this note, the physical and biological dose distributions of a typical therapeutic carbon-ion beam are calculated using the GEANT4 Monte Carlo simulation toolkit in comparison with those with the biological dose estimate system based on the one-dimensional beam model currently used in treatment planning. The results differed between the GEANT4 simulation and the one-dimensional beam model, indicating the physical limitations in the beam model. This study demonstrates that the Monte Carlo physics simulation technique can be applied to improve the accuracy of the biological dose distribution in treatment planning of heavy-ion radiotherapy.
Abstract: The commissioning of conformal radiotherapy system using heavy-ion beams at the Heavy Ion Medical Accelerator in Chiba (HIMAC) is described in detail. The system at HIMAC was upgraded for a clinical trial using a new technique: large spot uniform scanning with conformal layer stacking. The system was developed to localize the irradiation dose to the target volume more effectively than with the old system. With the present passive irradiation method using a ridge filter, a scatterer, a pair of wobbler magnets, and a multileaf collimator, the width of the spread-out Bragg peak (SOBP) in the radiation field could not be changed. With dynamic control of the beam-modifying devices during irradiation, a more conformal radiotherapy could be achieved. In order to safely perform treatments with this conformal therapy, the moving devices should be watched during irradiation and the synchronousness among the devices should be verified. This system, which has to be safe for patient irradiations, was constructed and tested for safety and for the quality of the dose localization realized. Through these commissioning tests, we were successfully able to prepare the conformal technique using layer stacking for patients. Subsequent to commissioning the technique has been applied to patients in clinical trials.
Abstract: We have developed a simple collimator model to improve the accuracy of penumbra behaviour in pencil-beam dose calculation for proton radiotherapy. In this model, transmission of particles through a three-dimensionally extended opening of a collimator is calculated in conjunction with phase-space distribution of the particles. Comparison of the dose distributions calculated using the new three-dimensional collimator model and the conventional two-dimensional model to lateral dose profiles experimentally measured with collimated proton beams showed the superiority of the new model over the conventional one.
Abstract: The idea of a computer-controlled range-compensating system for heavy charged particle radiotherapy, the multibar compensator, is proposed. By stacking multiple energy-absorbing layers along the beam, each of which has structure and behaviour similar to those of a multileaf collimator, variable range compensation will be achieved. The analysis of the conventional range compensators actually used for treatment concluded that the proposed system would not seriously degrade the treatment quality for the most cases, except for tumours in the head and neck region where 1 mm precision may be required. The system will even be able to coexist with the conventional range compensators to provide either method depending on clinical situations.
Abstract: The goal of radiotherapy is not only to apply a high radiation dose to a tumor, but also to avoid side effects in the surrounding healthy tissue. Therefore, it is important for carbon-ion treatment planning to calculate accurately the effects of the lateral penumbra. In this article, for wobbled beams under various irradiation conditions, we focus on the lateral penumbras at several aperture positions of one side leaf of the multileaf collimator. The penumbras predicted by an analytical penumbra calculation model were compared with the measured results. The results calculated by the model for various conditions agreed well with the experimental ones. In conclusion, we found that the analytical penumbra calculation model could predict accurately the measured results for wobbled beams and it was useful for carbon-ion treatment planning to apply the model.
Abstract: Hot- and cold-dose spots at a shallow depth in a target are formed by carbon ions passing through the bolus with sharp gradients. These spots are caused by sidescatter disequilibrium due to various multiple scattering effects in the different bolus thicknesses. When the dose calculation method by the broad beam algorithm (BBA) is used for treatment planning, these spots cannot be predicted, because the BBA neglects the multiple scattering effects in materials (rms error of 3.9%). On the other hand, since the dose calculation method by the pencil beam algorithm (PBA) takes into account the scattering effects, the results calculated by the PBA agreed better than the BBA with the measured hot- and cold-dose spots, having a rms error of 1.9%. Thus, dose calculation by the PBA improves the accuracy of dose prediction at the shallow depth. However, since dose distributions at deeper positions are affected by many light fragment particles generated by fragment reactions, the results calculated by the PBA disagree with the experimental ones. It is necessary that even the PBA accurately models behavior of fragment particles.
Abstract: A method to establish the relationship between CT number and effective density for therapeutic radiations is proposed. We approximated body tissues to mixtures of muscle, air, fat and bone. Consequently, the relationship can be calibrated only with a CT scan of their substitutes, for which we chose water, air, ethanol and potassium phosphate solution, respectively. With simple and specific corrections for non-equivalencies of the substitutes, a calibration accuracy of 1% will be achieved. We tested the calibration method with some biological materials to verify that the proposed method would offer the accuracy, simplicity and specificity required for a standard in radiotherapy treatment planning, in particular with heavy charged particles.
Abstract: A guideline for quality assurance of CT systems for radiotherapy treatment planning is proposed. Quantitative interpretation of CT number is very important especially for range calculation in heavy charged particle radiotherapy, for which we adopted the polybinary calibration method to correct variations among CT systems and scanning conditions. Practical procedures for commissioning and constancy testing are documented along with the methodologies against various sources of uncertainty. We propose this guideline for quality assurance of heavy charged particle radiotherapy.
Abstract: We have upgraded a heavy-ion radiotherapy treatment-planning system to adapt for the layer-stacking irradiation method, which is to conform a variable spread-out Bragg peak to a target volume by means of dynamic control of the conventional beam-modifying devices. The biophysical model, the beam-setup logic, and the dose-calculation algorithm implemented for the layer-stacking method are described and the expected clinical usability is discussed. The layer-stacking method was integrated in perfect accordance with the ongoing conventional treatments so that the established protocols, which are the clinically optimized dose fractionation schemes, will still be valid. On the other hand, a simulation study indicated a substantial improvement of dose distribution with the layer-stacking method though the significance may depend on the size, shape, and location of the tumor. The completed treatment system will provide an option for improved conformal radiotherapy without interfering with the conventional method and we expect a gradual expansion of the clinical cases applicable to the layer-stacking method.
Abstract: We report results from Experiment 871, performed at the BNL AGS, of a measurement of the branching ratio K(0)(L)-->&mgr;(+)&mgr;(-) with respect to the CP-violating mode K(0)(L)-->pi(+)pi(-). This experiment detected over 6200 candidate &mgr;(+)&mgr;(-) events, a factor of 6 more than that seen in all previous measurements combined. The resulting branching ratio gamma(K(0)(L)-->&mgr;(+)&mgr;(-))/gamma(K(0)(L)-->pi(+)pi(-)) = (3. 474+/-0.057)x10(-6) leads to a branching fraction B(K(0)(L)-->&mgr;(+)&mgr;(-)) = (7.18+/-0.17)x10(-9), which is consistent with the current world average, and reduces the uncertainty in this decay mode by a factor of 3.
