Search

Purpose: The measurement of depth dose profiles for range and energy verification of heavy ion beams is an important aspect of quality assurance procedures for heavy ion therapy facilities. The steep dose gradients in the Bragg peak region of these profiles require the use of detectors with high spatial resolution. The aim of this work is to characterize a one dimensional monolithic silicon detector array called the "serial Dose Magnifying Glass" (sDMG) as an independent ion beam energy and range verification system used for quality assurance conducted for ion beams used in heavy ion therapy., Methods: The sDMG detector consists of two linear arrays of 128 silicon sensitive volumes each with an effective size of 2mm × 50μm × 100μm fabricated on a p-type substrate at a pitch of 200 μm along a single axis of detection. The detector was characterized for beam energy and range verification by measuring the response of the detector when irradiated with a 290 MeV/u 12 C ion broad beam incident along the single axis of the detector embedded in a PMMA phantom. The energy of the 12 C ion beam incident on the detector and the residual energy of an ion beam incident on the phantom was determined from the measured Bragg peak position in the sDMG. Ad hoc Monte Carlo simulations of the experimental setup were also performed to give further insight into the detector response., Results: The relative response profiles along the single axis measured with the sDMG detector were found to have good agreement between experiment and simulation with the position of the Bragg peak determined to fall within 0.2 mm or 1.1% of the range in the detector for the two cases. The energy of the beam incident on the detector was found to vary less than 1% between experiment and simulation. The beam energy incident on the phantom was determined to be (280.9 ± 0.8) MeV/u from the experimental and (280.9 ± 0.2) MeV/u from the simulated profiles. These values coincide with the expected energy of 281 MeV/u., Conclusions: The sDMG detector response was studied experimentally and characterized using a Monte Carlo simulation. The sDMG detector was found to accurately determine the 12 C beam energy and is suited for fast energy and range verification quality assurance. It is proposed that the sDMG is also applicable for verification of treatment planning systems that rely on particle range., (© 2017 American Association of Physicists in Medicine.)

Background: Carbon-ion radiotherapy (CIRT) is a promising treatment for locally advanced non-small-cell lung cancer, especially for patients with inoperable lung cancer. Although the incidence of CIRT-induced radiation pneumonitis (RP) ≥ grade 2 ranges from 2.5 to 9.9%, the association between CIRT-induced RP and dosimetric parameters is not clear. Herein, we identified prognostic factors associated with symptomatic RP after CIRT for patients with non-small-cell lung cancer., Methods: Clinical results of 65 patients treated with CIRT between 2000 and 2015 at the National Institute of Radiological Sciences were retrospectively analyzed. Clinical stage II B disease (TNM classification) was the most common stage among the patients (45%). The median radiation dose was 72 Gy (68-76) relative biological effectiveness (RBE) in 16 fractions. In cases involving metastatic lymph nodes, prophylactic irradiation of mediastinal lymph nodes was performed at a median dose of 49.5 Gy (RBE). The median follow-up was 22 months., Results: Grade 2 and grade 3 RP occurred in 6 and 3 patients (9 and 5%), respectively. No patients developed grade 4 or 5 RP. Using univariate analysis, vital capacity as a percentage of predicted (%VC), forced expiratory volume in 1 s (FEV1), mean lung dose (MLD), volume of lung receiving ≥5 Gy (RBE) (V 5 ), V 10 , V 20 and V 30 were determined to be the significant predictive factors for ≥ grade 2 RP. The receiver operating characteristic (ROC) analysis revealed the cutoff values for %VC, FEV1, MLD, V 5 , V 10 , V 20 and V 30 for ≥ grade 2 RP, which were 86.9%, 1.16 L, 12.5 Gy (RBE), 28.8, 29.9, 20.1 and 15.0%, respectively. In addition, the multivariate analysis revealed that %VC <86.9% (odds ratio = 13.7; p = 0.0041) and V 30  ≥ 15% (odds ratio = 6.1; p = 0.0221) were significant risk factors., Conclusions: Our study demonstrated the risk factors for ≥ grade 2 RP after carbon-ion radiotherapy for patients with locally advanced lung cancer.

Silicon microdosimetry is a promising technology for heavy ion therapy (HIT) quality assurance, because of its sub-mm spatial resolution and capability to determine radiation effects at a cellular level in a mixed radiation field. A drawback of silicon is not being tissue-equivalent, thus the need to convert the detector response obtained in silicon to tissue. This paper presents a method for converting silicon microdosimetric spectra to tissue for a therapeutic 12 C beam, based on Monte Carlo simulations. The energy deposition spectra in a 10 μm sized silicon cylindrical sensitive volume (SV) were found to be equivalent to those measured in a tissue SV, with the same shape, but with dimensions scaled by a factor κ equal to 0.57 and 0.54 for muscle and water, respectively. A low energy correction factor was determined to account for the enhanced response in silicon at low energy depositions, produced by electrons. The concept of the mean path length [Formula: see text] to calculate the lineal energy was introduced as an alternative to the mean chord length [Formula: see text] because it was found that adopting Cauchy's formula for the [Formula: see text] was not appropriate for the radiation field typical of HIT as it is very directional. [Formula: see text] can be determined based on the peak of the lineal energy distribution produced by the incident carbon beam. Furthermore it was demonstrated that the thickness of the SV along the direction of the incident 12 C ion beam can be adopted as [Formula: see text]. The tissue equivalence conversion method and [Formula: see text] were adopted to determine the RBE 10 , calculated using a modified microdosimetric kinetic model, applied to the microdosimetric spectra resulting from the simulation study. Comparison of the RBE 10 along the Bragg peak to experimental TEPC measurements at HIMAC, NIRS, showed good agreement. Such agreement demonstrates the validity of the developed tissue equivalence correction factors and of the determination of [Formula: see text].

Background and Purpose: Radiation-induced cancer is a serious late effect that may follow radiotherapy. A considerable uncertainty is associated with carcinogenesis from photon-based treatment, and even less established when including relative biological effectiveness (RBE) for particle therapy. The aim of this work was therefore to estimate and in particular explore relative risks (RR) of secondary cancer (SC) following particle therapy as applied in treatment of prostate cancer., Material and Methods: RRs of radiation-induced SC in the bladder and rectum were estimated using a bell-shaped dose-response model incorporating RBE and fractionation effects. The risks from volumetric modulated arc therapy (VMAT) were compared to intensity-modulated proton therapy (IMPT) and scanning carbon ions for ten patients., Results: The mean estimated RR (95% CI) of SC for VMAT/C-ion was 1.31 (0.65-2.18) for the bladder and 0.58 (0.41-0.80) for the rectum. Corresponding values for VMAT/IMPT were 1.72 (1.06-2.37) and 1.10 (0.78-1.43). The radio-sensitivity parameter α had the strongest influence on the results with decreasing RR for increasing values of α., Conclusion: Based on the wide spread in RR between patients and variations across the included parameter values, the risk profiles of the rectum and bladder were not dramatically different for the investigated radiotherapy techniques., (Copyright © 2016 Elsevier Ireland Ltd. All rights reserved.)

Background and Purpose: In carbon ion radiotherapy (CIRT), the use of different relative biological effectiveness (RBE) models in the RBE-weighted dose (DRBE) calculation can lead to deviations in the physical dose (Dphy) delivered to the patient. Our aim is to reduce target Dphy deviations by converting prescription dose values., Material and Methods: Planning data of patients treated at the National Institute of Radiological Sciences (NIRS) were collected, with prescribed doses per fraction ranging from 3.6Gy (RBE) to 4.6Gy (RBE), according to the Japanese semi-empirical model. The Dphy was Monte Carlo (MC) re-calculated simulating the NIRS beamline. The local effect model (LEM)&#95;I was then applied to estimate DRBE. Target median DRBE ratios between MC+LEM&#95;I and NIRS plans determined correction factors for the conversion of prescription doses. Plans were re-optimized in a LEM&#95;I-based commercial system, prescribing the NIRS uncorrected and corrected DRBE., Results: The MC+LEM&#95;I target median DRBE was respectively 15% and 5% higher than the NIRS reference, for the lowest and highest dose levels. Uncorrected DRBE prescription resulted in significantly lower target Dphy in re-optimized plans, with respect to NIRS plans., Conclusions: Prescription dose conversion factors could minimize target physical dose variations due to the use of different radiobiological models in the calculation of CIRT RBE-weighted dose., (Copyright © 2016 Elsevier Ireland Ltd. All rights reserved.)

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., (Copyright © 2015 Elsevier Ireland Ltd. All rights reserved.)

We determined the relative biological effectiveness (RBE) and oxygen enhancement ratio (OER) of micronuclei (MN) formation in clamped (hypoxic) and non-clamped (normoxic) solid tumors in mice legs following exposure to X-rays and heavy ions. Single-cell suspensions (aerobic) of non-irradiated tumors were prepared in parallel and used directly to determine the radiation response for aerobic cells. Squamous cell carcinoma (SCCVII) cells were transplanted into the right hind legs of syngeneic C3H/He male mice. Irradiation doses with either X-rays or heavy ions at a dose-averaged LET (linear energy transfer) of 14-192keV/μm were delivered to 5-mm diameter tumors and aerobic single-cells in sample-tubes. After irradiation, the tumors were excised and trypsinized to observe MN in single-cells. The single-cell suspensions were used for MN formation assays. The RBE values increased with increasing LET. The maximum RBE values for the three different oxygen conditions; hypoxic tumor, normoxic tumor, and aerobic cells, were 8.18, 5.30, and 3.76 at an LET of 192keV/μm, respectively. After X-irradiation, the OERh/n values (hypoxic tumor/normoxic tumor) were lower than the OERh/a (hypoxic tumor/aerobic cells), and were 1.73 and 2.58, respectively. We found that the OER for the in vivo studies were smaller in comparison to that for the in vitro studies. Both of the OER values at 192keV/μm were small in comparison to those of the X-ray irradiated samples. The OERh/n and OERh/a values at 192keV/μm were 1.12 and 1.19, respectively. Our results suggest that high LET radiation has a large biological effect even if a solid tumor includes substantial numbers of hypoxic cells. To conclude, we found that the RBE values under each oxygen state for non-MN fraction increased with increasing LET and that the OER values for both tumors in vivo and cells in vitro decreased with increasing LET., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Background and Purpose: Carbon-ion radiotherapy uses spread-out Bragg peaks (SOBP) to produce uniform biological effects within a target volume. The relative biological effectiveness is determined by the in vitro cell kill after a single dose is employed to design the SOBP. A question remains as to whether biological effects for in vivo tissues after fractionated doses are also uniform within the SOBP., Material and Methods: Mouse foot skin was irradiated with fractionated doses of carbon ions at various linear energy transfer (LET) values. A new ridge filter was designed based on alpha and beta values for each LET to cause moderate skin reaction, and was studied concerning its uniformity., Results: The reciprocal total doses of intermediate-LET carbon ions and of reference gamma rays linearly increased with an increase of a dose per fraction in Fe-plots. As the single total dose of higher LET run off linearity, data obtained from 2 to 6 fractions were used to design a new ridge filter. The physical dose distribution of the new ridge filter was almost identical to, and indistinguishable from, the ridge filter designed based on the in vitro cell kill., Conclusions: The LET dependence of alpha is a principle of the biological factor to be used for designing spread-out Bragg peaks of carbon-ion radiotherapy., (Copyright © 2015 Elsevier Ireland Ltd. All rights reserved.)

At the National Institute of Radiological Sciences (NIRS), more than 8,000 patients have been treated for various tumors with carbon-ion (C-ion) radiotherapy in the past 20 years based on a radiobiologically defined clinical-dose system. Through clinical experience, including extensive dose escalation studies, optimum dose-fractionation protocols have been established for respective tumors, which may be considered as the standards in C-ion radiotherapy. Although the therapeutic appropriateness of the clinical-dose system has been widely demonstrated by clinical results, the system incorporates several oversimplifications such as dose-independent relative biological effectiveness (RBE), empirical nuclear fragmentation model, and use of dose-averaged linear energy transfer to represent the spectrum of particles. We took the opportunity to update the clinical-dose system at the time we started clinical treatment with pencil beam scanning, a new beam delivery method, in 2011. The requirements for the updated system were to correct the oversimplifications made in the original system, while harmonizing with the original system to maintain the established dose-fractionation protocols. In the updated system, the radiation quality of the therapeutic C-ion beam was derived with Monte Carlo simulations, and its biological effectiveness was predicted with a theoretical model. We selected the most used C-ion beam with αr = 0.764 Gy(-1) and β = 0.0615 Gy(-2) as reference radiation for RBE. The C-equivalent biological dose distribution is designed to allow the prescribed survival of tumor cells of the human salivary gland (HSG) in entire spread-out Bragg peak (SOBP) region, with consideration to the dose dependence of the RBE. This C-equivalent biological dose distribution is scaled to a clinical dose distribution to harmonize with our clinical experiences with C-ion radiotherapy. Treatment plans were made with the original and the updated clinical-dose systems, and both physical and clinical dose distributions were compared with regard to the prescribed dose level, beam energy, and SOBP width. Both systems provided uniform clinical dose distributions within the targets consistent with the prescriptions. The mean physical doses delivered to targets by the updated system agreed with the doses by the original system within ± 1.5% for all tested conditions. The updated system reflects the physical and biological characteristics of the therapeutic C-ion beam more accurately than the original system, while at the same time allowing the continued use of the dose-fractionation protocols established with the original system at NIRS.

To understand the biological response of normal cells to fractionated carbon beam irradiation, the effects of potentially lethal damage repair (PLDR) and sublethal damage repair (SLDR) were both taken into account in a linear-quadratic (LQ) model. The model was verified by the results of a fractionated cell survival experiment with normal human fibroblast cells. Cells were irradiated with 200-kV X-rays and monoenergetic carbon ion beams (290 MeV/u) at two irradiation depths, corresponding to linear energy transfers (LETs) of approximately 13 keV/μm and 75 keV/μm, respectively, at the Heavy Ion Medical Accelerator in Chiba of the National Institute of Radiological Sciences. When we only took into account the repair factor of PLDR, γ, which was derived from the delayed assay, the cell survival response to fractionated carbon ion irradiation was not fully explained in some cases. When both the effects of SLDR and PLDR were taken into account in the LQ model, the cell survival response was well reproduced. The model analysis suggested that PLDR occurs in any type of radiation. The γ factors ranged from 0.36-0.93. In addition, SLD was perfectly repaired during the fraction interval for the lower LET irradiations but remained at about 30% for the high-LET irradiation.

The aim of this study was to measure the RBE (relative biological effectiveness) and OER (oxygen enhancement ratio) for survival of cells within implanted solid tumors following exposure to 290MeV/nucleon carbon-ion beams or X-rays. Squamous cell carcinoma cells (SCCVII) were transplanted into the right hind legs of syngeneic C3H male mice. Irradiation with either carbon-ion beams with a 6-cm spread-out Bragg peak (SOBP, at 46 and 80keV/μm) or X-rays was delivered to 5-mm or less diameter tumors. We defined three different oxygen statuses of the irradiated cells. Hypoxic and normoxic conditions in tumors were produced by clamping or not clamping the leg to avoid blood flow. Furthermore, single-cell suspensions were prepared from non-irradiated tumors and directly used to determine the radiation response of aerobic cells. Single-cell suspensions (aerobic condition) were fully air-saturated. Single-cell suspensions were prepared from excised and trypsinized tumors, and were used for in vivo-in vitro colony formation assays to obtain cell survival curves. The RBE values increased with increasing LET in SOBP beams. The maximum RBE values in three different oxygen conditions; hypoxic tumor, normoxic tumor and aerobic cells, were 2.16, 1.76 and 1.66 at an LET of 80keV/μm, respectively. After X-ray irradiation the OERh/n values (hypoxic tumor/normoxic tumor) were lower than the OERh/a (hypoxic tumor/aerobic cells), and were 1.87±0.13 and 2.52±0.11, respectively. The OER values of carbon-ion irradiated samples were small in comparison to those of X-ray irradiated samples. However, no significant changes of the OER at proximal and distal positions within the SOBP carbon-ion beams were observed. To conclude, we found that the RBE values for cell survival increased with increasing LET and that the OER values changed little with increasing LET within the SOBP carbon-ion beams., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Heavy particle irradiation produces complex DNA double strand breaks (DSBs) which can arise from primary ionisation events within the particle trajectory. Additionally, secondary electrons, termed delta-electrons, which have a range of distributions can create low linear energy transfer (LET) damage within but also distant from the track. DNA damage by delta-electrons distant from the track has not previously been carefully characterised. Using imaging with deconvolution, we show that at 8 hours after exposure to Fe (∼200 keV/µm) ions, γH2AX foci forming at DSBs within the particle track are large and encompass multiple smaller and closely localised foci, which we designate as clustered γH2AX foci. These foci are repaired with slow kinetics by DNA non-homologous end-joining (NHEJ) in G1 phase with the magnitude of complexity diminishing with time. These clustered foci (containing 10 or more individual foci) represent a signature of DSBs caused by high LET heavy particle radiation. We also identified simple γH2AX foci distant from the track, which resemble those arising after X-ray exposure, which we attribute to low LET delta-electron induced DSBs. They are rapidly repaired by NHEJ. Clustered γH2AX foci induced by heavy particle radiation cause prolonged checkpoint arrest compared to simple γH2AX foci following X-irradiation. However, mitotic entry was observed when ∼10 clustered foci remain. Thus, cells can progress into mitosis with multiple clusters of DSBs following the traversal of a heavy particle.

The authors attempt to establish the relative biological effectiveness (RBE) calculation for designing therapeutic proton beams on the basis of microdosimetry. The tissue-equivalent proportional counter (TEPC) was used to measure microdosimetric lineal energy spectra for proton beams at various depths in a water phantom. An RBE-weighted absorbed dose is defined as an absorbed dose multiplied by an RBE for cell death of human salivary gland (HSG) tumor cells in this study. The RBE values were calculated by a modified microdosimetric kinetic model using the biological parameters for HSG tumor cells. The calculated RBE distributions showed a gradual increase to about 1cm short of a beam range and a steep increase around the beam range for both the mono-energetic and spread-out Bragg peak (SOBP) proton beams. The calculated RBE values were partially compared with a biological experiment in which the HSG tumor cells were irradiated by the SOBP beam except around the distal end. The RBE-weighted absorbed dose distribution for the SOBP beam was derived from the measured spectra for the mono-energetic beam by a mixing calculation, and it was confirmed that it agreed well with that directly derived from the microdosimetric spectra measured in the SOBP beam. The absorbed dose distributions to planarize the RBE-weighted absorbed dose were calculated in consideration of the RBE dependence on the prescribed absorbed dose and cellular radio-sensitivity. The results show that the microdosimetric measurement for the mono-energetic proton beam is also useful for designing RBE-weighted absorbed dose distributions for range-modulated proton beams.

We investigated the applicability of the repairable-conditionally repairable (RCR) model and the multi-target (MT) model to dose conversion in high-dose-per-fraction radiotherapy in comparison with the linear-quadratic (LQ) model. Cell survival data of V79 and EMT6 single cells receiving single doses of 2-12 Gy or 2 or 3 fractions of 4 or 5 Gy each, and that of V79 spheroids receiving single doses of 5-26 Gy or 2-5 fractions of 5-12 Gy, were analyzed. Single and fractionated doses to actually reduce cell survival to the same level were determined by a colony assay. Single doses used in the experiments and surviving fractions at the doses were substituted into equations of the RCR, MT and LQ models in the calculation software Mathematica, and each parameter coefficient was computed. Thereafter, using the coefficients and the three models, equivalent single doses for the hypofractionated doses were calculated. They were then compared with actually-determined equivalent single doses for the hypofractionated doses. The equivalent single doses calculated using the RCR, MT and LQ models tended to be lower than the actually determined equivalent single doses. The LQ model seemed to fit relatively well at doses of 5 Gy or less. At 6 Gy or higher doses, the RCR and MT models seemed to be more reliable than the LQ model. In hypofractionated stereotactic radiotherapy, the LQ model should not be used, and conversion models incorporating the concept of the RCR or MT models, such as the generalized linear-quadratic models, appear to be more suitable.

In carbon ion radiotherapy there is an urgent clinical need to develop objective tools for the conversion of relative biological effectiveness (RBE)-weighted doses based on different models. In this work we introduce a clinically oriented method to compare NIRS-based and LEM-based GyE systems, minimizing differences in physical dose distributions between treatment plans. Carbon ion plans were optimized on target volumes of cubic and spherical shapes, for RBE-weighted dose prescription levels ranging from 3.6 to 4.4 GyE. Plans were calculated for target sizes from 4 to 12 cm defining three beam geometries: single beam, opposed beam and orthogonal beam configurations. The two treatment planning systems currently employed in clinical practice were used, providing the NIRS-based and LEM-based GyE calculations. Physical dose distributions of NIRS-based and LEM-based treatment plans were compared. LEM-based prescription doses that minimize differences in physical dose distributions between the two systems were found. These doses were compared with the mean RBE-weighted dose obtained with a Monte Carlo code (FLUKA) interfaced with LEM I. In the investigated dose range, LEM-based RBE-weighted prescription doses, that minimize differences with NIRS plans, should be higher than NIRS reported prescription doses. The optimal dose depends on target size, shape and position, number of beams and dose level. The opposed beam configuration resulted in the smallest average prescription dose difference (0.45 ± 0.09 GyE). The second approach of recalculating NIRS RBE-weighted dose with a Monte Carlo code interfaced with LEM resulted in no significant difference with the results obtained from the planning study. The delivery of a voxel by voxel iso-effective plan, if different RBE models are employed, is not feasible; it is however possible to minimize differences in a treatment plan with the simple approach presented here. Dose prescription ultimately represents a clinical task under the responsibility of the radiation oncologist, the presented analysis intends to be a quantitative and objective way to assist the clinical decision.

Background and Purpose: Characterization of the out-of-field dose profile following irradiation of the target with a 3D treatment plan delivered with modern techniques., Methods: An anthropomorphic RANDO phantom was irradiated with a treatment plan designed for a simulated 5 × 2 × 5 cm(3) tumor volume located in the center of the head. The experiment was repeated with all most common radiation treatment types (photons, protons and carbon ions) and delivery techniques (Intensity Modulated Radiation Therapy, passive modulation and spot scanning). The measurements were performed with active diamond detector and passive thermoluminescence (TLD) detectors to investigate the out-of-field dose both inside and outside the phantom., Results: The highest out-of-field dose values both on the surface and inside the phantom were measured during the treatment with 25 MV photons. In the proximity of the Planned Target Volume (PTV), the lowest lateral dose profile was observed for passively modulated protons mainly because of the presence of the collimator in combination with the chosen volume shape. In the far out-of-field region (above 100mm from the PTV), passively modulated ions were characterized by a less pronounced dose fall-off in comparison with scanned beams. Overall, the treatment with scanned carbon ions delivered the lowest dose outside the target volume., Conclusions: For the selected PTV, the use of the collimator in proton therapy drastically reduced the dose deposited by ions or photons nearby the tumor. Scanning modulation represents the optimal technique for achieving the highest dose reduction far-out-of-field., (Copyright © 2012 Elsevier Ireland Ltd. All rights reserved.)

Purpose: Recent radiotherapy technologies including carbon-ion radiotherapy can improve the dose concentration in the target volume, thereby not only reducing side effects in organs at risk but also the secondary cancer risk within or near the irradiation field. However, secondary cancer risk in the low-dose region is considered to be non-negligible, especially for younger patients. To achieve a dose estimation of the whole body of each patient receiving carbon-ion radiotherapy, which is essential for risk assessment and epidemiological studies, Monte Carlo simulation plays an important role because the treatment planning system can provide dose distribution only in∕near the irradiation field and the measured data are limited. However, validation of Monte Carlo simulations is necessary. The primary purpose of this study was to establish a calculation method using the Monte Carlo code to estimate the dose and quality factor in the body and to validate the proposed method by comparison with experimental data. Furthermore, we show the distributions of dose equivalent in a phantom and identify the partial contribution of each radiation type., Methods: We proposed a calculation method based on a Monte Carlo simulation using the PHITS code to estimate absorbed dose, dose equivalent, and dose-averaged quality factor by using the Q(L)-L relationship based on the ICRP 60 recommendation. The values obtained by this method in modeling the passive beam line at the Heavy-Ion Medical Accelerator in Chiba were compared with our previously measured data., Results: It was shown that our calculation model can estimate the measured value within a factor of 2, which included not only the uncertainty of this calculation method but also those regarding the assumptions of the geometrical modeling and the PHITS code. Also, we showed the differences in the doses and the partial contributions of each radiation type between passive and active carbon-ion beams using this calculation method. These results indicated that it is essentially important to include the dose by secondary neutrons in the assessment of the secondary cancer risk of patients receiving carbon-ion radiotherapy with active as well as passive beams., Conclusions: We established a calculation method with a Monte Carlo simulation to estimate the distribution of dose equivalent in the body as a first step toward routine risk assessment and an epidemiological study of carbon-ion radiotherapy at NIRS. This method has the advantage of being verifiable by the measurement.

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., (© 2012 Institute of Physics and Engineering in Medicine)

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.

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.

We describe a method to calculate the relative biological effectiveness in mixed radiation fields of therapeutic ion beams based on the modified microdosimetric kinetic model (modified MKM). In addition, we show the procedure for integrating the modified MKM into a treatment planning system for a scanned carbon beam. With this procedure, the model is fully integrated into our research version of the treatment planning system. To account for the change in radiosensitivity of a cell line, we measured one of the three MKM parameters from a single survival curve of the current cells and used the parameter in biological optimization. Irradiation of human salivary gland tumor cells was performed with a scanned carbon beam in the Heavy Ion Medical Accelerator in Chiba (HIMAC), and we then compared the measured depth-survival curve with the modified MKM predicted survival curve. Good agreement between the two curves proves that the proposed method is a candidate for calculating the biological effects in treatment planning for ion irradiation.

Purpose: Successful results in carbon-ion and proton radiotherapies can extend patients' lives and thus present a treatment option for younger patients; however, the undesired exposure to normal tissues outside the treatment volume is a concern. Organ-specific information on the absorbed dose and the biological effectiveness in the patient is essential for assessing the risk, but experimental dose assessment has seldom been done. In this study, absorbed doses, quality factors, and dose equivalents in water phantom outside of the irradiation field were determined based on lineal energy distributions measured with a commercial tissue equivalent proportional counter (TEPC) at passive carbon-ion and proton radiotherapy facilities., Methods: Measurements at eight positions in the water phantom were carried out at the Heavy-Ion Medical Accelerator in Chiba of the National Institute of Radiological Sciences for 400 and 290 MeV/u carbon beams and at the National Cancer Center Hospital East for a 235 MeV proton beam., Results: The dose equivalent per treatment absorbed dose at the center of the range-modulated region H/Dt, decreased as the position became farther from the beam axis and farther from the phantom surface. The values of H/Dt ranged from 6.7 to 0.16 mSv/Gy for the 400 MeV/u carbon beam, from 1.3 to 0.055 mSv/Gy for the 290 MeV/u carbon beam, and from 4.7 to 0.24 mSv/GV for the 235 MeV proton beam. The values of the dose-averaged quality factor QD ranged from 2.4 to 4.6 for the 400 MeV/u beam, from 2.8 to 5.3 for the 290 MeV/u beam, and from 5.1 to 8.2 for the proton beam. The authors also observed differences in the distributions of H/Dt and QD between the carbon and proton beams., Conclusions: The authors experimentally obtained absorbed doses, dose-averaged quality factors, and dose equivalents in water phantom outside of the irradiation field in passive carbon-ion and proton radiotherapies with TEPC. These data are very useful for estimating the risk of secondary cancer after receiving passive radiotherapies and for verifying Monte Carlo calculations.

In the last few years, hospital-based facilities for carbon-ion radiotherapy are being constructed and proposed in Europe and Asia. During the next few years, several new facilities will be opened for carbon-ion radiotherapy in the world. These facilities in operation or under construction are categorized in two types by the beam shaping method used. One is the passive beam shaping method that is mainly improved and systematized for routine clinical use at HIMAC, Japan. The other method is active beam shaping which is also known as beam scanning adopted at GSI/HIT, Germany. In this paper an overview of some technical aspects for beam shaping is reported. The technique of passive beam shaping is established for stable clinical application and has clinical result of over 4000 patients in HIMAC. In contrast, clinical experience of active beam shaping is about 400 patients, and there is no clinical experience to respiratory moving target. A great advantage of the active beam shaping method is patient-specific collimator-less and compensator-less treatment. This may be an interesting potential for adaptive radiotherapy.

Purpose: Recent successful results in passive carbon-ion radiotherapy allow the patient to live for a longer time and allow younger patients to receive the radiotherapy. Undesired radiation exposure in normal tissues far from the target volume is considerably lower than that close to the treatment target, but it is considered to be non-negligible in the estimation of the secondary cancer risk. Therefore, it is very important to reduce the undesired secondary neutron exposure in passive carbon-ion radiotherapy without influencing the clinical beam. In this study, the source components in which the secondary neutrons are produced during passive carbon-ion radiotherapy were identified and the method to reduce the secondary neutron dose effectively based on the identification of the main sources without influencing the clinical beam was investigated., Methods: A Monte Carlo study with the PHITS code was performed by assuming the beamline at the Heavy-Ion Medical Accelerator in Chiba (HIMAC). At first, the authors investigated the main sources of secondary neutrons in passive carbon-ion radiotherapy. Next, they investigated the reduction in the neutron dose with various modifications of the beamline device that is the most dominant in the neutron production. Finally, they investigated the use of an additional shield for the patient., Results: It was shown that the main source is the secondary neutrons produced in the four-leaf collimator (FLC) used as a precollimator at HIAMC, of which contribution in the total neutron ambient dose equivalent is more than 70%. The investigations showed that the modification of the FLC can reduce the neutron dose at positions close to the beam axis by 70% and the FLC is very useful not only for the collimation of the primary beam but also the reduction in the secondary neutrons. Also, an additional shield for the patient is very effective to reduce the neutron dose at positions farther than 50 cm from the beam axis. Finally, they showed that the neutron dose can be reduced by approximately 70% at any position without influencing the primary beam used in treatment., Conclusions: This study was performed by assuming the HIMAC beamline; however, this study provides important information for reoptimizing the arrangement and the materials of beamline devices and designing a new facility for passive carbon-ion radiotherapy and probably passive proton radiotherapy.

Particle radiotherapy using proton and heavy ion beams has shown improved clinical results and is a promising cancer therapy which is expected to gradually spread in Japan. There are, however, no special regulations for radiotherapy treatment facilities. They have been operated under the same safety regulations as for a research facility using a research accelerator. Significantly high-energy radiation is necessary for particle radiotherapy compared with conventional radiation therapy. The treatment facility, therefore, should have a large accelerator, which is installed in a room with a thick shield wall. Data on radiation protection for such high energy medical facilities is fragmentary and insufficient. In this study, we examined the necessity of other regulations for the safe operation of medical facilities for particle radiotherapy. First, we measured activation levels of the therapeutic devices and of patients. Next the safety level of the medical facility was evaluated from the viewpoint of radiation protection. We have confirmed the facilities can be safely operated by present regulations given in the Law Concerning Prevention from Radiation Hazards due to Radiation Isotopes, etc. or the Law for Health Protection and Medical Care.

We developed a calculation technique to evaluate the biological dose distribution of heavy ion beams, and verified its reliability by comparison with experimental results. The calculation technique was developed by connecting two general-purpose Monte Carlo codes. In order to evaluate the radiation quality and biological effect, the microdosimetric kinetic model was adopted. We estimated the distribution of physical dose and radiation quality for the carbon broad beam with experiments and with Monte Carlo calculations. Relative biological effectiveness (RBE) was estimated from radiation quality, and biological dose could be calculated as the product of the physical dose and the RBE. Our calculations showed good agreement with experiments, not only for the physical dose, but also for the dose-averaged lineal energy as an expression of radiation quality and therefore biological dose. This finding indicates that our calculation tool will be useful to estimate biological dose distribution in the design of heavy ion radiation facilities or for quality assurance in treatment planning.

Secondary neutron ambient dose equivalents per the treatment absorbed dose in passive carbon-ion and proton radiotherapies were measured using a rem meter, WENDI-II at two carbon-ion radiotherapy facilities and four proton radiotherapy facilities in Japan. Our measured results showed that (1) neutron ambient dose equivalent in carbon-ion radiotherapy is lower than that in proton radiotherapy, and (2) the difference to the measured neutron ambient dose equivalents among the facilities is within a factor of 3 depending on the operational beam setting used at the facility and the arrangement of the beam line, regardless of the method for making a laterally uniform irradiation field: the double scattering method or the single-ring wobbling method. The reoptimization of the beam line in passive particle radiotherapy is an effective way to reduce the risk of secondary cancer because installing an adjustable precollimator and designing the beam line devices with consideration of their material, thickness and location, etc., can significantly reduce the neutron exposure. It was also found that the neutron ambient dose equivalent in passive particle radiotherapy is equal to or less than that in the photon radiotherapy. This result means that not only scanning particle radiotherapy but also passive particle radiotherapy can provide reduced exposure to normal tissues around the target volume without an accompanied increase in total body dose.

We investigated the earliest possible chromosome break and repair process in normal human fibroblasts irradiated with low and high LET (linear energy transfer) heavy ion radiation using the modified premature chromosome condensation (PCC) technique utilizing wortmannin (WM) during the fusion incubation period [M. Okada, S. Saito, R. Okayasu, Facilitated detection of chromosome break and repair at low levels of ionizing radiation by addition of wortmannin to G1-type PCC fusion incubation, Mutat. Res., 562 (2004) 11-17]. The initial numbers of breaks were approximately 10/cell/Gy in X-irradiated samples, followed by carbon (LET: 70 keV/microm), neon, and the number was around 5/cell/Gy in silicon (LET: 70 and 200 keV/microm) and iron (LET: 200 keV/microm) samples. If WM was not used, the initial numbers of breaks with silicon and iron were higher than those of X-rays. To quantify these data, we used initial repair ratio (IRR) defined as the number of G1 PCC breaks with WM divided by the number of breaks without WM. X-irradiation gave the maximum IRR ( approximately 2.0), while iron as well as silicon irradiation showed the minimum IRR ( approximately 1.0), suggesting almost no rejoining at the initial stage. Although there is a comparatively good correlation between the IRR value and the cell survival, the survival fraction with the repair data at 2 or 6h correlates better statistically. Our data indicate that high LET heavy ion irradiation induces a lower number of initial chromosome breaks with minimal repair when compared with low LET irradiation. These results at the chromosome level substantiate and extend the notion that high LET radiation produces complex-type DNA double strand breaks (DSBs).

Both the microdosimetric kinetic model (MKM) and the local effect model (LEM) can be used to calculate the surviving fraction of cells irradiated by high-energy ion beams. In this study, amorphous track structure models instead of the stochastic energy deposition are used for the MKM calculation, and it is found that the MKM calculation is useful for predicting the survival curves of the mammalian cells in vitro for (3)He-, (12)C- and (20)Ne-ion beams. The survival curves are also calculated by two different implementations of the LEM, which inherently used an amorphous track structure model. The results calculated in this manner show good agreement with the experimental results especially for the modified LEM. These results are compared to those calculated by the MKM. Comparison of the two models reveals that both models require three basic constituents: target geometry, photon survival curve and track structure, although the implementation of each model is significantly different. In the context of the amorphous track structure model, the difference between the MKM and LEM is primarily the result of different approaches calculating the biological effects of the extremely high local dose in the center of the ion track.

The aim of this work is to define a simplified semi-analytical beam transportation code that can calculate the spatial distribution of projectile fragments which are widely distributed in a patient's body during heavy-ion beam radiotherapy. In this code, we employed an elemental pencil beam model where the spatial distribution of radiation quality for an elemental beam is calculated and superposed according to the emittance ellipse of the narrow heavy-ion beam determined at the entrance of the target. The radiation quality for an elemental beam was calculated using Goldhaber's model of fragment distribution. The calculation results were compared with the experimental observations for a mono-energetic narrow (12)C beam measured at the secondary beam line in HIMAC. Despite its simplicity, the developed code could reproduce the experimental results well.

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.

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.

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.

Purpose: A retrospective analysis was made to examine appropriateness in the estimation of the biologic effectiveness of carbon-ion radiotherapy using resultant data from clinical trials at the heavy-ion medical accelerator complex (HIMAC) at the National Institute of Radiological Sciences in Chiba, Japan., Methods and Materials: At HIMAC, relative biologic effectiveness (RBE) values of therapeutic carbon beams were determined based on experimental results of cell responses, on values expected with the linear-quadratic model, and based on experiences with neutron therapy. We use fixed RBE values independent of dose levels, although this apparently contradicts radiobiologic observations. Our RBE system depends only on LET of the heavy-ion radiation fields. With this RBE system, over 2,000 patients have been treated by carbon beams. With data from these patients, the local control rate of non-small-cell lung cancer was analyzed to verify the clinical RBE of the carbon beam. The local control rate was compared with rates published by groups from Gunma University and Massachusetts General Hospital. Using a simplified tumor control probability (TCP) model, clinical RBE values were obtained for different levels of TCP., Results: For the 50% level of the clinical TCP, the RBE values nearly coincide with those for in vitro human salivary gland cell survival at 10%. For the higher levels of clinical TCP, the RBE values approach closer to those adapted in clinical trials at HIMAC.

The latest heavy ion therapy tends to require information about the spatial distribution of the quality of radiation in a patient's body in order to make the best use of any potential advantage of swift heavy ions for the therapeutic treatment of a tumour. The deflection of incident particles is described well by Molière's multiple-scattering theory of primary particles; however, the deflection of projectile fragments is not yet thoroughly understood. This paper reports on our investigation of the spatial distribution of fragments produced from a therapeutic carbon beam through nuclear reactions in thick water. A DeltaE-E counter telescope system, composed of a plastic scintillator, a gas-flow proportional counter and a BGO scintillator, was rotated around a water target in order to measure the spatial distribution of the radiation quality. The results revealed that the observed deflection of fragment particles exceeded the multiple scattering effect estimated by Molière's theory. However, the difference can be sufficiently accounted for by considering one term involved in the multiple-scattering formula; this term corresponds to a lateral 'kick' at the point of production of the fragment. This kick is successfully explained as a transfer of the intra-nucleus Fermi momentum of a projectile to the fragment; the extent of the kick obeys the expectation derived from the Goldhaber model.

The production of projectile fragments is one of the most important, but not yet perfectly understood, problems to be considered when planning for the utilization of high-energy heavy charged particles for radiotherapy. This paper reports our investigation of the fragments' fluence and linear energy transfer (LET) spectra produced from various incident ions using an experimental approach to reveal these physical qualities of the beams. Polymethyl methacrylate, as a substitute for the human body, was used as a target. A deltaE-E counter telescope with a plastic scintillator and a BGO scintillator made it possible to identify the species of fragments based on differences of various elements. By combining a gas-flow proportional counter with a counter telescope system, LET spectra as well as fluence spectra of the fragments were derived for each element down from the primary particles to hydrogen. Among them, the information on hydrogen and helium fragments was derived for the first time. The result revealed that the number of light fragments, such as hydrogen and helium, became larger than the number of primaries in the vicinity of the range end. However, the greater part of the dose delivered to a cell was still governed by the primaries. The calculated result of a simulation used for heavy-ion radiotherapy indicated room for improving the reaction model.

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.

Single event spectra of a clinical carbon beam have been measured by an ultra-miniature tissue-equivalent proportional counter (UMC). In order to cover the energy range of the Bragg peak, the incident energy of the carbon beam was degraded by aluminium plates. Single event spectra for carbon-events incident to the UMC were analysed and selected at several carbon energies using thin scintillation counters. It was found that the dose weighted lineal energy distributions have a doublet peak structure due to incident carbon beam and fragment contributions.

In heavy-ion radiotherapy, an accelerated beam is modified to realize a desired dose distribution in patients. The setup of the beam-modifying devices in the irradiation system is changed according to the patient, and it is important to check the depth dose distributions in the patient. In order to measure dose distributions realized by an irradiation system for heavy-ion radiotherapy, a multi-layer ionization chamber(MLIC) was developed. The MLIC consists of 64 ionization chambers, which are stacked mutually. The interval between each ionization chamber is about 4.1 mm water. There are signal and high voltage plates in the MILC, which are used as electrodes of the ionization chambers and phantom. Depth dose distribution from 5.09 mm to 261.92 mm water can be measured in about 30 seconds using this MLIC. Thus, it is possible to check beam quality in a short amount of time.

The feasibility of heavy ion computed tomography was investigated. A cylindrical polyethylene sample of 50 mm in diameter was used to establish the pencil beam scanning method and estimate the performance of the system. The sample was irradiated with carbon ions of 290 MeV/u and 400 MeV/u, and helium ions of 150 MeV/u. The residual beam energy was measured by a BGO scintillator in coincidence with two small plastic scintillators placed in front of and behind the sample to restrict the beam path and improve the spatial resolution in the CT reconstruction. The projection data were obtained by moving the sample only in the transverse direction because of its symmetrical structure. A Lucite sample with many holes of different diameters was used to demonstrate the spatial resolution of this CT system. From the reconstructed images the spatial resolution was estimated to be less than 2 mm under the conditions in this work. The electron density ratios of ethyl alcohol, water, and Lucite to polyethylene were obtained, which are in good agreement with the calculated values.

For exact determination of absorbed dose in heavy-ion irradiation fields which are used in radiation therapy and biological experiments, ionization chambers have been characterized with defined heavy-ion beams and correction factors. The LET (linear energy transfer) dependence of columnar recombination in a parallel-plate ionization chamber has been examined. Using 135 MeV/u carbon and neon beams, the ion collection efficiency was measured for several gases (air, carbon dioxide, argon and tissue-equivalent gas). 95 MeV/u argon beams and 90 MeV/u iron beams were also used for measurements of columnar recombination in air. As expected by Jaffe theory, the inverse of the ratio of the ionization charge to the saturated ionization charge had a linear relationship with the inverse of the electric field strength in the region below 0.002 V(-1) cm. The gradient of the line increases as the LET of the heavy ions increases. A strong LET dependence of the gradient was observed in air and carbon dioxide. The LET dependence was not observed in tissue-equivalent gas, nitrogen or argon. The exact depth-dose distribution of the heavy-ion beam was obtained by this correction of the initial recombination effect for the collected ionization charge. The columnar recombination in air was analysed using Jaffe theory; the obtained parameter b (a track radius) should be in the range between 0.001 cm and 0.005 cm, whereas the value obtained by Jaffe is 0.00179 cm. The value of the parameter b should increase as the LET of the heavy-ion beam increases in order to reproduce the experimental values of the initial recombination.

The precise conversion of CT numbers to their electron densities is essential in treatment planning for hadron therapy. Although some conversion methods have already been proposed, it is hard to check the conversion accuracy during practical therapy. We have estimated the CT numbers of real tissues by a calculational method established by Mustafa and Jackson. The relationship between the CT numbers and the electron densities was investigated for various body tissues as well as some tissue-equivalent materials used for a conversion to check the accuracy of the current conversion methods. The result indicates a slight disagreement at the high-CT-number region. A precise estimation of the multiple scattering, nuclear reaction and range straggling of incident particles has been considered as being important to realize higher-level conformal therapy in the future. The relationship between these parameters and the CT numbers was also investigated for tissues and water. The result shows that it is sufficiently practical to replace these parameters for real tissues with those for water by adjusting the density.