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.

Since the end of 1996, we have treated more than 160 patients at PSI using spot-scanned protons. The range of indications treated has been quite wide and includes, in the head region, base-of-skull sarcomas, low-grade gliomas, meningiomas, and para-nasal sinus tumors. In addition, we have treated bone sarcomas in the neck and trunk--mainly in the sacral area--as well as prostate cases and some soft tissue sarcomas. PTV volumes for our treated cases are in the range 20-4500 ml, indicating the flexibility of the spot scanning system for treating lesions of all types and sizes. The number of fields per applied plan ranges from between 1 and 4, with a mean of just under 3 beams per plan, and the number of fluence modulated Bragg peaks delivered per field has ranged from 200 to 45 000. With the current delivery rate of roughly 3000 Bragg peaks per minute, this translates into delivery times per field of between a few seconds to 20-25 min. Bragg peak weight analysis of these spots has shown that over all fields, only about 10% of delivered spots have a weight of more than 10% of the maximum in any given field, indicating that there is some scope for optimizing the number of spots delivered per field. Field specific dosimetry shows that these treatments can be delivered accurately and precisely to within +/-1 mm (1 SD) orthogonal to the field direction and to within 1.5 mm in range. With our current delivery system the mean widths of delivered pencil beams at the Bragg peak is about 8 mm (sigma) for all energies, indicating that this is an area where some improvements can be made. In addition, an analysis of the spot weights and energies of individual Bragg peaks shows a relatively broad spread of low and high weighted Bragg peaks over all energy steps, indicating that there is at best only a limited relationship between pencil beam weighting and depth of penetration. This latter observation may have some consequences when considering strategies for fast re-scanning on second generation scanning gantries.

In this paper, we report on the clinical application of fully automated three-dimensional intensity modulated proton therapy, as applied to a 34-year-old patient presenting with a thoracic chordoma. Due to the anatomically challenging position of the lesion, a three-field technique was adopted in which fields incident through the lungs and heart, as well as beams directed directly at the spinal cord, could be avoided. A homogeneous target dose and sparing of the spinal cord was achieved through field patching and computer optimization of the 3D fluence of each field. Sensitivity of the resultant plan to delivery and calculational errors was determined through both the assessment of the potential effects of range and patient setup errors, and by the application of Monte Carlo dose calculation methods. Ionization chamber profile measurements and 2D dosimetry using a scintillator/CCD camera arrangement were performed to verify the calculated fields in water. Modeling of a 10% overshoot of proton range showed that the maximum dose to the spinal cord remained unchanged, but setup error analysis showed that dose homogeneity in the target volume could be sensitive to offsets in the AP direction. No significant difference between the MC and analytic dose calculations was found and the measured dosimetry for all fields was accurate to 3% for all measured points. Over the course of the treatment, a setup accuracy of +/-4 mm (2 s.d.) could be achieved, with a mean offset in the AP direction of 0.1 mm. Inhalation/exhalation CT scans indicated that organ motion in the region of the target volume was negligible. We conclude that 3D IMPT plans can be applied clinically and safely without modification to our existing delivery system. However, analysis of the calculated intensity matrices should be performed to assess the practicality, or otherwise, of the plan.

No abstract