Experimental comparison of cylindrical and plane parallel ionization chambers for reference dosimetry in continuous and pulsed scanned proton beams.

Objective. In this experimental work we compared the determination of absorbed dose to water using four ionization chambers (ICs), a PTW-34045 Advanced Markus, a PTW-34001 Roos, an IBA-PPC05 and a PTW-30012 Farmer, irradiated under the same conditions in one continuous- and in two pulsed-scanned proton beams.Approach. The ICs were positioned at 2 cm depth in a water phantom in four square-field single-energy scanned-proton beams with nominal energies between 80 and 220 MeV and in the middle of 10 × 10 × 10 cm3dose cubes centered at 10 cm or 12.5 cm depth in water. The water-equivalent thickness (WET) of the entrance window and the effective point of measurement was considered when positioning the plane parallel (PP) ICs and the cylindrical ICs, respectively. To reduce uncertainties, all ICs were calibrated at the same primary standards laboratory. We used the beam quality (kQ) correction factors for the ICs under investigation from IAEA TRS-398, the newly calculated Monte Carlo (MC) values and the anticipated IAEA TRS-398 updated recommendations.Main results. Dose differences among the four ICs ranged between 1.5% and 3.7% using both the TRS-398 and the newly recommendedkQvalues. The spread among the chambers is reduced with the newlykQvalues. The largest differences were observed between the rest of the ICs and the IBA-PPC05 IC, obtaining lower dose with the IBA-PPC05.Significance. We provide experimental data comparing different types of chambers in different proton beam qualities. The observed dose differences between the ICs appear to be related to inconsistencies in the determination of thekQvalues. For PP ICs, MC studies account for the physical thickness of the entrance window rather than the WET. The additional energy loss that the wall material invokes is not negligible for the IBA-PPC05 and might partially explain the lowkQvalues determined for this IC. To resolve this inconsistency and to benchmark MC values,kQvalues measured using calorimetry are needed.

Quality assurance of scanned proton beams at different gantry angles using an ionization chamber array in a rotational phantom.

PURPOSE: To implement a single set-up monthly QA procedure for 9 different beam parameters at different gantry angles and evaluate its clinical implementation over a 12 month period. METHODS: We developed a QA procedure using an array detector (PTW Octavius 1500XDR) embedded in a rotational unit (PTW Octavius 4D) at our proton facility. With a single set-up we can monitor field central axis position, field symmetry, field size, flatness, penumbrae, output, spot size, spot position and range at different gantry angles (AAPM TG 224). The set-up is irradiated with homogenous 2D fields with dynamic aperture and spot patterns at five gantry angles. A modular top is used to check the range consistency. Absolute γ analysis were performed to compare measured dose distributions to calculated dose. All other parameters are directly extracted from the measurements. Additionally, the sensitivity of the set-up to small changes in beam parameters were compared to the Lynx detector (IBA). RESULTS: Over a 12 month period, output, symmetry, and flatness were within ± 2 %; FWHM, spot positions, penumbra widths, and central axis fields were within ± 1 mm. Range differences were all within 1/2 of the energy spacing (±0.6 MeV) relative to baseline. Most (2 %, 2 mm) γ-analysis showed agreement scores higher than 90 %. The sensitivity is comparable to the Lynx detector and measurement time is reduced by 40 %. CONCLUSION: The time-efficient monthly QA procedure that we developed can accurately be used to measure a large range of beam parameters at different gantry angles, within the TG 224 AAPM recommendations.

EURADOS REM-COUNTER INTERCOMPARISON AT MAASTRO PROTON THERAPY CENTRE: COMPARISON WITH LITERATURE DATA.

The Maastro Proton Therapy Centre is the first European facility housing the Mevion S250i Hyperscan synchrocyclotron. The proximity of the accelerator to the patient, the presence of an active pencil beam delivery system downstream of a passive energy degrader and the pulsed structure of the beam make the Mevion stray neutron field unique amongst proton therapy facilities. This paper reviews the results of a rem-counter intercomparison experiment promoted by the European Radiation Dosimetry Group at Maastro and compares them with those at other proton therapy facilities. The Maastro neutron H*(10) in the room (100-200 µSv/Gy at about 2 m from the isocentre) is in line with accelerators using purely passive or wobbling beam delivery modalities, even though Maastro shows a dose gradient peaked near the accelerator. Unlike synchrotron- and cyclotron-based facilities, the pulsed beam at Maastro requires the employment of rem-counters specifically designed to withstand pulsed neutron fields.

Treatment planning comparison in the PROTECT-trial randomising proton versus photon beam therapy in oesophageal cancer: Results from eight European centres.

PURPOSE: To compare dose distributions and robustness in treatment plans from eight European centres in preparation for the European randomized phase-III PROTECT-trial investigating the effect of proton therapy (PT) versus photon therapy (XT) for oesophageal cancer. MATERIALS AND METHODS: All centres optimized one PT and one XT nominal plan using delineated 4DCT scans for four patients receiving 50.4 Gy (RBE) in 28 fractions. Target volume receiving 95% of prescribed dose (V95%iCTVtotal) should be >99%. Robustness towards setup, range, and respiration was evaluated. The plans were recalculated on a surveillance 4DCT (sCT) acquired at fraction ten and robustness evaluation was performed to evaluate the effect of respiration and inter-fractional anatomical changes. RESULTS: All PT and XT plans complied with V95%iCTVtotal >99% for the nominal plan and V95%iCTVtotal >97% for all respiratory and robustness scenarios. Lung and heart dose varied considerably between centres for both modalities. The difference in mean lung dose and mean heart dose between each pair of XT and PT plans was in median [range] 4.8 Gy [1.1;7.6] and 8.4 Gy [1.9;24.5], respectively. Patients B and C showed large inter-fractional anatomical changes on sCT. For patient B, the minimum V95%iCTVtotal in the worst-case robustness scenario was 45% and 94% for XT and PT, respectively. For patient C, the minimum V95%iCTVtotal was 57% and 72% for XT and PT, respectively. Patient A and D showed minor inter-fractional changes and the minimum V95%iCTVtotal was >85%. CONCLUSION: Large variability in dose to the lungs and heart was observed for both modalities. Inter-fractional anatomical changes led to larger target dose deterioration for XT than PT plans.

Joint EURADOS WG9-WG11 rem-counter intercomparison in a Mevion S250i proton therapy facility with Hyperscan pulsed synchrocyclotron.

Objective Proton therapy is gaining popularity because of the improved dose delivery over conventional radiation therapy. The secondary dose to healthy tissues is dominated by secondary neutrons. Commercial rem-counters are valuable instruments for the on-line assessment of neutron ambient dose equivalent (H*(10)). In general, however, a priori knowledge of the type of facility and of the radiation field is required for the proper choice of any survey meter. The novel Mevion S250i Hyperscan synchrocyclotron mounts the accelerator directly on the gantry. It provides a scanned 227 MeV proton beam, delivered in pulses with a pulse width of 10 µs at 750 Hz frequency, which is afterwards degraded in energy by a range shifter modulator system. This environment is particularly challenging for commercial rem-counters; therefore, we tested the reliability of some of the most widespread rem-counters to understand their limits in the Mevion S250i stray neutron field. Approach This work, promoted by the European Radiation Dosimetry Group (EURADOS), describes a rem-counter intercomparison at the Maastro Proton Therapy centre in the Netherlands, which houses the novel Mevion S250i Hyperscan system. Several rem-counters were employed in the intercomparison (LUPIN, LINUS, WENDI-II, LB6411, NM2B-458, NM2B-495Pb), which included simulation of a patient treatment protocol employing a water tank phantom. The outcomes of the experiment were compared with models and data from the literature. Main results We found that only the LUPIN allowed for a correct assessment of H*(10) within a 20% uncertainty. All other rem-counters underestimated the reference H*(10) by factors from 2 to more than 10, depending on the detector model and on the neutron dose per pulse. In pulsed fields, the neutron dose per pulse is a fundamental parameter, while the average neutron dose rate is a secondary quantity. An average 150-200 µSv/GyRBE neutron H*(10) at various positions around the phantom and at distances between 186 cm and 300 cm from it was measured per unit therapeutic dose delivered to the target. Significance Our results are partially in line with results obtained at similar Mevion facilities employing passive energy modulation. Comparisons with facilities employing active energy modulation confirmed that the neutron H*(10) can increase up to more than a factor of 10 when passive energy modulation is employed. The challenging environment of the Mevion stray neutron field requires the use of specific rem-counters sensitive to high-energy neutrons (up to a few hundred MeV) and specifically designed to withstand pulsed neutron fields.

Visually guided inspiration breath-hold facilitated with nasal high flow therapy in locally advanced lung cancer.

BACKGROUND AND PURPOSE: Reducing breathing motion in radiotherapy (RT) is an attractive strategy to reduce margins and better spare normal tissues. The objective of this prospective study (NCT03729661) was to investigate the feasibility of irradiation of non-small cell lung cancer (NSCLC) with visually guided moderate deep inspiration breath-hold (IBH) using nasal high-flow therapy (NHFT). MATERIAL AND METHODS: Locally advanced NSCLC patients undergoing photon RT were given NHFT with heated humidified air (flow: 40 L/min with 80% oxygen) through a nasal cannula. IBH was monitored by optical surface tracking (OST) with visual feedback. At a training session, patients had to hold their breath as long as possible, without and with NHFT. For the daily cone beam CT (CBCT) and RT treatment in IBH, patients were instructed to keep their BH as long as it felt comfortable. OST was used to analyze stability and reproducibility of the BH, and CBCT to analyze daily tumor position. Subjective tolerance was measured with a questionnaire at 3 time points. RESULTS: Of 10 included patients, 9 were treated with RT. Seven (78%) completed the treatment with NHFT as planned. At the training session, the mean BH length without NHFT was 39 s (range 15-86 s), and with NHFT 78 s (range 29-223 s) (p = .005). NHFT prolonged the BH duration by a mean factor of 2.1 (range 1.1-3.9s). The mean overall stability and reproducibility were within 1 mm. Subjective tolerance was very good with the majority of patients having no or minor discomfort caused by the devices. The mean inter-fraction tumor position variability was 1.8 mm (-1.1-8.1 mm;SD 2.4 mm). CONCLUSION: NHFT for RT treatment of NSCLC in BH is feasible, well tolerated and significantly increases the breath-hold duration. Visually guided BH with OST is stable and reproducible. We therefore consider this an attractive patient-friendly approach to treat lung cancer patients with RT in BH.

A Monte Carlo based scatter removal method for non-isocentric cone-beam CT acquisitions using a deep convolutional autoencoder.

The primary cone-beam computed tomography (CBCT) imaging beam scatters inside the patient and produces a contaminating photon fluence that is registered by the detector. Scattered photons cause artifacts in the image reconstruction, and are partially responsible for the inferior image quality compared to diagnostic fan-beam CT. In this work, a deep convolutional autoencoder (DCAE) and projection-based scatter removal algorithm were constructed for the ImagingRingTM system on rails (IRr), which allows for non-isocentric acquisitions around virtual rotation centers with its independently rotatable source and detector arms. A Monte Carlo model was developed to simulate (i) a non-isocentric training dataset of ≈1200 projection pairs (primary + scatter) from 27 digital head-and-neck cancer patients around five different virtual rotation centers (DCAENONISO), and (ii) an isocentric dataset existing of ≈1200 projection pairs around the physical rotation center (DCAEISO). The scatter removal performance of both DCAE networks was investigated in two digital anthropomorphic phantom simulations and due to superior performance only the DCAENONISO was applied on eight real patient acquisitions. Measures for the quantitative error, the signal-to-noise ratio, and the similarity were evaluated for two simulated digital head-and-neck patients, and the contrast-to-noise ratio (CNR) was investigated between muscle and adipose tissue in the real patient image reconstructions. Image quality metrics were compared between the uncorrected data, the currently implemented heuristic scatter correction data, and the DCAE corrected image reconstruction. The DCAENONISO corrected image reconstructions of two digital patient simulations showed superior image quality metrics compared to the uncorrected and corrected image reconstructions using a heuristic scatter removal. The proposed DCAENONISO scatter correction in this study was successfully demonstrated in real non-isocentric patient CBCT acquisitions and achieved statistically significant higher CNRs compared to the uncorrected or the heuristic corrected image data. This paper presents for the first time a projection-based scatter removal algorithm for isocentric and non-isocentric CBCT imaging using a deep convolutional autoencoder trained on Monte Carlo composed datasets. The algorithm was successfully applied to real patient data.

Beam commissioning of the first compact proton therapy system with spot scanning and dynamic field collimation.

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.

Deriving the mean excitation energy map from dual-energy and proton computed tomography.

The mean excitation energy, I, is an essential quantity for proton treatment planning. This work investigated the feasibility of extracting the spatial distribution of I by combining two computed tomography (CT) modalities, dual-energy CT and proton CT, which provided the spatial distribution of the relative electron density and the stopping power relative to water, respectively. We provided the analytical derivation of I as well as its uncertainty. Results were validated on simulated X-ray and proton CT images of a digital anthropomorphic phantom. Accuracy was below 15% with a large uncertainty, which demonstrated the potential and limits of the technique.

Optimization of dual-energy CT acquisitions for proton therapy using projection-based decomposition.

PURPOSE: Dual-energy computed tomography (DECT) has been presented as a valid alternative to single-energy CT to reduce the uncertainty of the conversion of patient CT numbers to proton stopping power ratio (SPR) of tissues relative to water. The aim of this work was to optimize DECT acquisition protocols from simulations of X-ray images for the treatment planning of proton therapy using a projection-based dual-energy decomposition algorithm. METHODS: We have investigated the effect of various voltages and tin filtration combinations on the SPR map accuracy and precision, and the influence of the dose allocation between the low-energy (LE) and the high-energy (HE) acquisitions. For all spectra combinations, virtual CT projections of the Gammex phantom were simulated with a realistic energy-integrating detector response model. Two situations were simulated: an ideal case without noise (infinite dose) and a realistic situation with Poisson noise corresponding to a 20 mGy total central dose. To determine the optimal dose balance, the proportion of LE-dose with respect to the total dose was varied from 10% to 90% while keeping the central dose constant, for four dual-energy spectra. SPR images were derived using a two-step projection-based decomposition approach. The ranges of 70 MeV, 90 MeV, and 100 MeV proton beams onto the adult female (AF) reference computational phantom of the ICRP were analytically determined from the reconstructed SPR maps. RESULTS: The energy separation between the incident spectra had a strong impact on the SPR precision. Maximizing the incident energy gap reduced image noise. However, the energy gap was not a good metric to evaluate the accuracy of the SPR. In terms of SPR accuracy, a large variability of the optimal spectra was observed when studying each phantom material separately. The SPR accuracy was almost flat in the 30-70% LE-dose range, while the precision showed a minimum slightly shifted in favor of lower LE-dose. Photon noise in the SPR images (20 mGy dose) had lower impact on the proton range accuracy as comparable results were obtained for the noiseless situation (infinite dose). Root-mean-square range errors averaged over all irradiation angles associated to dual-energy imaging were comprised between 0.50 mm and 0.72 mm for the noiseless situation and between 0.51 mm and 0.77 mm for the realistic scenario. CONCLUSIONS: The impact of the dual-energy spectra and the dose allocation between energy levels on the SPR accuracy and precision determined through a projection-based dual-energy algorithm were evaluated to guide the choice of spectra for dual-energy CT for proton therapy. The dose balance between energy levels was not found to be sensitive for the SPR estimation. The optimal pair of dual-energy spectra was material dependent but on a heterogeneous anthropomorphic phantom, there was no significant difference in range accuracy and the choice of spectra could be driven by the precision, i.e., the energy gap.

Technical Note: Procedure for the calibration and validation of kilo-voltage cone-beam CT models.

PURPOSE: The aim of this work is to propose a general and simple procedure for the calibration and validation of kilo-voltage cone-beam CT (kV CBCT) models against experimental data. METHODS: The calibration and validation of the CT model is a two-step procedure: the source model then the detector model. The source is described by the direction dependent photon energy spectrum at each voltage while the detector is described by the pixel intensity value as a function of the direction and the energy of incident photons. The measurements for the source consist of a series of dose measurements in air performed at each voltage with varying filter thicknesses and materials in front of the x-ray tube. The measurements for the detector are acquisitions of projection images using the same filters and several tube voltages. The proposed procedure has been applied to calibrate and assess the accuracy of simple models of the source and the detector of three commercial kV CBCT units. If the CBCT system models had been calibrated differently, the current procedure would have been exclusively used to validate the models. Several high-purity attenuation filters of aluminum, copper, and silver combined with a dosimeter which is sensitive to the range of voltages of interest were used. A sensitivity analysis of the model has also been conducted for each parameter of the source and the detector models. RESULTS: Average deviations between experimental and theoretical dose values are below 1.5% after calibration for the three x-ray sources. The predicted energy deposited in the detector agrees with experimental data within 4% for all imaging systems. CONCLUSIONS: The authors developed and applied an experimental procedure to calibrate and validate any model of the source and the detector of a CBCT unit. The present protocol has been successfully applied to three x-ray imaging systems. The minimum requirements in terms of material and equipment would make its implementation suitable in most clinical environments.