A static beam delivery device for fast scanning proton arc-therapy.

Arc-therapy is a dose delivery technique regularly applied in photon radiation therapy, and is currently subject of great interest for proton therapy as well. In this technique, proton beams are aimed at a tumor from different continuous ranges of incident directions (so called 'arcs'). This technique can potentially yield a better dose conformity around the tumor and a very low dose in the surrounding healthy tissue. Currently, proton-arc therapy is performed by rotating a proton gantry around the patient, adapting the normally used dose-delivery method to the arc-specific motion of the gantry. Here we present first results from a feasibility study of the conceptual design of a new static fast beam delivery device/system for proton-arc therapy, which could be used instead of a gantry. In this novel concept, the incident angle of proton beams can be set rapidly by only changing field strengths of small magnets. This device eliminates the motion of the heavy gantry and related hardware. Therefore, a reduction of the total treatment time is expected. In the feasibility study presented here, we concentrate on the concept of the beam transport. Based on several simple, but realistic assumptions and approximations, proton tracking calculations were performed in a 3D magnetic field map, to calculate the beam transport in this device and to investigate and address several beam-optics challenges. We propose and simulate corresponding solutions and discuss their outcomes. To enable the implementation of some usually applied techniques in proton therapy, such as pencil beam scanning, energy modulation and beam shaping, we present and discuss our proposals. Here we present the concept of a new idea to perform fast proton arc-scanning and we report on first results of a feasibility study. Based on these results, we propose several options and next steps in the design.

Beam characterization and feasibility study for a small animal irradiation platform at clinical proton therapy facilities.

A deeper understanding of biological mechanisms to promote more efficient treatment strategies in proton therapy demands advances in preclinical radiation research. However this is often limited by insufficient availability of adequate infrastructures for precision image guided small animal proton irradiation. The project SIRMIO aims at filling this gap by developing a portable image-guided research platform for small animal irradiation, to be used at clinical facilities and allowing for a precision similar to a clinical treatment, when scaled down to the small animal size. This work investigates the achievable dosimetric properties of different lowest energy clinical proton therapy beams, manipulated by a dedicated portable beamline including active focusing after initial beam energy degradation and collimation. By measuring the lateral beam size in air close to the beam nozzle exit and the laterally integrated depth dose in water, an analytical beam model based on the beam parameters of the clinical beam at the Rinecker Proton Therapy Center was created for the lowest available clinical beam energy. The same approach was then applied to estimate the lowest energy beam model of different proton therapy facilities, Paul Scherrer Institute, Centre Antoine Lacassagne, Trento Proton Therapy Centre and the Danish Centre for Particle Therapy, based on their available beam commissioning data. This comparison indicated similar beam properties for all investigated sites, with emittance values of a few tens of mm·mrad. Finally, starting from these beam models, we simulated propagation through a novel beamline designed to manipulate the beam energy and size for precise small animal irradiation, and evaluated the resulting dosimetric properties in water. For all investigated initial clinical beams, similar dosimetric results suitable for small animal irradiation were found. This work supports the feasibility of the proposed SIRMIO beamline, promising suitable beam characteristics to allow for precise preclinical irradiation at clinical treatment facilities.

Uncertainty quantification analysis and optimization for proton therapy beam lines.

Since many years proton therapy is an effective treatment solution against deep-seated tumors. A precise quantification of sources of uncertainty in each proton therapy aspect (e.g. accelerator, beam lines, patient positioning, treatment planning) is of profound importance to increase the accuracy of the dose delivered to the patient. Together with Monte Carlo techniques, a new research field called Uncertainty Quantification (UQ) has been recently introduced to verify the robustness of the treatment planning. In this work we present the first application of UQ as a method to identify typical errors in the transport lines of a cyclotron-based proton therapy facility and analyze their impact on the properties of the therapeutic beams. We also demonstrate the potential of UQ methods in developing optimized beam optics solutions for high-dimensional problems. Sensitivity analysis and surrogate models offer a fast way to exclude unimportant parameters frcomplex optimization problems such as the design of a superconducting gantry performed at Paul Scherrer Institute in Switzerland.

Large energy acceptance gantry for proton therapy utilizing superconducting technology.

When using superconducting (SC) magnets in a gantry for proton therapy, the gantry will benefit from some reduction in size and a large reduction in weight. In this contribution we show an important additional advantage of SC magnets in proton therapy treatments. We present the design of a gantry with a SC bending section and achromatic beam optics with a very large beam momentum acceptance of [Formula: see text]15% (corresponding to about [Formula: see text]30% in the energy domain). Due to the related very large energy acceptance, approximately 70% of the treatments can be performed without changing the magnetic field for synchronization with energy modulation. In our design this is combined with a 2D lateral scanning system and a fast degrader mounted in the gantry, so that this gantry will be able to perform pencil beam scanning with very rapid energy variations at the patient, allowing a significant reduction of the irradiation time. We describe the iterative process we have applied to design the magnets and the beam transport, for which we have used different codes. COSY Infinity and OPAL have been used to design the beam transport optics and to track the particles in the magnetic fields, which are produced by the magnets designed in Opera. With beam optics calculations we have derived an optimal achromatic beam transport with the large momentum acceptance of the proton pencil beam and we show the agreement with particle tracking calculations in the 3D magnetic field map. A new cyclotron based facility with this gantry will have a significantly smaller footprint, since one can refrain from the standard degrader and energy selection system behind the cyclotron. In the treatments, this gantry will enable a very fast proton beam delivery sequence, which may be of advantage for treatments in moving tissue.

Proton range verification using a range probe: definition of concept and initial analysis.

Proton therapy can be particularly sensitive to changes or errors in range. Thus, methods for the in vivo measurement of range could be of great use to improve the quality and accuracy of proton-based radiotherapy. In this paper, we introduce the concept of the 'range probe'. This is a low-dose, high-energy proton pencil beam that would pass through a patient, and whose integral Bragg peak would be measured on the exit side using a multi-layer detector. We propose that by comparing the measured integral Bragg peak with that calculated based on the patient's planning CT, such a range probe could provide useful information about the accuracy of range calculations in vivo. To study the feasibility of this approach, a Monte Carlo-based study has been performed. Using a patient's planning CT, MC simulations (VMCpro) have been made for single pencil beams laterally traversing the head and stopping in a simulated range telescope behind the patient. Range probes have been calculated for different locations, and the residual range from the Bragg peak 'signal' in the range telescope has been assessed for different assumed detector thicknesses. The sensitivity of this approach to changes in CT values, calibration curve and positional shifts of the CT have been investigated. From our analysis, range resolutions of 1 mm may be possible with a detector thickness of 4 mm for homogeneous regions. Additionally, for heterogeneous regions, changes of the Bragg peak shape due to spatial shifts of the CT could be a sensitive measure for detecting patient set-up errors directly in the treatment position. The concept of the proton 'range probe' appears to be feasible for high-resolution range verification. We now want to test this concept experimentally using different possible range telescope detectors.

Geant4 simulations of proton beam transport through a carbon or beryllium degrader and following a beam line.

Monte Carlo simulations based on the Geant4 simulation toolkit were performed for the carbon wedge degrader used in the beam line at the Center of Proton Therapy of the Paul Scherrer Institute (PSI). The simulations are part of the beam line studies for the development and understanding of the GANTRY2 and OPTIS2 treatment facilities at PSI, but can also be applied to other beam lines. The simulated stopping power, momentum distributions at the degrader exit and beam line transmission have been compared to accurate benchmark measurements. Because the beam transport through magnetic elements is not easily modeled using Geant4a connection to the TURTLE beam line simulation program was made. After adjusting the mean ionization potential of the carbon degrader material from 78 eV to 95 eV, we found an accurate match between simulations and benchmark measurements, so that the simulation model could be validated. We found that the degrader does not completely erase the initial beam phase space even at low degraded beam energies. Using the validation results, we present a study of the usability of beryllium as a degrader material (mean ionization potential 63.7 eV). We found an improvement in the transmission of 30-45%, depending on the degraded beam energy, the higher value for the lower energies.

2D dosimetry in a proton beam with a scintillating GEM detector.

A two-dimensional position-sensitive dosimetry system based on a scintillating gas detector is being developed for pre-treatment verification of dose distributions in particle therapy. The dosimetry system consists of a chamber filled with an Ar/CF(4) scintillating gas mixture, inside which two gas electron multiplier (GEM) structures are mounted (Seravalli et al 2008b Med. Phys. Biol. 53 4651-65). Photons emitted by the excited Ar/CF(4) gas molecules during the gas multiplication in the GEM holes are detected by a mirror-lens-CCD camera system. The intensity distribution of the measured light spot is proportional to the 2D dose distribution. In this work, we report on the characterization of the scintillating GEM detector in terms of those properties that are of particular importance in relative dose measurements, e.g. response reproducibility, dose dependence, dose rate dependence, spatial and time response, field size dependence, response uniformity. The experiments were performed in a 150 MeV proton beam. We found that the detector response is very stable for measurements performed in succession (sigma = 0.6%) and its response reproducibility over 2 days is about 5%. The detector response was found to be linear with the dose in the range 0.05-19 Gy. No dose rate effects were observed between 1 and 16 Gy min(-1) at the shallow depth of a water phantom and 2 and 38 Gy min(-1) at the Bragg peak depth. No field size effects were observed in the range 120-3850 mm(2). A signal rise and fall time of 2 micros was recorded and a spatial response of

Characterization of a scintillating GEM detector with low energy x-rays.

A two-dimensional position-sensitive dosimetry system based on a scintillating gas detector is being developed with the aim of using it for pre-treatment verification of dose distributions in charged particle therapy. The dosimetry system consists of a chamber filled with an Ar/CF(4) scintillating gas mixture, inside which two cascaded gas electron multipliers (GEMs) are mounted. A GEM is a thin kapton foil with copper cladding structured with a regular pattern of sub-mm holes. In such a system, light quanta are emitted by the scintillating gas mixture during the electron avalanches in the GEM holes when radiation traverses the detector. The light intensity distribution is proportional to the energy deposited in the detector's sensitive volume by the beam. In the present work, we investigated the optimization of the scintillating GEM detector light yield. The light quanta are detected by means of a CCD camera or a photomultiplier tube coupled to a monochromator. The GEM charge signal is measured simultaneously. We have found that with 60 microm diameter double conical GEM holes, a brighter light signal and a higher electric signal are obtained than with 80 microm diameter holes. With an Ar + 8% CF(4) volume concentration, the highest voltage across the GEMs and the largest light and electric signals were reached. Moreover, we have found that the emission spectrum of Ar/CF(4) is independent of (1) the voltages applied across the GEMs, (2) the x-ray beam intensity and (3) the GEM hole diameter. On the other hand, the ratio of Ar to CF(4) peaks in the spectrum changes when the concentration of the latter gas is varied.

A scintillating gas detector for 2D dose measurements in clinical carbon beams.

A two-dimensional position sensitive dosimetry system based on a scintillating gas detector has been developed for pre-treatment verification of dose distributions in hadron therapy. The dosimetry system consists of a chamber filled with an Ar/CF4 scintillating gas mixture, inside which two cascaded gas electron multipliers (GEMs) are mounted. A GEM is a thin kapton foil with copper cladding structured with a regular pattern of sub-mm holes. The primary electrons, created in the detector's sensitive volume by the incoming beam, drift in an electric field towards the GEMs and undergo gas multiplication in the GEM holes. During this process, photons are emitted by the excited Ar/CF4 gas molecules and detected by a mirror-lens-CCD camera system. Since the amount of emitted light is proportional to the dose deposited in the sensitive volume of the detector by the incoming beam, the intensity distribution of the measured light spot is proportional to the 2D hadron dose distribution. For a measurement of a 3D dose distribution, the scintillating gas detector is mounted at the beam exit side of a water-bellows phantom, whose thickness can be varied in steps. In this work, the energy dependence of the output signal of the scintillating gas detector has been verified in a 250 MeV/u clinical 12C ion beam by means of a depth-dose curve measurement. The underestimation of the measured signal at the Bragg peak depth is only 9% with respect to an air-filled ionization chamber. This is much smaller than the underestimation found for a scintillating Gd2O2S:Tb ('Lanex') screen under the same measurement conditions (43%). Consequently, the scintillating gas detector is a promising device for verifying dose distributions in high LET beams, for example to check hadron therapy treatment plans which comprise beams with different energies.

Estimation of parameters of dose-volume models and their confidence limits.

Predictions of the normal-tissue complication probability (NTCP) for the ranking of treatment plans are based on fits of dose-volume models to clinical and/or experimental data. In the literature several different fit methods are used. In this work frequently used methods and techniques to fit NTCP models to dose response data for establishing dose-volume effects, are discussed. The techniques are tested for their usability with dose-volume data and NTCP models. Different methods to estimate the confidence intervals of the model parameters are part of this study. From a critical-volume (CV) model with biologically realistic parameters a primary dataset was generated, serving as the reference for this study and describable by the NTCP model. The CV model was fitted to this dataset. From the resulting parameters and the CV model, 1000 secondary datasets were generated by Monte Carlo simulation. All secondary datasets were fitted to obtain 1000 parameter sets of the CV model. Thus the 'real' spread in fit results due to statistical spreading in the data is obtained and has been compared with estimates of the confidence intervals obtained by different methods applied to the primary dataset. The confidence limits of the parameters of one dataset were estimated using the methods, employing the covariance matrix, the jackknife method and directly from the likelihood landscape. These results were compared with the spread of the parameters, obtained from the secondary parameter sets. For the estimation of confidence intervals on NTCP predictions, three methods were tested. Firstly, propagation of errors using the covariance matrix was used. Secondly, the meaning of the width of a bundle of curves that resulted from parameters that were within the one standard deviation region in the likelihood space was investigated. Thirdly, many parameter sets and their likelihood were used to create a likelihood-weighted probability distribution of the NTCP. It is concluded that for the type of dose response data used here, only a full likelihood analysis will produce reliable results. The often-used approximations, such as the usage of the covariance matrix, produce inconsistent confidence limits on both the parameter sets and the resulting NTCP values.

Techniques for precision irradiation of the lateral half of the rat cervical spinal cord using 150 MeV protons [corrected].

Techniques for high precision irradiation experiments with protons, to investigate the volume dependence of the tolerance dose of the rat cervical spinal cord are described. In the present study, 50% of the lateral cross section of the spinal cord was irradiated. The diameter of the cross section of this part of the rat spinal cord is at maximum 3.5 mm. Therefore, a dedicated procedure was developed to comply with the needs for a very high positioning accuracy and high spatial resolution dosimetry. By using 150 MeV protons a steep dose gradient (20-80% = 1 mm) in the centre of the spinal cord was achieved. This yields a good dose contrast between the left and right halves of the cord. A home-made digital x-ray imager with a pixel resolution of 0.18 mm/pixel was used for position verification of the spinal cord. A positioning accuracy of 0.09 mm was obtained by using information of multiple pixels. The average position stability during the irradiation was found to be 0.08 mm (1 SD) without significant systematic deviations. Profiles of the dose distribution were measured with a 2D dosimetry system consisting of a scintillating screen and a CCD camera. Dose volume histograms of the whole spinal cord as well as separately of the white and grey matters were calculated using MRI imaging of the cross section of the rat cervical spinal cord. From the irradiation of 20 animals a dose-response curve has been established. MRI showed radiation-induced damage at the high dose side of the spinal cord. Analysis of the preliminary dose-response data shows a significant dose-volume effect. With the described procedure and equipment it is possible to perform high precision irradiations on selected parts of the spinal cord.

Collimator scatter and 2D dosimetry in small proton beams.

Monte Carlo simulations have been performed to determine the influence of collimator-scattered protons from a 150 MeV proton beam on the dose distribution behind a collimator. Slit-shaped collimators with apertures between 2 and 20 mm have been simulated. The Monte Carlo code GEANT 3.21 has been validated against one-dimensional dose measurements with a scintillating screen, observed by a CCD camera. In order to account for the effects of the spatial response of the CCD/scintillator system, the line-spread function was determined by comparison with measurements made with a diamond detector. The line-spread function of the CCD/scintillator system is described by a Gaussian distribution with a standard deviation of 0.22 mm. The Monte Carlo simulations show that protons that hit the collimator on the entrance face and leave it through the wall of the aperture make the largest scatter contribution. Scatter on air is the major contribution to the extent of the penumbra. From the energy spectra it is derived that protons with a relative biological effectiveness greater than 1 cause at most 1% more damage in tissue than what would be expected from the physical dose.

Performance of a fluorescent screen and CCD camera as a two-dimensional dosimetry system for dynamic treatment techniques.

A two-dimensionally position sensitive dosimetry system has been tested for different dosimetric applications in a radiation therapy facility with a scanning proton beam. The system consists of a scintillating (fluorescent) screen, mounted at the beam-exit side of a phantom and it is observed by a charge coupled device (CCD) camera. The observed light distribution at the screen is equivalent to the two-dimensional (2D)-dose distribution at the screen position. It has been found that the dosimetric properties of the system, measured in a scanning proton beam, are equal to those measured in a proton beam broadened by a scattering system. Measurements of the transversal dose distribution of a single pencil beam are consistent with dose measurements as well as with dose calculations in clinically relevant fields made with multiple pencil beams. Measurements of inhomogeneous dose distributions have shown to be of sufficient accuracy to be suitable for the verification of dose calculation algorithms. The good sensitivity and sub-mm spatial resolution of the system allows for the detection of deviations of a few percent in dose from the expected (intended or calculated) dose distribution. Its dosimetric properties and the immediate availability of the data make this device a useful tool in the quality control of scanning proton beams.

Verification of the alignment of a therapeutic radiation beam relative to its patient positioner.

An easily-used system has been developed for routine measurements of the alignment of beams used for radiation therapy. The position of a beam of circular cross section is measured with respect to a steel sphere fixed to the patient positioning table and which should coincide with the isocenter. Since measurements can be done at all gantry angles (if one is available) and with all possible orientations of the patient table, the system is particularly suited for rapid and accurate measurements of gantry and/or couch isocentricity. Because it directly measures beam-to-positioner offset, the system provides an inclusive alignment verification of the total treatment system. The system has been developed for use with proton beams, but it could equally be used for alignment checks of an x-ray beam from a linear accelerator or other source. The measuring instrument consists of a scintillation screen viewed by a CCD camera, mounted on the gantry downstream of the sphere. The steel sphere is not large enough to stop protons of all energies of interest; however, it will always modify the energy and direction of protons which intersect it, creating a region of lower intensity (a "shadow") in the light spot created by the proton beam hitting the screen. The position of the shadow with respect to the light spot is a measure of the alignment of the system. An image-analysis algorithm has been developed for an automatic determination of the position of the shadow with respect to the light spot. The specifications and theoretical analysis of the system have been derived from Monte Carlo simulations, which are validated by measurements. We have demonstrated that the device detects beam misalignments with an accuracy (1 s.d.) of 0.05 mm, which is in agreement with the expected performance. This accuracy is more than sufficient to detect the maximum allowed misalignment of +/-0.5 mm.

Fast 2D phantom dosimetry for scanning proton beams.

A quality control system especially designed for dosimetry in scanning proton beams has been designed and tested. The system consists of a scintillating screen (Gd2O2S:Tb), mounted at the beam-exit side of a phantom, and observed by a low noise CCD camera with a long integration time. The purpose of the instrument is to make a fast and accurate two-dimensional image of the dose distribution at the screen position in the phantom. The linearity of the signal with the dose, the noise in the signal, the influence of the ionization density on the signal, and the influence of the field size on the signal have been investigated. The spatial resolution is 1.3 mm (1 s.d.), which is sufficiently smaller than typical penumbras in dose distributions. The measured yield depends linearly on the dose and agrees within 5% with the calculations. In the images a signal to noise ration (signal/1 s.d.) of 10(2) has been found, which is in the same order of magnitude as expected from the calculations. At locations in the dose distribution possessing a strong contribution of high ionization densities (i.e., in the Bragg peak), we found some quenching of the light output, which can be described well by existing models if the beam characteristics are known. For clinically used beam characteristics such as a Spread Out Bragg peak, there is at most 8% deviation from the NACP ionization chamber measurements. The conclusion is that this instrument is a useful tool for quick and reliable quality control of proton beams. The long integration-time capabilities of the system make it worthwhile to investigate its applicability in scanning proton beams and other dynamic treatment modalities.