We present a comprehensive analytical comparison of four types of proton imaging set-ups and, to this end, develop a mathematical framework to calculate the width of the uncertainty envelope around the most likely proton path depending on set-up geometry, detector properties, and proton beam parameters. As a figure of merit for the spatial resolution achievable with each set-up, we use the frequency [Formula: see text] at which the modular transfer function of a density step decreases below 10%. We verify the analytical results with Monte Carlo simulations. We find that set-ups which track the angle and position of individual protons in front of and behind the phantom would yield an average spatial resolution of 0.3-0.35 lp mm assuming realistic geometric parameters (i.e. 30-40 cm distance between detector and phantom, 15-20 cm phantom thickness). For set-ups combining pencil beam scanning with either a position sensitive detector, e.g. an x-ray flat panel, or with a position insensitive detector, e.g. a range telescope, we find an average spatial resolution of about 0.1 lp mm for an 8 mm FWHM beam spot size. The pixel information improves the spatial resolution by less than 10%. In both set-up types, performance can be significantly improved by reducing the pencil beam size down to 2 mm FWHM. In this case, the achievable spatial resolution reaches about 0.25 lp mm. Our results show that imaging set-ups combining double scattering with a pixel detector can provide sufficient spatial resolution only under very stringent conditions and are not ideally suited for computed tomography applications. We further propose a region-of-interest method for set-ups with a pixel detector to filter out protons which have undergone nuclear reactions and discuss the impact of tracker detector uncertainties on the most likely path.

Therapeutic proton and heavier ion beams generate prompt gamma photons that may escape from the patient. In principle, this allows for real-time, in situ monitoring of the treatment delivery, in particular, the hadron range within the patient, by imaging the emitted prompt gamma rays. Unfortunately, the neutrons simultaneously created with the prompt photons create a background that may obscure the prompt gamma signal. To enhance the accuracy of proton dose verification by prompt gamma imaging, we therefore propose a time-of-flight (TOF) technique to reject this neutron background, involving a shifting time window to account for the propagation of the protons through the patient. Time-resolved Monte Carlo simulations of the generation and transport of prompt gamma photons and neutrons upon irradiation of a PMMA phantom with 100, 150 and 200 MeV protons were performed using Geant4 (version 9.2.p02) and MCNPX (version 2.7.D). The influence of angular collimation and TOF selection on the prompt gamma and neutron longitudinal profiles is studied. Furthermore, the implications of the proton beam microstructure (characterized by the proton bunch width and repetition period) are investigated. The application of a shifting TOF window having a width of ΔTOF(z) = 1.0 ns appears to reduce the neutron background by more than 99%. Subsequent application of an energy threshold does not appear to sharpen the distal falloff of the prompt gamma profile but reduces the tail that is observed beyond the proton range. Investigations of the influence of the beam time structure show that TOF rejection of the neutron background is expected to be effective for typical therapeutic proton cyclotrons.

The aim of this work is to extend a widely used proton Monte Carlo tool, TOPAS, towards the modeling of relative biological effect (RBE) distributions in experimental arrangements as well as patients. TOPAS provides a software core which users configure by writing parameter files to, for instance, define application specific geometries and scoring conditions. Expert users may further extend TOPAS scoring capabilities by plugging in their own additional C++ code. This structure was utilized for the implementation of eight biophysical models suited to calculate proton RBE. As far as physics parameters are concerned, four of these models are based on the proton linear energy transfer (LET), while the others are based on DNA Double Strand Break (DSB) induction and the frequency-mean specific energy, lineal energy, or delta electron generated track structure. The biological input parameters for all models are typically inferred from fits of the models to radiobiological experiments. The model structures have been implemented in a coherent way within the TOPAS architecture. Their performance was validated against measured experimental data on proton RBE in a spread-out Bragg peak using V79 Chinese Hamster cells. This work is an important step in bringing biologically optimized treatment planning for proton therapy closer to the clinical practice as it will allow researchers to refine and compare pre-defined as well as user-defined models.

Ion beams exhibit a finite range and an inverted depth-dose profile, the Bragg peak. These favorable physical properties allow excellent tumor-dose conformality. However, they introduce sensitivity to range uncertainties. Although these uncertainties are typically taken into account in treatment planning, delivery of the intended dose to the patient has to be ensured daily to prevent underdosage of the tumor or overdosage of surrounding critical structures. Thus, imaging techniques play an increasingly important role for treatment planning and in situ monitoring in ion beam therapy. At the Heidelberg Ion Beam Therapy (HIT) center, a prototype detector system based on a stack of 61 ionization chambers has been assembled for the purpose of radiographic and tomographic imaging of transmitted energetic ions. Its applicability to ion-based transmission imaging was investigated experimentally. An extensive characterization of the set-up in terms of beam parameters and settings of the read-out electronics was performed. Overall, the findings of this work support the potential of an efficient experimental set-up as the range telescope equipped with high sensitivity and fast electronics to perform heavy ion radiography and tomography at HIT.

Prompt-gamma emission detection is a promising technique for hadrontherapy monitoring purposes. In this regard, obtaining prompt-gamma yields that can be used to develop monitoring systems based on this principle is of utmost importance since any camera design must cope with the available signal. Herein, a comprehensive study of the data from ten single-slit experiments is presented, five consisting in the irradiation of either PMMA or water targets with lower and higher energy carbon ions, and another five experiments using PMMA targets and proton beams. Analysis techniques such as background subtraction methods, geometrical normalization, and systematic uncertainty estimation were applied to the data in order to obtain absolute prompt-gamma yields in units of prompt-gamma counts per incident ion, unit of field of view, and unit of solid angle. At the entrance of a PMMA target, where the contribution of secondary nuclear reactions is negligible, prompt-gamma counts per incident ion, per millimetre and per steradian equal to (124 ± 0.7stat ± 30sys) × 10(-6) for 95 MeV u(-1) carbon ions, (79 ± 2stat ± 23sys) × 10(-6) for 310 MeV u(-1) carbon ions, and (16 ± 0.07stat ± 1sys) × 10(-6) for 160 MeV protons were found for prompt gammas with energies higher than 1 MeV. This shows a factor 5 between the yields of two different ions species with the same range in water (160 MeV protons and 310 MeV u(-1) carbon ions). The target composition was also found to influence the prompt-gamma yield since, for 300/310 MeV u(-1) carbon ions, a 42% greater yield ((112 ± 1stat ± 22sys) × 10(-6) counts ion(-1) mm(-1) sr(-1)) was obtained with a water target compared to a PMMA one.

Objectives: To describe the measurements and to present the results of the beam commissioning and the beam model validation of a compact, gantry-mounted, spot scanning proton accelerator system with dynamic layer-by-layer field collimation. Methods: We performed measurements of depth dose distributions in water, spot and scanned field size in air at different positions from the isocenter plane, spot position over the 20 × 20 cm2 scanned area, beam monitor calibration in terms of absorbed dose to water and specific field collimation measurements at different gantry angles to commission the system. To validate the beam model in the treatment planning system (TPS), we measured spot profiles in water at different depths, absolute dose in water of single energy layers of different field sizes and inversely optimised spread-out Bragg peaks (SOBP) under normal and oblique beam incidence, field size and penumbra in water of SOBPs, and patient treatment specific quality assurance in homogeneous and heterogeneous phantoms. Results: Energy range, spot size, spot position and dose output were consistent at all gantry angles with 0.3 mm, 0.4 mm, 0.6 mm and 0.5% maximum deviations, respectively. Uncollimated spot size (one sigma) in air with an air-gap of 10 cm ranged from 4.1 to 16.4 mm covering a range from 32.2 to 1.9 cm in water, respectively. Absolute dose measurements were within 3% when comparing TPS and experimental data. Gamma pass rates >98% and >96% at 3%/3 mm were obtained when performing 2D dose measurements in homogeneous and in heterogeneous media, respectively. Leaf position was within ±1 mm at all gantry angles and nozzle positions. Conclusions: Beam characterisation and machine commissioning results, and the exhaustive end-to-end tests performed to assess the proper functionality of the system, confirm that it is safe and accurate to treat patients. Advances in knowledge: This is the first paper addressing the beam commissioning and the beam validation of a compact, gantry-mounted, pencil beam scanning proton accelerator system with dynamic layer-by-layer multileaf collimation.

Ion beams offer an excellent tumor-dose conformality due to their inverted depth-dose profile and finite range in tissue, the Bragg peak (BP). However, they introduce sensitivity to range uncertainties. Imaging techniques play an increasingly important role in ion beam therapy to support precise diagnosis and identification of the target volume at the planning stage as well as to ensure the correspondence between the planning and treatment situation at the actual irradiation. For the purpose of improved treatment quality, ion-based radiographic images could be acquired at the treatment site before or during treatment and be employed to monitor the patient positioning and to check the patient-specific ion range. This work presents the initial experimental investigations carried out to address the feasibility of carbon ion radiography at the Heidelberg ion therapy center using a prototype range telescope set-up and an active raster scanning ion beam delivery system. Bragg curves are measured with a stack of ionization chambers (IC) synchronously to the beam delivery. The position of the BP is extracted from the data by locating the channel of maximum current signal for each delivered beam. Each BP is associated to the lateral and vertical positions of the scanned raster point extrapolated from the beam monitor system to build up a radiography. The radiographic images are converted into water equivalent thickness (WET) based on two calibrations of the detector. Radiographies of two phantoms of different complexities are reconstructed and their image quality is analyzed. A novel method proposed to increase the nominal range resolution of the IC stack is applied to the carbon ion radiography of an Alderson head phantom. Moreover, an x-ray digitally reconstructed radiography of the same anthropomorphic head phantom is converted in WET through the clinically used ion range calibration curve and compared with the carbon ion radiography based on a γ-index approach, yielding a good correspondence in terms of absolute WET within ±3%, 3 mm distance-to-agreement and, 87% passing ratio. Imaging artifacts at interfaces within the irradiated phantom due to the finite size of the beam, resulting in multiple maxima, are addressed. Overall, this work demonstrates the feasibility of the prototype range telescope to acquire ion-based transmission imaging with a resolution of up to 0.8 mm WET.

Ion beam therapy enables a highly accurate dose conformation delivery to the tumor due to the finite range of charged ions in matter (i.e. Bragg peak (BP)). Consequently, the dose profile is very sensitive to patients anatomical changes as well as minor mispositioning, and so it requires improved dose control techniques. Proton interaction vertex imaging (IVI) could offer an online range control in carbon ion therapy. In this paper, a statistical method was used to study the sensitivity of the IVI technique on experimental data obtained from the Heidelberg Ion-Beam Therapy Center. The vertices of secondary protons were reconstructed with pixelized silicon detectors. The statistical study used the [Formula: see text] test of the reconstructed vertex distributions for a given displacement of the BP position as a function of the impinging carbon ions. Different phantom configurations were used with or without bone equivalent tissue and air inserts. The inflection points in the fall-off region of the longitudinal vertex distribution were computed using different methods, while the relation with the BP position was established. In the present setup, the resolution of the BP position was about 4-5 mm in the homogeneous phantom under clinical conditions (10 incident carbon ions). Our results show that the IVI method could therefore monitor the BP position with a promising resolution in clinical conditions.