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.

A three-dimensional spot-scanning technique for radiotherapy with protons is being developed at the Paul Scherrer Institute. As part of the effort to optimize the design and ensure clinically useful dose distributions, a computer simulation of the dose deposition in the presence of respiratory motion was performed. Preliminary experiments have characterized the proton beam and the scanning procedure. Using these parameters, the computer program calculated the dose within a uniform volume of water in the presence of respiratory motion. Respiration amplitude, respiration period, respiration direction, number of fractions, size and position of the beamspots and rescanning multiplicity were systematically varied and the effect on the dose distribution determined. The dose uniformity is very dependent on the direction of the respiration relative to the three independent beam scanning directions. The dose uniformity decreases with increasing respiration amplitude, but has little response to changes in respiration frequency. Rescanning the volume, such as with fractionation, improves the dose uniformity roughly as the square root of the number of fractions. Broad, Gaussian beams result in better dose uniformity than narrow, sharply delineated ones, but produce slower dose fall-off at the edges of the scanned volume. Results of this work are being incorporated into the design of the system.

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