Shunsuke Yonai

Measurement of neutron ambient dose equivalent in passive carbon-ion and proton radiotherapies.

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.

The Fukushima Health Management Survey: estimation of external doses to residents in Fukushima Prefecture

The Fukushima Health Management Survey (including the Basic Survey for external dose estimation and four detailed surveys) was launched after the Fukushima Dai-ichi Nuclear Power Plant accident. The Basic Survey consists of a questionnaire that asks Fukushima Prefecture residents about their behavior in the first four months after the accident; and responses to the questionnaire have been returned from many residents. The individual external doses are estimated by using digitized behavior data and a computer program that included daily gamma ray dose rate maps drawn after the accident. The individual external doses of 421,394 residents for the first four months (excluding radiation workers) had a distribution as follows: 62.0%, <1 mSv; 94.0%, <2 mSv; 99.4%, <3 mSv. The arithmetic mean and maximum for the individual external doses were 0.8 and 25 mSv, respectively. While most dose estimation studies were based on typical scenarios of evacuation and time spent inside/outside, the Basic Survey estimated doses considering individually different personal behaviors. Thus, doses for some individuals who did not follow typical scenarios could be revealed. Even considering such extreme cases, the estimated external doses were generally low and no discernible increased incidence of radiation-related health effects is expected.

NIRS external dose estimation system for Fukushima residents after the Fukushima Dai-ichi NPP accident

The great east Japan earthquake and subsequent tsunamis caused Fukushima Dai-ichi Nuclear Power Plant (NPP) accident. National Institute of Radiological Sciences (NIRS) developed the external dose estimation system for Fukushima residents. The system is being used in the Fukushima health management survey. The doses can be obtained by superimposing the behavior data of the residents on the dose rate maps. For grasping the doses, 18 evacuation patterns of the residents were assumed by considering the actual evacuation information before using the survey data. The doses of the residents from the deliberate evacuation area were relatively higher than those from the area within 20 km radius. The estimated doses varied from around 1 to 6 mSv for the residents evacuated from the representative places in the deliberate evacuation area. The maximum dose in 18 evacuation patterns was estimated to be 19 mSv.

The Fukushima Nuclear Power Plant accident and exposures in the environment

The Great East Japan Earthquake has occurred on March 11, 2011, in the Tohoku District of Japan. Due to the earthquake, big tsunamis were induced, and they rushed to the Fukushima Nuclear Power Stations, causing severe accidents. Radioactive materials including I-131, Cs-137 and so on were emitted from the plant to the environment. The Japanese government, Fukushima prefectural government and other local governments have struggled against the accidents. The restricted area and deliberate evacuation area are set by the government, and the residents are evacuated. The dose rates in and around Fukushima Prefecture have been monitored by the governments and other involved organizations. Fukushima government has started the health management survey for all residents in Fukushima Prefecture including the questions on their activities for the estimations of their external doses.

Extended collimator model for pencil-beam dose calculation in proton radiotherapy

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.

Evaluation of beam wobbling methods for heavy-ion radiotherapy

The National Institute of Radiological Sciences (NIRS) has extensively studied carbon-ion radiotherapy at the Heavy-Ion Medical Accelerator in Chiba (HIMAC) with some positive outcomes, and has established its efficacy. Therefore, efforts to distribute the therapy to the general public should be made, for which it is essential to enable direct application of clinical and technological experiences obtained at NIRS. For widespread use, it is very important to reduce the cost through facility downsizing with minimal acceleration energy to deliver the HIMAC-equivalent clinical beams. For the beam delivery system, the requirement of miniaturization is translated to reduction in length while maintaining the clinically available field size and penetration range for range-modulated uniform broad beams of regular fields that are either circular or square for simplicity. In this paper, we evaluate the various wobbling methods including original improvements, especially for application to the compact facilities through the experimental and computational studies. The single-ring wobbling method used at HIMAC is the best one including a lot of experience at HIMAC but the residual range is a fatal problem in the case of a compact facility. On the other hand, uniform wobbling methods such as the spiral and zigzag wobbling methods are effective and suitable for a compact facility. Furthermore, these methods can be applied for treatment with passive range modulation including respiratory gated irradiation. In theory, the choice between the spiral and zigzag wobbling methods depends on the shape of the required irradiation field. However, we found that it is better to use the zigzag wobbling method with transformation of the wobbling pattern even when a circular uniform irradiation field is required, because it is difficult to maintain the stability of the wobbler magnet due to the rapid change of the wobbler current in the spiral wobbling method. The regulated wobbling method, which is our improvement, can well expand the uniform irradiation field and lead to reducing the power requirement of the wobbler magnets. Our evaluations showed that the regulated zigzag wobbling method is the most suitable method for use in currently designed compact carbon-therapy facilities.

Dynamic splitting of Gaussian pencil beams in heterogeneity-correction algorithms for radiotherapy with heavy charged particles

The pencil-beam algorithm is valid only when elementary Gaussian beams are small enough compared to the lateral heterogeneity of a medium, which is not always true in actual radiotherapy with protons and ions. This work addresses a solution for the problem. We found approximate self-similarity of Gaussian distributions, with which Gaussian beams can split into narrower and deflecting daughter beams when their sizes have overreached lateral heterogeneity in the beam-transport calculation. The effectiveness was assessed in a carbon-ion beam experiment in the presence of steep range compensation, where the splitting calculation reproduced a detour effect amounting to about 10% in dose or as large as the lateral particle disequilibrium effect. The efficiency was analyzed in calculations for carbon-ion and proton radiations with a heterogeneous phantom model, where the beam splitting increased computing times by factors of 4.7 and 3.2. The present method generally improves the accuracy of the pencil-beam algorithm without severe inefficiency. It will therefore be useful for treatment planning and potentially other demanding applications.

Field-size dependence of doses of therapeutic carbon beams.

To estimate the physical dose at the center of spread-out Bragg peaks (SOBP) for various conditions of the irradiation system, a semiempirical approach was applied. The dose at the center of the SOBP depends on the field size because of large-angle scattering particles in the water phantom. For a small field of 5 x 5 cm2, the dose was reduced to 99.2%, 97.5%, and 96.5% of the dose used for the open field in the case of 290, 350, and 400 MeV/n carbon beams, respectively. Based on the three-Gaussian form of the lateral dose distributions of the carbon pencil beam, which has previously been shown to be effective for describing scattered carbon beams, we reconstructed the dose distributions of the SOBP beam. The reconstructed lateral dose distribution reproduced the measured lateral dose distributions very well. The field-size dependencies calculated using the reconstructed lateral dose distribution of the therapeutic carbon beam agreed with the measured dose dependency very well. The reconstructed beam was also used for irregularly shaped fields. The resultant dose distribution agreed with the measured dose distribution. The reconstructed beams were found to be applicable to the treatment-planning system.

Computational modeling of beam-customization devices for heavy-charged-particle radiotherapy

A model for beam customization with collimators and a range-compensating filter based on the phase-space theory for beam transport is presented for dose distribution calculation in the treatment planning of radiotherapy with protons and heavier ions. Independent handling of pencil beams in conventional pencil-beam algorithms causes unphysical collimator-height dependence in the middle of large fields, which is resolved by the framework comprised of generation, transport, collimation, regeneration, range-compensation and edge-sharpening processes with a matrix of pencil beams. The model was verified to be consistent with measurement and analytic estimation at a submillimeter level in the penumbra of individual collimators with a combinational-collimated carbon-ion beam. The model computation is fast, accurate and readily applicable to pencil-beam algorithms in treatment planning with the capability of combinational collimation to make the best use of the beam-customization devices.

The grid-dose-spreading algorithm for dose distribution calculation in heavy charged particle radiotherapy

A new variant of the pencil-beam (PB) algorithm for dose distribution calculation for radiotherapy with protons and heavier ions, the grid-dose spreading (GDS) algorithm, is proposed. The GDS algorithm is intrinsically faster than conventional PB algorithms due to approximations in convolution integral, where physical calculations are decoupled from simple grid-to-grid energy transfer. It was effortlessly implemented to a carbon-ion radiotherapy treatment planning system to enable realistic beam blurring in the field, which was absent with the broad-beam (BB) algorithm. For a typical prostate treatment, the slowing factor of the GDS algorithm relative to the BB algorithm was 1.4, which is a great improvement over the conventional PB algorithms with a typical slowing factor of several tens. The GDS algorithm is mathematically equivalent to the PB algorithm for horizontal and vertical coplanar beams commonly used in carbon-ion radiotherapy while dose deformation within the size of the pristine spread occurs for angled beams, which was within 3 mm for a single 150-MeV proton pencil beam of 30 degrees incidence, and needs to be assessed against the clinical requirements and tolerances in practical situations.