Calibration method and performance of a time-of-flight detector to measure absolute beam energy in proton therapy.

BACKGROUND: The beam energy is one of the most significant parameters in particle therapy since it is directly correlated to the particles' penetration depth inside the patient. Nowadays, the range accuracy is guaranteed by offline routine quality control checks mainly performed with water phantoms, 2D detectors with PMMA wedges, or multi-layer ionization chambers. The latter feature low sensitivity, slow collection time, and response dependent on external parameters, which represent limiting factors for the quality controls of beams delivered with fast energy switching modalities, as foreseen in future treatments. In this context, a device based on solid-state detectors technology, able to perform a direct and absolute beam energy measurement, is proposed as a viable alternative for quality assurance measurements and beam commissioning, paving the way for online range monitoring and treatment verification. PURPOSE: This work follows the proof of concept of an energy monitoring system for clinical proton beams, based on Ultra Fast Silicon Detectors (featuring tenths of ps time resolution in 50 µm active thickness, and single particle detection capability) and time-of-flight techniques. An upgrade of such a system is presented here, together with the description of a dedicated self-calibration method, proving that this second prototype is able to assess the mean particles energy of a monoenergetic beam without any constraint on the beam temporal structure, neither any a priori knowledge of the beam energy for the calibration of the system. METHODS: A new detector geometry, consisting of sensors segmented in strips, has been designed and implemented in order to enhance the statistics of coincident protons, thus improving the accuracy of the measured time differences. The prototype was tested on the cyclotron proton beam of the Trento Protontherapy Center (TPC). In addition, a dedicated self-calibration method, exploiting the measurement of monoenergetic beams crossing the two telescope sensors for different flight distances, was introduced to remove the systematic uncertainties independently from any external reference. RESULTS: The novel calibration strategy was applied to the experimental data collected at TPC (Trento) and CNAO (Pavia). Deviations between measured and reference beam energies in the order of a few hundreds of keV with a maximum uncertainty of 0.5 MeV were found, in compliance with the clinically required water range accuracy of 1 mm. CONCLUSIONS: The presented version of the telescope system, minimally perturbative of the beam, relies on a few seconds of acquisition time to achieve the required clinical accuracy and therefore represents a feasible solution for beam commission, quality assurance checks, and online beam energy monitoring.

Proton therapy monitoring: spatiotemporal emission reconstruction with prompt gamma timing and implementation with PET detectors.

Objective. In this study we introduce spatiotemporal emission reconstruction prompt gamma timing (SER-PGT), a new method to directly reconstruct the prompt photon emission in the space and time domains inside the patient in proton therapy.Approach. SER-PGT is based on the numerical optimisation of a multidimensional likelihood function, followed by a post-processing of the results. The current approach relies on a specific implementation of the maximum-likelihood expectation maximisation algorithm. The robustness of the method is guaranteed by the complete absence of any information about the target composition in the algorithm.Main results. Accurate Monte Carlo simulations indicate a range resolution of about 0.5 cm (standard deviation) when considering 107primary protons impinging on an homogeneous phantom. Preliminary results on an anthropomorphic phantom are also reported.Significance. By showing the feasibility for the reconstruction of the primary particle range using PET detectors, this study provides significant basis for the development of an hybrid in-beam PET and prompt photon device.

Accuracy assessment of the CNAO dose delivery system in the initial period of clinical activity and impact of later improvements on delivered dose distributions.

PURPOSE: A retrospective analysis of the dose delivery system (DDS) performances of the initial clinical operation at CNAO (Centro Nazionale di Adroterapia Oncologica) is reported, and compared with the dose delivery accuracy following the implementation of a position feedback control. METHODS: Log files and raw data of the DDS were analyzed for every field of patients treated with protons and carbon ions between January 2012 and April 2013 (~3800 fields). To investigate the DDS accuracy, the spot positions and the number of particles per spot measured by the DDS and prescribed by the treatment planning system were compared for each field. The impact of deviations on dose distributions was studied by comparing, through the gamma-index method, 2 three-dimensional (3D) physical dose maps (one for prescribed, one for measured data), generated by a validated dose computation software. The maximum gamma and the percentage of points with gamma ≤ 1 (passing volume) were studied as a function of the treatment day, and correlated with the deviations from the prescription in the measured number of particles and spot positions. Finally, delivered dose distributions of same treatment plans were compared before and after the implementation of a feedback algorithm for the correction of small position deviations, to study the effect on the delivery quality. A double comparison of prescribed and measured 3D maps, before and after feedback implementation, is reported and studied for a representative treatment delivered in 2012, redelivered on a polymethyl methacrylate (PMMA) block in 2018. RESULTS: Systematic deviations of spot positions, mainly due to beam lateral offsets, were always found within 1.5 mm, with the exception of the initial clinical period. The number of particles was very stable, as possible deviations are exclusively related to the quantization error in the conversion from monitor counts to particles. For the chosen representative patient treatment, the gamma-index evaluation of prescribed and measured dose maps, before and after feedback implementation, showed a higher variability of maximum gamma for the 2012 irradiation, with respect to the reirradiation of 2018. However, the 2012 passing volume is >99.8% for the sum of all fields, which is comparable to the value of 2018, with the exception of one day with 98.2% passing volume, probably related to an instability of the accelerating system. CONCLUSIONS: A detailed retrospective analysis of the DDS performances in the initial period of CNAO clinical activity is reported. The spot position deviations are referable to beam lateral offset fluctuations, while almost no deviation was found in the number of particles. The impact of deviations on dose distributions showed that the position feedback implementation and the increased beam control capability acquired after the first years of clinical experience led to an evident improvement in the DDS stability, evaluated in terms of gamma-index as a measure of the impact on dose distributions. However, the clinical effect of the maximum gamma variability found in the 2012 representative irradiation is mitigated by averaging along the number of fractions, and the high percentage of passing volumes confirmed the accuracy of the delivery even before the feedback implementation.

Dosimetry in Radiosynoviorthesis: 90Y VS. 153Sm.

Although there are several radionuclides suitable for radiosynoviorthesis (RSO), not all of them can irradiate deeper synovium. Yttrium-90 (Y) is the beta radionuclide with more penetration range; therefore, it is predominantly used to treat knees. The aim of this paper is to highlight several dosimetry concepts to compare Y and Sm, also discussing the feasibility of implementing a dose planning methodology for both in RSO. The MCNPX Monte Carlo nuclear code version 2.6 was used for calculating S-values from which the activity to be injected into the joint was obtained. This activity is considered sufficient to deliver a 100-Gy absorbed dose in 1 mm of synovial tissue. The simulated mathematical model consisted of a system formed by several cylindrical slabs of 1-mm thickness, aligned consecutively. The different areas of the cylinder base simulate several synovial membrane sizes. The effective treatment range for each radionuclide was also calculated. Quantification of the synovial joint features (synovial thickness and synovial surface) by diagnostic imaging, such as magnetic resonance (MRI) combined with a Monte Carlo simulation, can be used to achieve a treatment planning strategy in RSO with the available radionuclides.

A new approach for radiosynoviorthesis: A dose-optimized planning method based on Monte Carlo simulation and synovial measurement using 3D slicer and MRI.

PURPOSE: Recently, there has been a growing interest in a methodology for dose planning in radiosynoviorthesis to substitute fixed activity. Clinical practice based on fixed activity frequently does not embrace radiopharmaceutical dose optimization in patients. The aim of this paper is to propose and discuss a dose planning methodology considering the radiological findings of interest obtained by three-dimensional magnetic resonance imaging combined with Monte Carlo simulation in radiosynoviorthesis treatment applied to hemophilic arthropathy. METHOD: The parameters analyzed were: surface area of the synovial membrane (synovial size), synovial thickness and joint effusion obtained by 3D MRI of nine knees from nine patients on a SIEMENS AVANTO 1.5 T scanner using a knee coil. The 3D Slicer software performed both the semiautomatic segmentation and quantitation of these radiological findings. A Lucite phantom 3D MRI validated the quantitation methodology. The study used Monte Carlo N-Particle eXtended code version 2.6 for calculating the S-values required to set up the injected activity to deliver a 100 Gy absorbed dose at a determined synovial thickness. The radionuclides assessed were: 90Y, 32P, 188Re, 186Re, 153Sm, and 177Lu, and the present study shows their effective treatment ranges. RESULT: The quantitation methodology was successfully tested, with an error below 5% for different materials. S-values calculated could provide data on the activity to be injected into the joint, considering no extra-articular leakage from joint cavity. Calculation of effective treatment range could assist with the therapeutic decision, with an optimized protocol for dose prescription in RSO. CONCLUSION: Using 3D Slicer software, this study focused on segmentation and quantitation of radiological features such as joint effusion, synovial size, and thickness, all obtained by 3D MRI in patients' knees with hemophilic arthropathy. The combination of synovial size and thickness with the parameters obtained by Monte Carlo simulation such as effective treatment range and S-value, from which is calculated the injected activity, could be used for treatment planning in RSO. Data from this methodology could be a potential aid to clinical decision making by selecting the most suitable radionuclide; justifying the procedure, fractioning the dose, and the calculated injected activity for children and adolescents, considering both the synovial size and thickness.