Search

Microdosimetry is a particularly powerful method to estimate the relative biological effectiveness (RBE) of any mixed radiation field. This is particularly convenient for therapeutic heavy ion therapy (HIT) beams, referring to ions larger than protons, where the RBE of the beam can vary significantly along the Bragg curve. Additionally, due to the sharp dose gradients at the end of the Bragg peak (BP), or spread out BP, to make accurate measurements and estimations of the biological properties of a beam a high spatial resolution is required, less than a millimetre. This requirement makes silicon microdosimetry particularly attractive due to the thicknesses of the sensitive volumes commonly being  ∼10 [Formula: see text]m or less. Monte Carlo (MC) codes are widely used to study the complex mixed HIT radiation field as well as to model the response of novel microdosimeter detectors when irradiated with HIT beams. Therefore it is essential to validate MC codes against experimental measurements. This work compares measurements performed with a silicon microdosimeter in mono-energetic [Formula: see text], [Formula: see text] and [Formula: see text] ion beams of therapeutic energies, against simulation results calculated with the Geant4 toolkit. Experimental and simulation results were compared in terms of microdosimetric spectra (dose lineal energy, [Formula: see text]), the dose mean lineal energy, y  D and the RBE 10 , as estimated by the microdosimetric kinetic model (MKM). Overall Geant4 showed reasonable agreement with experimental measurements. Before the distal edge of the BP, simulation and experiment agreed within  ∼10% for y  D and  ∼2% for RBE 10 . Downstream of the BP less agreement was observed between simulation and experiment, particularly for the [Formula: see text] and [Formula: see text] beams. Simulation results downstream of the BP had lower values of y  D and RBE 10 compared to the experiment due to a higher contribution from lighter fragments compared to heavier fragments.

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].

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.

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.

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)

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.

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.

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.

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.

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.