Computer tomographic (CT) scans are used to correct for tissue inhomogeneities in radiotherapy treatment planning. In order to guarantee a precise treatment, it is important to obtain the relationship between CT Hounsfield units and electron densities (or proton stopping powers for proton radiotherapy), which is the basic input for radiotherapy planning systems which consider tissue heterogeneities. A method is described to determine improved CT calibrations for biological tissue (a stoichiometric calibration) based on measurements using tissue equivalent materials. The precision of this stoichiometric calibration and the more usual tissue substitute calibration is determined by a comparison of calculated proton radiographic images based on these calibrations and measured radiographs of a biological sample. It has been found that the stoichiometric calibration is more precise than the tissue substitute calibration.

In this paper we present the pencil beam dose model used for treatment planning at the PSI proton gantry, the only system presently applying proton therapy with a beam scanning technique. The scope of the paper is to give a general overview on the various components of the dose model, on the related measurements and on the practical parametrization of the results. The physical model estimates from first physical principles absolute dose normalized to the number of incident protons. The proton beam flux is measured in practice by plane-parallel ionization chambers (ICs) normalized to protons via Faraday-cup measurements. It is therefore possible to predict and deliver absolute dose directly from this model without other means. The dose predicted in this way agrees very well with the results obtained with ICs calibrated in a cobalt beam. Emphasis is given in this paper to the characterization of nuclear interaction effects, which play a significant role in the model and are the major source of uncertainty in the direct estimation of the absolute dose. Nuclear interactions attenuate the primary proton flux, they modify the shape of the depth-dose curve and produce a faint beam halo of secondary dose around the primary proton pencil beam in water. A very simple beam halo model has been developed and used at PSI to eliminate the systematic dependences of the dose observed as a function of the size of the target volume. We show typical results for the relative (using a CCD system) and absolute (using calibrated ICs) dosimetry, routinely applied for the verification of patient plans. With the dose model including the nuclear beam halo we can predict quite precisely the dose directly from treatment planning without renormalization measurements, independently of the dose, shape and size of the dose fields. This applies also to the complex non-homogeneous dose distributions required for the delivery of range-intensity-modulated proton therapy, a novel therapy technique developed at PSI.

The new proton therapy facility is being assembled at the Paul Scherrer Institute (PSI). The beam delivered by the PSI sector cyclotron can be split and brought into a new hall where it is degraded from 590 MeV down to an energy in the range of 85-270 MeV. A new beam line following the degrader is used to clean the low-energetic beam in phase space and momentum band. The analyzed beam is then injected into a compact isocentric gantry, where it is applied to the patient using a new dynamic treatment modality, the so-called spot-scanning technique. This technique will permit full three-dimensional conformation of the dose to the target volume to be realized in a routine way without the need for individualized patient hardware like collimators and compensators. By combining the scanning of the focused pencil beam within the beam optics of the gantry and by mounting the patient table eccentrically on the gantry, the diameter of the rotating structure has been reduced to only 4 m. In the article the degrees of freedom available on the gantry to apply the beam to the patient (with two rotations for head treatments) are also discussed. The devices for the positioning of the patient on the gantry (x rays and proton radiography) and outside the treatment room (the patient transporter system and the modified mechanics of the computer tomograph unit) are briefly presented. The status of the facility and first experimental results are introduced for later reference.

The precision in proton radiotherapy treatment planning depends on the accuracy of the information used to calculate the stopping power properties of the tissues in the patient's body. This information is obtained from computed tomography (CT) images using a calibration curve to convert CT Hounsfield units into relative proton stopping power values. The validity of a stoichiometric method to create the calibration curve has been verified by measuring pairs of Hounsfield units and stopping power values for animal tissue samples. It was found that the agreement between measurement and calibration curve is better than 1% if beam hardening effects in the acquisition of the CT images can be neglected. The influence of beam hardening effects on the quantitative reading of the CT measurements is discussed and an estimation for the overall range precision of proton beams is given. It is expected that the range of protons in the human body can be controlled to better than +/-1.1% of the water equivalent range in soft tissue and +/-1.8% in bone, which translates into a range precision of about 1-3 mm in typical treatment situations.

The gantry for proton radiotherapy at the Paul Scherrer Institute (PSI) is designed specifically for the spot-scanning technique. Use of this technique to its full potential requires dose calculation algorithms which are capable of precisely simulating each scanned beam individually. Different specialized analytical dose calculations have been developed, which attempt to model the effects of density heterogeneities in the patient's body on the dose. Their accuracy has been evaluated by a comparison with Monte Carlo calculated dose distributions in the case of a simple geometrical density interface parallel to the beam and typical anatomical situations. A specialized ray casting model which takes range dilution effects (broadening of the spectrum of proton ranges) into account has been found to produce results of good accuracy. This algorithm can easily be implemented in the iterative optimization procedure used for the calculation of the optimal contribution of each individual scanned pencil beam. In most cases an elemental pencil beam dose calculation has been found to be most accurate. Due to the long computing time, this model is currently used only after the optimization procedure as an alternative method of calculating the dose.