A diamond detector based dosimetric system for instantaneous dose rate measurements in FLASH electron beams.

Objective.A reliable determination of the instantaneous dose rate (I-DR) delivered in FLASH radiotherapy treatments is believed to be crucial to assess the so-called FLASH effect in preclinical and biological studies. At present, no detectors nor real-time procedures are available to do that in ultra high dose rate (UH-DR) electron beams, typically consisting ofµs pulses characterized by I-DRs of the order of MGy/s. A dosimetric system is proposed possibly overcoming the above reported limitation, based on the recently developed flashDiamond (fD) detector (model 60025, PTW-Freiburg, Germany).Approach.A dosimetric system is proposed, based on a flashDiamond detector prototype, properly modified and adapted for very fast signal transmission. It was used in combination with a fast transimpedance amplifier and a digital oscilloscope to record the temporal traces of the pulses delivered by an ElectronFlash linac (SIT S.p.A., Italy). The proposed dosimetric systems was investigated in terms of the temporal characteristics of its response and the capability to measure the absolute delivered dose and instantaneous dose rate (I-DR). A 'standard' flashDiamond was also investigated and its response compared with the one of the specifically designed prototype.Main results. Temporal traces recorded in several UH-DR irradiation conditions showed very good signal to noise ratios and rise and decay times of the order of a few tens ns, faster than the ones obtained by the current transformer embedded in the linac head. By analyzing such signals, a calibration coefficient was derived for the fD prototype and found to be in agreement within 1% with the one obtained under reference60Co irradiation. I-DRs as high as about 2 MGy s-1were detected without any undesired saturation effect. Absolute dose per pulse values extracted by integrating the I-DR signals were found to be linear up to at least 7.13 Gy and in very good agreement with the ones obtained by connecting the fD to a UNIDOS electrometer (PTW-Freiburg, Germany). A good short term reproducibility of the linac output was observed, characterized by a pulse-to-pulse variation coefficient of 0.9%. Negligible differences were observed when replacing the fD prototype with a standard one, with the only exception of a somewhat slower response time for the latter detector type.Significance.The proposed fD-based system was demonstrated to be a suitable tool for a thorough characterization of UH-DR beams, providing accurate and reliable time resolved I-DR measurements from which absolute dose values can be straightforwardly derived.

A new calculation method for the free electron fraction of an ionization chamber in the ultra-high-dose-per-pulse regimen.

The free electron fraction is the fraction of electrons, produced inside the cavity of an ionization chamber after irradiation, which does not bind to gas molecules and thereby reaches the electrode as free electrons. It is a fundamental quantity to describe the recombination processes of an ionization chamber, as it generates a gap of positive charges compared to negative ones, which certainly will not undergo recombination. The free electron fraction depends on the specific chamber geometry, the polarizing applied voltage and the gas thermodynamic properties. Therefore, it is necessary to evaluate such fraction in an accurate and easy way for any measurement condition. In this paper, a simple and direct method for evaluating the free electron fraction of ionization chambers is proposed. We first model the capture process of the electrons produced inside an ionization chamber after the beam pulse; then we present a method to evaluate the free electron fraction based on simple measurements of collected charge, by varying the applied voltage. Finally, the results obtained using an Advanced Markus chamber irradiated with a Flash Radiotherapy dedicated research Linac (ElectronFlash) to estimate the free electron fraction are presented. The proposed method allows the use of a conventional ionization chamber for measurements in ultra-high-dose-per-pulse (UHDP) conditions, up to values of dose-per-pulse at which the perturbation of the electric field due to the generated charge can be considered negligible.

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO.

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam's Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

A new solution for UHDP and UHDR (Flash) measurements: Theory and conceptual design of ALLS chamber.

Ultra-High dose-per-pulse regimens (UHDP), necessary to trigger the "FLASH" effect, still pose serious challenges to dosimetry. Dosimetry plays a crucial role, both to significantly improve the accuracy of the radiobiological experiments necessary to fully understand the mechanisms underlying the effect and its dependencies on the beam parameters, and to be able to translate such effect into clinical practice. The standard ionization chamber in UHDP region is significantly affected by the effects of the electric field generated by the enormous density of charges produced by the dose pulse. This work describes the theory and the conceptual design of a gas chamber (the ALLS chamber) which overcomes the above-mentioned problems.

Localization of anatomical changes in patients during proton therapy with in-beam PET monitoring: A voxel-based morphometry approach exploiting Monte Carlo simulations.

PURPOSE: In-beam positron emission tomography (PET) is one of the modalities that can be used for in vivo noninvasive treatment monitoring in proton therapy. Although PET monitoring has been frequently applied for this purpose, there is still no straightforward method to translate the information obtained from the PET images into easy-to-interpret information for clinical personnel. The purpose of this work is to propose a statistical method for analyzing in-beam PET monitoring images that can be used to locate, quantify, and visualize regions with possible morphological changes occurring over the course of treatment. METHODS: We selected a patient treated for squamous cell carcinoma (SCC) with proton therapy, to perform multiple Monte Carlo (MC) simulations of the expected PET signal at the start of treatment, and to study how the PET signal may change along the treatment course due to morphological changes. We performed voxel-wise two-tailed statistical tests of the simulated PET images, resembling the voxel-based morphometry (VBM) method commonly used in neuroimaging data analysis, to locate regions with significant morphological changes and to quantify the change. RESULTS: The VBM resembling method has been successfully applied to the simulated in-beam PET images, despite the fact that such images suffer from image artifacts and limited statistics. Three dimensional probability maps were obtained, that allowed to identify interfractional morphological changes and to visualize them superimposed on the computed tomography (CT) scan. In particular, the characteristic color patterns resulting from the two-tailed statistical tests lend themselves to trigger alarms in case of morphological changes along the course of treatment. CONCLUSIONS: The statistical method presented in this work is a promising method to apply to PET monitoring data to reveal interfractional morphological changes in patients, occurring over the course of treatment. Based on simulated in-beam PET treatment monitoring images, we showed that with our method it was possible to correctly identify the regions that changed. Moreover we could quantify the changes, and visualize them superimposed on the CT scan. The proposed method can possibly help clinical personnel in the replanning procedure in adaptive proton therapy treatments.

Normalized glandular dose coefficients for digital breast tomosynthesis systems with a homogeneous breast model.

This work aims at calculating and releasing tabulated values of dose conversion coefficients, DgNDBT, for mean glandular dose (MGD) estimates in digital breast tomosynthesis (DBT). The DgNDBT coefficients are proposed as unique conversion coefficients for MGD estimates, in place of dose conversion coefficients in mammography (DgNDM or c, g, s triad as proposed in worldwide quality assurance protocols) used together with the T correction factor. DgNDBT is the MGD per unit incident air kerma measured at the breast surface for a 0° projection and the entire tube load used for the scan. The dataset of polyenergetic DgNDBT coefficients was derived via a Monte Carlo software based on the Geant4 toolkit. Dose coefficients were calculated for a grid of values of breast characteristics (breast thickness in the range 20-90 mm and glandular fraction by mass of 1%, 25%, 50%, 75%, 100%) and the simulated geometries, scan protocols, irradiation geometries and typical spectral qualities replicated those of six commercial DBT systems (GE SenoClaire, Hologic Selenia Dimensions, GE Senographe Pristina, Fujifilm Amulet Innovality, Siemens Mammomat Inspiration and IMS Giotto Class). For given breast characteristics, target/filter combination, tube voltage and half value layer (HVL), two spectra with two HVL values have been simulated in order to permit MGD estimates from experimental HVL values via mathematical interpolation from tabulated values. The adopted breast model assumes homogenous composition of glandular and adipose tissues; it includes a 1.45 mm thick skin envelope in place of the 4-5 mm envelope commonly adopted in dosimetry protocols. The simulation code was validated versus AAPM Task group 195 Monte Carlo reference data sets (absolute differences not higher than 1.1%) and by comparison to relative dosimetry measurements with radiochromic film in a PMMA test object (differences within the maximum experimental uncertainty of 11%). The calculated coefficients show maximum relative deviations of -17.6% and +6.1% from those provided by the DBT dose coefficients adopted in the EUREF protocol and of 1.5%, on average, from data in the AAPM TG223 report. A spreadsheet is provided for interpolating the tabulated DgNDBT coefficients for arbitrary values of HVL, compressed breast thickness and glandular fraction, in the corresponding investigated ranges, for each DBT unit modeled in this work.