We describe a novel code system, TRiP, dedicated to the planning of radiotherapy with energetic ions, in particular 12C. The software is designed to cooperate with three-dimensional active dose shaping devices like the GSI raster scan system. This unique beam delivery system allows us to select any combination from a list of 253 individual beam energies, 7 different beam spot sizes and 15 intensity levels. The software includes a beam model adapted to and verified for carbon ions. Inverse planning techniques are implemented in order to obtain a uniform target dose distribution from clinical input data, i.e. CT images and patient contours. This implies the automatic generation of intensity modulated fields of heavy ions with as many as 40000 raster points, where each point corresponds to a specific beam position, energy and particle fluence. This set of data is directly passed to the beam delivery and control system. The treatment planning code has been in clinical use since the start of the GSI pilot project in December 1997. Forty-eight patients have been successfully planned and treated.

The ripple filter was designed to broaden the Bragg maximum of carbon beams for the raster-scan technique, a special type of tumour-conformal ion beam treatment. In this technique the target volume is divided into individual layers that are treated sequentially by varying the energy from the accelerator stepwise. Because the unmodified Bragg maximum has a small half-width, below 1 mm for small energies (< 160 MeV u(-1)), homogeneous irradiation at small penetration depths of 2-6 cm can only be obtained by using a large number of energy steps. If the energy step is too large, ripples are produced in the superimposed depth dose distribution. The ripple filter widens a Bragg peak to a Gaussian peak with a half-width of more than 2 mm. This helps to smooth the extended Bragg peak and to reduce the number of energy steps required by a factor of two to three, leading to significantly shorter overall irradiation times and a higher particle fluence per layer. The ripple filter consists of a 2 mm thick Plexiglass (PMMA) plate with a periodic structure of fine grooves. It can be mounted 60 cm upstream of the patient as a stationary device, because the fine structure of the grooves is completely washed out by the lateral scattering of the beam.

The purpose of this work was to design and manufacture a 3D range-modulator for scanned particle therapy. The modulator is intended to create a highly conformal dose distribution with only one fixed energy, simultaneously reducing considerably the treatment time. As a proof of concept, a 3D range-modulator was developed for a spherical target volume with a diameter of 5 cm, placed at a depth of 25 cm in a water phantom. It consists of a large number of thin pins with a well-defined shape and different lengths to modulate the necessary shift of the Bragg peak. The 3D range-modulator was manufactured with a rapid prototyping technique. The FLUKA Monte Carlo package was used to simulate the modulating effect of the 3D range-modulator and the resulting dose distribution. For that purpose, a special user routine was implemented to handle its complex geometrical contour. Additionally, FLUKA was extended with the capability of intensity modulated scanning. To validate the simulation results, dose measurements were carried out at the Heidelberg Ion Beam Therapy Center with a 400.41 MeV/u C beam. The high resolution dosimetric measurements show a good agreement between simulated and measured dose distributions. Irradiation of the monoenergetic raster plan took 3 s, which is approximately 20 times shorter than a comparable plan with 16 different energies. The combination of only one energy and a 3D range-modulator leads to a tremendous decrease in irradiation time. 'Interplay effects', typical for moving targets and pencil beam scanning, can be immensely reduced or disappear completely, making the delivery of a homogeneous dose to moving targets more reliable. Combining high dose conformity, very good homogeneity and extremely short irradiation times, the 3D range-modulator is considered to become a clinically applicable method for very fast treatment of lung tumours.

At present, beam ion beam therapy has started to spread worldwide. In Europe and Asia, combined carbon/proton facilities are favored, but in the US, only proton centers are under construction. This development is partially due to the different funding procedures and partially due to the more complex physical and especially biologic features of the heavy ions. In this article, the basic properties of both ions are presented, and their features for therapy are outlined. This refers to the dose conformity, the general precision of the treatment, and the ability to monitor via in-beam positron emission tomography the ions range inside the patient. Then the very complex biologic features are treated, and, finally, the treatment plans are compared.

Sub-millimetre-sized heterogeneities such as lung parenchyma cause Bragg peak degradation which can lead to an underdose of the tumor and an overdose of healthy tissue when not accounted for in treatment planning. Since commonly used treatment-planning CTs do not resolve the fine structure of lungs, this degradation can hardly be considered. We present a mathematical model capable of predicting and describing Bragg peak degradation due to a lung-equivalent geometry consisting of sub-millimetre voxels filled with either lung tissue or air. The material characteristic 'modulation power' is introduced to quantify the Bragg peak degradation. A strategy was developed to transfer the modulating effects of such fine structures to rougher structures such as 2 mm thick CT voxels, which is the resolution of typically used CTs. This is done by using the modulation power to derive a density distribution applicable to these voxels. By replacing the previously used sub-millimetre voxels by 2 mm thick voxels filled with lung tissue and modulating the lung tissue's density in each voxel individually, we were able to reproduce the Bragg peak degradation. Hence a solution is found to include Bragg curve degradation due to lung-equivalent materials in Monte Carlo-based treatment-planning systems.